Logarithmic corrected Einstein-Gauss-Bonnet inflation compatible with GW170817

In this paper we investigate the inflationary phenomenology of an Einstein-Gauss-Bonnet theory with the extension of a logarithmic modified $f(R)$ gravity, compatible with the GW170817 event. The main idea of our work is to study different results for an almost linear Ricci scalar through logarithmic corrections and examine whether such model is viable. First of all, the theoretical framework under slow-roll evolution of the scalar field is presented and also developed the formalism of the constant-roll evolution making predictions for the non- Gaussianities of the models is developed , since the constant-roll evolution is known to enhance non-Gaussianities. As shown, the non-Gaussianities are of the order $\mathcal{O}\sim(10^{-1})$. Furthermore, the slow-roll indices and the observational indices of inflation, are calculated for several models of interest. As demonstrated, the phenomenological viability of the models at hand is achieved for a wide range of the free parameters and the logarithmic term has a minor contribution to numerical calculations, as expected.


I. INTRODUCTION
The recent years, one of the most remarkable achievements in theoretical and observational Cosmology is, without a doubt, the detection of the gravitational waves about one hundred years after Einstein's predictions. In 2017, the LIGO-Virgo detectors observed a gravitational wave coming from the merging of two neutron stars [1]. The interesting fact of the observation was that the gravitational wave arrived almost equal with the gamma-ray burst. This means that, the speed of the gravitational wave c T is approximately equal to the speed of light, namely c 2 T = 1 in Natural Units. This constraint imposes stringent conditions on modified theories of gravity that may describe successfully the nature on such scales. Many generalized theories of gravity provide viable models compatible with this astrophysical event, see Ref. [2].
One attractive class of theories, in order to approach the inflationary era of the universe or other astrophysical phenomena is the Einstein-Gauss-Bonnet theories . These theories involve the Gauss-Bonnet term, more specifically, quadratic terms of the Ricci scalar and of the Riemann and Ricci tensors, in the context of Einstein's gravity. Our motivation to utilize this class of gravitational theories in order to describe the inflationary era of the universe originate from the string property. In detail, the whole theory is string-corrected canonical scalar field theory minimally coupled to gravity. In our previous work [36], we proved that this theory can be rectified in view of the GW170817 event, by simply setting the gravitational speed wave equal to unity and expressing all the physical quantities in terms of the scalar field.
In this paper we shall extent our previous works [36,46], presenting a modified theory of gravity f (R) with a logarithmic term, based on Ref. [48], also modified gravity with logarithmic terms is presented in Ref. [49]. First of all, we present in detail the theoretical framework of the evolution of the scalar field under slow-roll and constant-roll approximations in the background of Friedmann-Roberton-Walker spacetime. As mentioned before, we express all the involved physical quantities in terms of functions of the scalar field and their higher derivatives with respect to the scalar field. By imposing the slow-roll conditions into the gravitational equations of motion, the slow-roll indices and the observational indices have quite simple and elegant final expressions. Afterwards, we consider the constant roll evolution of the scalar field, by imposing the conditionφ = βHφ and we repeat the same process. Moreover, the amount of Non-Gaussianities is predicted and as expected is not significantly altered in comparison to the R gravity. Finally, we shall examine the compatibility of the theoretical framework with the latest Plack Data see, Ref [47] in the context of the two possible ways of the scalar field evolution. In the last section, viable models for the inflationary era are constructed by considering coupling scalar field functions ξ(φ) with compatible numerical values with the observations.
Before we begin our analysis, it is necessary to explain the reason why the theory needs massless gravitons in order to be consistent with the GW170817 event. In nature the interactions between particles and fields are achieved through the propagators of the fundamental forces. The graviton is the progagator of the gravitational waves produced either in the early universe or from astrophysical events. From the perspective of the Elementary Particle Physics the graviton must be massless during the inflationary and post-inflationary era. Despite the appealing string property of the Einstein-Gauss-Bonnet theory there is a serious drawback. The theory predicts that the primordial tensor perturbations are incopatible with the GW170817 event. As we have proved in our previous works [32] it is possible to get a massless graviton under certain constraints. Only if the coupling scalar function ξ(φ) satisfies the differential equationξ − Hξ = 0, we can demand c 2 T = 1. Hence, we can obtain a viable model for the inflationary era in the context of Einstein-Gauss-Bonnet gravity.

