We study a metric cubic gravity theory considering odd-parity modes of linear inhomogeneous perturbations on a spatially homogeneous Bianchi type I manifold close to the isotropic de Sitter spacetime. We show that in the regime of small anisotropy, the theory possesses new degrees of freedom compared to General Relativity, whose kinetic energy vanishes in the limit of exact isotropy. From the mass dispersion relation we show that such theory always possesses at least one ghost mode as well as a very short-time-scale (compared to the Hubble time) classical tachyonic (or ghost-tachyonic) instability. In order to confirm our analytic analysis, we also solve the equations of motion numerically and we find that this instability is developed well before a single e-fold of the scale factor. This shows that this gravity theory, as it is, cannot be used to construct viable cosmological models.

We study the phenomenology of a class of minimally modified gravity theories called $f(\mathcal{H})$ theories, in which the usual general relativistic Hamiltonian constraint is replaced by a free function of it. After reviewing the construction of the theory and a consistent matter coupling, we analyze the dynamics of cosmology at the levels of both background and perturbations, and present a concrete example of the theory with a $3$-parameter family of the function $f$. Finally, we compare this example model to Planck data as well as some later-time probes, showing that such a realization of $f(\mathcal{H})$ theories fits the data significantly better than the standard $\Lambda$CDM model, in particular by modifying gravity at intermediate redshifts, $z\simeq743$.

We propose a modified gravity theory that propagates only two local gravitational degrees of freedom and that does not have an Einstein frame. According to the classification in JCAP 01 (2019) 017 [arXiv:1810.01047 [gr-qc]], this is a type-II minimally modified gravity theory. The theory is characterized by the gravitational constant $G_{\rm N}$ and a function $V(\phi)$ of a non-dynamical auxiliary field $\phi$ that plays the role of dark energy. Once one fixes a homogeneous and isotropic cosmological background, the form of $V(\phi)$ is determined and the theory no longer possesses a free parameter or a free function, besides $G_{\rm N}$. For $V'(\phi) = 0$ the theory reduces to general relativity (GR) with $G_N$ being the Newton's constant and $V=const.$ being the cosmological constant. For $V'(\phi) \ne 0$, it is shown that gravity behaves differently from GR but that GR with $G_{\rm N}$ being the Newton's constant is recovered for weak gravity at distance and time scales sufficiently shorter than the scale associated with $V(\phi)$. Therefore this theory provides the simplest framework of cosmology in which deviations from GR can be tested by observational data.

We study a coupled dark energy scenario in which a massive vector field $A_{\mu}$ with broken $U(1)$ gauge symmetry interacts with the four-velocity $u_c^{\mu}$ of cold dark matter (CDM) through the scalar product $Z=-u_c^{\mu} A_{\mu}$. This new coupling corresponds to the momentum transfer, so that the background vector and CDM continuity equations do not have explicit interacting terms analogous to the energy exchange. Hence the observational preference of uncoupled generalized Proca theories over the $\Lambda$CDM model can be still maintained at the background level. Meanwhile, the same coupling strongly affects the evolution of cosmological perturbations. While the effective sound speed of CDM vanishes, the propagation speed and no-ghost condition of a longitudinal scalar of $A_{\mu}$ and the CDM no-ghost condition are subject to nontrivial modifications by the $Z$ dependence in the Lagrangian. We propose a concrete dark energy model and show that the gravitational interaction on scales relevant to the linear growth of large-scale structures can be smaller than the Newton constant at low redshifts. This leads to the suppression of growth rates of both CDM and total matter density perturbations, so our model allows an interesting possibility for reducing the tension of matter density contrast $\sigma_8$ between high- and low-redshift measurements.

We study the possibility to obtain cosmological late-time acceleration from a geometry changing with the scale, in particular, in the so-called multi-fractional theories with q-derivatives and with weighted derivatives. In the theory with q-derivatives, the luminosity distance is the same as in general relativity and, therefore, geometry cannot act as dark energy. In the theory with weighted derivatives, geometry alone is able to sustain a late-time acceleration phase without fine tuning, while being compatible with structure-formation and big-bang nucleosynthesis bounds. This suggests to extend the theory, in a natural way, from just small-scale to also large-scale modifications of gravity. Surprisingly, the Hausdorff dimension of spacetime is constrained to be close to the topological dimension 4. Promoting this finding to a principle, we conclude that present-day acceleration can be regarded as the effect of a new restoration law for spacetime geometry.

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 investigate the cosmological viability of the generalized proca theory. We first implement the background and linear perturbation equations of motion in the Boltzmann code and then study the constraints on the parameters of the generalized proca theory after running MCMC against the cosmological data set. With Planck + HST data, we obtain the constraint $h=0.7334_{-0.0269}^{+0.0246}$, which indicates that the tension between early universe and late time universe within this theory is removed. By adding other late-time data sets (BAO, RSD, etc.) we show that the tension is reduced, as the 2$\sigma$ allowed region for $h$ in Proca, $h=0.7041_{-0.0087}^{+0.0094}$, overlaps with the 2$\sigma$ region of the HST data.

We provide two derivations of the baryonic equations that can be straightforwardly implemented in existing Einstein--Boltzmann solvers. One of the derivations begins with an action principle, while the other exploits the conservation of the stress-energy tensor. While our result is manifestly covariant and satisfies the Bianchi identities, we point out that this is not the case for the implementation of the seminal work by Ma and Bertschinger and in the existing Boltzmann codes. We also study the tight coupling approximation up to the second order without choosing any gauge using the covariant full baryon equations. We implement the improved baryon equations in a Boltzmann code and investigate the change in the estimate of cosmological parameters by performing an MCMC analysis. With the covariantly correct baryon equations of motion, we find 1% deviation for the best fit values of the cosmological parameters that should be taken into account. While in this paper, we study the Lambda-CDM model only, our baryon equations can be easily implemented in other models and various modified gravity theories.

The gravitational-wave event GW170817, together with the electromagnetic counterpart, shows that the speed of tensor perturbations $c_T$ on the cosmological background is very close to that of light $c$ for the redshift $z<0.009$. In generalized Proca theories, the Lagrangians compatible with the condition $c_T=c$ are constrained to be derivative interactions up to cubic order, besides those corresponding to intrinsic vector modes. We place observational constraints on a dark energy model in cubic-order generalized Proca theories with intrinsic vector modes by running the Markov chain Monte Carlo (MCMC) code. We use the cross-correlation data of the integrated Sachs-Wolfe (ISW) signal and galaxy distributions in addition to the data sets of cosmic microwave background, baryon acoustic oscillations, type Ia supernovae, local measurements of the Hubble expansion rate, and redshift-space distortions. We show that, unlike cubic-order scalar-tensor theories, the existence of intrinsic vector modes allows the possibility for evading the ISW-galaxy anticorrelation incompatible with the current observational data. As a result, we find that the dark energy model in cubic-order generalized Proca theories exhibits a better fit to the data than the cosmological constant, even by including the ISW-galaxy correlation data in the MCMC analysis.

Several modified cosmological models exist, which also try to address the tensions between data and predictions of the $\Lambda$-CDM model. Galileon models are particular scalar tensor theories that represent one such possibilities. While it is commonly understood that there may be inconsistencies between predictions of some Galileon models and observations, in particular concerning ISW-galaxy cross-correlations, there is no proof yet that these models are completely ruled out. Indeed, by using a specific background in the generalized covariant Galileon theory known as the the tracker solution, here we show that, after imposing all standard theoretical stability constraints, it is still possible to identify a region in the parameter space of the model that allows for positive ISW-galaxy cross-correlations. By a physical interpretation in terms of a chi-square analysis, we confirm the expectation that in this viable region the predictions of generalized covariant Galileon theory on the tracker solution background have higher likelihood when they approach the physics of the $\Lambda$-CDM model.

