pith. sign in

arxiv: 2606.26339 · v1 · pith:Z4CIZ2XYnew · submitted 2026-06-24 · 🌀 gr-qc · astro-ph.HE

Time-domain framework for the Teukolsky equation with a particle source using comoving hyperboloidal coordinates

Aditya Vaswani , Leor Barack , Oliver Long , Rodrigo Panosso Macedo This is my paper

Pith reviewed 2026-06-26 01:13 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HE
keywords Teukolsky equationtime domainhyperboloidal coordinatesself-forceextreme mass ratioblack hole perturbationnumerical methodspoint particle source
0
0 comments X

The pith

Comoving hyperboloidal coordinates enable stable time-domain solutions of the Teukolsky equation with particle sources.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper introduces a numerical method for integrating the Teukolsky equation in time domain when a point particle is present as a source. By using comoving hyperboloidal coordinates that are compactified in space, the scheme avoids the unphysical growing modes that appear in other approaches. The comoving nature makes it simpler to apply the necessary jump conditions across the particle's path. Tests for scalar fields on circular and scattering orbits, and for the gravitational perturbation case, show long-term stability. This opens the door to computing the complete gravitational self-force for extreme mass ratio scattering scenarios.

Core claim

The scheme based on comoving, spatially compactified hyperboloidal coordinates for the 1+1 dimensional Teukolsky equation with a point-particle source evades nonphysical growing modes and simplifies jump conditions on the particle worldline, as shown in scalar field tests for circular and scattering geodesics and in a stable implementation for the s = -2 case.

What carries the argument

Comoving hyperboloidal coordinates that are spatially compactified, which compactify the spatial domain and align with the particle motion to simplify jumps and ensure stability.

If this is right

  • The method allows long-term stable evolutions without growing modes for the gravitational Teukolsky equation.
  • Jump conditions at the particle location are simplified by the comoving frame.
  • It provides a path to calculating the full gravitational self-force in extreme-mass-ratio scattering.
  • Performance is demonstrated for both circular and scattering geodesic sources in the scalar case.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • This approach may generalize to other perturbation equations in black hole spacetimes.
  • It could facilitate waveform modeling for extreme mass ratio events in gravitational wave astronomy.
  • Comparisons with noncompactified methods highlight the role of compactification in stability.

Load-bearing premise

The assumption that comoving hyperboloidal coordinates can be implemented to both simplify jump conditions and maintain long-term stability for the s = -2 Teukolsky equation without introducing other numerical artifacts.

What would settle it

A long-term numerical evolution of the s = -2 Teukolsky equation with a scattering particle source that develops growing nonphysical modes would falsify the stability claim.

Figures

Figures reproduced from arXiv: 2606.26339 by Aditya Vaswani, Leor Barack, Oliver Long, Rodrigo Panosso Macedo.