II. THEORETICAL FRAMEWORK OF LOGARITHMIC CORRECTED f (R) GRAVITY
We begin our analysis by considering the following gravitational theory of a scalar field φ, since all the information about the universe in the era of inflation is encoded in it. Let us assume that the action is defined as, where, g is the determinant of the metric tensor, κ = 1 MP is the gravitational constant while M P denotes the reduced Planck mass, V (φ) is the scalar potential and ξ(φ) signifies the Gauss-Bonnet coupling scalar function. We assume a modified theory of gravity f (R) where, f (R) = R 1+δ with δ ≪ 1. Considering a power-law model, becomes abundantly clear that for small deviations from unity, one could easily obtain logarithmic corrections to gravity via Taylor expansion, such terms are known for describing quantum corrections. Finally, the Gauss-Bonnet term G is given by the expression G = R 2 − 4R αβ R αβ + R αβγδ R αβγδ , with R αβ and R αβγδ being the Ricci and Riemann tensor respectively. Furthermore, the line-element is assumed to have the Friedmann-Robertson-Walker form, where a(t) is the scale factor of the Universe and the metric tensor has the form of g µν = diag(−1, a(t) 2 , a(t) 2 , a(t) 2 ). The effective Lagrangian of inflation is not specified by the data at present time. Thus, although the inflationary era is a classical era of our Universe, which is described by a four dimensional spacetime, it still is possible that the quantum era may have a direct imprint on the effective Lagrangian of inflation. Therefore, the two most simple corrections of the inflationary effective Lagrangian may be provided by higher curvature terms , like f (R) gravity corrections, and Einstein-Gauss-Bonnet corrections. As long as the metric is flat, the Ricci scalar and the Gauss-Bonnet term are topological invariant and can be written as R = 12H 2 + 6Ḣ, G = 24H 2 (Ḣ + H 2 ) respectively. H is Hubble's parameter and in addition, the "dot" denotes differentiation with respect to the cosmic time. We expand the modified gravitational function f (R) as follows, where α is a constant with mass dimensions [m] −2 for consistency. It is expected that the logarithmic term ln(αR) has minor contribution to the equations of motion because represents quantum corrections. Differentiating Eq. (3) with respect to the Ricci scalar gives, Implementing the variation principle with respect to the metric tensor and the scalar field in Eq. (1) generates the field equations of gravity and the continuity equation of the scalar field. By splitting the field equations in time and space components, the gravitational equations of motion are then derived which read, φ + 3Hφ + V ′ + ξ ′ G = 0.
As proved in a recent work of ours [36], certain additional constraints on the gravitational wave speed need to be imposed so as to achieve compatibility with recent striking observations from GW170817. Gravitational waves are perturbations in the metric which travel through spacetime with the speed of light. The gravitational wave speed in natural units for Einstein-Gauss-Bonnet theories has the form, where Q f = 16(ξ − Hξ) and Q t = F κ 2 − 8ξH, are auxiliary functions depending on the scalar field and the Ricci scalar. Compatibility can be achieved by equating the velocity of gravitational waves with unity, or making it infinitesimally close to unity. In other words, we demand Q f = 0. The constraint leads to an ordinary differential equationξ = Hξ. Now we will solve this equation in terms of the derivatives of scalar field. Assuming thatξ = ξ ′φ and d dt =φ d dφ the constraint equation has the form, Considering the approximation Eq. (9) can be solved easily with respect to the derivative of the scalar field, In order to study the inflationary era of the Universe it is necessary to solve analytically the system of equations of motion. It is obvious that, this system is very difficult to study analytically. Thus, we assume the slow-roll approximations during inflation. Mathematically speaking, the following conditions are assumed to hold true, thus, the equations of motion can be simplified greatly. Hence, after imposing the constraint of the gravitational wave and considering the slow-roll approximations the equations of motion have the following elegant forms, 3Hφ However, even with the slow-roll approximations holding true, the system of differential equations still remains intricate and cannot be solved. Further approximations are needed in order to derive the inflationary phenomenology, so we neglect string corrections themselves. This is a reasonable assumption since even though the Gauss-Bonnet scalar coupling function is seemingly neglected, it participates indirectly from the gravitational wave condition. Also, in many cases, string corrections are proven to be subleading. Moreover, under slow-roll assumptions the Ricci scalar is written as, R ≃ 12H 2 . Recalling that f (R) ≃ R + δRln(αR), one obtains elegant simplifications and functional expressions for the equations of motion. The first and the second derivatives of the function F are respectively, H 2 ). The last two terms in the first equation of motion are quite smaller in order of magnitude than the scalar potential of the field hence, as we will prove in the fourth section numerically, these terms can be neglected. The same approximation can be applied in the second equation of motion for the last term of the right hand side. Thus, the final simplified equations of motion are, where the parameter δ can be discarded for |δ| ≪ 1. In the next section we shall prove that the discarded terms are quite smaller in order of magnitude. The dynamics of inflation can be described by six parameters named the slow-roll indices, defined as follows [3,36], where the auxiliary functions are defined as, in detail, the auxiliary functions and the slow-roll indices can be written as follows, where, the index ǫ 4 is omitted due to the perplexed expression. The scalar potential with an unspecified scalar coupling function ξ(φ) can be written as, with V 0 being the amplitude of the scalar potential with mass dimensions [m] 4 . In order to examine the validity of a model, the results which the model produces must be confronted to the recent Planck observational data [47]. In the following model, we shall derive the values for the quantities, namely the spectral index of primordial curvature perturbations n S , the tensor-to-scalar-ratio r and finally, the tensor spectral index n S [3,36]. These quantities are connected with the slow-roll indices introduced previously, as shown below, where c A the sound wave velocity defined as, Based on the latest Planck observational data [47] the spectral index of primordial curvature perturbations is n S = 0.9649 ± 0.0042 and the tensor-to-scalar-ratio r must be r < 0.064. Our goal now is to evaluate the observational indices during the first horizon crossing. However, instead of using wavenumbers, we shall use the values of the scalar potential during the initial stage of inflation. Taking it as an input, we can obtain the actual values of the observational quantities. We can do so by firstly evaluating the final value of the scalar field. This value can be derived by equating slow-roll index ǫ 1 in equation (25) to unity. Consequently, the initial value can be evaluated from the e-foldings number, defined as N = where the difference t f − t i signifies the duration of the inflationary era. Recalling the definition ofφ in Eq. (11), one finds that the proper relation from which the initial value of the scalar field can be derived is, From this equation, as well as equation (25), it is obvious that choosing an appropriate coupling function, is the key in order to simplify the results.

