Jacobi stability of circular orbits around conformally invariant Weyl gravity black holes

This paper investigates the timelike circular geodesics of spherically symmetric Weyl black holes by analyzing their effective potentials and assessing the Jacobi and Lyapunov stability of these orbits to provide new insights into the stability properties of Weyl conformal gravity and the influence of its free parameters.

The Gist

The Big Picture: A New Kind of Gravity

Imagine that for the last 100 years, we've been driving a car (General Relativity) that works perfectly on city streets (our Solar System). But when we try to drive it across the vast, empty deserts of the universe (galaxies), the car starts to sputter. We need "Dark Matter" (a mysterious, invisible fuel) to keep it moving, but nobody has ever seen this fuel.

Enter Weyl Gravity. Think of this as a new, experimental engine design proposed back in the 1920s. It suggests that gravity isn't just about the shape of space, but also about how that shape scales (like zooming in or out on a map). This paper asks: If we swap our old engine for this new Weyl engine, do black holes still behave safely?

The Main Characters: The Black Hole and the Test Particle

The authors are studying a specific type of black hole solution in this new gravity theory.

• The Black Hole: It's not just a simple sphere like in Einstein's theory. It has a few extra "knobs" or dials (parameters named $\beta$, $\gamma$, and $k$) that can be turned. These knobs change how the gravity feels as you get closer or further away.
• The Test Particle: Imagine a tiny, fearless astronaut on a spaceship trying to orbit this black hole. We want to know: Can they stay in a perfect circle without crashing or flying off into space?

The Concept: The "Cosmic Roller Coaster"

To understand if the astronaut is safe, the authors look at something called the Effective Potential.

• The Metaphor: Imagine the gravitational field is a giant, 3D landscape. The astronaut is a marble rolling on this landscape.
• The Valley: If the marble rolls into a deep valley, it gets stuck at the bottom. This is a stable orbit. If you nudge the marble, it wobbles a bit but stays in the valley.
• The Hilltop: If the marble is balanced on the very peak of a hill, it's an unstable orbit. The slightest nudge sends it tumbling down the other side.
• The Flat Plain: If the landscape is flat, the marble just rolls away.

The authors mapped out this landscape for the Weyl black hole. They found that depending on how you turn those "knobs" ($\beta, \gamma, k$), the shape of the valley changes. Sometimes, the valley is deep and safe; other times, it's a precarious peak.

The Two Ways to Check Safety

The paper uses two different "safety inspectors" to check if the astronaut's orbit is stable.

1. The "Linear" Inspector (Lyapunov Stability)

This is the standard, common-sense check.

• The Analogy: You gently push the marble. Does it wobble and come back to the center, or does it roll away?
• The Math: They look at the slope of the hill right where the marble is sitting. If the slope curves upward (a valley), it's safe. If it curves downward (a hilltop), it's dangerous.

2. The "Geometric" Inspector (Jacobi Stability)

This is the fancy, high-tech check introduced in the paper.

• The Analogy: Instead of just pushing one marble, imagine you have a whole swarm of marbles orbiting side-by-side. Jacobi stability asks: Do these marbles stay close together, or do they scatter apart like a flock of birds?
• The Math: This uses a complex geometric theory (KCC theory) to measure how the "fabric" of space itself bends the paths of these neighboring marbles. It's like checking if the road itself is twisting in a way that forces cars apart, even if the road looks flat.

The Big Discovery: They Agree!

In many complex physics systems, these two inspectors often disagree. One might say, "It's stable!" while the other says, "It's chaotic!"

The main result of this paper is that for Weyl Black Holes, they agree perfectly.

• If the "Linear" inspector says the orbit is a safe valley, the "Geometric" inspector also says the swarm of marbles stays together.
• If the Linear inspector says it's a dangerous hill, the Geometric inspector says the swarm will scatter.

This is a huge relief for physicists. It means the math is consistent. The "knobs" on the black hole ($\gamma$ and $k$) change the size of the safe zone, but they don't break the rules of stability.

The "Innermost" Limit

The authors also calculated the Innermost Stable Circular Orbit (ISCO).