Figure 1
Figure 1. Figure 1: FIG. 1. A sample scattering orbit with [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2. An illustration of our hyperboloidal coordinate grid, [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4. Same as in Fig. 3, but now using a scattering geodesic [PITH_FULL_IMAGE:figures/full_fig_p009_4.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3. Decay of the [PITH_FULL_IMAGE:figures/full_fig_p009_3.png] view at source ↗
Figure 5
Figure 5. Figure 5: FIG. 5. Magnitude of the scalar field along the particle’s orbit [PITH_FULL_IMAGE:figures/full_fig_p010_5.png] view at source ↗
Figure 1
Figure 1. Figure 1: Here we cannot leverage any simplifications when [PITH_FULL_IMAGE:figures/full_fig_p010_1.png] view at source ↗
Figure 6
Figure 6. Figure 6: FIG. 6. The [PITH_FULL_IMAGE:figures/full_fig_p011_6.png] view at source ↗
Figure 8
Figure 8. Figure 8: FIG. 8. Relative difference between the hyperboloidal and uv [PITH_FULL_IMAGE:figures/full_fig_p012_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: FIG. 9. The [PITH_FULL_IMAGE:figures/full_fig_p012_9.png] view at source ↗
Figure 10
Figure 10. Figure 10: FIG. 10. Long-term stability of the hyperboloidal scheme (in [PITH_FULL_IMAGE:figures/full_fig_p014_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: FIG. 11. Same as Fig. 10, here for the [PITH_FULL_IMAGE:figures/full_fig_p015_11.png] view at source ↗
Figure 12
Figure 12. Figure 12: FIG. 12. The [PITH_FULL_IMAGE:figures/full_fig_p015_12.png] view at source ↗
read the original abstract

We present a scheme and implementation code for time-domain integration of the Teukolsky equation in 1+1 dimensions with a point-particle source, based on comoving, spatially compactified hyperboloidal coordinates. We demonstrate that the scheme evades the problem of nonphysical growing modes that plague some numerical evolution schemes without compactification. Our use of comoving coordinates greatly simplifies the application of jump conditions on the particle's worldline. We develop our method and test its performance for a scalar field on a Schwarzschild background, first for a circular geodesic orbit source and then for a scattering geodesic orbit. We then present a test implementation of the method for the $s = -2$ Teukolsky equation, illustrating its long-term stability and absence of growing-mode behavior through comparison with similar results using noncompactified characteristic coordinates. Our method paves the way to calculations of the full gravitational self-force in extreme-mass-ratio scattering.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

0 major / 2 minor

Summary. The manuscript presents a time-domain numerical scheme for integrating the Teukolsky equation sourced by a point particle in 1+1 dimensions, using comoving spatially compactified hyperboloidal coordinates. The method is developed and tested first for a scalar field on Schwarzschild with circular and scattering geodesic sources, then extended to the s = -2 gravitational case. It claims to evade nonphysical growing modes that affect some non-compactified schemes, to simplify jump conditions at the particle worldline via the comoving choice, and to demonstrate long-term stability through explicit comparisons, thereby enabling future gravitational self-force calculations for extreme-mass-ratio scattering.

Significance. If the reported stability and absence of growing modes hold under the presented tests, the work provides a concrete, implementable framework that addresses a known obstacle in time-domain self-force computations for scattering trajectories. The explicit numerical demonstrations for both scalar and s = -2 cases, together with the code implementation, constitute a practical contribution that could be directly extended to full gravitational self-force calculations.

minor comments (2)
  1. [Section describing jump conditions] The manuscript states that the comoving coordinates 'greatly simplify' the jump conditions, but does not quantify the reduction in complexity (e.g., number of non-zero jump terms before and after the coordinate choice) in the relevant section describing the source implementation.
  2. [Figures showing time-domain evolution] Figure captions for the stability comparisons should explicitly state the evolution duration in units of M and the grid resolution used, to allow direct assessment of the 'long-term' claim.

Simulated Author's Rebuttal

0 responses · 0 unresolved

We thank the referee for their positive summary, significance assessment, and recommendation of minor revision. No specific major comments were provided in the report, so we have no individual points requiring point-by-point response or manuscript changes at this stage. We remain available for any additional clarifications or minor edits the editor may request.

Circularity Check

0 steps flagged

No significant circularity in numerical method development

full rationale

The paper presents a numerical time-domain integration scheme for the Teukolsky equation using comoving hyperboloidal coordinates, with direct tests for scalar (circular/scattering geodesics) and s=-2 cases demonstrating stability and absence of growing modes via comparison to non-compactified methods. No load-bearing derivation steps reduce by construction to fitted parameters, self-definitions, or self-citation chains; claims rest on explicit implementation and numerical verification rather than tautological reductions.

Axiom & Free-Parameter Ledger

0 free parameters · 1 axioms · 0 invented entities

The central claim rests on standard properties of hyperboloidal foliations and coordinate transformations in general relativity, plus the numerical stability of the chosen discretization; no free parameters or invented entities are indicated in the abstract.

axioms (1)
  • domain assumption Hyperboloidal coordinates allow spatial compactification at future null infinity without introducing instabilities for wave equations on black hole backgrounds.
    Invoked to explain evasion of growing modes.

pith-pipeline@v0.9.1-grok · 5706 in / 1258 out tokens · 22580 ms · 2026-06-26T01:13:48.819795+00:00 · methodology

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

57 extracted references · 2 canonical work pages

  1. [1]

    height” and “shift

    The total scattering angle is given by ∆φp =φ p(χ∞)−φ p(−χ∞)−π.(14) A section of a sample scattering orbit with(v∞, b) = (0.2,21M)is depicted in Figure 1. This will be our reference scattering geodesic for various numerical tests in this work. 4 FIG. 1. A sample scattering orbit with(v∞, b) = (0.2,21M). The black disc represents the central Schwarzschild ...

    1944

  2. [2]

    comoving

    has been a roadblock in the program to calculate the gravitational self-force in extreme-mass-ratio scattering. In this work we have confirmed the expectation, expressed in Ref. [19], that the problem can be effectively mitigated through the use of spatially compactified hyperboloidal coordinates. This we have illustrated here with an explicit numerical i...

    arXiv

  3. [3]

    Damour, Phys

    T. Damour, Phys. Rev. D97, 044038 (2018)

    2018

  4. [4]

    Damour, Phys

    T. Damour, Phys. Rev. D102, 024060 (2020)

    2020

  5. [5]

    Damour, Phys

    T. Damour, Phys. Rev. D94, 104015 (2016)

    2016

  6. [6]

    Kälin and R

    G. Kälin and R. Porto, J. High Energ. Phys. (2020), https://doi.org/10.1007/JHEP01(2020)072

    work page doi:10.1007/jhep01(2020)072 2020

  7. [7]

    G.KälinandR.A.Porto,JournalofHighEnergyPhysics 2020, 120 (2020), arXiv:1911.09130 [hep-th]

    arXiv 2020

  8. [8]

    D. A. Kosower, B. Maybee, and D. O’Connell, Journal of High Energy Physics2019, 137 (2019)

    2019

  9. [9]

    Z. Bern, A. Luna, R. Roiban, C.-H. Shen, and M. Zeng, Phys. Rev. D104, 065014 (2021)

    2021

  10. [10]

    Cheung, I

    C. Cheung, I. Z. Rothstein, and M. P. Solon, Phys. Rev. Lett.121, 251101 (2018)

    2018

  11. [11]

    Barack and A

    L. Barack and A. Pound, Rep. Prog. Phys (2019), 10.1088/1361-6633/aae552

    work page doi:10.1088/1361-6633/aae552 2019

  12. [12]

    Black hole perturbation theory and gravitational self-force,

    A. Pound and B. Wardell, “Black hole perturbation theory and gravitational self-force,” inHandbook of Gravitational Wave Astronomy, edited by C. Bambi, S. Katsanevas, and K. D. Kokkotas (Springer Singapore, Singapore, 2020) pp. 1–119

    2020

  13. [13]

    Warburton, (2025), arXiv:2512.02274 [gr-qc]

    N. Warburton, (2025), arXiv:2512.02274 [gr-qc]

    Pith/arXiv arXiv 2025

  14. [14]

    Hopper, Phys

    S. Hopper, Phys. Rev. D97, 064007 (2018)

    2018

  15. [15]

    Hopper and V

    S. Hopper and V. Cardoso, Phys. Rev. D97, 044031 (2018)

    2018

  16. [16]

    Barack and O

    L. Barack and O. Long, Physical Review D106, 104031 (2022), arXiv:2209.03740 [astro-ph, physics:gr-qc]

    arXiv 2022

  17. [17]

    Barack, Z

    L. Barack, Z. Bern, E. Herrmann, O. Long, J. Parra- Martinez, R. Roiban, M. S. Ruf, C.-H. Shen, M. P. Solon, F. Teng, and M. Zeng, Physical Review D108, 024025 (2023), arXiv:2304.09200 [hep-th]

    arXiv 2023

  18. [18]

    Whittall and L

    C. Whittall and L. Barack, Phys. Rev. D108, 064017 (2023), arXiv:2305.09724 [gr-qc]

    arXiv 2023

  19. [19]

    Whittall, L

    C. Whittall, L. Barack, and O. Long, Phys. Rev. D112, 124045 (2025), arXiv:2509.19439 [gr-qc]

    arXiv 2025

  20. [20]

    O. Long, C. Whittall, and L. Barack, Phys. Rev. D110, 044039 (2024), arXiv:2406.08363 [gr-qc]

    arXiv 2024

  21. [21]

    O.LongandL.Barack,Phys.Rev.D104,024014(2021), arXiv:2105.05630 [gr-qc]. 18

    arXiv 2021

  22. [22]

    Barack and N

    L. Barack and N. Sago, Phys. Rev. D75, 064021 (2007), arXiv:gr-qc/0701069

    Pith/arXiv arXiv 2007

  23. [23]

    Barack and N

    L. Barack and N. Sago, Phys. Rev. D81, 084021 (2010), arXiv:1002.2386 [gr-qc]

    Pith/arXiv arXiv 2010

  24. [24]

    O’Toole, A

    C. O’Toole, A. Ottewill, and B. Wardell, Phys. Rev. D 103, 124022 (2021)

    2021

  25. [25]

    O’Toole,Green Function Methods in Black Hole Spacetimes, Ph.D

    C. O’Toole,Green Function Methods in Black Hole Spacetimes, Ph.D. thesis, University College Dublin, School of Mathematical Sciences (2022)

    2022

  26. [26]

    Vishal, Private communication

    M. Vishal, Private communication

  27. [27]

    S. D. Bishoyi, S. E. Field, and S. R. Lau, (2026), arXiv:2604.22734 [gr-qc]

    Pith/arXiv arXiv 2026

  28. [28]

    Zenginoglu, D

    A. Zenginoglu, D. Nunez, and S. Husa, Class. Quant. Grav.26, 035009 (2009), arXiv:0810.1929 [gr-qc]

    Pith/arXiv arXiv 2009

  29. [29]

    Harms, S

    E. Harms, S. Bernuzzi, and B. Brügmann, Class. Quant. Grav.30, 115013 (2013), arXiv:1301.1591 [gr-qc]

    Pith/arXiv arXiv 2013

  30. [30]

    Panosso Macedo and M

    R. Panosso Macedo and M. Ansorg, J. Comput. Phys. 276, 357 (2014), arXiv:1402.7343 [physics.comp-ph]

    Pith/arXiv arXiv 2014

  31. [31]

    Csukás, I

    K. Csukás, I. Rácz, and G. Z. Tóth, Phys. Rev. D100, 104025 (2019), arXiv:1905.09082 [gr-qc]

    arXiv 2019

  32. [32]

    Csukás and I

    K. Csukás and I. Rácz, Phys. Rev. D103, 084035 (2021), arXiv:2101.05530 [gr-qc]

    arXiv 2021

  33. [33]

    Harms, S

    E. Harms, S. Bernuzzi, A. Nagar, and A. Zenginoglu, Class. Quant. Grav.31, 245004 (2014), arXiv:1406.5983 [gr-qc]

    Pith/arXiv arXiv 2014

  34. [34]

    M. F. O’Boyle and C. Markakis, (2023), arXiv:2308.02385 [math.NA]

    arXiv 2023

  35. [35]

    L. J. G. Da Silva, R. Panosso Macedo, J. E. Thompson, J. A. V. Kroon, L. Durkan, and O. Long, (2023), arXiv:2306.13153 [gr-qc]

    arXiv 2023

  36. [36]

    Vishal, S

    M. Vishal, S. E. Field, K. Rink, S. Gottlieb, and G. Khanna, Phys. Rev. D110, 104009 (2024), arXiv:2307.01349 [gr-qc]

    arXiv 2024

  37. [37]

    L. J. G. Da Silva, (2024), arXiv:2401.08758 [gr-qc]

    arXiv 2024

  38. [38]

    Vishal, S

    M. Vishal, S. E. Field, S. Gottlieb, and J. Ryan, Gen. Rel. Grav.57, 104 (2025), arXiv:2503.11523 [gr-qc]

    arXiv 2025

  39. [39]

    A. Roy, L. Küchler, A. Pound, and R. Panosso Macedo, (2025), arXiv:2510.11793 [gr-qc]

    Pith/arXiv arXiv 2025

  40. [40]

    Panosso Macedo, B

    R. Panosso Macedo, B. Leather, N. Warburton, B. Wardell, and A. Zenginoğlu, Phys. Rev. D105, 104033 (2022), arXiv:2202.01794 [gr-qc]

    arXiv 2022

  41. [41]

    Leather, Gen

    B. Leather, Gen. Rel. Grav.57, 112 (2025), arXiv:2411.14976 [gr-qc]

    arXiv 2025

  42. [42]

    Zenginoğlu, Gen

    A. Zenginoğlu, Gen. Rel. Grav.57, 75 (2025), arXiv:2502.08581 [gr-qc]

    arXiv 2025

  43. [43]

    Pound, C

    A. Pound, C. Merlin, and L. Barack, Phys. Rev. D89, 024009 (2014), arXiv:1310.1513 [gr-qc]

    Pith/arXiv arXiv 2014

  44. [44]

    Zenginoglu, Class

    A. Zenginoglu, Class. Quant. Grav.25, 145002 (2008), arXiv:0712.4333 [gr-qc]

    Pith/arXiv arXiv 2008

  45. [45]

    Panosso Macedo, Phil

    R. Panosso Macedo, Phil. Trans. Roy. Soc. Lond. A382, 20230046 (2024), arXiv:2307.15735 [gr-qc]

    arXiv 2024

  46. [46]

    Ansorg and R

    M. Ansorg and R. Panosso Macedo, Phys. Rev. D93, 124016 (2016), arXiv:1604.02261 [gr-qc]

    Pith/arXiv arXiv 2016

  47. [47]

    Panosso Macedo, J

    R. Panosso Macedo, J. L. Jaramillo, and M. Ansorg, Phys. Rev. D98, 124005 (2018), arXiv:1809.02837 [gr- qc]

    Pith/arXiv arXiv 2018

  48. [48]

    Panosso Macedo, Class

    R. Panosso Macedo, Class. Quant. Grav.37, 065019 (2020), arXiv:1910.13452 [gr-qc]

    arXiv 2020

  49. [49]

    W. H. Press, S. A. Teukolsky, W. T. Vetterling, and B. P. Flannery,Numerical recipes 3rd edition(Cambridge University Press, 2007)

    2007

  50. [50]

    R. H. Price, Phys. Rev. D5, 2419 (1972)

    1972

  51. [51]

    R. H. Price, Phys. Rev. D5, 2439 (1972)

    1972

  52. [52]

    Thornburg, B

    J. Thornburg, B. Wardell, and M. van de Meent, Phys. Rev. Res.2, 013365 (2020), arXiv:1906.06791 [gr-qc]

    arXiv 2020

  53. [53]

    Nasipak, T

    Z. Nasipak, T. Osburn, and C. R. Evans, Phys. Rev. D 100, 064008 (2019), arXiv:1905.13237 [gr-qc]

    arXiv 2019

  54. [54]

    Merlin, A

    C. Merlin, A. Ori, L. Barack, A. Pound, and M. van de Meent, Phys. Rev. D94, 104066 (2016), arXiv:1609.01227 [gr-qc]

    Pith/arXiv arXiv 2016

  55. [55]

    Barack and P

    L. Barack and P. Giudice, Phys. Rev. D95, 104033 (2017), arXiv:1702.04204 [gr-qc]

    Pith/arXiv arXiv 2017

  56. [56]

    Kelly, A

    C. Kelly, A. Vaswani, L. Barack, and A. Pound, In progress

  57. [57]

    Black Hole Perturbation Toolkit,

    “Black Hole Perturbation Toolkit,” (bhptoolkit.org)