III. CONSTANT-ROLL EVOLUTION OF THE SCALAR FIELD IN LOGARITHMIC CORRECTED f (R) GRAVITY AND PRIMORDIAL NON-GAUSSIANITIES
In the following section we shall analyse the theoretical framework of the evolution of the scalar field under constantroll conditionφ = βHφ, where β is defined as the constant-roll parameter to be specified later. After the specific assumption and neglecting string corrections the equations of motion modified as follows, as a result, the equations of motion are, In the context of the constant-roll evolution of the scalar field, the slow-roll indices for an arbitrary coupling scalar function ξ(φ) are given by the following expressions, where again as before, the index ǫ 4 is omitted due to the lengthy expression. The auxiliary functions is obvious that can be found easily. The e-folding number can be written as, When the scalar field evolves with a constant-rate of roll enhances the non-Gaussianities features. Until now, the perturbations in the Cosmic Microwave Background (CMB) are described perfectly as Gaussian distributions, since no practical evidence is found pointing out a non-Gaussian pattern in the CMB. It is possible though, not evident for the moment, that in the following years the observations may reveal a non-Gaussian pattern in the CMB primordial power spectrum. In this section we shall discuss how to evaluate the non-Gaussianities quantitatively in the context of the GW170817-compatible Einstein-Gauss-Bonnet gravity, using the formalism and notation of [50]. Even though a f (R) logarithmic gravity is assumed, we consider the same equations as the [50]. We first define the following quantities, Here, we shall implement a different formula for the sound wave speed, however equivalent to the previous, which is based on these newly defined quantities for convenience and reads, Recalling equations (34), (36) and (37), one finds that the aforementioned auxiliary terms have the following forms for an unspecified scalar coupling function ξ(φ), In certain examples, we shall demonstrate that by choosing appropriately the coupling function, the quantity ǫ s is simplified greatly, and it shall also coincide with ǫ 1 , along with δ X . No matter the form of the sound wave velocity, the derivative c A is very complex, hence for obvious reasons its analytic expression is omitted. These forms are very useful due to the fact that the power spectra P S of the primordial curvature perturbations and the equilateral momentum approximation term f eq N L can be derived from such terms. These quantities are defined as, f eq N L ≃ 55 36 ǫ s + 5 12 n + 10 3 δ ξ .
In the following we shall appropriately specify the value of the term f eq N L during the first horizon crossing, to see what the constant-roll condition brings along. The evaluation shall be performed by using the values of the free parameters in such a way so that the viability of the observational indices of inflation is achieved according to the 2018 Planck data [47].

