Late time cosmology in f(R,G)- gravity with interacting fluids

Cosmological models are obtained in a $f(R)$ modified gravity with a coupled Gauss-Bonnet (GB) terms in the gravitational action. The dynamical role of the GB terms is explored with a coupled dilaton field in two different cases (I) $f(R)= R + \gamma R^2- \lambda \left( \frac{R}{3m_s^2} \right)^{\delta}$ where $\gamma$, $\lambda$ and $\delta$ are arbitrary constants and (II) $f(R)=R$ and estimate the constraints on the model parameters. In the first case we choose GB terms coupled with a free scalar field in the presence of interacting fluid and in the second case GB terms coupled with scalar field in a self interacting potential to compare the observed universe. The evolutionary scenario of the universe is obtained adopting a numerical technique as the field equations are highly non-linear. Defining a new density parameter $\Omega_{H}$, a ratio of the dark energy density to the present energy density of the non-relativistic matter, we look for a late accelerating universe. The state finder parameters $\Omega_{H}$, deceleration parameter ($q$), jerk parameter ($j$) are plotted. It is noted that a non-singular universe with oscillating cosmological parameters for a given strength of interactions is admitted in Model-I. The gravitational coupling constant $\lambda$ is playing an important role. The Lagrangian density of $f(R)$ is found to dominate over the GB terms when oscillating phase of dark energy arises. In Model-II, we do not find oscillation of the cosmological parameters as the universe evolves. In the presence of interaction the energy from radiation sector of matter cannot flow to the other two sectors of fluid. The range of values of the strengths of interaction of the fluids are estimated for a stable universe assuming the primordial gravitational wave speed equal to unity.