Skip to content
DraconiticResearch
Search Sources API
203,203 studies31 sourcesUpdated 39 min ago
PreprintarXiv licencearXiv

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

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.

Versions

Related

  1. Parametrized beyond-Teukolsky framework in the time domain
    Ciro De Simone, Sebastian H. Völkel, Kostas D. Kokkotas +1
    arXivPreprintarXiv licenceFull textHigh-energy astrophysicsGravitation & general relativityarXiv
  2. Reduced magnetohydrodynamics with electron inertia in pseudodipole coordinates
    Khilav Majmudar, Robert L. Lysak
    arXivPreprintCC BY 4.0Full textSolar & stellar physicsComputational physicsarXiv
  3. Natural coordinates for constrained correlation functions: Partial autocorrelations and the geometry of positive power spectra
    Thomas Erben
    arXivPreprintCC BY 4.0Full textInstrumentation & methodsMathematical physicsarXiv
  4. The Equations of Reduced Magnetohydrodynamics in Dipole Coordinates
    Ian M. DesJardin, Christopher Bard, Natalia Y. Buzulukova +1
    arXivPreprintCC BY-NC-ND 4.0Full textSolar & stellar physicsPlasma physicsarXiv
  5. Analysis of wave processes using beam-driven Langmuir/Z-mode waveforms generated in particle-in-cell simulations
    F J Polanco-Rodríguez, C. Krafft, P. Savoini
    Astronomy & Astrophysics - A&APeer-reviewedCC BY 4.0Full text1 citedSolar & stellar physicsHAL
BibTeX
@misc{vaswani2026timedomain,
  title = {{Time-domain framework for the Teukolsky equation with a particle source using comoving hyperboloidal coordinates}},
  author = {Aditya Vaswani and Leor Barack and Oliver Long and Rodrigo Panosso Macedo},
  howpublished = {arXiv},
  year = {2026},
  month = oct,
  doi = {10.48550/arxiv.2606.26339},
  eprint = {2606.26339},
  archivePrefix = {arXiv},
  url = {https://arxiv.org/abs/2606.26339},
}
RIS
TY  - UNPB
TI  - Time-domain framework for the Teukolsky equation with a particle source using comoving hyperboloidal coordinates
AU  - Aditya Vaswani
AU  - Leor Barack
AU  - Oliver Long
AU  - Rodrigo Panosso Macedo
PY  - 2026
DA  - 2026/10/26
PB  - arXiv
DO  - 10.48550/arxiv.2606.26339
UR  - https://arxiv.org/abs/2606.26339
AB  - 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.
ER  -