Analytical solutions of bound timelike geodesic orbits in effective-one-body frame

Basic Information

• Paper ID: 2511.06684
• Title: Analytical solutions of bound timelike geodesic orbits in effective-one-body frame
• Authors: Chen Zhang, Wen-Biao Han
• Classification: gr-qc (General Relativity and Quantum Cosmology)
• Publication Date: November 11, 2025
• Paper Link: https://arxiv.org/abs/2511.06684

Abstract

This paper derives approximate analytical solutions for bound timelike geodesic orbits in the effective-one-body (EOB) framework in the extreme mass ratio limit. The analytical solutions use Mino time λ as the independent variable and are expressed in terms of elliptic integrals. Since Mino time decouples the r and θ motions, the authors provide explicit expressions for three orbital frequencies Ωr, Ωθ, Ωϕ (via Fourier series expansion). Using these analytical expressions, one can perform Fourier expansion in Mino time λ for any function expressed in coordinates (r, θ, φ). In particular, the observer time t is decomposed into Mino time λ, and frequency-domain descriptions are constructed from λ-Fourier expansions and expansions of t. These analytical expressions are simple to implement and applicable to gravitational wave calculations for extreme mass ratio inspirals (EMRIs) using frequency-domain Teukolsky equations.

Research Background and Motivation

Research Problem

The core problem addressed in this paper is: How to obtain analytical solutions for bound geodesic orbits in the EOB framework with mass ratio corrections, particularly analytical expressions for orbital frequencies?

Problem Significance

1. Gravitational wave astronomy requirements: Space-based gravitational wave detectors (LISA, Taiji, Tianqin) will observe extreme mass ratio inspiral systems (EMRIs)—the process of stellar-mass compact objects spiraling into supermassive black holes (10⁵-10⁷M⊙)
2. Waveform template construction: Accurate EMRI gravitational wave templates are crucial for data analysis and parameter estimation
3. Computational efficiency: Fast and accurate geodesic calculations form the foundation for studying EMRI gravitational waves and constructing efficient templates

Limitations of Existing Methods

1. Test particle approximation: Traditional Kerr geodesic studies neglect the mass of the small body, unsuitable for real EMRI systems
2. Numerical integration methods: Although the authors' previous work obtained geometric solutions in the EOB framework, the r and θ components remain coupled, requiring numerical integration to compute frequencies with lower efficiency
3. Lack of analytical expressions with mass ratio corrections: This is the first derivation of semi-analytical fundamental frequency expressions including mass ratio corrections

Research Motivation

Combine the Mino time method (decoupling r and θ motions) with EOB formalism (including mass ratio corrections) to derive analytical orbit solutions with high computational efficiency and sufficient accuracy, providing foundational tools for computing EMRI gravitational waves using the Teukolsky equation.

Core Contributions

1. First derivation of analytical geodesic solutions in EOB framework: For the first time, analytical expressions (in elliptic integral form) for bound orbits are obtained in the deformed Kerr metric with mass ratio corrections
2. Analytical frequency expressions: First derivation of semi-analytical fundamental frequency expressions (Ωr, Ωθ, Ωϕ) including mass ratio corrections
3. Fourier expansion framework: Establishes a complete Fourier expansion framework in Mino time λ, enabling frequency-domain decomposition of arbitrary orbital functions fr(t), θ(t)
4. Computational efficiency improvement: The analytical method is at least one order of magnitude faster than numerical integration while maintaining sufficient accuracy (relative error < 10⁻²ν, where ν is the mass ratio)
5. Foundation for EOB-Teukolsky waveform algorithm: These analytical expressions can be directly applied to frequency-domain Teukolsky equations for generating EMRI gravitational wave waveforms

Detailed Methodology

Task Definition

Input: EMRI system parameters

• Black hole mass ratio ν = m₁m₂/M² (extreme mass ratio limit)
• Effective Kerr parameter a = |S_Kerr|/M
• Orbital parameters: semi-latus rectum p, eccentricity e, minimum inclination angle θ_min

Output:

• Orbital solutions r(λ), θ(λ), φ(λ), t(λ) (λ is Mino time)
• Fundamental frequencies Ωr, Ωθ, Ωϕ (relative to observer time t)
• Fourier expansion coefficients