We study cosmology in a class of minimally modified gravity (MMG) with two local gravitational degrees of freedom. We classify modified gravity theories into type-I and type-II: theories of type-I have an Einstein frame and can be recast by change of variables as general relativity (GR) with a non-minimal matter coupling, while theories of type-II have no Einstein frame. Considering a canonical transformation of the lapse, the 3-dimensional induced metric and their conjugate momenta we generate type-I MMG. We then show that phenomenological deviations from GR, such as the speed of gravitational waves $c_T$ and the effective gravitational constant for scalar perturbations $G_{\rm eff}$, are characterized by two functions of an auxiliary variable. We study the phenomenology of several models all having $c_T=1$. We obtain a scenario with $c_T=1$ in which the effective equation-of-state parameter of dark energy is different from $-1$ even though the cosmic acceleration is caused by a bare cosmological constant, and we find that it is possible to reconstruct the theory on choosing a selected time-evolution for the effective dark energy component.

The minimal theory of massive gravity (MTMG) has two branches of stable cosmological solutions: a self-accelerating branch, which, except for the mass of tensor modes has exactly the same behavior of linear perturbations as $\Lambda$CDM in general relativity (GR), and a normal branch with nontrivial behavior. We explore the influence of the integrated Sachs-Wolfe-galaxy correlation constraints on the normal branch of MTMG, which, in its simplest implementation, has one free parameter more than $\Lambda$CDM in GR (or the self-accelerating branch of MTMG): $\theta$. This parameter is related to the graviton mass and only affects the behavior of the cosmological linear perturbation dynamics. Using 2d-mass and SDSS data, we check which values of $\theta$ lead to a positive or negative cross-correlation. We find that positive cross-correlation is achieved for a large parameter-space interval. Within this allowed region of parameter space, we perform a $\chi^2$ analysis in terms of the parameter $\theta$, while keeping the other background parameters fixed to the best-fit values of Planck. We then infer that the normal branch of MTMG fits the data well in a nontrivial portion of the parameter space, and future experiments should be able to distinguish such a model from $\Lambda$CDM in GR (or the self-accelerating branch of MTMG).

The minimal theory of quasidilaton massive gravity with or without a Horndeski-type kinetic term for the quasidilaton field propagates only three physical modes: the two massive tensor polarizations and one scalar mode. This reduced number of degrees of freedom is realized by a Lorentz symmetry violation at cosmological scales and the presence of appropriate constraints that remove unwanted modes. Vacuum cosmological solutions have been considered in a previous work, and it has been shown that the late-time de Sitter attractor is stable under inhomogeneous perturbations. In this work, we explore the stability of cosmological solutions in the presence of matter fields. Assuming for simplicity that the quasidilaton scalar is on an attractor at the level of the background, we derive stability conditions in the subhorizon limit, and find the scalar sound speeds, as well as the modification with respect to general relativity to the gravitational potential in the quasistatic approximation. We also find that the speed limit of gravitational waves coincides with the speed of light for any homogeneous and isotropic cosmological background, on or away from the attractor.

The recently proposed chameleonic extension of bigravity theory, by including a scalar field dependence in the graviton potential, avoids several fine-tunings found to be necessary in usual massive bigravity. In particular it ensures that the Higuchi bound is satisfied at all scales, that no Vainshtein mechanism is needed to satisfy solar-system experiments, and that the strong coupling scale is always above the scale of cosmological interest all the way up to the early universe. This paper extends the previous work by presenting a stable example of cosmology in the chameleon bigravity model. We find a set of initial conditions and parameters such that the derived stability conditions on general flat Friedmann background are satisfied at all times. The evolution goes through radiation dominated, matter dominated, and de Sitter eras. We argue that the parameter space allowing for such a stable evolution may be large enough to encompass an observationally viable evolution. We also argue that our model satisfies all known constraints due to gravitational wave observations so far and thus can be considered as a unique testing ground of gravitational wave phenomenologies in bimetric theories of gravity.

This article extends bimetric formulations of massive gravity to make the mass of the graviton to depend on its environment. This minimal extension offers a novel way to reconcile massive gravity with local tests of general relativity without invoking the Vainshtein mechanism. On cosmological scales, it is argued that the model is stable and that it circumvents the Higuchi bound, hence relaxing the constraints on the parameter space. Moreover, with this extension the strong coupling scale is also environmentally dependent in such a way that it is kept sufficiently higher than the expansion rate all the way up to the very early universe, while the present graviton mass is low enough to be phenomenologically interesting. In this sense the extended bigravity theory serves as a partial UV completion of the standard bigravity theory. This extension is very generic and robust and a simple specific example is described.

We present results from a controlled numerical experiment investigating the effect of stellar density gas on the coalescence of binary black holes (BBHs) and the resulting gravitational waves (GWs). This investigation is motivated by the proposed stellar core fragmentation scenario for BBH formation and the associated possibility of an electromagnetic counterpart to a BBH GW event. We employ full numerical relativity coupled with general-relativistic hydrodynamics and set up a $30 + 30 M_\odot$ BBH (motivated by GW150914) inside gas with realistic stellar densities. Our results show that at densities $\rho \gtrsim 10^6 - 10^7 \, \mathrm{g \, cm}^{-3}$ dynamical friction between the BHs and gas changes the coalescence dynamics and the GW signal in an unmistakable way. We show that for GW150914, LIGO observations conclusively rule out BBH coalescence inside stellar gas of $\rho \gtrsim 10^7 \, \mathrm{g\,cm}^{-3}$. Typical densities in the collapsing cores of massive stars are in excess of this density. This excludes the fragmentation scenario for the formation of GW150914.

The minimal theory of quasidilaton massive gravity allows for a stable self-accelerating de Sitter solution in a wide range of parameters. On the other hand, in order for the theory to be compatible with local gravity tests, the fifth force due to the quasidilaton scalar needs to be screened at local scales. The present paper thus extends the theory by inclusion of a cubic Horndeski term in a way that (i) respects the quasidilaton global symmetry, that (ii) maintains the physical degrees of freedom in the theory being three, that (iii) can accommodate the Vainshtein screening mechanism and that still (iv) allows for a stable self-accelerating de Sitter solution. After adding the Horndeski term (and a k-essence type nonlinear kinetic term as well) to the precursor action, we switch to the Hamiltonian language and find a complete set of independent constraints. We then construct the minimal theory with three physical degrees of freedom by carefully adding a pair of constraints to the total Hamiltonian of the precursor theory. Switching back to the Lagrangian language, we study cosmological solutions and their stability in the minimal theory. In particular, we show that a self-accelerating de Sitter solution is stable for a wide range of parameters. Furthermore, as in the minimal theory of massive gravity, the propagation speed of the massive gravitational waves in the high momentum limit precisely agrees with the speed of light.

