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

This paper investigates stellar dynamics within $f(G, B)$ gravity, demonstrating that the theory yields second-order field equations free of ghosts and utilizing dynamical systems analysis to reveal that the autonomous isotropy equation possesses generally stable invariant submanifolds.

The Gist

Imagine the universe as a giant, complex machine. For over a century, our best manual for how this machine works has been General Relativity, written by Albert Einstein. It explains how gravity bends space and time, much like a heavy bowling ball warping a trampoline. This theory works perfectly for most things, but it hits a wall when we try to explain the very beginning of the universe or the very end of a star's life. It also struggles to explain why the universe is expanding faster and faster without inventing invisible "dark energy."

To fix these gaps, scientists are trying to write "new manuals" or modified gravity theories. This paper is about one specific new manual called $f(G, B)$ gravity.

Here is a simple breakdown of what the authors did, using everyday analogies.

1. The "Split" in the Recipe

In standard gravity, the "recipe" for how space bends involves a single ingredient called the Ricci Scalar (let's call it the "Curvature Score").

The authors of this paper decided to take that single ingredient and split it into two parts, like separating a cake batter into the cake (the main body) and the icing (the boundary).

• The Bulk ($G$): This is the main cake. It does the heavy lifting of creating gravity.
• The Boundary ($B$): This is the icing. In their theory, the icing doesn't actually change how the cake tastes (the dynamics); it just sits on the edge.

Why do this?
In many new gravity theories, the math gets so messy that it produces "ghosts"—mathematical errors that make the theory physically impossible (like predicting negative energy). By ignoring the "icing" (the boundary term) and only focusing on the "cake" (the bulk term), they managed to create a theory that is mathematically stable and doesn't have these ghosts.

2. The Goal: Modeling Stars

The authors wanted to see if this new theory could describe stars, specifically the dense, compact ones like neutron stars.

Usually, describing a star is like trying to solve a puzzle where you have to figure out the pressure inside the star at every single point. In standard physics, this is a very hard puzzle. In this new theory, the authors found something magical: The puzzle simplified itself.

They discovered that the equation governing the star's shape (the "isotropy equation," which ensures pressure is equal in all directions) could be rewritten in a way that didn't depend on the specific size of the star. It became a "self-contained" system.

3. The "Dynamical Systems" Approach: A Map of Possibilities

Instead of trying to solve the puzzle to find one single, perfect answer (which is often impossible), the authors used a technique called Dynamical Systems.

The Analogy: A Hiker on a Mountain
Imagine the possible shapes of a star are a vast, foggy mountain landscape.

• The Peaks and Valleys: These are "Fixed Points." They represent specific, stable shapes a star could take.
• The Rivers: These are "Trajectories." They show how a star evolves from one shape to another.
• The Goal: Instead of trying to walk every single path, the authors wanted to find the rivers that naturally flow toward the stable valleys. If a river flows toward a valley, it means that shape is stable and likely to exist in nature.

They created a "map" (called a Phase Portrait) of this landscape.

• They found that the landscape has invariant submanifolds. Think of these as highways or rivers that the system naturally wants to follow.
• Their analysis showed that many of these highways are attractors. If a star starts slightly off the highway, the laws of this new gravity will gently nudge it back onto the path. This suggests that stars formed under these rules would naturally settle into very specific, stable shapes.

4. The Vacuum Surprise: Two Roads to Nowhere

Before looking at stars, they looked at empty space (a vacuum). In Einstein's theory, empty space around a star has one unique shape (the Schwarzschild solution).

In this new theory, they found two possible shapes for empty space:

1. The Flat Road: A completely flat, boring space (like standard empty space).
2. The Curved Road with a Hole: A space that is curved but has a "singularity" (a point where the math breaks down, like a hole in the fabric of space).

The Catch: The "Curved Road" has a hole in it. However, the authors realized that if you build a star big enough, this hole sits outside the star, in the empty space beyond. So, the star itself is safe and stable, even though the empty space around it has a weird, singular feature. This is a new kind of "vacuum" that doesn't exist in Einstein's universe.

5. What Does This Mean for Us?

This paper is a theoretical "proof of concept." It shows that:

• It's possible to have a modified gravity theory that is mathematically clean (no ghosts).
• It's possible to describe stars within this theory.
• The stars behave predictably. The "attractor" behavior means that even if we don't know the exact details of the star's interior, we know it will likely settle into a stable, self-similar shape where the gravity and pressure balance out in a specific way.

The Bottom Line:
The authors didn't just find a new equation; they found a new map for how stars might behave if gravity works slightly differently than Einstein thought. They showed that in this new universe, stars have "natural highways" they tend to follow, leading to stable, compact objects.

The next step (which they mention they are working on) is to take these mathematical highways and match them with real-world data (like the mass and size of actual neutron stars) to see if this new theory fits our universe better than the old one.

Technical Summary

1. Problem Statement

The paper addresses the challenge of modeling stellar interiors within modified gravity theories, specifically the f(G, B) gravity framework proposed by Böhmer and Jensko.

• Context: Standard General Relativity (GR) faces issues regarding dark energy and renormalizability. Modified gravity theories (like $f(R)$ or $f(G)$) often introduce higher-order derivatives, leading to "ghost" instabilities (Ostrogradsky instability).
• The f(G, B) Proposal: This theory splits the Ricci scalar $R$ into a bulk term $G$ (quadratic in connection coefficients) and a boundary term $B$. By constructing the action using only the bulk term $G$ (and ignoring $B$ for dynamics), the theory yields field equations of second order, thereby avoiding ghost modes.
• The Gap: While $f(G, B)$ has been studied in cosmology, its application to stellar structure (compact objects like neutron stars) is unexplored. Furthermore, finding exact solutions for stellar interiors in modified gravity is notoriously difficult due to the complexity of the field equations.
• Objective: To derive the field equations for a pure quadratic $f(G, B)$ model, analyze vacuum solutions, and employ a dynamical systems approach to understand the stability and behavior of isotropic stellar models without relying solely on finding exact analytical solutions.

2. Methodology

The authors utilize a multi-stage analytical approach:

1. Formalism Derivation:

   • They start with a static, spherically symmetric metric in isotropic coordinates: $ds^2 = -e^{\nu(r)}dt^2 + e^{\mu(r)}(dx^2 + dy^2 + dz^2)$.
   • They derive the general field equations for an arbitrary function $f(G, B)$.
   • They focus on the pure quadratic case: $f(G) = \epsilon G^2$ (setting the linear Einstein term $\alpha G$ to zero to isolate the new physics).

2. Vacuum Analysis:

   • They set the energy-momentum tensor to zero ($\rho = p = 0$) to find vacuum solutions.
   • Unlike GR (where Birkhoff's theorem guarantees a unique Schwarzschild solution), they solve the resulting algebraic constraints to identify all possible vacuum branches.

3. Dynamical Systems Transformation:

   • The pressure isotropy condition ($p_r = p_t$) yields a complex, non-linear second-order differential equation (the "master equation").
   • Scale-Invariant Substitution: To convert this into an autonomous system, they introduce variables:
     • $U = r\mu'$ and $V = r\nu'$ (where prime denotes $d/dr$).
     • $X = \ln r$.
   • This transformation exploits the homogeneity of the equation, removing explicit $r$-dependence.

4. Phase Space Analysis:

   • Since the master equation is a single constraint on two variables ($\dot{U}, \dot{V}$), the authors introduce a gauge choice (specifically the minimal-norm projection) to split the constraint into a system of two first-order autonomous ODEs: $\dot{U} = F(U, V)$ and $\dot{V} = G(U, V)$.
   • They perform a phase plane analysis to locate fixed points and invariant submanifolds (curves where the system is at rest).
   • Stability Analysis: They calculate the Jacobian matrix to determine the eigenvalues, classifying the invariant submanifolds as stable (attracting) or unstable (repelling).
   • Poincaré Compactification: They map the infinite phase space to a finite disk to analyze the behavior of solutions at infinity.

3. Key Contributions and Results

A. Vacuum Solutions

The study reveals that the pure quadratic $f(G, B)$ theory admits two distinct vacuum branches, contrasting with the unique Schwarzschild solution in GR:

1. Flat Spatial Slices: A solution where $\nu$ is constant and $e^\mu \propto (r+A)^2$. The Kretschmann scalar vanishes, indicating a locally flat (Minkowski) spacetime in non-standard coordinates.
2. Curved Spacetime with Singularity: A solution where $e^\nu \propto (r+A)^4$ and $e^\mu \propto (r+A)^{-2}$. This spacetime is curved and possesses a power-law curvature singularity at $r = -A$. Crucially, this singularity is naked (not hidden by a horizon) but can be physically excised if the stellar radius is chosen such that the singularity lies outside the physical domain.

B. The Autonomous Isotropy Equation

A major theoretical achievement is the reduction of the pressure isotropy equation to an autonomous system via the variables $U$ and $V$.

• In standard GR, the isotropy equation does not typically admit such a scale-invariant autonomous form in these variables.
• This reduction allows for a global analysis of the solution space, revealing that the system admits curves of equilibrium states (invariant submanifolds) rather than just isolated fixed points.

C. Stability and Invariant Submanifolds

The phase portrait analysis identifies three primary invariant submanifolds:

1. $U = 0$ (The V-axis): Corresponds to $\mu' = 0$ (constant spatial conformal factor).
   • Stability: Unstable for $V < 8$, stable (attracting) for $V > 8$.
2. $U + 2V = 0$: Corresponds to the relation $\mu' = -2\nu'$ found in the second vacuum branch.
   • Stability: Unstable for $U > -16/3$, stable (attracting) for $U < -16/3$.
3. The Conic $Q(U, V) = 0$: A shifted hyperbola defined by $3U^2 + 4UV + 8U - V^2 + 8V = 0$.
   • Stability: Segments in the third quadrant ($U<0, V<0$) are transversely attracting.

D. Physical Interpretation

• Attractor Behavior: The fact that these submanifolds are attracting implies that physically reasonable stellar models (even those not exactly on the fixed curves) will asymptotically approach these geometric profiles.
• Metric Potentials: The attracting nature of the curves suggests that metric potentials $\mu$ and $\nu$ generally exhibit an inverse fall-off with radius ($\mu' \sim 1/r, \nu' \sim 1/r$).
• Stellar Structure: Trajectories attracted to the conic in the third quadrant correspond to decreasing pressure profiles, consistent with realistic stellar interiors.

4. Significance

• Novelty in Stellar Modeling: This is the first application of the Böhmer–Jensko $f(G, B)$ framework to stellar interiors, moving beyond its previous cosmological applications.
• Ghost-Free Higher-Order Gravity: It demonstrates that a theory with quadratic curvature terms can remain second-order and free of ghosts, providing a viable alternative to GR for compact objects.
• Methodological Advance: The successful reduction of the stellar isotropy equation to an autonomous dynamical system is rare in relativistic astrophysics. It provides a robust method to study stability and solution behavior without needing exact analytical solutions, which are often impossible to find in modified gravity.
• Rich Vacuum Structure: The discovery of two vacuum branches (one flat, one singular) challenges the uniqueness of the Schwarzschild solution and highlights new gravitational phenomenology accessible in $f(G, B)$ theories.
• Future Directions: The paper lays the groundwork for constructing complete stellar models by matching these interior dynamical trajectories to the derived vacuum exteriors using realistic equations of state (e.g., linear barotropic or polytropic), a task the authors note is currently in progress.

In summary, the paper establishes that pure quadratic $f(G, B)$ gravity offers a mathematically consistent and physically rich framework for stellar modeling, characterized by unique vacuum geometries and stable attractor solutions that guide the structure of compact stars.