Quadrupolar tidal effects destroy the integrability of black hole geodesics: analytic proof and numerical evidence of chaos

Technical Summary: Quadrupolar Tidal Effects and the Loss of Integrability in Kerr Geodesics

Problem Statement
In General Relativity, the motion of a test mass around a rotating (Kerr) black hole is described by geodesics. Due to the spacetime's symmetries, specifically the existence of a Killing–Stäckel tensor, these geodesics form an integrable Hamiltonian system with four constants of motion: energy ($E$), axial angular momentum ($L_z$), dynamical mass ($\mu$), and the Carter constant ($Q_0$). This integrability is the foundational pillar for modeling Extreme Mass-Ratio Inspirals (EMRIs), the primary sources for future space-based gravitational wave detectors like LISA. It allows for the use of action-angle variables, fundamental frequencies, and two-timescale expansions to generate waveforms.

However, real compact bodies are not structureless test masses; their internal structure couples to the background curvature. While linear-in-spin corrections preserve integrability (via the Rüdiger invariant), the effect of a tidally induced quadrupole on a non-spinning body remains an open question. This paper investigates whether the leading-order tidal dynamics in a Kerr background retains the crucial property of integrability. Specifically, it asks: Do the equations of motion remain Hamiltonian? Do they admit four independent constants of motion in involution?

Methodology
The authors employ a dual approach combining rigorous analytic proofs with numerical diagnostics.

1. Hamiltonian Formulation:

   • The authors first establish a covariant Hamiltonian formulation for a non-spinning particle with a tidally induced quadrupole in an arbitrary background spacetime.
   • They demonstrate that the Mathisson-Papapetrou-Tulczyjew-Dixon (MPTD) equations, truncated at quadrupolar order, are generated by a Hamiltonian $H = H_0 + H_1$ on an 8-dimensional phase space.
   • $H_0$ is the standard geodesic Hamiltonian, while $H_1$ represents the tidal correction. Crucially, they define a conserved mass $\mu_T$ (conserved to leading order) and show that $H = -\mu_T^2/2$.

2. Analytic Non-Integrability Proof:

   • Specializing to the Kerr spacetime, the authors derive closed-form expressions for the electric ($E_2$) and magnetic ($B_2$) tidal scalars. A key intermediate result is showing that these scalars depend on the particle's momentum only through the geodesic Carter constant $Q_0$ and the Weyl scalar $\Psi$.
   • They attempt to construct a deformed Carter constant $Q = Q_0 + Q_1$ (where $Q_1$ is linear in tidal effects) that remains conserved under the full Hamiltonian.
   • Using the cohomological equation $\{Q_1, H_0\} = -\{Q_0, H_1\}$, they reduce the general polynomial Ansatz for $Q_1$ to a unique form involving three unknown scalar functions of the coordinates $(r, \theta)$.
   • They derive necessary integrability conditions (Schwarz's theorem applied to the resulting partial differential equations) and demonstrate that these conditions are violated for generic tidal couplings and any non-zero Kerr spin ($a \neq 0$). The only exception is the Schwarzschild limit ($a=0$), where spherical symmetry preserves integrability.

3. Numerical Diagnostics:

   • To confirm the analytic result, the authors perform numerical integrations of the tidally perturbed equations of motion.
   • They utilize three standard chaos diagnostics:
     • Poincaré Sections: Visualizing the phase space structure by recording turning points.
     • Lyapunov Exponents: Measuring the exponential divergence of neighboring orbits.
     • Escape-Time Maps: Analyzing the fractal dependence of plunge times on initial conditions.
   • A scaling analysis is performed to ensure the observed chaos is a genuine first-order effect in the tidal coupling $\epsilon$, distinguishing it from numerical artifacts or higher-order ($O(\epsilon^2)$) effects.

Key Contributions and Results

• Hamiltonian Structure: The paper provides a covariant Hamiltonian formulation for tidally induced quadrupole dynamics, valid in arbitrary backgrounds, identifying the conserved mass $\mu_T$ and the specific form of the tidal Hamiltonian.
• New Tidal Scalar Formulae: The authors derive novel closed-form expressions for the tidal invariants $E_2$ and $B_2$ in Kerr spacetime, expressing them solely in terms of the Weyl scalar and the normalized Carter constant.
• Proof of Non-Integrability: The central result is an analytic proof that no polynomial-in-momenta deformation of the Carter constant exists for a non-spinning body with a generic tidally induced quadrupole in a Kerr background ($a \neq 0$). Consequently, the leading-order tidal dynamics is non-integrable.
• Numerical Confirmation: The numerical evidence robustly supports the analytic proof:
  • Poincaré sections reveal the destruction of invariant tori, the emergence of stochastic layers, and resonant island chains.
  • Positive Lyapunov exponents confirm the exponential divergence of trajectories in the chaotic sea.
  • Escape-time maps exhibit fractal basin boundaries, indicating extreme sensitivity to initial conditions.
• Order of Effects: The scaling analysis confirms that the observed chaotic behavior (specifically the variation in turning points) scales linearly with the tidal coupling ($\Delta r \propto \epsilon$), proving that the loss of integrability is a leading-order effect, not a sub-leading artifact.

Significance and Claims
The paper claims that the integrability of Kerr geodesics is fragile: while spin-induced quadrupoles can preserve integrability (conditionally on the body being a black hole with a specific deformability parameter), tidally induced quadrupoles generically destroy integrability.

The authors frame this as a dynamical characterization of black holes: at the quadrupolar order, only black holes (which have vanishing tidal Love numbers, i.e., $c_E = c_B = 0$) preserve the integrable structure. Any other compact body, possessing non-zero tidal deformability, introduces non-integrable dynamics.

Regarding astrophysical implications, the paper remains modest. It acknowledges that for realistic EMRIs, the tidal coupling is extremely small, and KAM theory suggests that most invariant tori will survive as slightly deformed structures, with chaos confined to thin stochastic layers. Therefore, the immediate impact on waveform modeling for LISA may be limited. However, the loss of a fourth constant of motion fundamentally alters the mathematical tractability of the problem, potentially affecting resonance crossings and the precise definition of plunge times. The work establishes a rigorous baseline for understanding how internal structure breaks the symmetries that make black hole perturbation theory solvable.