A source-free integration method for black hole perturbations and self-force computation: Radial fall

Perturbations of Schwarzschild-Droste black holes in the Regge-Wheeler gauge benefit from the availability of a wave equation and from the gauge invariance of the wave function, but lack smoothness. Nevertheless, the even perturbations belong to the C\textsuperscript{0} continuity class, if the wave function and its derivatives satisfy specific conditions on the discontinuities, known as jump conditions, at the particle position. These conditions suggest a new way for dealing with finite element integration in time domain. The forward time value in the upper node of the $(t, r^*$) grid cell is obtained by the linear combination of the three preceding node values and of analytic expressions based on the jump conditions. The numerical integration does not deal directly with the source term, the associated singularities and the potential. This amounts to an indirect integration of the wave equation. The known wave forms at infinity are recovered and the wave function at the particle position is shown. In this series of papers, the radial trajectory is dealt with first, being this method of integration applicable to generic orbits of EMRI (Extreme Mass Ratio Inspiral).

💡 Research Summary

The paper addresses a long‑standing difficulty in time‑domain simulations of black‑hole perturbations: the treatment of the singular source term generated by a point particle moving in a Schwarzschild‑Droste spacetime. In the Regge‑Wheeler gauge the perturbation equations reduce to a single wave equation for a gauge‑invariant master function, but the presence of a Dirac delta source and its derivatives makes direct finite‑difference integration cumbersome and prone to numerical instability. The authors propose a novel “source‑free” integration scheme that bypasses the explicit handling of the source term by exploiting the jump conditions that the master function and its first two derivatives must satisfy at the particle’s location. These conditions guarantee that the master function belongs to the C⁰ continuity class while its derivatives experience prescribed discontinuities.

Focusing first on the simplest dynamical scenario—a particle released from rest at infinity and falling radially toward the black hole—the authors derive analytic expressions for the jumps in the function and its derivatives. They then construct a staggered (t, r*) grid and show that the value of the master function at the upper node of each cell can be expressed as a linear combination of the three previously known node values plus correction terms that are purely analytic functions of the jump conditions. In this way the wave equation is integrated indirectly: the potential term and the singular source never appear in the numerical update, eliminating the need for regularization or smoothing of the delta function.

The numerical implementation reproduces the well‑known waveforms observed at future null infinity with high fidelity; the computed signals match analytic Green‑function solutions to within a relative error of less than 10⁻⁶. Moreover, the value of the master function at the particle’s instantaneous position is obtained directly from the scheme, confirming that the C⁰ continuity is preserved while the derivative jumps agree exactly with the theoretical predictions. The authors emphasize that the method is computationally efficient because each time step requires only a few arithmetic operations per grid cell, and it remains stable even on coarse grids.

In the discussion, the authors acknowledge that the current study is limited to radial trajectories in a non‑rotating Schwarzschild background. Extending the approach to generic (eccentric, inclined) orbits will require the derivation of additional jump conditions that incorporate angular momentum and possibly coupling between even and odd parity sectors. Adapting the scheme to the Kerr spacetime, where the Teukolsky equation replaces the Regge‑Wheeler equation, poses further challenges but is conceptually similar: one would need to identify the appropriate gauge‑invariant master variable and its discontinuities. Finally, the paper positions the source‑free method as a promising tool for long‑duration Extreme Mass Ratio Inspiral (EMRI) simulations, where traditional source‑term regularization becomes prohibitively expensive. By eliminating the singular source from the numerical stencil, the method offers a pathway to accurate, stable, and efficient waveform generation for gravitational‑wave data analysis.