We introduce a quasidilaton scalar field to the minimal theory of massive gravity with the Minkowski fiducial metric, in such a way that the quasidilaton global symmetry is maintained and that the theory admits a stable self-accelerating de Sitter solution. We start with a precursor theory that contains three propagating gravitational degrees of freedom without a quasidilaton scalar and introduce St\"uckelberg fields to covariantize its action. This makes it possible for us to formulate the quasidilaton global symmetry that mixes the St\"uckelberg fields and the quasidilaton scalar field. By the Hamiltonian analysis we confirm that the precursor theory with the quasidilaton scalar contains four degrees of freedom, three from the precursor massive gravity and one from the quasidilaton scalar. We further remove one propagating degree of freedom to construct the minimal quasidilaton theory with three propagating degrees of freedom, corresponding to two polarizations of gravitational waves from the minimal theory of massive gravity and one scalar from the quasidilaton field, by carefully introducing two additional constraints to the system in the Hamiltonian language. Switching to the Lagrangian language, we find self-accelerating de Sitter solutions in the minimal quasidilaton theory and analyze their stability. It is found that the self-accelerating de Sitter solution is stable in a wide range of parameters.

We investigate the minimal theory of massive gravity (MTMG) recently introduced. After reviewing the original construction based on its Hamiltonian in the vielbein formalism, we reformulate it in terms of its Lagrangian in both the vielbein and the metric formalisms. It then becomes obvious that, unlike previous attempts in the literature of Lorentz-violating massive gravity, not only the potential but also the kinetic structure of the action is modified from the de Rham-Gabadadze-Tolley (dRGT) massive gravity theory. We confirm that the number of physical degrees of freedom in MTMG is two at fully nonlinear level. This proves the absence of various possible pathologies such as superluminality, acausality and strong coupling. Afterwards, we discuss the phenomenology of MTMG in the presence of a dust fluid. We find that on a flat homogeneous and isotropic background we have two branches. One of them (self-accelerating branch) naturally leads to acceleration without the genuine cosmological constant or dark energy. For this branch both the scalar and the vector modes behave exactly as in general relativity (GR). The phenomenology of this branch differs from GR in the tensor modes sector, as the tensor modes acquire a non-zero mass. Hence, MTMG serves as a stable nonlinear completion of the self-accelerating cosmological solution found originally in dRGT theory. The other branch (normal branch) has a dynamics which depends on the time-dependent fiducial metric. For the normal branch, the scalar mode sector, even though as in GR only one scalar mode is present (due to the dust fluid), differs from the one in GR, and, in general, structure formation will follow a different phenomenology. The tensor modes will be massive, whereas the vector modes, for both branches, will have the same phenomenology as in GR.

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.

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.

The standard $\Lambda$-CDM predicts a growth of structures which tends to be higher than the values of redshift space distortion (RSD) measurements, if the cosmological parameters are fixed by the CMB data. In this paper we point out that this discrepancy can be resolved/understood if we assume that the graviton has a small but non-zero mass. In the context of the Minimal Theory of Massive Gravity (MTMG), due to infrared Lorentz violations measurable only at present cosmological scales, the graviton acquires a mass without being haunted by unwanted extra degrees of freedom. While the so-called self-accelerating branch of cosmological solutions in MTMG has the same phenomenology for the background as well as the scalar- and vector-type linear perturbations as the $\Lambda$-CDM in General Relativity (GR), it is possible to choose another branch so that the background is the same as that in GR but the evolution of matter perturbations gets modified by the graviton mass. On studying the fit of such modified dynamics to the above-mentioned RSD measurements, we find that the $\Lambda$-CDM model is less probable than MTMG by two orders of magnitude. With the help of the cross correlation between the integrated Sachs-Wolfe (ISW) effect and the large scale structure (LSS), the data also pin-down the graviton mass squared around $\mu^2\approx - (3\times 10^{-33} \rm{eV})^2$, which is consistent with the latest bound $|\mu^2|<(1.2\times 10^{-22} \rm{eV})^2$ set by the recent LIGO observation.

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.

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.

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.

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.

In Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, it is known that the conical singularity arises at the center of a spherically symmetric body ($r=0$) in the case where the parameter $\alpha_{{\rm H}4}$ characterizing the deviation from the Horndeski Lagrangian $L_4$ approaches a non-zero constant as $r \to 0$. We derive spherically symmetric solutions around the center in full GLPV theories and show that the GLPV Lagrangian $L_5$ does not modify the divergent property of the Ricci scalar $R$ induced by the non-zero $\alpha_{{\rm H}4}$. Provided that $\alpha_{{\rm H}4}=0$, curvature scalar quantities can remain finite at $r=0$ even in the presence of $L_5$ beyond the Horndeski domain. For the theories in which the scalar field $\phi$ is directly coupled to $R$, we also obtain spherically symmetric solutions inside/outside the body to study whether the fifth force mediated by $\phi$ can be screened by non-linear field self-interactions. We find that there is one specific model of GLPV theories in which the effect of $L_5$ vanishes in the equations of motion. We also show that, depending on the sign of a $L_5$-dependent term in the field equation, the model can be compatible with solar-system constraints under the Vainshtein mechanism or it is plagued by the problem of a divergence of the field derivative in high-density regions.

We propose a dark energy model where a scalar field $\phi$ has nonlinear self-interactions in the presence of a dilatonic coupling with the Ricci scalar. This belongs to a sub-class of theories beyond Horndeski, which accommodates covariant Galileons and Brans-Dicke theories as specific cases. We derive conditions under which the scalar sound speed squared $c_{\rm s}^2$ is positive from the radiation era to today. Since $c_{\rm s}^2$ remains to be smaller than the order of 1, the deviation from Horndeski theories does not cause heavy oscillations of gauge-invariant gravitational potentials. In this case, the evolution of matter perturbations at low redshifts is similar to that in the coupled dark energy scenario with an enhanced gravitational interaction. On the spherically symmetric background with a matter source, the existence of field self-interactions suppresses the propagation of fifth force inside a Vainshtein radius. We estimate an allowed parameter space in which the model can be compatible with solar system constraints while satisfying conditions for the cosmological viability of background and perturbations.

In a class of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, we derive both vacuum and interior Schwarzschild solutions under the condition that the derivatives of a scalar field $\phi$ with respect to the radius $r$ vanish. If the parameter $\alpha_{\rm H}$ characterizing the deviation from Horndeski theories approaches a non-zero constant at the center of a spherically symmetric body, we find that the conical singularity arises at $r=0$ with the Ricci scalar given by $R=-2\alpha_{\rm H}/r^2$. This originates from violation of the geometrical structure of four-dimensional curvature quantities. The conical singularity can disappear for the models in which the parameter $\alpha_{\rm H}$ vanishes in the limit that $r \to 0$. We propose explicit models without the conical singularity by properly designing the classical Lagrangian in such a way that the main contribution to $\alpha_{\rm H}$ comes from the field derivative $\phi'(r)$ around $r=0$. We show that the extension of covariant Galileons with a diatonic coupling allows for the recovery of general relativistic behavior inside a so-called Vainshtein radius. In this case, both the propagation of a fifth force and the deviation from Horndeski theories are suppressed outside a compact body in such a way that the model is compatible with local gravity experiments inside the solar system.

We propose a new theory of massive gravity with only two propagating degrees of freedom. After defining the theory in the unitary gauge in the vielbein language, we shall perform a Hamiltonian analysis to count the number of physical degrees of freedom, and then study some phenomenologies. While the homogeneous and isotropic background cosmology and the tensor linear perturbations around it are described by exactly the same equations as those in the de Rham-Gabadadze-Tolley (dRGT) massive gravity, the scalar and vector gravitational degrees of freedom are absent in the new theory at the fully nonlinear level. Hence the new theory provides a stable nonlinear completion of the self-accelerating cosmological solution that was originally found in the dRGT theory.

