Cosmological phase-space analysis of $f(G)$-theories of gravity

The impact of topological terms that modify the Hilbert-Einstein action is here explored by virtue of a further $f(G)$ contribution. In particular, we investigate the phase-space stability and critical points of an equivalent scalar field representation that makes use of a massive field, whose potential is function of the topological correction. To do so, we introduce to the gravitational Action Integral a Lagrange multiplier and model the modified Friedmann equations by virtue of new non-dimensional variables. We single out dimensionless variables that permit \emph{a priori} the Hubble rate change of sign, enabling regions in which the Hubble parameter either vanishes or becomes negative. In this respect, we thus analyze the various possibilities associated with a universe characterized by such topological contributions and find the attractors, saddle points and regions of stability, in general. The overall analysis is carried out considering the exponential potential first and then shifting to more complicated cases, where the underlying variables do not simplify. We compute the eigenvalues of our coupled differential equations and, accordingly, the stability of the system, in both a spatially-flat and non-flat universe. Quite remarkably, regardless of the spatial curvature, we show that a stable de Sitter-like phase that can model current time appears only a small fraction of the entire phase-space, suggesting that the model under exam is unlikely in describing the whole universe dynamics, i.e., the topological terms appear disfavored in framing the entire evolution of the universe.