We study the inflationary phenomenology of a non-minimally coupled Einstein Gauss-Bonnet gravity theory, in the presence of a scalar potential, under the condition that the gravitational wave speed of the primordial gravitational waves is equal to unity, that is $c_T^2=1$, in natural units. The equations of motion, which are derived directly from the gravitational action, form a system of differential equations with respect to Hubble's parameter and the inflaton field which are very complicated and cannot be solved analytically, even in the minimal coupling case. In this paper, we present a variety of different approximations which could be used, along with the constraint $c_T^2=1$, in order to produce an inflationary phenomenology compatible with recent observations. All the different approaches are able to lead to viable results if the model coupling functions obey simple relations, however, different approaches contain different approximations which must be obeyed during the first horizon crossing, in order for the model to be rendered correct. Models which may lead to a non-viable phenomenology are presented as well in order to understand better the inner framework of this theory. Furthermore, since the velocity of the gravitational waves is set equal to $c_T^2=1$, as stated by the striking event of GW170817 recently, the non-minimal coupling function, the Gauss-Bonnet scalar coupling and the scalar potential are related to each other. Here, we shall assume no particular form of the scalar potential and we choose freely the scalar functions coupled to the Ricci scalar and the Gauss-Bonnet invariant. Certain models are also studied in order to assess the phenomenological validity of the theory, but we need to note that all approximations must hold true in order for a particular model to be valid.

After the recent GW170817 event of the two neutron stars merging, many string corrected cosmological theories confronted the non-viability peril. This was due to the fact that most of these theories produce massive gravitons primordially. Among these theories were the ones containing a non-minimal kinetic coupling correction term in the Lagrangian, which belong to a subclass of Horndeski theories. In this work we demonstrate how these theories may be revived and we show how these theories can produce primordial gravitational waves with speed $c_T^2=1$ in natural units, thus complying with the GW170817 event. As we show, if the gravitational action of an Einstein-Gauss-Bonnet theory also contains a kinetic coupling of the form $\sim \xi(\phi) G^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$, the condition of having primordial massless gravitons, or equivalently primordial gravitational waves with speed $c_T^2=1$ in natural units, results to certain conditions on the scalar field dependent coupling function of the Gauss-Bonnet term, which is also the non-minimal coupling of the kinetic coupling. We extensively study the phenomenological implications of such a theory focusing on the inflationary era, by only assuming slow-roll dynamics for the scalar field. Accordingly, we briefly study the case that the scalar field evolves in a constant-roll way. By using some illustrative examples, we demonstrate that the viability of the theoretical framework at hand may easily be achieved. Also, theories containing terms of the form $\sim \xi(\phi)\Box\phi g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi$ and $\sim \xi(\phi)\left(g^{\mu\nu}\partial_\mu\phi\partial_\nu\phi\right)^2$ also lead to the same gravitational wave speed as the theory we shall study in this paper, so this covers a larger class of Horndeski theories.

In this work we shall study the late-time behavior of $k$-Essence $f(R)$ gravity without scalar potential, in the presence of matter and radiation perfect fluids. We quantify the late-time study by using the statefinder function $Y_H(z)=\frac{\rho_{DE}}{\rho_m^{(0)}}$, which is a function of the redshift and of the Hubble rate. By appropriately rewriting the Friedmann equation in terms of the redshift and of the function $Y_H(z)$, we numerically solve it using appropriate initial conditions, and we critically examine the effects of the $k$-Essence higher order kinetic terms. As we demonstrate, the effect of the higher order scalar field kinetic terms on the late-time dynamics is radical, since the dark energy oscillations are absent, and in addition, the cosmological physical quantities are compatible with the latest Planck data and also the model is almost indistinguishable from the $\Lambda$ Cold Dark Matter model. This is in contrast to the standard $f(R)$ gravity case, where the oscillations are present. Furthermore, by choosing a different set of values of two of the free parameters of the model, and specifically the coefficient of the higher order kinetic term and of the exponent of $R^{\delta}$ appearing in the $f(R)$ action, we demonstrate that it is possible to obtain $\rho_{DE}<0$ for redshifts $z\sim 2-3.8$, which complies phenomenologically with, and seems to explain, the observational data for the same redshifts, and also to obtain a viable cosmological evolution at $z\sim 0$, at least when the dark energy equation of state parameter and the dark energy density parameters are considered.

In this work we introduce a new theoretical framework for Einstein-Gauss-Bonnet theories of gravity, which results to particularly elegant, functionally simple and transparent gravitational equations of motion, slow-roll indices and the corresponding observational indices. The main requirement is that the Einstein-Gauss-Bonnet theory has to be compatible with the GW170817 event, so the gravitational wave speed $c_T^2$ is required to be $c_T^2\simeq 1$ in natural units. This assumption was also made in a previous work of ours, but in this work we express all the related quantities as functions of the scalar field. The constraint $c_T^2\simeq 1$ restricts the functional form of the scalar Gauss-Bonnet coupling function $\xi(\phi)$ and of the scalar potential $V(\phi)$, which must satisfy a differential equation. However, by also assuming that the slow-roll conditions hold true, the resulting equations of motion and the slow-roll indices acquire particularly simple forms, and also the relation that yields the $e$-foldings number is $N=\int_{\phi_i}^{\phi_f}\xi''/\xi'd \phi$, a fact that enables us to perform particularly simple calculations in order to study the inflationary phenomenological implications of several models. As it proves, the models we presented are compatible with the observational data, and also satisfy all the assumptions made during the process of extracting the gravitational equations of motion. More interestingly, we also investigated the phenomenological implications of an additional condition $\xi'/\xi''\ll 1$, which is motivated by the slow-roll conditions that are imposed on the scalar field evolution and on the Hubble rate, in which case the study is easier. Our approach opens a new window in viable Einstein-Gauss-Bonnet theories of gravity.