The recent observation of the the gravitational wave event GW170817 and of its electromagnetic counterpart GRB170817A, from a binary neutron star merger, has established that the speed of gravitational waves deviates from the speed of light by less than one part in $10^{15}$. As a consequence, many extensions of General Relativity are inevitably ruled out. Among these we find the most relevant sectors of Horndeski gravity. In its original formulation, mimetic gravity is able to mimic cosmological dark matter, has tensorial perturbations that travel exactly at the speed of light but has vanishing scalar perturbations and this fact persists if we combine mimetic with Horndeski gravity. In this work, we show that implementing the mimetic gravity action with higher-order terms that break the Horndeski structure yields a cosmological model that satisfies the constraint on the speed of gravitational waves and mimics both dark energy and dark matter with a non-vanishing speed of sound. In this way, we are able to reproduce the $\Lambda$CDM cosmological model without introducing particle cold dark matter.

We consider a special class of vacuum $F(R)$-modified gravity models. The form of their Lagrangian is such that the field equations are trivially satisfied when the Ricci scalar is constant. There are many interesting $F(R)$-models for inflation and dark energy that fall in this class. However, little is known outside the domain of cosmology therefore we aim to explore the class of solutions that are static and spherically symmetric. After some general considerations, we investigate in more detail black hole solutions, traversable wormhole metrics and, finally, configurations that can match the anomalous rotation curves of galaxies.

In this paper we investigate how to describe in a unified way a constant-roll inflationary era with a dark energy era, by using the theoretical framework of $F(R)$ gravity. To this end, we introduce some classes of appropriately chosen $F(R)$ gravity models, and we examine in detail how the unification of early and late-time acceleration eras can be achieved. We study in detail the inflationary era, and as we demonstrate it is possible to achieve a viable inflationary era, for which the spectral index of primordial curvature perturbations and the scalar-to-tensor ratio can be compatible with the latest observational data. In addition, the graceful exit issue is briefly discussed for a class of models. Finally, we discuss the dark energy oscillations issue, and we investigate which model from one of the classes we introduced, can produce oscillations with the smallest amplitude.

We consider higher-derivative quantum gravity where renormalization group improved effective action beyond one-loop approximation is derived. Using this effective action, the quantum-corrected FRW equations are analyzed. De Sitter universe solution is found. It is demonstrated that such de Sitter inflationary universe is instable. The slow-roll inflationary parameters are calculated. The contribution of renormalization group improved Gauss-Bonnet term to quantum-corrected FRW equations as well as to instability of de Sitter universe is estimated. It is demonstrated that in this case the spectral index and tensor-to-scalar ratio are consistent with Planck data.

We study inflation induced by (power-low) scalar curvature corrections to General Relativity. The class of inflationary scalar potentials $V(\sigma)\sim\exp[n\,\sigma]$, $n$ general parameter, is investigated in the Einsein frame and the corresponding actions in the Jordan frame are derived. We found the conditions for which these potentials are able to reproduce viable inflation according with the last cosmological data and lead to large scalar curvature corrections which emerge only at a mass scale larger than the Planck mass. Cosmological constant may appear or be set equal to zero in the Jordan frame action without changing the behaviour of the model during inflation. Moreover, polynomial corrections to General Relativity are analyzed in detail. When de Sitter space-time emerges as an exact solution of the models, it is necessary to use perturbative equations in the Jordan framework to study their dynamics during the inflation. In this case, we demonstrate that the Ricci scalar decreases after a correct amount of inflation, making the models consistent with the observable evolution of the universe.

We explore conformal-anomaly driven inflation in $F(R)$ gravity without invoking the scalar-tensor representation. We derive the stress-energy tensor of the quantum anomaly in the flat homogeneous and isotropic universe. We investigate a suitable toy model of exponential gravity plus the quantum contribution due to the conformal anomaly, which leads to the de Sitter solution. It is shown that in $F(R)$ gravity model, the curvature perturbations with its enough amplitude consistent with the observations are generated during inflation. We also evaluate the number of $e$-folds at the inflationary stage and the spectral index $n_\mathrm{s}$ of scalar modes of the curvature perturbations by analogy with scalar tensor theories, and compare them with the observational data. As a result, it is found that the Ricci scalar decreases during inflation and the standard evolution history of the universe is recovered at the small curvature regime. Furthermore, it is demonstrated that in our model, the tensor-to-scalar ratio of the curvature perturbations can be a finite value within the $68\%\,\mathrm{CL}$ error of the very recent result found by the BICEP2 experiment.

