Caustic echoes from a Schwarzschild black hole

We present the first numerical construction of the scalar Schwarzschild Green function in the time-domain, which reveals several universal features of wave propagation in black hole spacetimes. We demonstrate the trapping of energy near the photon sphere and confirm its exponential decay. The trapped wavefront propagates through caustics resulting in echoes that propagate to infinity. The arrival times and the decay rate of these caustic echoes are consistent with propagation along null geodesics and the large l-limit of quasinormal modes. We show that the four-fold singularity structure of the retarded Green function is due to the well-known action of a Hilbert transform on the trapped wavefront at caustics. A two-fold cycle is obtained for degenerate source-observer configurations along the caustic line, where the energy amplification increases with an inverse power of the scale of the source. Finally, we discuss the tail piece of the solution due to propagation within the light cone, up to and including null infinity, and argue that, even with ideal instruments, only a finite number of echoes can be observed. Putting these pieces together, we provide a heuristic expression that approximates the Green function with a few free parameters. Accurate calculations and approximations of the Green function are the most general way of solving for wave propagation in curved spacetimes and should be useful in a variety of studies such as the computation of the self-force on a particle.

💡 Research Summary

The authors present the first time‑domain numerical construction of the scalar Green function on a Schwarzschild background, revealing a rich structure of wave propagation that had previously been inferred only from analytic approximations. By evolving an initially compact scalar pulse with a high‑resolution finite‑difference scheme, they demonstrate that a substantial fraction of the wave energy becomes trapped near the photon sphere (r = 3M). The trapped front circles the photon sphere repeatedly, following null geodesics, and each time it passes through a caustic (the focal line of the photon sphere) it undergoes a Hilbert‑transform phase shift of π/2. This phase shift produces a four‑fold singularity cycle in the retarded Green function: a Dirac delta, a principal‑value 1/t term, a derivative of a delta, and another principal‑value term. Each cycle corresponds to a distinct “echo” that propagates outward to infinity.

The arrival times of the echoes follow tₙ ≈ 2π n · 3M, i.e. the light‑travel time for n complete orbits around the photon sphere. Their amplitudes decay exponentially as e^{‑γ n}, where γ equals the instability exponent of the photon sphere (≈0.0889/M). This decay matches the imaginary part of the large‑ℓ quasinormal‑mode frequencies ω_{ℓn}=ℓ Ω_c − i (n+½) |λ|, confirming that the echoes are the time‑domain manifestation of the high‑ℓ QNM spectrum.

In the special degenerate configuration where source and observer lie on the caustic line, the singularity cycle reduces to a two‑fold pattern, and the echo amplitude scales as ε^{‑1/2} with the source’s spatial scale ε. This inverse‑power amplification arises from the strong focusing of the wavefront at the caustic.

Beyond the discrete echoes, the authors analyse the “tail” component that lives inside the light cone. The tail originates from waves that never reach the photon sphere and instead disperse throughout the interior region. Its late‑time behaviour follows a power law t^{‑2ℓ‑3}, consistent with classic Price‑law results, and it contributes a slowly decaying background that would be observable only with unrealistically sensitive detectors.

Importantly, the paper argues that even with ideal instruments only a finite number of echoes (typically three to five) can be detected, because each subsequent echo is exponentially weaker and increasingly delayed. To make the results practically useful, the authors propose a compact heuristic expression for the full Green function that involves a small set of parameters: echo amplitudes, decay rates, arrival times, and tail strength. This formula reproduces the numerical Green function to high accuracy and can be employed in self‑force calculations, waveform modelling, and any problem that requires an accurate propagator in curved spacetime.

Overall, the work provides a concrete, numerically validated picture of how scalar (and by extension, electromagnetic or gravitational) waves behave near a Schwarzschild black hole: trapping at the photon sphere, periodic caustic‑induced phase shifts, exponential damping linked to quasinormal modes, and a residual power‑law tail. These insights deepen our theoretical understanding and furnish a practical tool for a wide range of relativistic astrophysics applications.