An open-source numerical tool for rational orbits and gravitational radiation in static spherically symmetric spacetimes
Pith reviewed 2026-06-26 13:59 UTC · model grok-4.3
The pith
A numerical framework computes rational orbits and gravitational radiation in any user-supplied static spherically symmetric spacetime.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
The paper establishes a computational framework that, given only a user-defined covariant metric for a static spherically symmetric spacetime, efficiently calculates rational orbits of timelike particles and the corresponding gravitational wave polarization states and characteristic strains, with validation on the Schwarzschild metric and an application showing detectability of waves from a specific extreme mass ratio inspiral system.
What carries the argument
Parallelized numerical integration and wave extraction routines that accept an arbitrary covariant metric as sole input to compute geodesic orbits and linearized gravitational radiation.
Load-bearing premise
The numerical integration and wave extraction routines correctly reproduce geodesic motion and linearized gravitational radiation for arbitrary user-supplied metrics without introducing integration artifacts or gauge-dependent errors.
What would settle it
A direct comparison of the code's computed orbital frequencies or wave strains against known analytic results for a specific rational orbit in Schwarzschild spacetime; disagreement beyond numerical tolerance would falsify correctness.
read the original abstract
Timelike orbits constitute a crucial probe for exploring the intrinsic properties of curved spacetimes, and the carried gravitational radiation signals provide a direct window into strong field gravity. In this paper, we develop a versatile computational framework based on Mathematica and the OpenMP parallel architecture to simulate the rational orbits of timelike particles and their gravitational radiation in static spherically symmetric spacetimes. Specifically, requiring only the user defined covariant metric, this numerical tool can efficiently calculate rational orbits across various configurations, as well as the corresponding gravitational wave polarization states and characteristic strains. The package presented here offers a highly efficient and comprehensive one-stop solution for investigating the properties of curved spacetimes and their potential observational signatures. To demonstrate the reliability and capability of our code, we apply it to the Schwarzschild spacetime as a test case, illustrating the functionality of the code across several key aspects, including the effective potential, stable orbital regions, rational and irrational orbits, and gravitational wave signals. Furthermore, we show that the gravitational waves emitted by an extreme mass ratio inspiral system composed of an intermediate mass black hole and the Galactic Center supermassive black hole have the potential to be identified by future space detectors.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript presents an open-source numerical framework in Mathematica with OpenMP parallelization to compute rational timelike orbits and associated gravitational wave polarizations and strains in static spherically symmetric spacetimes, requiring only a user-supplied covariant metric. The tool is validated on the Schwarzschild metric through calculations of effective potentials, stable orbital regions, rational and irrational orbits, and GW signals, and is applied to an extreme-mass-ratio inspiral consisting of an intermediate-mass black hole and Sgr A* to argue that the emitted waves have detection potential for future space-based detectors.
Significance. If the numerical implementation is accurate, the open-source release of a flexible, parallelized one-stop tool for arbitrary SSS metrics represents a useful contribution for exploring geodesic motion and linearized gravitational radiation beyond standard cases. The parallel architecture and requirement of only the metric as input are practical strengths that support reproducibility. The EMRI application provides a concrete example of observational relevance for the Galactic Center.
major comments (2)
- [demonstration section on Schwarzschild] Schwarzschild validation (abstract and demonstration section): the description of successful tests on effective potentials, orbits, and GW signals provides no quantitative error metrics, convergence checks, or comparisons to independent codes or analytic results. This is load-bearing for the central claim that the routines correctly handle arbitrary user-defined SSS metrics without integration artifacts.
- [EMRI section] EMRI application (abstract): the claim that the gravitational waves from the IMBH-Sgr A* system have detection potential for future detectors rests on numerical output whose details, error estimates, and direct comparison to detector noise curves are not quantified.
minor comments (1)
- The abstract would benefit from an explicit link to the code repository or GitHub page to facilitate immediate access by readers.
Simulated Author's Rebuttal
We thank the referee for the careful review and constructive feedback on our manuscript. We address each major comment below and commit to revisions that strengthen the quantitative validation.
read point-by-point responses
-
Referee: [demonstration section on Schwarzschild] Schwarzschild validation (abstract and demonstration section): the description of successful tests on effective potentials, orbits, and GW signals provides no quantitative error metrics, convergence checks, or comparisons to independent codes or analytic results. This is load-bearing for the central claim that the routines correctly handle arbitrary user-defined SSS metrics without integration artifacts.
Authors: We acknowledge the absence of explicit quantitative metrics in the current demonstration section. The revised manuscript will add: (i) relative errors in the conserved energy and angular momentum along integrated orbits (targeting <10^{-8} or better), (ii) direct comparisons of orbital frequencies, radial periods, and periastron precession rates against the known analytic Schwarzschild expressions, (iii) convergence tests varying the integration step size and tolerance, and (iv) side-by-side comparison of computed GW polarizations and strains with published results for circular and eccentric Schwarzschild orbits. These additions will be inserted into the demonstration section to support the claim of reliable handling of arbitrary SSS metrics. revision: yes
-
Referee: [EMRI section] EMRI application (abstract): the claim that the gravitational waves from the IMBH-Sgr A* system have detection potential for future detectors rests on numerical output whose details, error estimates, and direct comparison to detector noise curves are not quantified.
Authors: We agree that the EMRI claim requires more quantitative support. The revision will expand the application section to report: the precise masses (IMBH and Sgr A*), orbital parameters (semi-major axis, eccentricity), integration error estimates on the strain time series, the characteristic strain spectrum h_c(f), and direct overlay comparisons against the noise curves of LISA, TianQin, and Taiji (including approximate SNR values where the signal exceeds the noise). This will make the detection-potential statement evidence-based. revision: yes
Circularity Check
No significant circularity
full rationale
The paper describes a numerical implementation (Mathematica + OpenMP) that takes a user-supplied metric as input and computes orbits plus linearized GWs. Validation is performed against known Schwarzschild results; the EMRI detectability statement is a qualitative extrapolation from those outputs. No equation reduces to a fitted parameter renamed as a prediction, no self-citation supplies a load-bearing uniqueness theorem, and no ansatz is smuggled in. The derivation chain is therefore self-contained against external benchmarks.
Axiom & Free-Parameter Ledger
axioms (2)
- standard math Timelike geodesics in a static spherically symmetric metric obey the standard Euler-Lagrange equations derived from the line element.
- domain assumption Gravitational wave polarization states and strains can be extracted from the quadrupole or higher moments of the orbiting particle in the weak-field limit at large distance.
Reference graph
Works this paper leans on
-
[1]
LIGO Scientific Collaboration and Virgo Collaboration,Observation of Gravitational Waves from a Binary Black Hole Merger,Phys. Rev. Lett.116(2016) 061102 [1602.03837]
Pith/arXiv arXiv 2016
-
[2]
Event Horizon Telescope Collaboration,First M87 Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole,Astrophys. J. Lett.875(2019) L1 [1906.11238]
Pith/arXiv arXiv 2019
-
[3]
Event Horizon Telescope Collaboration,First Sagittarius A* Event Horizon Telescope Results. I. The Shadow of the Supermassive Black Hole in the Center of the Milky Way,Astrophys. J. Lett.930(2022) L12
2022
- [4]
-
[5]
M. Takahashi and H. Koyama,Chaotic Motion of Charged Particles in an Electromagnetic Field Surrounding a Rotating Black Hole,Astrophys. J.693(2009) 472 [0807.0277]
Pith/arXiv arXiv 2009
-
[6]
Z. Stuchlík and M. Kološ,Acceleration of the charged particles due to chaotic scattering in the combined black hole gravitational field and asymptotically uniform magnetic field,Eur. Phys. J. C76(2016) 32 [1511.02936]
Pith/arXiv arXiv 2016
-
[7]
R. Pánis, M. Kološ and Z. Stuchlík,Determination of chaotic behaviour in time series generated by charged particle motion around magnetized Schwarzschild black holes,Eur. Phys. J. C79(2019) 479 [1905.01186]
Pith/arXiv arXiv 2019
-
[8]
D. Li and X. Wu,Chaotic motion of neutral and charged particles in a magnetized Ernst-Schwarzschild spacetime,Eur. Phys. J. Plus134(2019) 96 [1803.02119]
Pith/arXiv arXiv 2019
-
[9]
S. Hu, X. Wu and E. Liang,Construction of a Second-order Six-dimensional Hamiltonian-conserving Scheme,Astrophys. J. Suppl. S.257(2021) 40 [2108.09882]
arXiv 2021
-
[10]
W. Sun, Y. Wang, F. Liu and X. Wu,Applying explicit symplectic integrator to study chaos of charged particles around magnetized Kerr black hole,Eur. Phys. J. C81(2021) 785 [2109.02295]
arXiv 2021
-
[11]
Huang, G
Z. Huang, G. Huang and A. Hu,Application of Explicit Symplectic Integrators in a Magnetized Deformed Schwarzschild Black Spacetime,Astrophys. J.925(2022) 158. – 13 –
2022
- [12]
-
[13]
Z. Xu, D. Ma, W. Cao and K. Li,Chaotic motion of charged test particles in a Kerr-MOG black hole with explicit symplectic algorithms,Eur. Phys. J. C85(2025) 770 [2412.06122]
arXiv 2025
-
[14]
O. Kopáček and V. Karas,Inducing Chaos by Breaking Axial Symmetry in a Black Hole Magnetosphere,Astrophys. J.787(2014) 117 [1404.5495]
Pith/arXiv arXiv 2014
-
[15]
T. Jacobson and T.P. Sotiriou,Spinning Black Holes as Particle Accelerators,Phys. Rev. Lett. 104(2010) 021101 [0911.3363]
Pith/arXiv arXiv 2010
-
[16]
S.-W. Wei, Y.-X. Liu, H. Guo and C.-E. Fu,Charged spinning black holes as particle accelerators,Phys. Rev. D82(2010) 103005 [1006.1056]
Pith/arXiv arXiv 2010
-
[17]
M. Guo and S. Gao,Kerr black holes as accelerators of spinning test particles,Phys. Rev. D93 (2016) 084025 [1602.08679]
Pith/arXiv arXiv 2016
-
[18]
J. Levin and G. Perez-Giz,A periodic table for black hole orbits,Phys. Rev. D77(2008) 103005 [0802.0459]
Pith/arXiv arXiv 2008
-
[19]
M. Alloqulov, T. Xamidov, S. Shaymatov and B. Ahmedov,Gravitational waveforms from periodic orbits around a Schwarzschild black hole embedded in a Dehnen-type dark matter halo, Eur. Phys. J. C85(2025) 798 [2504.05236]
arXiv 2025
-
[20]
V. Misra and J. Levin,Rational Orbits around Charged Black Holes,Phys. Rev. D82(2010) 083001 [1007.2699]
Pith/arXiv arXiv 2010
-
[21]
S.-W. Wei, J. Yang and Y.-X. Liu,Geodesics and periodic orbits in Kehagias-Sfetsos black holes in deformed Hořava-Lifshitz gravity,Phys. Rev. D99(2019) 104016 [1904.03129]
Pith/arXiv arXiv 2019
-
[22]
C. Liu, C. Ding and J. Jing,Periodic orbits around Kerr Sen black holes,Commun. Theor. Phys.71(2019) 1461 [1804.05883]
arXiv 2019
-
[23]
Deng,Geodesics and periodic orbits around quantum-corrected black holes,Phys
X.-M. Deng,Geodesics and periodic orbits around quantum-corrected black holes,Phys. Dark Univ.30(2020) 100629
2020
-
[24]
Deng,Periodic orbits around brane-world black holes,Eur
X.-M. Deng,Periodic orbits around brane-world black holes,Eur. Phys. J. C80(2020) 489
2020
-
[25]
Lin and X.-M
H.-Y. Lin and X.-M. Deng,Rational orbits around 4DEinstein–Lovelock black holes,Phys. Dark Univ.31(2021) 100745
2021
-
[26]
R. Wang, F. Gao and H. Chen,Periodic orbits around a static spherically symmetric black hole surrounded by quintessence,Ann. Phys.447(2022) 169167
2022
-
[27]
Gao and X.-M
B. Gao and X.-M. Deng,Bound orbits around Bardeen black holes,Ann. Phys.418(2020) 168194
2020
-
[28]
Lin and X.-M
H.-Y. Lin and X.-M. Deng,Precessing and periodic orbits around Lee–Wick black holes,Eur. Phys. J. Plus137(2022) 176
2022
-
[29]
Zhang and Y
J. Zhang and Y. Xie,Probing a black-bounce-Reissner–Nordström spacetime with precessing and periodic motion,Eur. Phys. J. C82(2022) 854
2022
-
[30]
Gao and X.-M
B. Gao and X.-M. Deng,Bound orbits around modified Hayward black holes,Mod. Phys. Lett. A36(2021) 2150237
2021
-
[31]
Lin and X.-M
H.-Y. Lin and X.-M. Deng,Bound Orbits and Epicyclic Motions around Renormalization Group Improved Schwarzschild Black Holes,Universe8(2022) 278
2022
-
[32]
Lin and X.-M
H.-Y. Lin and X.-M. Deng,Precessing and periodic orbits around hairy black holes in Horndeski’s Theory,Eur. Phys. J. C83(2023) 311. – 14 –
2023
- [33]
-
[34]
Y.-K. Lim and Z.C. Yeo,Energies and angular momenta of periodic Schwarzschild geodesics, Phys. Rev. D109(2024) 024037 [2401.13894]
arXiv 2024
-
[35]
A. Bragado and G.J. Olmo,Periodic equatorial orbits in a black bounce scenario,Class. Quant. Grav.42(2025) 185004 [2506.15600]
arXiv 2025
-
[36]
S. Shaymatov, B. Ahmedov, M. De Laurentis, M. Jamil, Q. Wu, A. Wang et al.,On the Parameters of the Spherically Symmetric Parameterized Rezzolla–Zhidenko Spacetime through Solar System Tests, the Orbit of the S2 Star about Sgr A*, and Quasiperiodic Oscillations, Astrophys. J.959(2023) 6 [2307.10804]
arXiv 2023
-
[37]
T. Xamidov, S. Shaymatov, B. Ahmedov and T. Zhu,Probing quantum corrected black hole through astrophysical tests with the orbit of S2 star and quasiperiodic oscillations,J. Cosmol. Astropart. Phys.01(2026) 044 [2503.06750]
arXiv 2026
-
[38]
Z.-Y. Tu, T. Zhu and A. Wang,Periodic orbits and their gravitational wave radiations in a polymer black hole in loop quantum gravity,Phys. Rev. D108(2023) 024035 [2304.14160]
arXiv 2023
-
[39]
L. Meng, Z. Xu and M. Tang,Bound orbits and gravitational wave radiation around the hairy black hole,Eur. Phys. J. C85(2025) 306 [2411.01858]
arXiv 2025
- [40]
-
[41]
C.-H. Wang, X.-C. Meng, Y.-P. Zhang, T. Zhu and S.-W. Wei,Equatorial periodic orbits and gravitational waveforms in a black hole free of Cauchy horizon,J. Cosmol. Astropart. Phys.07 (2025) 021 [2502.08994]
arXiv 2025
-
[42]
S. Lu, H.-J. Lin, T. Zhu, Y.-X. Liu and X. Zhang,Gravitational radiations from periodic orbits around a black hole in the effective field theory extension of general relativity,Eur. Phys. J. C 86(2026) 283 [2512.11911]
Pith/arXiv arXiv 2026
-
[43]
S. Lu and T. Zhu,Gravitational radiations from periodic orbits around Einstein-Æther black holes,Phys. Dark Univ.50(2025) 102141 [2505.00294]
Pith/arXiv arXiv 2025
- [44]
- [45]
-
[46]
W. Deng, S. Long, Q. Tan and J. Jing,Gravitational waveforms from periodic orbits around a charged black hole with scalar hair*,Chin. Phys. C50(2026) 035103 [2510.24468]
arXiv 2026
-
[47]
M. Alloqulov, S. Shaymatov, B. Ahmedov and T. Zhu,Gravitational waveforms from periodic orbits around a dyonic ModMax black hole,Eur. Phys. J. C86(2026) 259 [2511.15237]
arXiv 2026
-
[48]
S. Yang, Y.-P. Zhang, T. Zhu, L. Zhao and Y.-X. Liu,Gravitational waveforms from periodic orbits around a quantum-corrected black hole,J. Cosmol. Astropart. Phys.01(2025) 091 [2407.00283]
arXiv 2025
-
[49]
M. Alloqulov, S. Shaymatov, B. Ahmedov and T. Zhu,Regular black hole’s impact on the gravitational waveforms from periodic orbits,Eur. Phys. J. C86(2026) 117 [2508.05245]
arXiv 2026
-
[50]
L. Zhao, M. Tang and Z. Xu,Periodic orbits and gravitational wave radiation in short hair black hole spacetimes for an extreme mass ratio system,Eur. Phys. J. C85(2025) 36 [2411.01979]. – 15 –
arXiv 2025
-
[51]
O. Shabbir, M. Jamil and M. Azreg-Aïnou,Periodic orbits and their gravitational wave radiations around the Schwarzschild-MOG black hole,Phys. Dark Univ.47(2025) 101816 [2501.04367]
arXiv 2025
-
[52]
S. Haroon and T. Zhu,Periodic orbits and their gravitational wave radiations in a black hole with a dark matter halo,Phys. Rev. D112(2025) 044046 [2502.09171]
arXiv 2025
-
[53]
C. Zhang and T. Zhu,Periodic orbits and their gravitational wave radiations inγ-metric,Eur. Phys. J. C86(2026) 350 [2511.14080]
Pith/arXiv arXiv 2026
-
[54]
E.L.B. Junior, J.T.S.S. Junior, F.S.N. Lobo, M.E. Rodrigues, D. Rubiera-Garcia, L.F.D. da Silva et al.,Periodical orbits and waveforms with spontaneous Lorentz symmetry-breaking in Kalb–Ramond gravity,Eur. Phys. J. C85(2025) 557 [2412.00769]
arXiv 2025
-
[55]
S. Zare, T. Zhu, L.M. Nieto, S. Lu and H. Hassanabadi,Probing regular black holes with sub-Planckian curvature through periodic orbits and their gravitational wave radiation,J. Cosmol. Astropart. Phys.01(2026) 059 [2510.05166]
arXiv 2026
-
[56]
Huang,Probing holonomy corrected Schwarzschild black holes with precessing and periodic orbits,Phys
L. Huang,Probing holonomy corrected Schwarzschild black holes with precessing and periodic orbits,Phys. Rev. D111(2025) 084038
2025
-
[57]
O. Shabbir, A. Shermatov, B. Majeed, T. Zahra, M. Jamil and J. Rayimbaev,Probing the nature of Einstein-nonlinear Maxwell-Yukawa black hole through gravitational wave forms from periodic orbits and quasiperiodic oscillations,Nucl. Phys. B1026(2026) 117432 [2601.02904]
Pith/arXiv arXiv 2026
-
[58]
Z. Hua, Z.-T. He, J.-Q. Lai, J. Jiao and Y. Tian,Taxonomy of periodic orbits and gravitational waves in a non-rotating Destounis-Suvorov-Kokkotas black hole spacetime,Phys. Lett. B876 (2026) 140402 [2601.00550]
Pith/arXiv arXiv 2026
-
[59]
Y.-Z. Li and X.-M. Kuang,The bound orbits and gravitational waveforms of timelike particles around renormalization group improved Kerr black holes,Eur. Phys. J. C86(2026) 261 [2509.07333]
arXiv 2026
-
[60]
Huang and X.-M
L. Huang and X.-M. Deng,Can a particle’s motion distinguish scale-dependent Planck stars from renormalization group improved Schwarzschild black holes?,Phys. Rev. D109(2024) 124005
2024
-
[61]
L. Huang,Distinguishing scale-dependent Planck stars from renormalization group improved Schwarzschild black holes by gravitational waves,J. Cosmol. Astropart. Phys.10(2025) 067 [2506.17909]
arXiv 2025
-
[62]
Jiang, M
H. Jiang, M. Alloqulov, Q. Wu, S. Shaymatov and T. Zhu,Periodic orbits and plasma effects on gravitational weak lensing by self-dual black hole in loop quantum gravity,Phys. Dark Univ. 46(2024) 101627
2024
- [63]
-
[64]
A. Kumar, A. Mohammadi and S.G. Ghosh,Probing Gravitational Wave Signatures from Periodic Orbits of Regular Black Holes in Asymptotically Safe Gravity,2605.07187
-
[65]
T. Robson, N.J. Cornish and C. Liu,The construction and use of LISA sensitivity curves, Class. Quant. Grav.36(2019) 105011 [1803.01944]. [66]LISA Sciencecollaboration,Science of the LISA mission: A Summary for the European Strategy for Particle Physics,2507.05130
Pith/arXiv arXiv 2019
-
[66]
J. Luo, L.-S. Chen, H.-Z. Duan, Y.-G. Gong, S. Hu, J. Ji et al.,TianQin: a space-borne gravitational wave detector,Classical and Quantum Gravity33(2016) 035010 [1512.02076]
Pith/arXiv arXiv 2016
- [67]
-
[68]
Li et al.,Gravitational wave astronomy with TianQin,Rept
E.-K. Li et al.,Gravitational wave astronomy with TianQin,Rept. Prog. Phys.88(2025) 056901 [2409.19665]
arXiv 2025
-
[69]
W.-H. Ruan, Z.-K. Guo, R.-G. Cai and Y.-Z. Zhang,Taiji program: Gravitational-wave sources,Int. J. Mod. Phys. A35(2020) 2050075 [1807.09495]
arXiv 2020
-
[70]
Liu et al.,Recent Development of the Laser Interferometer for Taiji Space Gravitational Wave Detection,Research9(2026) 1252
H. Liu et al.,Recent Development of the Laser Interferometer for Taiji Space Gravitational Wave Detection,Research9(2026) 1252
2026
-
[71]
S. Babak, H. Fang, J.R. Gair, K. Glampedakis and S.A. Hughes,“Kludge” gravitational waveforms for a test-body orbiting a Kerr black hole,Phys. Rev. D75(2007) 024005 [gr-qc/0607007]
Pith/arXiv arXiv 2007
-
[72]
L.S. Finn and K.S. Thorne,Gravitational waves from a compact star in a circular, inspiral orbit, in the equatorial plane of a massive, spinning black hole, as observed by LISA,Phys. Rev. D62(2000) 124021 [gr-qc/0007074]
Pith/arXiv arXiv 2000
-
[73]
S. Tan, C. Jiang, D. Li, S. Hu, C. Deng and W. Lin,Gravitational emissions and light curves of quasi-periodic orbits in Schwarzschild spacetime embedded in a Dehnen-type dark matter halo, 2604.13832. – 17 –
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.