One century after its formulation, Einstein's general relativity has made remarkable predictions and turned out to be compatible with all experimental tests. Most of these tests probe the theory in the weak-field regime, and there are theoretical and experimental reasons to believe that general relativity should be modified when gravitational fields are strong and spacetime curvature is large. The best astrophysical laboratories to probe strong-field gravity are black holes and neutron stars, whether isolated or in binary systems. We review the motivations to consider extensions of general relativity. We present a (necessarily incomplete) catalog of modified theories of gravity for which strong-field predictions have been computed and contrasted to Einstein's theory, and we summarize our current understanding of the structure and dynamics of compact objects in these theories. We discuss current bounds on modified gravity from binary pulsar and cosmological observations, and we highlight the potential of future gravitational wave measurements to inform us on the behavior of gravity in the strong-field regime.

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 approach of the effective field theory of modified gravity, we derive the equations of motion for linear perturbations in the presence of a barotropic perfect fluid on the flat isotropic cosmological background. In a simple version of Gleyzes-Langlois-Piazza-Vernizzi (GLPV) theories, which is the minimum extension of Horndeski theories, we show that a slight deviation of the tensor propagation speed squared $c_{\rm t}^2$ from 1 generally leads to the large modification to the propagation speed squared $c_{\rm s}^2$ of a scalar degree of freedom $\phi$. This problem persists whenever the kinetic energy $\rho_X$ of the field $\phi$ is much smaller than the background energy density $\rho_m$, which is the case for most of dark energy models in the asymptotic past. Since the scaling solution characterized by the constant ratio $\rho_X/\rho_m$ is one way out for avoiding such a problem, we study the evolution of perturbations for a scaling dark energy model in the framework of GLPV theories in the Jordan frame. Provided the oscillating mode of scalar perturbations is fine-tuned so that it is initially suppressed, the anisotropic parameter $\eta=-\Phi/\Psi$ between the two gravitational potentials $\Psi$ and $\Phi$ significantly deviates from 1 for $c_{\rm t}^2$ away from 1. For other general initial conditions, the deviation of $c_{\rm t}^2$ from 1 gives rise to the large oscillation of $\Psi$ with the frequency related to $c_{\rm s}^2$. In both cases, the model can leave distinct imprints for the observations of CMB and weak lensing.

In the approach of the effective field theory of modified gravity, we derive the second-order action and the equation of motion for tensor perturbations on the flat isotropic cosmological background. This analysis accommodates a wide range of gravitational theories including Horndeski theories, its generalization, and the theories with spatial derivatives higher than second order (e.g., Horava-Lifshitz gravity). We obtain the inflationary power spectrum of tensor modes by taking into account corrections induced by higher-order spatial derivatives and slow-roll corrections to the de Sitter background. We also show that the leading-order spectrum in concrete modified gravitational theories can be mapped on to that in General Relativity under a disformal transformation. Our general formula will be useful to constrain inflationary models from the future precise measurement of the B-mode polarization in the cosmic microwave background.

We consider a class of non-local gravity theories where the Lagrangian is a function of powers of the inverse d'Alembertian operator acting on the Ricci scalar. We take an approach in which the non-local Lagrangian is made local by introducing auxiliary scalar fields, and study the degrees of freedom of the localized Lagrangian. We find that among the auxiliary scalar fields introduced, some of them are always ghost-like. That is, in the Einstein frame they develop a negative kinetic term. Because of this, except for a particular case already known in the literature, in general, it is not clear how to quantize these models and how to interpret this theory in the light of standard field theory.

The general k-essence Lagrangian for the existence of cosmological scaling solutions is derived in the presence of multiple scalar fields coupled to a barotropic perfect fluid. In addition to the scaling fixed point associated with the dynamics during the radiation and matter eras, we also obtain a scalar-field dominated solution relevant to dark energy and discuss the stability of them in the two-field set-up. We apply our general results to a model of two canonical fields with coupled exponential potentials arising in string theory. Depending on model parameters and initial conditions, we show that the scaling matter-dominated epochs followed by an attractor with cosmic acceleration can be realized with/without the couplings to scalar fields. The different types of scaling solutions can be distinguished from each other by the evolution of the dark energy equation of state from high-redshifts to today.

We study cosmological perturbations in bimetric theory with two fluids each of which is coupled to one of the two metrics. Focusing on a healthy branch of background solutions, we clarify the stability of the cosmological perturbations. For this purpose, we extend the condition for the absence of the so-called Higuchi ghost, and show that the condition is guaranteed to be satisfied on the healthy branch. We also calculate the squared propagation speeds of perturbations and derive the conditions for the absence of the gradient instability. To avoid the gradient instability, we find that the model parameters are weakly constrained.

Recently the first example of a unitary theory of Lorentz-invariant massive gravity allowing for stable self-accelerating de Sitter solutions was found, extending the quasidilaton theory. In this paper we further generalize this new action for the quasidilaton field by introducing general Lagrangian terms which are consistent with the quasidilaton symmetry while leading to second order equations of motion. We find that the structure of the theory, compared to the simplest stable example, does not change on introducing these new terms.

We present a review of cosmological solutions in non-linear massive gravity, focusing on the stability of perturbations. Although homogeneous and isotropic solutions have been found, these are now known to suffer from either Higuchi ghost or a new non-linear ghost instability. We discuss two approaches to alleviate this issue. By relaxing the symmetry of the background by e.g. breaking isotropy in the hidden sector, it is possible to accommodate a stable cosmological solution. Alternatively, extending the theory to allow for new dynamical degrees of freedom can also remove the conditions which lead to the instability. As examples for this case, we study the stability of self-accelerating solutions in the quasi-dilatonic extension and generic cosmological solutions in the varying mass extension. While the quasi-dilaton case turns out to be unstable, the varying mass case allows stable regimes of parameters. Viable self-accelerating solutions in the varying mass theory yet remain to be found.

