Pith. Sign in

REVIEW

Linear stability matrix and Lyapunov exponents for stationary spacetime

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Xu Ye, Shao-Wen Wei

High Energy Physics - Theory · Submitted 14 Oct 2026 · 39 references

Pith's verdict: Extends Lyapunov exponent formalism to spherical/non-equatorial orbits in stationary spacetimes with a Killing tensor, but the stress-test time-parameter issue is real and undermines the claimed agreement with prior results.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 15, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 4 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.19418 (version 1)
Pith IDHEK4O77O, Pith's permanent record of this paper and its review
CitationShao-Wen Wei, Xu Ye. "Linear stability matrix and Lyapunov exponents for stationary spacetime." arXiv:2610.19418 (2026). https://pith.science/paper/HEK4O77O
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

Block-diagonal matrix yields Lyapunov exponents for stationary holes

The reading

The paper establishes that for a general stationary spacetime, the linear stability matrix describing geodesic perturbations has rank 4, requiring a four-dimensional phase space spanned by $(r,p_r,\theta,p_\theta)$. When the spacetime admits a rank-2 Killing tensor—such as in the Kerr family—the matrix decouples into two $2\times2$ blocks, one for radial perturbations and one for polar perturbations. This yields analytical Lyapunov exponents given by simple formulas involving second derivatives of the effective potentials $R(r)$ and $\Theta(\theta)$ and the metric functions. The result generalizes previous work that was restricted to equatorial circular orbits and two-dimensional phase space. By imposing appropriate constraints, the authors obtain explicit Lyapunov exponents for circular and spherical (non-equatorial) orbits, and demonstrate consistency with classical examples like the harmonic oscillator and Kepler motion, as well as with known results for Schwarzschild and Kerr black holes.

What carries the argument

The block-diagonal linear stability matrix $M_{LP}$ (equation 2.60) together with the effective radial and polar potentials $R(r)$ and $\Theta(\theta)$ defined via $\Sigma^2 \dot{r}^2 = R(r)$ and $\Sigma^2 \dot{\theta}^2 = \Theta(\theta)$, where $\Sigma = r^2 + a^2\cos^2\theta$. The rank-2 Killing tensor ensures separability of the Hamilton–Jacobi equation and decouples the $r$ and $\theta$ dynamics, allowing the matrix elements to be expressed purely in terms of $R$, $\Theta$, $\Delta_r$, $\Delta_\theta$, and their derivatives.

Extended reading notes

Core claim

The central discovery is that the linear stability matrix for geodesic perturbations in a stationary spacetime, when expressed in the basis $(\delta r, \delta p_r, \delta\theta, \delta p_\theta)$, can be reduced to block-diagonal form (equation 2.60) if a rank-2 Killing tensor exists. The diagonal blocks correspond to coupled radial and polar perturbations, with elements depending on the effective potentials $R(r)$ and $\Theta(\theta)$ and the metric components $\Delta_r(r)$ and $\Delta_\theta(\theta)$. The Lyapunov exponents are then given by equation 2.62, which explicitly separates radial and polar contributions. For circular orbits (constant radius and equatorial plane) and spherical orbits (constant radius but varying polar angle), the exponents simplify to expressions involving $R''(r_c)/2$ and $\Theta''(\theta_0)/2$ for the radial and polar sectors respectively. This provides a general method to compute stability properties for trajectories beyond the equatorial plane, a previously unresolved problem.

Load-bearing premise

The spacetime must positions a rank-2 Killing tensor so that the Hamilton–Jacobi equation separates and the linear stability matrix becomes block diagonal.

Editorial extensions