Theoretical Framework

1. EOB Formalism

The EOB method maps the two-body problem to motion of a test particle in an effective metric. The effective metric is a deformed Kerr metric:

$g_{\alpha\beta}^{\text{eff}} P^\alpha P^\beta = \frac{1}{r^2 + a^2\cos^2\theta}\left[\Delta_r P_r^2 + P_\theta^2 + \frac{1}{\sin^2\theta}(P_\phi + a\sin^2\theta P_t)^2 - \frac{1}{\Delta_t}((r^2+a^2)P_t + aP_\phi)^2\right]$

where the metric potential functions include mass ratio corrections (3PN accuracy): $A(u) = 1-2u+2\nu u^3 + \left(\frac{94}{3}-\frac{41\pi^2}{32}\right)\nu u^4$$D^{-1}(u) = 1 + 6\nu u^2 + 2\nu u^3(26-3\nu)$

2. Mino Time Parametrization

The key advantage of Mino time is decoupling r and θ motions: $d\lambda = \frac{g_{\theta\theta}}{E/M(g_{t\phi}\hat{L}_z - g_{tt}\hat{H}_{\text{eff}})}dt$

The geodesic equations become: $\left(\frac{dr}{d\lambda}\right)^2 = (1+6\nu u^2 + 2\nu u^3(26-3\nu))\mathcal{R}(r)$$\left(\frac{d\cos\theta}{d\lambda}\right)^2 = \Theta(\cos\theta)$

Core Method

Step 1: Quartic Polynomial Approximation of Radial Motion

The key innovation is approximating the radial potential function R(r) as a quartic polynomial R̄(r):

$\mathcal{R}(r) \approx \bar{\mathcal{R}}(r) = (1-\hat{H}_{\text{eff}}^2)(r_1-r)(r-r_2)(r-r_3)(r-r_4)$

where r₁ (apocenter) and r₂ (pericenter) are zeros of R(r), and r₃ and r₄ are determined by matching conditions: $\int_{r_2}^{r_1}\frac{dr}{\sqrt{\mathcal{R}(r)}} = \int_{r_2}^{r_1}\frac{dr}{\sqrt{\bar{\mathcal{R}}(r)}}$

This introduces an adjustable parameter C_R to compensate for approximation errors. As shown in Figure 1, even for extreme parameters near the last stable orbit (LSO), the integral error is less than 0.16ν.

Step 2: Elliptic Integral Representation

Using the quartic polynomial form, the radial and polar motion periods can be expressed as elliptic integrals:

Radial period: $\Lambda_r = \frac{4}{\sqrt{(1-\hat{H}_{\text{eff}}^2)(r_1-r_3)(r_2-r_4)}}\left[K(k_r) - 3\nu\frac{Z_2|_{\phi=\pi/2}}{r_2^2} - 26\nu\frac{Z_3|_{\phi=\pi/2}}{r_3^2}\right]$

Polar period: $\Lambda_\theta = \frac{4K(k_\theta)}{\sqrt{\hat{Q}/z_-}}$

where K is the complete elliptic integral of the first kind, and the moduli are: $k_r = \frac{r_1-r_2}{r_1-r_3}\frac{r_3-r_4}{r_2-r_4}, \quad k_\theta = \frac{z_-}{z_+}$

Step 3: Jacobi Elliptic Function Representation of Coordinates

For the test particle limit (ν→0), coordinates can be exactly expressed as Jacobi elliptic functions: $r'(w_r) = \frac{r_3(r_1-r_2)\text{sn}^2(\phi_r, k_r) - r_2(r_1-r_3)}{(r_1-r_2)\text{sn}^2(\phi_r, k_r) - (r_1-r_3)}$$\cos\theta(w_\theta) = \sqrt{z_-}\text{sn}(\phi_\theta, k_\theta)$

For cases with mass ratio corrections, corrections are obtained through second-order Taylor expansion: $F_1 + F_2(y_r - x) + F_3(y_r - x)^2 = 0$

where x = sn(φ_r, k_r) is the test particle solution and y_r is the corrected solution.

Step 4: Fourier Expansion of Frequencies