IV. TESTING THE THEORETICAL FRAMEWORK WITH THE LATEST OBSERVATIONAL DATA
In this section we shall present explicitly examples of GW170817 compatible Einstein-Gauss-Bonnet models, with the f (R) extension, that can yield a phenomenologically viable inflationary era. First of all, the Gauss-Bonnet coupling scalar function ξ(φ) must be defined properly, aiming for a simple ratio ξ ′ /ξ ′′ so as to facilitate our study.
A. Model compatible with the latest Planck Data under the slow-roll assumption Consider the coupling scalar function ξ(φ) being equal to, where λ 1 , γ 1 and n are dimensionless constants to be specified later while x is an auxiliary integration variable. This particular coupling scalar function is chosen due to the simplicity of the ratio ξ ′ /ξ ′′ , specifically even though the coupling itself seems superfluous. Nevertheless, it has been proven that it is a viable candidate for a Gauss-Bonnet scalar coupling function for the minimally coupled case so it is interesting to examine the impact of logarithmic corrections. One can specify the scalar potential from the equation (18), which has the following form, where V 1 is the integration constant. Let us now proceed with the evaluation of the slow-roll indices, It is obvious that the first three slow-roll indices have quite simple expressions in contrast to the next three due to the derivatives of the scalar potential and the derivatives of the coupling scalar function. Again ǫ 4 is omitted due to the length expression. Now, we can determine the final value of the field when the inflationary era ends, by setting the the first slow-roll index Eq. (56) into unity, The initial value of the scalar field can be calculated from the e-folds number see, equation (33) where we assumed that the number of e-folds is N ≃ 60. Considering the following values for the free parameters in Natural Units, namely, κ 2 = 1, (N, n, γ 1 , V 1 , λ 1 , δ, α) = (60, −2, −5, 1, 1, 0.001, 1) then, the spectral index of primordial curvature perturbations n S , the tensor-to-scalar-ratio r and finally, the tensor spectral index n T are respectively n S = 0.966381, r = 0.000014 and n T = −1.7 × 10 −6 , which are acceptable according to the latest Planck Data. Moreover, we found that the the initial value of the field is φ i = 0.236704 and the final numerical value of the field is φ f = 2.41827 which means that with the passage of time, based on the continuity, the field increases until the inflationary era ends. In addition, the numerical values of the slow-roll indices are ǫ 1 = 8.79 × 10 −7 , ǫ 2 = 0.016, ǫ 3 = ǫ 4 = ǫ 5 = −8.7 × 10 −10 , and ǫ 6 = −1.75 × 10 −9 , where it becomes apparent that all of them are subleading compared to the second slow-roll index.