• The Metaphor: Imagine a "No Parking" zone around the black hole. If you get too close, no amount of steering can keep you in a circle; you will inevitably crash.
• The Finding: In our standard Einstein universe, this "No Parking" line is at a specific distance (6 times the black hole's radius). In the Weyl universe, this line moves depending on the extra "knobs." If you turn the knobs just right, you can get closer to the black hole safely, or you might have to stay further away.

Why Does This Matter?

This paper is like a stress test for a new theory of gravity.

1. It validates the theory: It shows that Weyl gravity produces black holes that behave logically and predictably, just like the ones we see in Einstein's theory.
2. It offers a new tool: By proving that these two stability methods match, the authors give future scientists a reliable way to test other weird gravity theories. If a theory fails this stability test, it's probably wrong.
3. It helps us understand the Universe: If we can measure how stars orbit black holes in distant galaxies, we might be able to tell if the universe is running on Einstein's engine or Weyl's engine.

Summary

The authors took a complex, alternative theory of gravity (Weyl), built a model of a black hole within it, and asked: "Is it safe to orbit here?" They used two different mathematical methods to check. They found that yes, it is safe, and both methods agree on exactly where the safe zones are. This gives us confidence that this alternative theory is a viable candidate for explaining the mysteries of the universe without needing invisible "Dark Matter."

Technical Summary

1. Problem Statement

The paper addresses the stability of timelike circular geodesics (orbits of massive particles) around static, spherically symmetric black holes in Conformally Invariant Weyl Gravity. While General Relativity (GR) successfully describes local gravitational phenomena, it faces challenges at cosmological scales (e.g., dark energy, dark matter, Hubble tension). Weyl gravity, a fourth-order theory based on the conformal invariance of the Weyl tensor, offers an alternative framework that can explain galactic rotation curves without dark matter.

The specific problem investigated is whether the circular orbits in this modified gravity theory are stable. The authors aim to determine if the stability criteria differ from standard GR and to compare two distinct stability analysis methods: Lyapunov stability (linear stability) and Jacobi stability (geometric stability based on the Kosambi-Cartan-Chern theory).

2. Methodology

A. Theoretical Framework: Weyl Gravity

The authors utilize the exact vacuum solution for a static, spherically symmetric black hole derived by Mannheim and Kazanas. The metric is given by:
$$ds^2 = -B(r)dt^2 + \frac{dr^2}{B(r)} + r^2(d\theta^2 + \sin^2\theta d\phi^2)$$
where the lapse function $B(r)$ is:
$$B(r) = 1 - \frac{\beta(2-3\beta\gamma)}{r} - 3\beta\gamma + \gamma r - kr^2$$
Here, $\beta$ is related to the mass, $\gamma$ measures the deviation from GR, and $k$ is related to the cosmological curvature scale. The authors focus on the regime where $\gamma$ is small and $k$ is negative.

B. Geodesic Dynamics

1. Lagrangian Formalism: The authors derive the equations of motion for a massive test particle using the Euler-Lagrange formalism.
2. Effective Potential: By restricting motion to the equatorial plane ($\theta = \pi/2$) and utilizing conserved energy ($E$) and angular momentum ($L$), they derive the effective potential $V_{eff}(r)$:
   $$V_{eff}(r) = \frac{1}{2}B(r)\left(\frac{L^2}{r^2} + 1\right)$$
3. Circular Orbits: Circular orbits are identified by the conditions $V'_{eff}(r) = 0$ and $V_{eff}(r) = E^2$. The innermost stable circular orbit (ISCO) is found by solving $V'_{eff}(r) = 0$ and $V''_{eff}(r) = 0$ simultaneously.

C. Stability Analysis

The paper employs two complementary approaches:

1. Lyapunov Stability (Linear):

   • The system is converted into a first-order dynamical system in phase space $(r, p)$.
   • The Jacobian matrix $J$ of the system is computed at equilibrium points.
   • Stability is determined by the eigenvalues of $J$. If the real parts are negative (or purely imaginary for centers in conservative systems), the orbit is stable. Specifically, for circular orbits, stability depends on the sign of $V''_{eff}(r)$.

2. Jacobi Stability (Geometric/KCC Theory):

   • Based on the Kosambi-Cartan-Chern (KCC) theory, which treats the system of second-order ODEs as geodesics on a manifold.
   • The stability is determined by the deviation curvature tensor (second KCC invariant), denoted as $P^i_j$.
   • For the specific 2D system analyzed, the condition for Jacobi stability reduces to the sign of the discriminant $\Delta = (\text{tr} J)^2 - 4\det J$.
   • Criterion: The system is Jacobi stable if $\Delta < 0$ (equivalent to the deviation curvature being negative).

3. Key Contributions

• Derivation of ISCO in Weyl Gravity: The authors derive a quintic equation (Eq. 66) that determines the radius of the Innermost Stable Circular Orbit (ISCO) as a function of the Weyl parameters ($\tilde{\gamma}, \tilde{k}$). They show that in the limit where Weyl parameters vanish, the result recovers the standard Schwarzschild ISCO ($r = 6\beta$).
• Equivalence of Stability Criteria: A primary theoretical contribution is the rigorous demonstration that for circular timelike geodesics in this specific spacetime, Jacobi stability and Lyapunov stability are equivalent.
  • They prove that the condition $V''_{eff}(r) > 0$ (Lyapunov stable) is mathematically identical to the condition $\Delta < 0$ (Jacobi stable).
  • This contrasts with general dynamical systems where these two stability notions often diverge.
• Parameter Space Analysis: The paper analyzes how the stability regions change with the Weyl parameters. It confirms that circular orbits only exist for radii $r > 3\beta$ (outside the photon sphere) and investigates the existence of multiple circular orbits (stable and unstable) depending on the angular momentum $L$.

4. Results

• Orbital Existence: Circular orbits exist only for $r > 3\beta$. As $r \to 3\beta$, the required energy and angular momentum diverge.
• Stability Regions:
  • Unstable Orbits: Correspond to local maxima of the effective potential ($V''_{eff} < 0$). These are Jacobi unstable.
  • Stable Orbits: Correspond to local minima of the effective potential ($V''_{eff} > 0$). These are both Lyapunov and Jacobi stable.
  • ISCO: The transition point where $V''_{eff} = 0$ marks the boundary between stable and unstable regions.
• Numerical Verification: Using specific parameter values ($\beta=0.1, \gamma=0.1, k=-0.045$), the authors plot the effective potential and phase portraits.
  • For low angular momentum ($L=0.2$), no circular orbits exist.
  • For critical angular momentum ($L \approx 0.3793$), a single marginally stable orbit (cusp) exists.
  • For higher angular momentum ($L=0.5$), two orbits exist: an unstable one at $r \approx 0.353$ and a stable one at $r \approx 0.962$.
• Consistency: In all tested cases, the regions identified as stable by linearization (Lyapunov) perfectly match the regions identified as stable by the geometric deviation tensor (Jacobi).

5. Significance

• Validation of Geometric Methods: The paper validates the use of KCC theory and Jacobi stability as robust tools for analyzing gravitational systems, showing they yield consistent physical predictions with standard linear stability analysis in the context of black hole orbits.
• Insight into Weyl Gravity: The results provide a deeper understanding of the dynamical properties of Weyl black holes. The fact that the stability criteria remain equivalent to GR-like behavior (despite the fourth-order field equations) suggests that the conformal modifications primarily shift the location of stability boundaries (ISCO) rather than fundamentally altering the nature of orbital stability.
• Observational Relevance: By characterizing the ISCO and stability regions, the study offers potential observational signatures. Deviations in the ISCO radius or the stability of accretion disk orbits could, in principle, distinguish Weyl gravity from General Relativity in future high-precision astrophysical observations (e.g., via X-ray spectroscopy of accretion disks).
• Generalizability: The methodology presented can be extended to other modified gravity theories and black hole solutions, offering a unified framework for assessing the robustness of dynamical systems in curved spacetime.

In conclusion, the paper successfully bridges geometric stability theory (KCC) with classical dynamical systems analysis in the context of conformal gravity, confirming that for circular orbits around Weyl black holes, the two stability concepts are not just related but equivalent.