Time and azimuthal angle evolution can be expanded as Fourier series: $\frac{dt}{d\lambda} = \Gamma + \sum_{k\neq0}T_{k,0}e^{-ikw_r} + \sum_{n\neq0}T_{0,n}e^{-inw_\theta}$$\frac{d\phi}{d\lambda} = \Upsilon_\phi + \sum_{k\neq0}\Phi_{k,0}e^{-ikw_r} + \sum_{n\neq0}\Phi_{0,n}e^{-inw_\theta}$

Fundamental frequency relations: $\Omega_r = \frac{\Upsilon_r}{\Gamma}, \quad \Omega_\theta = \frac{\Upsilon_\theta}{\Gamma}, \quad \Omega_\phi = \frac{\Upsilon_\phi}{\Gamma}$

Technical Innovations

1. Quartic polynomial approximation strategy: Cleverly maintains physical accuracy by preserving orbital turning points (r₁, r₂) and introducing an adjustable parameter C_R while achieving analytical solvability
2. Hybrid analytical-numerical method: Taylor expansion corrections to test particle solutions avoid complete numerical solving, significantly improving efficiency
3. Error control: Systematically analyzes errors introduced by various approximations, demonstrating that for EMRI systems (ν ≪ 1), total errors remain within acceptable ranges
4. Frequency-domain framework construction: Completely establishes Fourier transformation framework from Mino time to observer time, providing direct input for frequency-domain Teukolsky equations

Experimental Setup

Verification Method

The paper validates analytical expressions by comparison with numerical integration methods.

Parameter Space

Seven representative orbital parameter sets are tested, covering different physical regions:

• Effective Kerr parameter: a/M = 0, 0.3, 0.6, 0.9
• Semi-latus rectum: p/M = 6, 8
• Eccentricity: e = 0.6, 0.9
• Minimum inclination angle: θ_min = π/4, π/3
• Mass ratio: ν = 10⁻³

These parameters include:

• Moderate orbits (e=0.6, p=8)
• Extreme orbits near LSO (e=0.9, p=6, a=0.9)
• Various spin and inclination combinations

Evaluation Metrics

Relative frequency difference: $\frac{\Delta\Omega}{\Omega} = \frac{\Omega_{\text{ana}} - \Omega_{\text{num}}}{\Omega_{\text{num}}}$

Normalized frequency shift: $\frac{\Delta\Omega}{\Omega\nu}$

used to assess accuracy of mass ratio corrections.

Comparison Methods

• Baseline method: Numerical integration method from the authors' previous work 46
• Test particle limit: Analytical Kerr geodesic solutions from Fujita and Hikida 15

Implementation Details

• Elliptic integrals computed using standard numerical libraries
• Fourier coefficients calculated via numerical integration (one-dimensional)
• Computational speed at least one order of magnitude faster than numerical methods

Experimental Results

Main Results

Table I presents frequency comparison results for seven orbital parameter sets:

Key Findings:

1. High accuracy: Relative errors ΔΩ/Ω < 10⁻²ν in all cases
2. Mass ratio correction consistency: Normalized frequency shifts ΔΩ/(Ων) are essentially independent of the mass ratio itself, ranging within 10⁻³ to 10⁻²
3. Extreme parameter stability: Even for extreme orbits near LSO (e=0.9, p=6, a=0.9), accuracy is maintained

Error Analysis

Figure 1 shows how integral error varies with radial coordinate:

• Error is maximum near pericenter r₂
• Error asymptotically approaches zero as r increases
• Even for extreme parameters (orange line, r_LSO=2.9M), error < 0.16ν

Computational Efficiency

• Analytical method: Frequency computation at least one order of magnitude faster than numerical integration
• Reason: Avoids iterative solving of coupled differential equations, requiring only elliptic integral calculations and one-dimensional Fourier transforms

Experimental Findings

1. Approximation validity verification: The quartic polynomial approximation R(r)→R̄(r) maintains high accuracy across the entire parameter space
2. Separability of mass ratio corrections: Mass ratio corrections primarily enter through metric potential functions A(u), D(u), and ω̃_fd, allowing systematic incorporation into the analytical framework
3. Completeness of frequency-domain representation: Through Mino time Fourier expansion, orbital dynamics can be completely described, providing foundation for frequency-domain gravitational wave calculations

Related Work

Kerr Geodesic Research

