From inflation to reheating and their dynamical stability analysis in Gauss-Bonnet gravity

We investigate the inflation and reheating phenomenology in scalar-Einstein-Gauss-Bonnet theory of gravity where a scalar field non-minimally couples with the Gauss-Bonnet (GB) curvature term. Regarding the inflationary phenomenology, we find -- (1) the inflation starts with a quasi de-Sitter phase and has an exit at a finite e-fold, (2) the scalar and tensor perturbations prove to be ghost free and do not suffer from gradient instability, (3) the curvature perturbation amplitude as well as its tilt and the tensor-to-scalar ratio turn out to be simultaneously compatible with the recent Planck data for suitable values of the parameters. After the inflation ends, the scalar field starts to decay to radiation with a constant decay width. For our considered scalar potential and the GB coupling function, the model results to an analytic power law solution of the Hubble parameter and a logarithmic solution of the scalar field during the reheating era, where the exponent of the Hubble parameter determines the effective EoS parameter ($w_\mathrm{eff}$) during the same. The stability of such reheating dynamics is examined by dynamical analysis which ensures that $w_\mathrm{eff}$ can go beyond unity and reach up-to the maximum value of $\mathrm{max}(w_\mathrm{eff}) = 1.56$. The scenario with $w_\mathrm{eff} > 1$ proves to be purely due to the presence of the GB coupling function, which in turn may have important consequences on enhancing the primordial gravitational waves' amplitude observed today. The inflationary e-fold number gets further constrained by the input of the reheating stage. We finally construct the complete forms of scalar potential ($V(\phi)$) and the GB coupling ($\xi(\phi)$) function that smoothly transits from inflation to reheating, and numerically solve the Hubble parameter and the scalar field for such complete forms of $V(\phi)$ and $\xi(\phi)$.