We place observational constraints on slow-variation single-field inflationary models by carrying out the cosmological Monte Carlo simulation with the recent data of Planck combined with the WMAP large-angle polarization, baryon acoustic oscillations, and ACT/SPT temperature data. Our analysis covers a wide variety of models with second-order equations of motion-- including potential-driven slow-roll inflation, non-minimally coupled models, running kinetic couplings, Brans-Dicke theories, potential-driven Galileon inflation, field-derivative couplings to the Einstein tensor, and k-inflation. In the presence of running kinetic exponential couplings, covariant Galileon terms, and field-derivative couplings, the tensor-to-scalar ratio of the self-coupling potential gets smaller relative to that in standard slow-roll inflation, but the models lie outside the 68 % CL observational contour. We also show that k-inflation models can be tightly constrained by adding the bounds from the scalar non-Gaussianities. The small-field inflationary models with asymptotic flat Einstein-frame potentials in the regime phi >> M_{pl} generally fit the data very well. These include the models such as Kahler-moduli inflation, non-minimally coupled Higgs inflation, and inflation in Brans-Dicke theories in the presence of the potential V(phi)=3M^2 (phi-M_{pl})^2/4 with the Brans-Dicke parameter omega_{BD}<O(1) (which covers the Starobinsky's model f(R)=R+R^2/(6M^2) as a special case).

We present the first example of a unitary theory of Lorentz-invariant massive gravity, with all degrees of freedom propagating on a strictly homogeneous and isotropic, self-accelerating de Sitter background. The theory is a simple extension of the quasidilaton theory, respecting the symmetry of the original theory but allowing for a new type of coupling between the massive graviton and the quasidilaton scalar.

We investigate nonlinear stability of two classes of cosmological solutions in massive gravity: isotropic Friedmann-Lemaitre-Robertson-Walker (FLRW) solutions and anisotropic FLRW solutions. For this purpose we construct the linear cosmological perturbation theory around axisymmetric Bianchi type--I backgrounds. We then expand the background around the two classes of solutions, which are fixed points of the background evolution equation, and analyze linear perturbations on top of it. This provides a consistent truncation of nonlinear perturbations around these fixed point solutions and allows us to analyze nonlinear stability in a simple way. In particular, it is shown that isotropic FLRW solutions exhibit nonlinear ghost instability. On the other hand, anisotropic FLRW solutions are shown to be ghost-free for a range of parameters and initial conditions.

We discuss graviton oscillations based on the ghost free bi-gravity theory. We point out that this theory possesses a natural cosmological background solution which is very close to the case of general relativity. Furthermore, interesting parameter range of the graviton mass, which can be explored by the observations of gravitational waves, is not at all excluded by the constraint from the solar system tests. Therefore the graviton oscillation with possible inverse chirp signal would be an interesting scientific target of KAGRA, adv LIGO, adv Virgo and GEO.

For two types of quintessence models having thawing and tracking properties, there exist analytic solutions for the dark energy equation of state w expressed in terms of several free parameters. We put observational bounds on the parameters in such scenarios by using the recent data of Supernovae type Ia, Cosmic Microwave Background, and Baryon Acoustic Oscillations. The observational constraints are quite different depending on whether or not the recent BAO data from BOSS are taken into account. With the BOSS data the upper bounds of today's values of w (=w_0) in thawing models is very close to -1, whereas without this data the values of w_0 away from -1 can be still allowed. The tracker equation of state w_{(0)} during the matter era is constrained to be w_{(0)}<-0.949 at 95 % confidence level even without the BOSS data, so that the tracker models with w away from -1 are severely disfavored. We also study observational constraints on scaling models in which w starts to evolve from 0 in the deep matter era and show that the transition to the equation of state close to w=-1 needs to occur at an early cosmological epoch. In the three classes of quintessence models studied in this paper, the past evolution of the Hubble parameters in the best-fit models shows only less than the 2.5 % difference compared to the LCDM.

In the Horndeski's most general scalar-tensor theories, we derive the three-point correlation function of scalar non-Gaussianities generated during single-field inflation in the presence of slow-variation corrections to the leading-order term. Unlike previous works, the resulting bispectrum is valid for any shape of non-Gaussianities. In the squeezed limit, for example, this gives rise to the same consistency relation as that derived by Maldacena in standard single-field slow-roll inflation. We estimate the shape close to the squeezed one at which the effect of the term inversely proportional to the scalar propagation speed squared begins to contribute to the bispectrum. We also show that the leading-order bispectrum can be expressed by the linear combination of two convenient bases whose shapes are highly correlated with equilateral and orthogonal types respectively. We present concrete models in which the orthogonal and enfolded shapes can dominate over the equilateral one.

The red-shift space distortions in the galaxy power spectrum can be used to measure the growth rate of matter density perturbations delta_m. For dynamical dark energy models in General Relativity we provide a convenient analytic formula of f(z) sigma_8(z) written as a function of the redshift z, where f=d ln delta_m/d ln a (a is the cosmological scale factor) and sigma_8 is the rms amplitude of over-density at the scale 8 h^{-1} Mpc. Our formula can be applied to the models of imperfect fluids, quintessence, and k-essence, provided that the dark energy equation of state w does not vary significantly and that the sound speed is not much smaller than 1. We also place observational constraints on dark energy models of constant w and tracking quintessence from the recent data of red-shift space distortions.

We argue that all homogeneous and isotropic solutions in nonlinear massive gravity are unstable. For this purpose, we study the propagating modes on a Bianchi type--I manifold. We analyze their kinetic terms and dispersion relations as the background manifold approaches the homogeneous and isotropic limit. We show that in this limit, at least one ghost always exists and that its frequency tends to vanish for large scales, meaning that it cannot be integrated out from the low energy effective theory. This ghost mode is interpreted as a leading nonlinear perturbation around a homogeneous and isotropic background.

We consider the 3-form field, which has been considered as a candidate for realizing inflation, coupled to a scalar field which models the relativistic matter particles produced during the reheating epoch. We have investigated the stability conditions for this theory and found that introducing such a coupling does not lead to any ghosts or Laplacian instabilities. We have also investigated the reheating temperature and the production of particles due to parametric resonances. We have found that this process is more efficient in this theory compared to the result of the standard-scalar-field inflationary scenario.

We consider the Eddington-Born-Infeld (EBI) model here without assuming any cosmological constant. The EBI scalar field is supposed to play a role of both dark matter and dark energy. Different eras in cosmology are reconstructed for the model. A comparison is drawn with $\Lambda$CDM model using Supernova Ia, WMAP7 and BAO data. It seems that the EBI field in this form does not give good fit to observational data in comparison to the $\Lambda$CDM model.

We investigate circumstances under which one can generalize Horndeski's most general scalar-tensor theory of gravity. Specifically we demonstrate that a nonlinear combination of purely kinetic gravity terms can give rise to an accelerating universe without the addition of extra propagating degrees of freedom on cosmological backgrounds, and exhibit self tuning to bring a large cosmological constant under control. This nonlinear approach leads to new properties that may be instructive for exploring the behaviors of gravity.

We consider the minimally coupled 3-form field which has been considered as a candidate to realize inflation. We have studied the conditions to avoid ghosts and Laplacian instabilities and found that some classes of potentials, e.g. the Mexican-hat one, will in general be unstable. We then propose other classes of potentials which are instead free from any instability, drive a long-enough slow-roll regime followed by an oscillatory epoch, and as a consequence, can provide successful inflation. Finally, we also provide stable potentials which lead to a small enough propagation speed for the scalar perturbations, giving a possibility for these models to produce non-Gaussianities.

We place observational constraints on models with the late-time cosmic acceleration based on a number of parametrizations allowing fast transitions for the equation of state of dark energy. In addition to the model of Linder and Huterer where the dark energy equation of state $w$ monotonically grows or decreases in time, we propose two new parametrizations in which $w$ has an extremum. We carry out the likelihood analysis with the three parametrizations by using the observational data of supernovae type Ia, cosmic microwave background, and baryon acoustic oscillations. Although the transient cosmic acceleration models with fast transitions can give rise to the total chi square smaller than that in the $\Lambda$-Cold-Dark-Matter ($\Lambda$CDM) model, these models are not favored over $\Lambda$CDM when one uses the Akaike information criterion which penalizes the extra degrees of freedom present in the parametrizations.

In second-order scalar-tensor theories we study how the Vainshtein mechanism works in a spherically symmetric background with a matter source. In the presence of the field coupling $F(\phi)=e^{-2Q\phi}$ with the Ricci scalar $R$ we generally derive the Vainshtein radius within which the General Relativistic behavior is recovered even for the coupling $Q$ of the order of unity. Our analysis covers the models such as the extended Galileon and Brans-Dicke theories with a dilatonic field self-interaction. We show that, if these models are responsible for the cosmic acceleration today, the corrections to gravitational potentials are generally small enough to be compatible with local gravity constraints.

The extended Galileon models possess tracker solutions with de Sitter attractors along which the dark energy equation of state is constant during the matter-dominated epoch, i.e. w_DE = -1-s, where s is a positive constant. Even with this phantom equation of state there are viable parameter spaces in which the ghosts and Laplacian instabilities are absent. Using the observational data of the supernovae type Ia, the cosmic microwave background (CMB), and baryon acoustic oscillations, we place constraints on the tracker solutions at the background level and find that the parameter s is constrained to be s=0.034 (-0.034,+0.327) (95% CL) in the flat Universe. In order to break the degeneracy between the models we also study the evolution of cosmological density perturbations relevant to the large-scale structure (LSS) and the Integrated-Sachs-Wolfe (ISW) effect in CMB. We show that, depending on the model parameters, the LSS and the ISW effect is either positively or negatively correlated. It is then possible to constrain viable parameter spaces further from the observational data of the ISW-LSS cross-correlation as well as from the matter power spectrum.

In the Horndeski's most general scalar-tensor theories the equations of scalar density perturbations are derived in the presence of non-relativistic matter minimally coupled to gravity. Under a quasi-static approximation on sub-horizon scales we obtain the effective gravitational coupling $G_{eff}$ associated with the growth rate of matter perturbations as well as the effective gravitational potential $\Phi_{eff}$ relevant to the deviation of light rays. We then apply our formulas to a number of modified gravitational models of dark energy--such as those based on f(R) theories, Brans-Dicke theories, kinetic gravity braidings, covariant Galileons, and field derivative couplings with the Einstein tensor. Our results are useful to test the large-distance modification of gravity from the future high-precision observations of large-scale structure, weak lensing, and cosmic microwave background.

In the Horndeski's most general scalar-tensor theories with second-order field equations, we derive the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar, tensor, and vector perturbations in the presence of two perfect fluids on the flat Friedmann-Lemaitre-Robertson-Walker (FLRW) background. Our general results are useful for the construction of theoretically consistent models of dark energy. We apply our formulas to extended Galileon models in which a tracker solution with an equation of state smaller than -1 is present. We clarify the allowed parameter space in which the ghosts and Laplacian instabilities are absent and we numerically confirm that such models are indeed cosmologically viable.

We study the cosmology of a generalized Galileon field $\phi$ with five covariant Lagrangians in which $\phi$ is replaced by general scalar functions $f_{i}(\phi)$ (i=1,...,5). For these theories, the equations of motion remain at second-order in time derivatives. We restrict the functional forms of $f_{i}(\phi)$ from the demand to obtain de Sitter solutions responsible for dark energy. There are two possible choices for power-law functions $f_{i}(\phi)$, depending on whether the coupling $F(\phi)$ with the Ricci scalar $R$ is independent of $\phi$ or depends on $\phi$. The former corresponds to the covariant Galileon theory that respects the Galilean symmetry in the Minkowski space-time. For generalized Galileon theories we derive the conditions for the avoidance of ghosts and Laplacian instabilities associated with scalar and tensor perturbations as well as the condition for the stability of de Sitter solutions. We also carry out detailed analytic and numerical study for the cosmological dynamics in those theories.

For very general scalar-field theories in which the equations of motion are at second-order, we evaluate the three-point correlation function of primordial scalar perturbations generated during inflation. We show that the shape of non-Gaussianities is well approximated by the equilateral type. The equilateral non-linear parameter f_NL^equil is derived on the quasi de Sitter background where the slow-variation parameters are much smaller than unity. We apply our formula for f_NL^equil to a number of single-field models of inflation--such as k-inflation, k-inflation with Galileon terms, potential-driven Galileon inflation, nonminimal coupling models (including field-derivative coupling models), and Gauss-Bonnet gravity.

We study chaotic inflation in the context of modified gravitational theories. Our analysis covers models based on (i) a field coupling $\omega(\phi)$ with the kinetic energy $X$ and a nonmimimal coupling $\zeta \phi^{2} R/2$ with a Ricci scalar $R$, (ii) Brans-Dicke (BD) theories, (iii) Gauss-Bonnet (GB) gravity, and (iv) gravity with a Galileon correction. Dilatonic coupling with the kinetic energy and/or negative nonminimal coupling are shown to lead to compatibility with observations of the Cosmic Microwave Background (CMB) temperature anisotropies for the self-coupling inflaton potential $V(\phi)=\lambda \phi^{4}/4$. BD theory with a quadratic inflaton potential, which covers Starobinsky's $f(R)$ model $f(R)=R+R^{2}/(6M^{2})$ with the BD parameter $\omega_{BD}=0$, gives rise to a smaller tensor-to-scalar ratio for decreasing $\omega_{BD}$. In the presence of a GB term coupled to the field $\phi$, we express the scalar/tensor spectral indices $n_{s}$ and $n_{t}$ as well as the tensor-to-scalar ratio $r$ in terms of two slow-roll parameters and place bounds on the strength of the GB coupling from the joint data analysis of WMAP 7yr combined with other observations. We also study the Galileon-like self-interaction $\Phi(\phi) X \square\phi$ with exponential coupling $\Phi(\phi) \propto e^{\mu\phi}$. Using a CMB likelihood analysis we put bounds on the strength of the Galileon coupling and show that the self coupling potential can in fact be made compatible with observations in the presence of the exponential coupling with $\mu>0$.

We derive the equation of matter density perturbations on sub-horizon scales around a flat Friedmann-Lema\^\i tre-Robertson-Walker background for the general Lagrangian density $f(R,\GB)$ that is a function of a Ricci scalar $R$ and a Gauss-Bonnet term $\GB$. We find that the effective gravitational constant generically scales as distance squared at small distances. The effect of this diminishing of the gravitational constant might be important in the gravitational dynamics of cosmic objects such as galaxies, which can be in principle tested by observations. We also provide the general expressions for the effective anisotropic stress, which is useful to constrain modified gravity models from observations of large-scale structure and weak lensing. We also find that there is a special class of theories which evade this unusual behaviour and that the condition to belong to this special class is exactly the same as the one for not having super-luminal modes with propagation speed proportional to their wavenumber.

We compute the three-point correlation function of primordial scalar density perturbations in a general single-field inflationary scenario, where a scalar field phi has a direct coupling with the Ricci scalar R and the Gauss-Bonnet term GB. Our analysis also covers the models in which the Lagrangian includes a function non-linear in the field kinetic energy X=-(nabla phi)^2/2, and a Galileon-type field self-interaction G(phi, X)*(Box phi), where G is a function of phi and X. We provide a general analytic formula for the equilateral non-Gaussianity parameter f_{NL}^{equil} associated with the bispectrum of curvature perturbations. A quasi de Sitter approximation in terms of slow-variation parameters allows us to derive a simplified form of f_{NL}^{equil} convenient to constrain various inflation models observationally. If the propagation speed of the scalar perturbations is much smaller than the speed of light, the Gauss-Bonnet term as well as the Galileon-type field self-interaction can give rise to large non-Gaussianities testable in future observations. We also show that, in Brans-Dicke theory with a field potential (including f(R) gravity), f_{NL}^{equil} is of the order of slow-roll parameters as in standard inflation driven by a minimally coupled scalar field.

We study the evolution of matter density perturbations in Galileon cosmology where the late-time cosmic acceleration can be realized by a field kinetic energy. We obtain full perturbation equations at linear order in the presence of five covariant Lagrangians ${cal L}_i$ ($i=1,...,5$) satisfying a Galilean symmetry in the flat space-time. The equations for a matter perturbation as well as an effective gravitational potential are derived under a quasi-static approximation on sub-horizon scales. This approximation can reproduce full numerical solutions with high accuracy for the wavelengths relevant to large-scale structures. For the model parameters constrained by the background expansion history of the Universe the growth rate of matter perturbations is larger than that in the LCDM model, with the growth index $gamma$ today typically smaller than 0.4. We also find that, even on very large scales associated with the Integrated-Sachs-Wolfe (ISW) effect in Cosmic Microwave Background (CMB) temperature anisotropies, the effective gravitational potential exhibits a temporal growth during the transition from the matter era to the epoch of cosmic acceleration. These properties are useful to distinguish the Galileon model from the LCDM in future high-precision observations.

We study linear metric perturbations around a spherically symmetric static spacetime for general f(R,G) theories, where R is the Ricci scalar and G is the Gauss-Bonnet term. We find that unless the determinant of the Hessian of f(R,G) is zero, even-type perturbations have a ghost for any multi-pole mode. In order for these theories to be plausible alternatives to General Relativity, the theory should satisfy the condition that the ghost is massive enough to effectively decouple from the other fields. We study the requirement on the form of f(R,G) which satisfies this condition. We also classify the number of propagating modes both for the odd-type and the even-type perturbations and derive the propagation speeds for each mode.

We study the cosmology of a covariant Galileon field with five covariant Lagrangians and confront this theory with the most recent cosmological probes: the type Ia supernovae data (Constitution and Union2 sets), cosmic microwave background (WMAP7) and the baryon acoustic oscillations (SDSS7). In the Galileon cosmology with a late-time de Sitter attractor, there is a tracker that attracts solutions with different initial conditions to a common trajectory. Including the cosmic curvature K, we place observational constraints on two distinct cases: (i) the tracker, and (ii) the generic solutions to the equations of motion. We find that the tracker solution can be consistent with the individual observational data, but it is disfavored by the combined data analysis. The generic solutions fare quite well when a non-zero curvature parameter is taken into account, but the Akaike and Bayesian information criteria show that they are not particularly favored over the LCDM model.

We study the cosmology of a covariant scalar field respecting a Galilean symmetry in flat space-time. We show the existence of a tracker solution that finally approaches a de Sitter fixed point responsible for cosmic acceleration today. The viable region of model parameters is clarified by deriving conditions under which ghosts and Laplacian instabilities of scalar and tensor perturbations are absent. The field equation of state exhibits a peculiar phantom-like behavior along the tracker, which allows a possibility to observationally distinguish the Galileon gravity from the Lambda-CDM model.

In Brans-Dicke theory a non-linear self interaction of a scalar field allows a possibility of realizing the late-time cosmic acceleration, while recovering the General Relativistic behavior at early cosmological epochs. We extend this to more general modified gravitational theories in which a de Sitter solution for dark energy exists without using a field potential. We derive a condition for the stability of the de Sitter point and study the background cosmological dynamics of such theories. We also restrict the allowed region of model parameters from the demand for the avoidance of ghosts and instabilities. A peculiar evolution of the field propagation speed allows us to distinguish those theories from the LCDM model.

The study of linear perturbation theory for general functions of the Ricci and Gauss-Bonnet scalars is done over an empty anisotropic universe, i.e. the Kasner-type background, in order to show that an anisotropic background in general has ghost degrees of freedom, which are absent on Friedmann-Lemaitre-Robertson-Walker (FLRW) backgrounds. The study of the scalar perturbation reveals that on this background the number of independent propagating degrees of freedom is four and reduces to three on FLRW backgrounds, as one mode becomes highly massive to decouple from the physical spectrum. When this mode remains physical, there is inevitably a ghost mode.

Over the past decade, f(R) theories have been extensively studied as one of the simplest modifications to General Relativity. In this article we review various applications of f(R) theories to cosmology and gravity - such as inflation, dark energy, local gravity constraints, cosmological perturbations, and spherically symmetric solutions in weak and strong gravitational backgrounds. We present a number of ways to distinguish those theories from General Relativity observationally and experimentally. We also discuss the extension to other modified gravity theories such as Brans-Dicke theory and Gauss-Bonnet gravity, and address models that can satisfy both cosmological and local gravity constraints.

We derive the equations of linear cosmological perturbations for the general Lagrangian density $f (R,\phi, X)/2+L_c$, where $R$ is a Ricci scalar, $\phi$ is a scalar field, and $X=-(\nabla \phi)^2/2$ is a field kinetic energy. We take into account a nonlinear self-interaction term $L_c$ recently studied in the context of "Galileon" cosmology, which keeps the field equations at second order. Taking into account a scalar-field mass explicitly, the equations of matter density perturbations and gravitational potentials are obtained under a quasi-static approximation on sub-horizon scales. We also derive conditions for the avoidance of ghosts and Laplacian instabilities associated with propagation speeds. Our analysis includes most of modified gravity models of dark energy proposed in literature and thus it is convenient to test the viability of such models from both theoretical and observational points of view.

In order to classify modified gravity models according to their physical properties, we analyze the cosmological linear perturbations for f(R,G) theories (R being the Ricci scalar and G, the Gauss-Bonnet term) with a minimally coupled perfect fluid. For the scalar type perturbations, we identify in general six degrees of freedom. We find that two of these physical modes obey the same dispersion relation as the one for a non-relativistic de Broglie wave. This means that spacetime is either highly unstable or its fluctuations undergo a scale-dependent super-luminal propagation. Two other modes correspond to the degrees of freedom of the perfect fluid, and propagate with the sound speed of such a fluid. The remaining two modes correspond to the entropy and temperature perturbations of the perfect fluid, and completely decouple from the other modes for a barotropic equation of state. We then provide a concise condition on f(R,G) theories, that both f(R) and R+f(G) do fulfill, to avoid the de Broglie type dispersion relation. For the vector type perturbation, we find that the perturbations decay in time. For the tensor type perturbation, the perturbations can be either super-luminal or sub-luminal, depending on the model. No-ghost conditions are also obtained for each type of perturbation.

We describe the linear cosmological perturbations of a perfect fluid at the level of an action, providing thus an alternative to the standard approach based only on the equations of motion. This action is suited not only to perfect fluids with a barotropic equation of state, but also to those for which the pressure depends on two thermodynamical variables. By quantizing the system we find that (1) some perturbation fields exhibit a noncommutativity quite analogous to the one observed for a charged particle moving in a strong magnetic field, (2) local curvature and pressure perturbations cannot be measured simultaneously, (3) ghosts appear if the null energy condition is violated.

We study the evolution of cosmological perturbations in f(G) gravity, where the Lagrangian is the sum of a Ricci scalar R and an arbitrary function f in terms of a Gauss-Bonnet term G. We derive the equations for perturbations assuming matter to be described by a perfect fluid with a constant equation of state w. We show that density perturbations in perfect fluids exhibit negative instabilities during both the radiation and the matter domination, irrespective of the form of f(G). This growth of perturbations gets stronger on smaller scales, which is difficult to be compatible with the observed galaxy spectrum unless the deviation from General Relativity is very small. Thus f(G) cosmological models are effectively ruled out from this Ultra-Violet instability, even though they can be compatible with the late-time cosmic acceleration and local gravity constraints.

We study the propagation of the scalar modes around a Friedmann-Lemaitre-Robertson-Walker universe for general modifications of gravity in the presence of a real scalar field. In general, there will be two propagating scalar perturbation fields, which will have in total four degrees of freedom. Two of these degrees will have a superluminal propagation--with k-dependent speed of propagation--whereas the other two will travel with the speed of light. Therefore, the scalar degrees of freedom do not modify the general feature of modified gravity models: the appearance of modes whose frequency depends on the second power of the modulus of the wave vector. Constraints are given and special cases are discussed.

We discuss solar system constraints on f(G) gravity models, where f is a function of the Gauss-Bonnet term G. We focus on cosmologically viable f(G) models that can be responsible for late-time cosmic acceleration. These models generally give rise to corrections of the form epsilon*(r/rs)^p to the vacuum Schwarzschild solution, where epsilon = H^2 rs^2 << 1, rs is the Schwarzschild radius of Sun, and H is the Hubble parameter today. We generally estimate the strength of modifications to General Relativity in order to confront models with a number of experiments such as the deflection of light and the perihelion shift. We show that cosmologically viable f(G) models can be consistent with solar system constraints for a wide range of model parameters.

We have found for the general class of Modified Gravity Models f(R,G) a new instability which can arise in vacuum for the scalar modes of the cosmological perturbations if the background is not de Sitter. In particular, the short-wavelength modes, if stable, in general have a group velocity which depends linearly in k, the wave number. Therefore these modes will be in general superluminal. We have also discussed the condition for which in general these scalar modes will be ghost-like. There is a subclass of these models, defined out of properties of the function f(R,G) and to which the f(R) and f(G) models belong, which however does not have this feature.

We derive exact solutions of the seven-dimensional Einstein-Maxwell equations for a spacetime exhibiting Poincare invariance along four-dimensions and spherical symmetry in the extra-dimensions. Such topology generically arises in the context of braneworld models. Our solutions generalise previous results on Ricci-flat spacetimes admitting the two-sphere and are shown to include wormhole configurations. A regular coordinate system suitable to describe the whole spacetime is singled-out and we discuss the physical relevance of the derived solutions.

We derive conditions under which f(G) gravity models, whose Lagrangian densities f are written in terms of a Gauss-Bonnet term G, are cosmologically viable. The most crucial condition to be satisfied is that f_GG, the second derivative of f with respect to G, must be positive, which is required to ensure the stability of a late-time de-Sitter solution as well as the existence of standard radiation/matter dominated epochs. We present a number of explicit f(G) models in which a cosmic acceleration is followed by the matter era. We find that the equation of state of dark energy can cross the phantom divide before reaching the present Universe. The viable models have asymptotic behavior f_GG goes to +0 when |G| goes to infinity, in which case a rapid oscillation of perturbations occurs unless such an oscillating degree of freedom is suppressed relative to a homogeneous mode in the early universe. We also introduce an iterative method to avoid numerical instabilities associated with a large mass of the oscillating mode.

We propose a regular classical field theory realisation of the Dvali-Gabadadze-Porrati mechanism by considering our universe to be the four-dimensional core of a seven dimensional 't Hooft-Polyakov hypermonopole. We show the existence of metastable gravitons trapped in the core. Their mass spectrum is discrete, positive definite, and computed for various values of the field coupling constants: the resulting Newton gravity law is seven-dimensional at small and large distances but can be made four-dimensional on intermediate length scales. There is no need of a cosmological constant in the bulk, the spacetime is asymptotically flat and of infinite volume in the extra-dimensions. Confinement is achieved through the local positive curvature of the extra-dimensions induced by the monopole-forming fields and for natural values of the coupling constants of order unity.

A class of Modified Gravity Models, consisting of inverse powers of linear combination of quadratic curvature invariants, is studied in the full parameter space. We find that singularity-free cosmological solutions, interpolating between an almost-Friedmann universe at Big Bang Nucleosynthesis and an accelerating universe today, exist only in a restricted parameter space. Furthermore, for all parameters of the models, there is an unstable scalar mode of the gravitational field. Therefore we conclude that this class of Modified Gravity Models is not viable.

Modified gravity provides a possible explanation for the currently observed cosmic accelaration. In this paper, we study general classes of modified gravity models. The Einstein-Hilbert action is modified by using general functions of the Ricci and the Gauss-Bonnet scalars, both in the metric and in the Palatini formalisms. We do not use an explicit form for the functions, but a general form with a valid Taylor expansion up to second order about redshift zero in the Riemann-scalars. The coefficients of this expansion are then reconstructed via the cosmic expansion history measured using current cosmological observations. These are the quantities of interest for theoretical considerations relating to ghosts and instabilities. We find that current data provide interesting constraints on the coefficients. The next-generation dark energy surveys should shrink the allowed parameter space for modifed gravity models quite dramatically.

We reconsider constraints on Brans-Dicke theories arising from the requirement of successful Big Bang Nucleosynthesis. Such constraints typically arise by imposing that the universe be radiation-dominated at early times, and therefore restricting the contribution that a Brans-Dicke scalar could make to the energy budget of the universe. However, in this paper we show how the dynamics of the Brans-Dicke scalar itself can mimic a portion of the radiation contribution, thereby allowing successful nucleosynthesis with a sizable contribution to the total cosmic energy density. This possibility significantly relaxes the existing bounds on Brans-Dicke fields, and opens the door to new possibilities for early universe cosmology. The necessary fine tunings required by such an arrangement are identified and discussed.

We consider Modified Gravity models involving inverse powers of fourth-order curvature invariants. Using these models' equivalence to the theory of a scalar field coupled to a linear combination of the invariants, we investigate the properties of the propagating modes. Even in the case for which the fourth derivative terms in the field equations vanish, we find that the second derivative terms can give rise to ghosts, instabilities, and superluminal propagation speeds. We establish the conditions which the theories must satisfy in order to avoid these problems in Friedmann backgrounds, and show that the late-time attractor solutions are generically afflicted by superluminally propagating tensor or scalar modes.

We investigate the stability against inhomogeneous perturbations and the appearance of ghost modes in Gauss-Bonnet gravitational theories with a non-minimally coupled scalar field, which can be regarded as either the dilaton or a compactification modulus in the context of string theory. Through cosmological linear perturbations we extract four no-ghost and two sub-luminal constraint equations, written in terms of background quantities, which must be satisfied for consistency. We also argue that, for a general action with quadratic Riemann invariants, homogeneous and inhomogeneous perturbations are, in general, inequivalent, and that attractors in the phase space can have ghosts. These results are then generalized to a two-field configuration. Single-field models as candidates for dark energy are explored numerically and severe bounds on the parameter space of initial conditions are placed. A number of cases proposed in the literature are tested and most of them are found to be unstable or observationally unviable.

We consider general curvature-invariant modifications of the Einstein-Hilbert action that become important only in regions of extremely low space-time curvature. We investigate the far future evolution of the universe in such models, examining the possibilities for cosmic acceleration and other ultimate destinies. The models generically possess de Sitter space as an unstable solution and exhibit an interesting set of attractor solutions which, in some cases, provide alternatives to dark energy models.

Dark energy candidates for which the equation-of-state parameter w is less than -1 violate the dominant energy condition, and are typically unstable. In scalar-tensor theories of gravity, however, the expansion of the universe can mimic the behavior of general relativity with w<-1 dark energy, without violating any energy conditions. We examine whether this possibility is phenomenologically viable by studying Brans-Dicke models and characterizing both the naturalness of the models themselves, and additional observational constraints from limits on the time-dependence of Newton's constant. We find that only highly contrived models would lead observers to measure w<-1.

The simplest explanation for early time acceleration (inflation) and the late time acceleration indicated by recent data is that they have a common origin. We investigate another generic cosmological implication of this possiblity, that the baryon asymmetry of the universe may be generated in such models. We identify several novel features of baryogenesis in such a universe, in which a rolling scalar field is always part of the cosmological energy budget. We also propose a concrete mechanism by which the baryon asymmetry of the universe may be generated in this context. We analyze the generic properties of and constraints on these cosmologies, and then demonstrate explicitly how a complete cosmology may develop in some specific classes of models.