B. Model with power-law coupling under the constant-roll approximation
In this subsection it is considered a much simpler Gauss-Bonnet coupling scalar function ξ(φ) that follows a powerlaw form, where λ 2 and m are dimensionless constants to be specified later. This is a very appealing function since the ratio ξ ′ ξ ′′ which appears in our calculations is greatly simplified, since, This model was also studied in, Ref. [46]. Since the scalar coupling function is specified, one can find the scalar potential from Eq. (15), where V 2 is an integration constant. For this model the constant-roll indices are given by the following equations, where ǫ 4 is not written analytically due to the perplexed form. Despite the fact that, the indices ǫ 4 to ǫ 6 are quite perplexed, the indices ǫ 1 to ǫ 3 have very simple forms. By setting the index ǫ 1 to unity, the final value of the field when inflationary era ends is equal to, while, the initial value of the field is, which is expected at least for ǫ 5 and ǫ 6 when F is more dominant than string corrections. Moreover, we mention that the initial and final value of the scalar field are φ i ≃ 0.002 and φ f = 10.0707 which indicates an increase in the scalar field. In Fig. 2 we plot the spectral index of primordial curvature perturbations n S (left) and the tensor-to-scalar ratio r (right) depending on parameters β and m ranging [0.001, 0.009] and [8,12] respectively. Moreover, we make predictions for the amount of non-Gaussianities in the primordial power spectrum of the curvature perturbations. The auxiliary terms of the specific model are defined as, From equation (52), the expected value of f eq N L , for the exact same set of parameters we used to obtain the viability of the model with the Planck data, is f eq N L = 0.117024 which is also an accepted value and may explain why non Gaussianities have yet to be observed. Finally, the parameters used to derive such values are δ ξ ≃ 2.1 × 10 −22 , ǫ s = 4.80305 × 10 −8 and η = 0.280857 which means that one of them is in fact dominant. These results imply that ǫ s = ǫ 1 . At this point, it is important to mention that, the observed quantities n S and r experience different changes after the alteration of the free parameters. Specifically, the primordial spectral index is affected only by the constant-roll parameter β while, the tensor-to-scalar ratio, is affected by the exponent m of the coupling scalar function along with the constant-roll parameter with the first being more dominant factor. This can easily be observed in Fig. 2 where one sees that the spectral index of scalar perturbations is depicted by a simple plot resembling vertical lines. In addition, while the term f eq N L is independent of parameter λ 2 , it can be enhanced by decreasing the exponent m but such a change leads to a subsequent decrease in the tensor-to-scalar ratio. For instance, numerically speaking, choosing m = 2 leads to f eq N L = 0.819167, n S = 0.966 and the effective value of the tensor-to-scalar ratio is 0, since, r ∼ O(10 −51 ).
Finally, we examine the validity of the approximations which were made during this model at the first horizon crossing. Firstly, the slow-roll approximations in Eq. (12) hold true sinceḢ ∼ O(10 −8 ) compared to H 2 ∼ O(10 −1 ) and 1 2φ 2 ∼ O(10 −8 ) compared to V ∼ O(10 1 ) are negligible. Also, the terms which were omitted in equations (13) and (14) are of of the order (in reduced Planck units) 24ξH 3 ∼ O(10 −21 ) while 16ξHḢ ∼ O(10 −29 ), which explains why these terms, compared to the scalar potential and the kinetic term, can be neglected and V ′ ∼ O(10 −4 ) whereas ξ ′ G ∼ O(10 −18 ) which explains the neglect of string term. Lastly, the following two terms in equation (13) F R−f

V. CONCLUSIONS
In this work it is presented an alteration of our previous work [36], by using a modified f (R) logarithmic gravity instead of Einstein's gravity, consistent with the GW170817 event. We focused our analysis on the inflationary era of the universe, by considering that the scalar field evolves with either under slow-roll assumptions or with a constant rate of roll. In both of these cases, the slow-roll indices and the observational quantities of inflation were evaluated along with the predicted amount of Non-Gaussianities, in the case of constant-roll evolution of the scalar field. After the theoretical framework we confronted two models with the observational data coming from the Planck 2018 collaboration, considering logarithmic corrections. As shown, the resulting inflationary phenomenology can be compatible with the latest Planck data, for a wide range of the free parameters of the theory, after neglecting string corrections, which is shown to be a reasonable approach. In our analysis we demonstrated that all the assumptions made were satisfied for all the models examined, and for the values of the free parameters that yield inflationary viability with respect to the latest Planck data. In the constant-roll case we also investigated the amount of non-Gaussianities that are predicted from the model, by calculating the nonlinear term f eq N L in the equilateral momentum approximation. Interestingly enough, we demonstrated that the amount of non-Gaussianities is quite small. Finally, we performed an analytic approximation in the differential equation that connects the scalar field potential and the scalar coupling function, and we examined the phenomenology of inflation in this case too. As evinced, the models can also be compatible with the Planck 2018 too, even with a logarithmic modified gravity f (R) with an almost linear Ricci scalar through logarithmic corrections.