Particles with precessing spin in Kerr spacetime: analytic solutions for eccentric orbits and homoclinic motion near the equatorial plane

Basic Information

• Paper ID: 2510.09597
• Title: Particles with precessing spin in Kerr spacetime: analytic solutions for eccentric orbits and homoclinic motion near the equatorial plane
• Author: Gabriel Andres Piovano (Université Libre de Bruxelles & Université de Mons)
• Classification: gr-qc (General Relativity and Quantum Cosmology)
• Submission Date: October 15, 2025 (v2)
• Paper Link: https://arxiv.org/abs/2510.09597

Abstract

This paper presents a family of analytic solutions for near-equatorial motion of test particles with precessing spin in Kerr spacetime. The author solves the equations of motion for small bodies at linear order in spin, encompassing both periodic and homoclinic orbits. At zeroth order, particles move along equatorial geodesics. When spin-curvature forces are introduced, geodesic corrections arise, leading to orbital plane precession for general spin orientations. The study derives eccentric orbit solutions for general reference geodesics (termed "spin gauges") using Legendre elliptic integrals and Jacobi elliptic functions. The analytic solutions show perfect agreement with numerical orbits from Drummond and Hughes as well as Piovano et al. Furthermore, analytic solutions for homoclinic orbits of spinning particles in Kerr spacetime are presented for the first time, along with spin corrections to the separatrix location. Homoclinic orbits are described in closed form using elementary functions. Finally, a new parametrization called the "fixed eccentricity spin gauge" is introduced, which is the unique spin gauge where periodic orbit corrections remain finite at the geodesic separatrix.

Research Background and Motivation

Scientific Background

1. Gravitational wave detection requirements: Future space-based gravitational wave detectors such as LISA, TianQin, and Taiji will observe gravitational wave sources inaccessible to ground-based detectors, particularly extreme mass ratio inspirals (EMRIs).
2. Importance of EMRI systems: EMRIs consist of stellar-mass compact objects captured by supermassive black holes (SMBHs), with mass ratios q = μ/M ∼ 10⁻⁵-10⁻⁷. The small body completes 1/q orbits in the strong gravitational field, producing rich gravitational wave signals.
3. Necessity of precise modeling: EMRI signals can persist in the detector band for years, requiring extremely high-precision waveform models to avoid systematic biases in parameter estimation.

Theoretical Challenges

1. Multi-timescale expansion framework: At the first post-adiabatic (1PA) order, one must include conservative first-order self-force effects, dissipative second-order effects, and secondary body spin effects.
2. Complexity of spin dynamics: Post-geodesic effects from spin-curvature forces contribute to the post-adiabatic phase at linear spin order. While numerical methods have advanced, analytic solutions remain lacking.
3. Singularities near the separatrix: Existing methods exhibit divergences near the geodesic separatrix, preventing continuous treatment of the transition from stable orbits to plunge.

Core Contributions

1. First complete family of analytic solutions: Provides complete analytic solutions for near-equatorial motion of spinning particles in Kerr spacetime, including both periodic and homoclinic orbits.
2. Novel spin gauge method: Introduces the "fixed eccentricity gauge" (FE gauge), the unique parametrization maintaining finite corrections at the geodesic separatrix.
3. First analytic description of homoclinic orbits: Presents the first closed-form solutions for homoclinic orbits of spinning particles in Kerr spacetime using elementary functions.
4. Spin corrections to separatrix location: Provides analytic expressions for linear spin corrections to the separatrix location.
5. Numerical verification: Analytic solutions show perfect agreement with existing numerical results, validating the methodology.

Detailed Methodology

Theoretical Framework

Mathisson-Papapetrou-Dixon Equations

Spinning particle dynamics are described by the MPD equations:

Dpμ/dτ = -1/2 Rμνρσ vν Sρσ
DSμν/dτ = 2p[μvν]

where the first term on the right represents the spin-curvature force, the primary post-geodesic effect from spin.

Linearization Treatment

Using the small parameter expansion S ∼ qχ ≪ 1, the MPD equations are linearized to:

Dgvμg/dτ = 0
Dgδvμ/dτ = -1/2 Rμνρσ vνg sρσ  
Dgssμ/dτ = 0

Near-equatorial Motion Decomposition

Near-equatorial motion decomposes into:

1. In-plane motion: Planar motion constrained to the equatorial plane
2. Vertical oscillations: Small oscillations perpendicular to the orbital plane, causing slow orbital plane precession

Spin Gauge Theory

Three Primary Gauges

1. Fixed constant gauge (FC): δE = δLz = 0
2. Fixed turning point gauge (FT): δr₁ = δr₂ = 0
3. Fixed eccentricity gauge (FE): Newly proposed gauge ensuring δr₂ = δr₃

Advantages of FE Gauge

The FE gauge eliminates singular terms near the separatrix through the condition δr₂ = δr₃, ensuring:

• Periodic orbit corrections remain finite at the separatrix
• Continuous transition to homoclinic orbits
• Smooth reduction to the last stable orbit in appropriate limits

Solution Methods

Radial Potential Function

The spin-corrected radial potential is a fifth-order polynomial:

R₅(r) = Σ(i=0 to 5) cᵢ rⁱ

The structure of its roots determines the orbit type.

Elliptic Function Solutions

Periodic orbit solutions are expressed using Legendre elliptic integrals and Jacobi elliptic functions:

δr(rg) = √Rg(rg) [(...) I^per(rg) + Σ δrᵢ I^per_1/(r-rᵢ)(rg) + ...]

Elementary Function Solutions for Homoclinic Orbits