We investigate the future evolution of Nariai black hole which is extremal limit of Schwarzschild-de Sitter one in modified gravity. The perturbations equations around Nariai black hole are derived in static and cosmological patches for general $F(R)$-gravity. The analytical and numerical study of several realistic $F(R)$-models shows the occurence of rich variety of scenarios: instabilities, celebrated Hawking evaporation and anti-evaporation of black hole. The realization of specific scenario depends on the model under consideration. It is remarkable that the presence of such primordial black holes at current universe may indicate towards the modified gravity which supports the anti-evaporation as preferrable model. As some generalization we extend the study of Nariai black hole evolution to modified Gauss-Bonnet gravity. The corresponding perturbations equations turn out to be much more complicated than in the case of $F(R)$ gravity. For specific example of modified Gauss-Bonnet gravity we demonstrate that Nariai solution maybe stable.

A generic feature of viable exponential $F(R)$-gravity is investigated. An additional modification to stabilize the effective dark energy oscillations during matter era is proposed and applied to two viable models. An analysis on the future evolution of the universe is performed. Furthermore, a unified model for early and late-time acceleration is proposed and studied.

A generic feature of viable $F(R)$ gravity is investigated: It is demonstrated that during the matter dominated era the large frequency oscillations of the effective dark energy may influence the behavior of higher derivatives of the Hubble parameter with the risk to produce some singular unphysical solutions at high redshift. This behavior is explicitly analyzed for realistic $F(R)$ models, in particular, exponential gravity and a power form model. To stabilize such oscillations, we consider the additional modification of the models via a correction term which does not destroy the viability properties. A detailed analysis on the future evolution of the universe and the evolution history of the growth index of the matter density perturbations are performed. Furthermore, we explore two applications of exponential gravity to the inflationary scenario. We show how it is possible to obtain different numbers of $e$-folds during the early-time acceleration by making different choices of the model parameters in the presence of ultrarelativistic matter, which destabilizes inflation and eventually leads to the exit from the inflationary stage. We execute the numerical analysis of inflation in two viable exponential gravity models. It is proved that at the end of the inflation, the effective energy density and curvature of the universe decrease and thus a unified description between inflation and the $\Lambda$CDM-like dark energy dominated era can be realized.

Oscillations of the $F(R)$ dark energy around the phantom divide line, $\omega_{DE}=-1$, both during the matter era and also in the de Sitter epoch are investigated. The analysis during the de Sitter epoch is revisited by expanding the modified equations of motion around the de Sitter solution. Then, during the matter epoch, the time dependence of the dark energy perturbations is discussed by using two different local expansions. For high values of the red shift, the matter epoch is a stable point of the theory, giving the possibility to expand the $F(R)$-functions in terms of the dark energy perturbations. In the late-time matter era, the realistic case is considered where dark energy tends to a constant. The results obtained are confirmed by precise numerical computation on a specific model of exponential gravity. A novel and very detailed discussion is provided on the critical points in the matter era and on the relation of the oscillations with possible singularities.

A theory of exponential modified gravity which explains both early-time inflation and late-time acceleration, in a unified way, is proposed. The theory successfully passes the local tests and fulfills the cosmological bounds and, remarkably, the corresponding inflationary era is proven to be unstable. Numerical investigation of its late-time evolution leads to the conclusion that the corresponding dark energy epoch is not distinguishable from the one for the $\Lambda$CDM model. Several versions of this exponential gravity, sharing similar properties, are formulated. It is also shown that this theory is non-singular, being protected against the formation of finite-time future singularities. As a result, the corresponding future universe evolution asymptotically tends, in a smooth way, to de Sitter space, which turns out to be the final attractor of the system.

We study all four types of finite-time future singularities emerging in late-time accelerating (effective quintessence/phantom) era from $\mathcal{F}(R,G)$-gravity, where $R$ and $G$ are the Ricci scalar and the Gauss-Bonnet invariant, respectively. As an explicit example of $\mathcal{F}(R,G)$-gravity, we also investigate modified Gauss-Bonnet gravity, so-called $F(G)$-gravity. In particular, we reconstruct the $F(G)$-gravity and $\mathcal{F}(R,G)$-gravity models where accelerating cosmologies realizing the finite-time future singularities emerge. Furthermore, we discuss a possible way to cure the finite-time future singularities in $F(G)$-gravity and $\mathcal{F}(R,G)$-gravity by taking into account higher-order curvature corrections. The example of non-singular realistic modified Gauss-Bonnet gravity is presented. It turns out that adding such non-singular modified gravity to singular Dark Energy makes the combined theory to be non-singular one as well.

A general approach to viable modified $f(R)$ gravity is developed in both the Jordan and the Einstein frames. A class of exponential, realistic modified gravities is introduced and investigated with care. Special focus is made on step-class models, most promising from the phenomenological viewpoint and which provide a natural way to classify all viable modified gravities. One- and two-steps models are explicitly considered, but the analysis is extensible to $N$-step models. Both inflation in the early universe and the onset of recent accelerated expansion arise in these models in a natural, unified way. Moreover, it is demonstrated that models in this category easily pass all local tests, including stability of spherical body solution, non-violation of Newton's law, and generation of a very heavy positive mass for the additional scalar degree of freedom.