1. Chandrasekhar 11: Classical review systematically summarizing geodesic motion in black hole spacetimes
2. Schmidt 12: Uses Hamilton-Jacobi formalism to derive fundamental frequencies, but does not resolve r-θ coupling
3. Mino 13: Introduces Mino time parameter, first achieving decoupling of r and θ motions
4. Drasco & Hughes 14: Combines Schmidt method with Mino time, deriving Fourier component calculation methods
5. Fujita & Hikida 15: Derives analytical expressions for bound orbits in Kerr spacetime (test particle)

EOB Formalism

1. Buonanno & Damour 16,17: Propose EOB method, mapping two-body problem to effective one-body problem
2. Barausse & Buonanno 50: Construct deformed Kerr metric, making EOB Hamiltonian exactly reproduce dynamics of spinning test particles in Kerr spacetime
3. Extreme mass ratio EOB research 26-38: Nagar et al. systematically study EOB waveforms and fluxes in extreme mass ratio limits

Gravitational Self-Force Research

1. Damour et al. 39-45: Analytically compute gravitational self-force terms, providing high-order corrections to EOB formalism
2. Numerical self-force calculations: Provide benchmarks for EOB parameter calibration

Unique Contributions of This Paper

This paper is the first to combine Mino time methods with EOB formalism, deriving semi-analytical orbit solutions and frequency expressions with mass ratio corrections, filling a gap in this field.

Conclusions and Discussion

Main Conclusions

1. Successful derivation of analytical solutions: For the first time, complete analytical expressions (in elliptic integral form) for bound orbits are obtained in EOB deformed Kerr spacetime
2. Efficient frequency calculation: Analytical expressions for fundamental frequencies Ωr, Ωθ, Ωϕ are at least one order of magnitude faster than numerical methods while meeting EMRI application accuracy requirements
3. Complete Fourier expansion framework: Establishes comprehensive frequency-domain description framework enabling Mino time Fourier expansion of arbitrary orbital functions
4. Mass ratio corrections incorporated: Systematically incorporates mass ratio corrections into analytical framework with errors controlled at O(10⁻²)ν level

Limitations

1. Effective spin of small body neglected: Temporarily omits effective spin of small body (~µa/M), introducing errors at least two orders lower than mass ratio relative to Hamiltonian
2. Approximations introduced:
   • Quartic polynomial approximation R(r)→R̄(r)
   • Approximation of metric denominator terms (∆_t∧_t)
   • 3PN truncation (not using 5PN logarithmic resummation potential)
3. Applicable range: Primarily for extreme mass ratio systems (ν ≪ 1); moderate mass ratios may require higher-order corrections
4. Radiation reaction not included: Currently considers only conservative orbital dynamics, excluding radiation reaction effects

Future Directions

1. Complete EMRI waveform generation: Apply analytical solutions to frequency-domain Teukolsky equations to generate EMRI gravitational wave waveforms with mass ratio corrections
2. EOB-Teukolsky algorithm refinement: Combine orbital evolution and radiation reaction to complete full EOB-Teukolsky EMRI waveform algorithm
3. Higher-order corrections: Incorporate effective spin, 5PN potential functions, and other higher-order corrections
4. Parameter estimation applications: Utilize efficient analytical expressions for LISA data analysis and parameter estimation studies
5. Self-force effects: Combine gravitational self-force calculation results to further improve EOB model accuracy

In-Depth Evaluation

Strengths

1. Important theoretical breakthrough: First analytical orbit solutions with mass ratio corrections in EOB framework represent significant progress in this field
2. Strong methodological innovation:
   • Quartic polynomial approximation strategy is clever and effective
   • Taylor expansion correction method balances accuracy and efficiency
   • Systematic error analysis ensures method reliability
3. High practical value:
   • Computational efficiency improved by over one order of magnitude
   • Directly applicable to frequency-domain Teukolsky equations
   • Provides tools for data analysis by LISA and other detectors
4. Comprehensive verification:
   • Systematic comparison with numerical methods
   • Covers extensive parameter space
   • Tests extreme orbital parameters
5. Clear presentation:
   • Detailed mathematical derivations
   • Clear physical pictures
   • Systematic result presentation

Weaknesses