In the limit r₃g → r₂g, elliptic functions degenerate to elementary functions:

sin(ξ) = tanh(1/2 Yr λ)
F(ξ|1) = 1/2 Yr λ

Experimental Setup

Numerical Verification Methods

1. Comparison benchmarks: Comparison with numerical orbits from Drummond-Hughes and Piovano et al.
2. Code implementation: Using the "KerrGeodesic" package from the Black Hole Perturbation Toolkit
3. Parameter ranges: Covering various black hole spin parameters a, orbital eccentricities eg, and spin parameters χ

Test Cases

• Zoom-whirl orbits: Intermediate mass ratio systems with q = 10⁻²
• Near-separatrix orbits: Orbits with pg = p*g + 1/10 near the separatrix
• Homoclinic orbits: Separatrix orbits with varying eccentricities

Experimental Results

Main Verification Results

1. Perfect numerical agreement: Analytic solutions for all three spin gauges (FT, FC, FE) show perfect agreement with numerical results
2. Separatrix behavior:
   • FT and FC gauges diverge near the separatrix
   • Only FE gauge remains finite, with continuous transition to homoclinic orbits
3. Special case verification:
   • Schwarzschild limit: δp*(0, eg) = -2√2(1+eg)/(3+e*g)
   • ISCO corrections: Consistent with Jefremov et al. results
   • Extreme Kerr: Separatrix corrections vanish

Orbital Characteristics Analysis

Periodic Orbits

• Spin corrections to radial, temporal, and azimuthal orbits appear as linear combinations of elliptic functions
• Correction terms include oscillatory and long-term evolution components
• Frequency corrections require Hadamard finite part regularization

Homoclinic Orbits

• Radial and azimuthal orbits exhibit logarithmic divergence near periapsis
• Spin corrections maintain the same divergence structure
• Frequencies remain finite at the separatrix: Υtg = Tg(r₂g), Υϕg = Φg(r₂g)

Physical Insights

Orbital Plane Precession

The polar coordinate equation of motion:

d²δz/dλ² + Υ²zg δz = -3s⊥ cos ψp/(r²) Lzrd √(L²zrd + r²g)

shows that orbital plane precession is driven by the perpendicular spin component χ⊥.

Spin Precession Phase

The evolution equation for precession phase ψp(rg) is separable with analytic expressions for both periodic and homoclinic orbits.

Comparison with Related Work

Limitations of Existing Methods

1. Witzany's virtual geodesic method: While providing analytic solutions, the deformed Mino time is not separable
2. Numerical methods: Drummond-Hughes frequency-domain methods are efficient but lack analytic insights
3. Hackmann et al.'s hyperelliptic function solutions: Applicable only to aligned spin with complex expressions

Advantages of This Work

1. Physical transparency: Solutions expressed as physical geodesics plus spin corrections
2. Complete separability: All integrals expressible in standard special functions
3. Continuity: FE gauge ensures smooth transitions from periodic to homoclinic orbits

Application Prospects

Gravitational Wave Astrophysics

1. EMRI waveform modeling: Provides analytic foundation for dual-timescale frameworks
2. Transition to plunge modeling: Homoclinic orbit solutions are particularly important for IMRI systems
3. Parameter estimation accuracy: Analytic expressions facilitate rapid waveform generation

Theoretical Extensions

1. Quadrupole moment effects: Can be extended to include spin-induced quadrupole moments
2. Higher-order spin effects: Provides foundation for second-order spin terms
3. Other spacetimes: Methods generalizable to Reissner-Nordström and other spacetimes

Conclusions and Discussion

Main Conclusions

1. Successfully constructs a complete analytic theory for near-equatorial motion of spinning particles in Kerr spacetime
2. FE gauge resolves singularities near the separatrix
3. Presents the first closed-form analytic solutions for homoclinic orbits
4. All results rigorously verified through numerical comparison

Limitations

1. Near-equatorial restriction: Applicable only to orbits with small inclination angles
2. Linear spin approximation: High-spin systems require inclusion of higher-order terms
3. Test particle approximation: Neglects self-force effects

Future Directions

1. Quasi-spherical orbit extension: Generalize methods to orbits with arbitrary inclination
2. Plunge trajectory study: Construct analytic solutions for bound-to-plunge orbits
3. Numerical-analytic mapping: Establish precise correspondence with existing numerical methods

In-Depth Evaluation

Strengths

1. Theoretical rigor: Mathematically rigorous derivations with clear physical interpretation
2. Methodological innovation: Introduction of FE gauge resolves long-standing technical issues
3. High practical value: Provides important tools for gravitational wave data analysis
4. Comprehensive verification: Comparison with multiple independent numerical results validates correctness

Technical Highlights

1. Clever application of elliptic functions: Transforms complex gravitational dynamics into standard elliptic integrals
2. Regularization techniques: Uses Hadamard finite part to handle divergent integrals
3. Continuity treatment: FE gauge ensures smooth transitions between different orbit types

Potential Impact

1. Academic influence: Provides new theoretical tools for spinning particle dynamics research
2. Application value: Directly serves scientific objectives of future gravitational wave detectors like LISA
3. Methodological contribution: Analytic methods generalizable to other complex gravitational systems

Applicable Scenarios

1. EMRI system modeling: Particularly suitable for near-equatorial, moderate eccentricity systems
2. Theoretical research: Provides analytic insights for understanding spin-orbit coupling
3. Numerical method validation: Serves as benchmark test for numerical codes

This work holds significant importance in gravitational wave astrophysics, providing powerful analytic tools for understanding and modeling the complex dynamics of spinning binary systems, with important implications for future gravitational wave detection and parameter estimation.