Dynamical systems approach to stellar modelling in $f(G, B)$ gravity

TL;DR

This paper introduces a dynamical systems approach to stellar modeling within $f(G, B)$ gravity, analyzing stability and geometric properties of solutions with novel boundary term splitting.

Contribution

It proposes a new method using dynamical systems to analyze stellar structures in $f(G, B)$ gravity, emphasizing boundary terms and stability of solutions.

Findings

Invariant submanifolds are generally stable.

The autonomous isotropy equation allows phase portrait analysis.

Two possible vacuum metrics emerge from the theory.

Abstract

The novel proposal to invoke the split of the Ricci scalar into bulk and boundary terms in the gravitational action, opens up a new avenue of investigation into stellar dynamics. The Lagrangian contains functional forms of the bulk term while the boundary term do not contribute to the dynamics. The advantage of the proposition is that the stellar structure equations are up to order two, thus the theory is not haunted by ghosts. We obtain explicitly the defining equations for the thermodynamical variables and the geometry for the pure quadratic case, since the linear case amounts to general relativity. In trying to establish the vacuum geometry associated with the theory it turns out that two possible metrics emerge through the vanishing of the energy-momentum tensor. Next, we analyse the isotropy equation and make the observation that it is autonomous. It is rare that this happens in astrophysical modelling. This behaviour prompted the use of dynamical systems to understand the stability properties of fixed points or invariant submanifolds. It was necessary to choose a gauge in order to split the autonomous equation into a system from which we could plot a phase portrait and deduce the stability of solution trajectories. We find that the invariant submanifolds were generally stable with nearby paths approaching them.