1. Insufficient systematic evaluation of approximations:
   • While specific error estimates are provided, analysis of interactions between different approximations is limited
   • Lacks systematic study of error variation with parameters
2. Unclear physical meaning of parameter adjustment:
   • Physical significance of adjustable parameter C_R is not sufficiently clear
   • Variation patterns of C_R across different orbital parameters not deeply discussed
3. Limited comparison with other methods:
   • Only compared with authors' previous numerical method
   • Lacks comparison with other semi-analytical methods (if they exist)
4. Limited application demonstration:
   • No actual gravitational waveform calculation shown
   • Lacks direct connection to LISA observational capabilities
5. Insufficient discussion of limiting cases:
   • Behavior near LSO requires more detailed analysis
   • Applicability in high-spin, high-eccentricity limits needs further discussion

Impact

1. Academic contribution:
   • Provides important theoretical tools for EMRI gravitational wave research
   • Advances application of EOB methods in extreme mass ratio systems
   • May become standard method in this field
2. Practical value:
   • Directly serves space-based gravitational wave detectors (LISA, Taiji, Tianqin)
   • Enables fast waveform template generation
   • Improves efficiency of EMRI parameter estimation
3. Reproducibility:
   • Detailed mathematical derivations facilitate implementation
   • Key formulas provided in explicit form
   • Code availability not mentioned
4. Potential for follow-up research:
   • Provides foundation for evolution studies including radiation reaction
   • Extensible to higher-order mass ratio corrections
   • Methods potentially applicable to other spacetimes

Applicable Scenarios

1. Primary applications:
   • Fast EMRI gravitational waveform generation
   • Frequency-domain Teukolsky equation solving
   • Rapid orbital frequency calculation
   • Parameter space exploration and template library construction
2. Applicable conditions:
   • Mass ratio ν ≪ 1 (extreme mass ratio)
   • Bound orbits (excluding scattering orbits)
   • Conservative dynamics (short timescales, neglecting radiation reaction)
   • Non-spinning small body
3. Inapplicable scenarios:
   • Moderate or comparable mass ratio systems
   • Applications requiring extreme precision (< 10⁻³ν)
   • Late inspiral stages (strong radiation reaction)
   • Spinning small bodies with significant spin

Technical Highlights

1. Skillful application of elliptic integral techniques: Transforms complex geodesic equations into standard elliptic integral forms
2. Application of Jacobi elliptic functions: Elegantly represents periodic orbital motion
3. Perturbation expansion strategy: Performs Taylor expansion on test particle solutions, avoiding complexity of solving from scratch
4. Frequency-domain methods: Fully exploits advantages of Mino time, establishing complete frequency-domain framework

Key References

13 Y. Mino, "Perturbative approach to an orbital evolution around a supermassive black hole," Phys. Rev. D 67, 084027 (2003)

14 S. Drasco and S. A. Hughes, "Rotating black hole orbit functionals in the frequency domain," Phys. Rev. D 69, 044015 (2004)

15 R. Fujita, W. Hikida, and H. Tagoshi, "An efficient numerical method for computing gravitational waves induced by a particle moving on eccentric inclined orbits around a Kerr black hole," Prog. Theor. Phys. 121, 843 (2009)

16 A. Buonanno and T. Damour, "Effective one-body approach to general relativistic two-body dynamics," Phys. Rev. D 59, 084006 (1999)

46 C. Zhang, W.-B. Han, X.-Y. Zhong, and G. Wang, "Geometrized effective-one-body formalism for extreme-mass-ratio limits: Generic orbits," Phys. Rev. D 104, 024050 (2021)

50 E. Barausse and A. Buonanno, "Improved effective-one-body hamiltonian for spinning black-hole binaries," Phys. Rev. D 81, 084024 (2010)

Overall Assessment: This is a high-quality theoretical physics paper achieving substantial progress on important problems in gravitational wave astronomy. The methods are innovative, derivations rigorous, and verification comprehensive, possessing significant academic and practical value. While certain limitations exist (such as approximation applicability ranges and lack of actual waveform demonstrations), these do not diminish the overall quality. The paper will likely have positive impact on EMRI gravitational wave research. It is recommended that follow-up work promptly demonstrate actual gravitational waveform calculations and establish more direct connections with observational data analysis requirements.