If this is right

  • The analytical formulas for Lyapunov exponents become directly applicable to any stationary black-hole spacetime that admits a rank-2 Killing tensor (e.g., Kerr–Newman–(AdS) without further calculation.
  • The method enables stability classification of spherical photon orbits, which are critical for black-hole shadow studies and for connecting quasinormal modes to geodesic properties in the eikonal limit.
  • The unified framework bridges classical mechanics (harmonic oscillator, Kepler) and relativistic systems, demonstrating that the same phase-space reduction and block-diagonalization procedure works across dynamical systems with mixed phase-space dimension.
  • The results provide a systematic way to compute Lyapunov exponents for non-equatorial orbits, which were previously inaccessible, and can be used to test chaotic behavior near black holes.
  • The explicit forms for $\lambda_{1,2}$ and $\lambda_{3,4}$ in terms of $R''$ and $\Theta''$ may lead to a geometric interpretation of stability via curvature of the effective potentials.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The requirement of a rank-2 Killing tensors is necessary for the block-diagonal structure; in generic stationary spacetimes without such a tensor, the linear stability matrix remains coupled and the Lyapunov exponents cannot be expressed transparently, so the framework's applicability is limited to the class of spacetimes with separable Hamilton–Jacobi equations.
  • Although the authors explicitly treat Schwarzschild and Kerr, the metric ansatz (2.27) also covers Kerr–Newman, (AdS) Kerr, and other Plebański–Demiański family solutions, so the results immediately extend to those spacetimes, potentially including cases with cosmological constant or electromagnetic charge.
  • The paper's classical examples (harmonic oscilator, invert oscilator, Kepler) serve as pedagogical checks but also suggest that the same Lyapunov exponent formula (2.62) might be transferable to any Hamiltonian system with two decoupled degrees of freedom and a conserved quantity analogous to the Carter constant.
Original abstract

Lyapunov exponents characterize the average rate at which nearby trajectories in phase space converge or diverge and therefore provide a useful criterion for dynamical stability. In this work, we develop a general framework for analyzing Lyapunov exponents in stationary spacetimes and extend previous studies beyond circular orbits. For a generic stationary spacetime described by $g_{\mu\nu}(r,\theta)$, we show that the phase space relevant to linear perturbations is spanned by $(r,p_r,\theta,p_\theta)$ and the corresponding linear stability matrix has rank four. When the spacetime admits a rank-2 Killing tensor, the linear stability matrix can be reduced to a block-diagonal form, allowing the Lyapunov exponents to be obtained analytically. By imposing the appropriate constraint conditions, the formalism further yields the linear stability matrices and Lyapunov exponents for circular and spherical orbits. The validity of the method is illustrated by several representative examples, including the (inverted) harmonic oscillator, planar Kepler motion, and particle motion in Schwarzschild and Kerr spacetimes. These results provide a unified framework for investigating the stability properties of particle orbits in both classical and relativistic dynamical systems.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-15: 1 request, 5,741 input tokens.

What the abstract itself states, as probabilities: claims first 0.07 · states limitation 0.12 · human participants 0.01 · controlled comparison 0.02 · randomized assignment 0.01 · states quantitative result 0.30 · compares with prior results 0.60 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Lyapunov exponents characterize the average rate at which nearby trajectories in phase space converge or diverge and therefore provide a useful criterion for dynamical stability. In this work, we develop a general framework for analyzing Lyapunov exponents in stationary spacetimes and extend previous studies beyond circular orbits. For a generic stationary spacetime described by $g_{\mu\nu}(r,\theta)$, we show that the phase space relevant to linear perturbations is spanned by $(r,p_r,\theta,p_\theta)$ and the corresponding linear stability matrix has rank four. When the spacetime admits a rank-2 Killing tensor, the linear stability matrix can be reduced to a block-diagonal form, allowing the Lyapunov exponents to be obtained analytically. By imposing the appropriate constraint conditions, the formalism further yields the linear stability matrices and Lyapunov exponents for circular and spherical orbits. The validity of the method is illustrated by several representative examples, including the (inverted) harmonic oscillator, planar Kepler motion, and particle motion in Schwarzschild and Kerr spacetimes. These results provide a unified framework for investigating the stability properties of particle orbits in both classical and relativistic dynamical systems.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Lyapunov exponents characterize the average rate at which nearby trajectories in phase space converge or diverge and therefore provide a useful criterion for dynamical stability. background 1.00 0.04 0.03 asserted 0.99
s1 In this work, we develop a general framework for analyzing Lyapunov exponents in stationary spacetimes and extend previous studies beyond circular orbits. contribution 1.00 0.94 0.12 asserted 1.00
s2 For a generic stationary spacetime described by $g_{\mu\nu}(r,\theta)$, we show that the phase space relevant to linear perturbations is spanned by $(r,p_r,\theta,p_\theta)$ and the corresponding linear stability matrix has rank four. result 0.54 0.93 0.19 asserted 1.00
s3 When the spacetime admits a rank-2 Killing tensor, the linear stability matrix can be reduced to a block-diagonal form, allowing the Lyapunov exponents to be obtained analytically. result 0.53 0.87 0.17 asserted 0.87
s4 By imposing the appropriate constraint conditions, the formalism further yields the linear stability matrices and Lyapunov exponents for circular and spherical orbits. result 0.50 0.93 0.13 asserted 0.99
s5 The validity of the method is illustrated by several representative examples, including the (inverted) harmonic oscillator, planar Kepler motion, and particle motion in Schwarzschild and Kerr spacetimes. method 0.35 0.88 0.08 asserted 1.00
s6 These results provide a unified framework for investigating the stability properties of particle orbits in both classical and relativistic dynamical systems. interpretation 0.99 0.59 0.60 asserted 0.99

What kind of paper the classifier reports this is: method or tool 0.88 · evidence from proof or derivation 1.00.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Trajectories of the one-dimensional harmonic oscillator. The curve x0(t) represents the background solution, while x1(t) = x0(t) + δx(t) denotes the perturbed trajectory. 3.2 One-dimensional inverted harmonic oscillator We next consider another example: the one-dimensional inverted harmonic oscillator. The corresponding Hamiltonian is H(px, x) = T + V = p 2 x 2m − 1 2 mω2x 2 . (3.12) – 14 –
Figure 2. Trajectories of the one-dimensional inverted harmonic oscillator. The curve x0(t) rep￾resents the background solution, while x1(t) = x0(t) + δx(t) denotes the perturbed trajectory. 3.3 Two-dimensional planar Kepler motion The final classical example is a two-dimensional circular orbit in Newtonian gravity. The corresponding Lagrangian is L = T − V = 1 2 m  r˙ 2 + r 2ϕ˙2  + GMm r . (3.23) Here, the gravitational potential is V (r) = − GMm r , and the kinetic energy is expressed in terms of the polar coordinates (r, ϕ). The canonical momenta conjugate to r and ϕ are pr =
Figure 3. Trajectories of the two-dimensional planar Kepler motion. The curve r0(t) represents the background solution, while r1(t) = r0(t) + δr(t) denotes the perturbed trajectory. 4 Relativistic example 1: Schwarzschild black holes In section 3, we investigated several classical examples. We now extend our analysis to relativistic systems, beginning with the Schwarzschild black hole, whose spacetime line element is given by ds2 = −∆r Σ [dt − χdϕ] 2 + sin2 θ∆θ Σ
Figure 4. Timelike circular orbit and its perturbed counterpart in Schwarzschild spacetime. (a) Radial motion with r0 = 8 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ0 = π/2 and θ1(τ ) = θ0 + δθ(τ ). The initial perturbations are chosen as δr(0) = 0.3 and δθ(0) = 0.04. The trajectories r1(τ ) and θ1(τ ) can be projected onto the (x, y) and (x, z) planes using the following spherical coordinate transformation: x = r sin θ cos ϕ, y = r sin θ sin ϕ, z = r cos θ. (4.31) The corresponding projected diagrams are shown in figure 5. Obviously, the circular orbits are stable from the linear…
Figure 5. Projections of a timelike circular orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Schwarzschild spacetime. Left: (x, y) plane. Right: (x, z) plane. – 21 –
Figure 6. Null circular orbit and its perturbed counterpart in Schwarzschild spacetime. (a) Radial motion with r0 = 3 and r1(τ ) = r0 +δr(τ ). (b) Polar motion with θ0 = π/2 and θ1(τ ) = θ0 +δθ(τ ). The initial perturbations are chosen as δr(0) = −0.1 and δθ(0) = 0.01. The projections of the trajectories r1(τ ) and θ1(τ ) onto the (x, y) and (x, z) planes are shown in figure 7. It is clear that the radial motion is unstable, while the θ motion is stable, which is consistent with teh Lyapunov exponents obtained in Eq. (4.34). – 22 –
Figure 7. Projections of a null circular orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Schwarzschild spacetime. Left: (x, y) plane. Right: (x, z) plane. 4.2 Lyapunov exponents for spherical orbits We assume that the spherical orbits are located at r = rs, while θ oscillates between θmin and θmax. The radial conditions for spherical orbits are R(rs) = 0, R′ (rs) = 0. (4.35) The first condition implies that pr = p R(rs)/∆r(rs) = 0 at r = rs. Moreover, at the turning points θmin and θmax, one has Θ(θmin) = Θ(θmax) = 0. Then the Carter constant is obtained as Cb= L 2 csc2…
Figure 8. Timelike spherical orbit and its perturbed counterpart in Schwarzschild spacetime. (a) Radial motion with r0 = 8 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ(τ ) ∈ [ π 3 , 2π 3 ] and θ1(τ ) = θ(τ ) + δθ(τ ). The initial perturbations are chosen as δr(0) = 0.1 and δθ(0) = 0.002. r0(τ) r1(τ) -5 0 5 -5 0 5 x y r0(τ) r1(τ) -5 0 5 -5 0 5 x z
Figure 9. Projections of a timelike spherical orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Schwarzschild spacetime. Left: (x, y) plane. Right: (x, z) plane. Note that, in the null case, the ratio L/E is invariant under reparametrizations of the affine parameter. Using the above expressions for the angular momentum L and orbital radius rs, we obtain the Lyapunov exponents for null spherical orbits λ1 = −λ2 = −3 √ 3Em, (4.51) λ3 = −λ4 = 3i √ 3Emp cos(2θ) + 2 csc2 (θ) sin(θ0). (4.52) For the null spherical orbit at rs = 3m in Schwarzschild spacetime, λ1 and λ2 are real,
Figure 10. Null spherical orbit and its perturbed counterpart in Schwarzschild spacetime. (a) Radial motion with r0 = 3 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ(τ ) ∈ [ π 3 , 2π 3 ] and θ1(τ ) = θ(τ ) + δθ(τ ). The initial perturbations are chosen as δr(0) = 0.01 and δθ(0) = 0.001. r0(τ) r1(τ) -4 -2 0 2 4 -4 -2 0 2 4 x y r0(τ) r1(τ) -3 -2 -1 0 1 2 3 -2 -1 0 1 2 x z
Figure 11. Projections of a timelike spherical orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Schwarzschild spacetime. Left: (x, y) plane. Right: (x, z) plane. 5 Relativistic example 2: Kerr black holes In the previous section 4, we discussed the Schwarzschild black hole, which describes a static and spherically symmetric spacetime. In this section, we turn our attention to the – 26 –
Figure 12. Timelike circular orbit and its perturbed counterpart in Kerr spacetime. (a) Radial motion with r0 = 6 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ0 = π 2 and θ1(τ ) = θ0 + δθ(τ ). The initial perturbations are chosen as δr(0) = 0.2 and δθ(0) = 0.04. r0(τ) r1(τ) -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 x y r0(τ) r1(τ) -6 -4 -2 0 2 4 6 -3 -2 -1 0 1 2 3 x z
Figure 13. Projections of a timelike circular orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Kerr spacetime. Left: (x, y) plane. Right: (x, z) plane. 5.1.2 Null case For the null circular orbits with m0 = 0 in Kerr spacetime, the radial effective potential from (5.4) is R(r) = E
Figure 14. Null circular orbit and its perturbed counterpart in Kerr spacetime. (a) Radial motion with r0 = 1.8110 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ0 = π 2 and θ1(τ ) = θ0 + δθ(τ ). The initial perturbations are chosen as δr(0) = −0.01 and δθ(0) = 0.005. – 31 –
Figure 15. Projections of a null circular orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Kerr spacetime. Left: (x, y) plane. Right: (x, z) plane. 5.2 Lyapunov exponents for spherical orbits We now investigate spherical orbits in the Kerr spacetime. Such orbits are characterized by a constant radial coordinate r = rs, whereas the polar angle θ varies within the range θmin ≤ θ ≤ θmax. The conditions for the radial and polar motions can be expressed as R(rs) = 0, R′ (rs) = 0, (5.46) Θ(θmin) = Θ(θmax) = 0. (5.47) From R(rs) = 0 and Θ(θmin) = 0, it follows that pr = p…
Figure 16. Timelike spherical orbit and its perturbed counterpart in Kerr spacetime. (a) Radial motion with r0 = 6 and r1(τ ) = r0 + δr(τ ). (b) Polar motion with θ(τ ) ∈ [ π 3 , 2π 3 ] and θ1(τ ) = θ(τ ) + δθ(τ ). The initial perturbations are chosen as δr(0) = 0.04 and δθ(0) = 0.005. r0(τ) r1(τ) -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 x y r0(τ) r1(τ) -6 -4 -2 0 2 4 6 -6 -4 -2 0 2 4 6 x z
Figure 17. Projections of a timelike spherical orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Kerr spacetime. Left: (x, y) plane. Right: (x, z) plane. From R(rs) = 0, the relation between L and E can be obtained as L1 = 4aEmrs − E csc(θ0) p a 2 − 2mrs + r 2 s
Figure 18. Projections of a null spherical orbit and its perturbed orbit onto the (x, y) and (x, z) planes in Kerr spacetime. Left: (x, y) plane. Right: (x, z) plane. Here, the null spherical orbit is described by r0 = 1.8832 and θ(τ ) ∈ [ π 3 , 2π 3 ]. The perturbed orbit is described by r1(τ ) = r0 +δr(τ ) and θ1(τ ) = θ(τ ) + δθ(τ ). 6 Conclusion In the present work, we studied the linear stability matrix and Lyapunov exponents in general stationary spacetimes. Our main results are summarized as follows. – 36 –

Figures reproduced from arXiv:2610.19418 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.