Efficient and stable computation of gravitational-wave fluxes from generic Kerr orbits via a unified Heun-function framework
Abstract
Modeling extreme-mass-ratio inspirals hinges on the accurate and efficient computation of gravitational-wave fluxes from generic Kerr orbits. Conventional frequency-domain techniques are often limited by costly auxiliary parameter searches and numerical instabilities in the strong-field or high-frequency regimes. We address these challenges by reformulating both the angular and radial Teukolsky equations in terms of confluent Heun functions. Employing a hybrid analytic continuation algorithm to compute the connection coefficients eliminates the dependence on auxiliary parameters, directly yielding globally convergent solutions and scattering amplitudes. To resolve the highly oscillatory source integrands for generic orbits, we implement an adaptive bi-power mapping quadrature. Comprehensive benchmarks under standard double-precision arithmetic demonstrate that, for the total radiative flux summed over 168 low-order modes, our method achieves relative errors of order 10−11, with computational costs typically reduced by factors of 3--13 compared to the state-of-the-art GeneralizedSasakiNakamura. jl and pybhpt packages. Notably, for highly oscillatory high-order modes, our framework achieves a speedup of up to 60 times compared to specialized oscillatory integrators like GeneralizedSasakiNakamura. jl. These demonstrated gains in precision and efficiency establish the framework as a robust tool for strong-field perturbation theory, providing the numerical foundation for high-order self-force calculations and rapid, high-precision waveform generation.
Versions
Related
- Gravitational waves from core collapse of rotating very-massive stars: 3D numerical relativity computation
- Enabling gravitational-wave astronomy with spin-precessing black holes on generic orbits
- Efficient quantum Gibbs sampling of stabilizer codes using hybrid computation
- Efficient Computation of QKD Key Rates without Semidefinite Programming
- Efficient computation of the asymptotics of extensive-rank HCIZ integrals
BibTeX
@misc{chen2026efficient,
title = {{Efficient and stable computation of gravitational-wave fluxes from generic Kerr orbits via a unified Heun-function framework}},
author = {Changkai Chen and Zhoujian Cao and Jiliang Jing},
howpublished = {arXiv},
year = {2026},
month = sep,
doi = {10.48550/arxiv.2605.09250},
eprint = {2605.09250},
archivePrefix = {arXiv},
url = {https://arxiv.org/abs/2605.09250},
note = {CC BY 4.0},
}
RIS
TY - UNPB
TI - Efficient and stable computation of gravitational-wave fluxes from generic Kerr orbits via a unified Heun-function framework
AU - Changkai Chen
AU - Zhoujian Cao
AU - Jiliang Jing
PY - 2026
DA - 2026/09/04
PB - arXiv
DO - 10.48550/arxiv.2605.09250
UR - https://arxiv.org/abs/2605.09250
AB - Modeling extreme-mass-ratio inspirals hinges on the accurate and efficient computation of gravitational-wave fluxes from generic Kerr orbits. Conventional frequency-domain techniques are often limited by costly auxiliary parameter searches and numerical instabilities in the strong-field or high-frequency regimes. We address these challenges by reformulating both the angular and radial Teukolsky equations in terms of confluent Heun functions. Employing a hybrid analytic continuation algorithm to compute the connection coefficients eliminates the dependence on auxiliary parameters, directly yielding globally convergent solutions and scattering amplitudes. To resolve the highly oscillatory source integrands for generic orbits, we implement an adaptive bi-power mapping quadrature. Comprehensive benchmarks under standard double-precision arithmetic demonstrate that, for the total radiative flux summed over 168 low-order modes, our method achieves relative errors of order 10^−11, with computational costs typically reduced by factors of 3--13 compared to the state-of-the-art GeneralizedSasakiNakamura. jl and pybhpt packages. Notably, for highly oscillatory high-order modes, our framework achieves a speedup of up to 60 times compared to specialized oscillatory integrators like GeneralizedSasakiNakamura. jl. These demonstrated gains in precision and efficiency establish the framework as a robust tool for strong-field perturbation theory, providing the numerical foundation for high-order self-force calculations and rapid, high-precision waveform generation.
ER -