Highly accurate and efficient self-force computations using time-domain methods: Error estimates, validation, and optimization
by Thornburg, Jonathan
27 pages, 12 eps figures (10 of them color, but all are viewable ok in black-and-white), uses RevTeX 4.1
If a small “particle” of mass $latex \mu M$ (with $latex \mu \ll 1$) orbits a Schwarzschild or Kerr black hole of mass $latex M$, the particle is subject to an $latex \O(\mu)$ radiation-reaction “self-force”. Here I argue that it’s valuable to compute this self-force highly accurately (relative error of $latex \ltsim 10^{-6}$) and efficiently, and I describe techniques for doing this and for obtaining and validating error estimates for the computation. I use an adaptive-mesh-refinement (AMR) time-domain numerical integration of the perturbation equations in the Barack-Ori mode-sum regularization formalism; this is efficient, yet allows easy generalization to arbitrary particle orbits. I focus on the model problem of a scalar particle in a circular geodesic orbit in Schwarzschild spacetime.
The mode-sum formalism gives the self-force as an infinite sum of regularized spherical-harmonic modes $latex \sum_{\ell=0}^\infty F_{\ell,\reg}$, with $latex F_{\ell,\reg}$ (and an “internal” error estimate) computed numerically for $latex \ell \ltsim 30$ and estimated for larger~$latex \ell$ by fitting an asymptotic “tail” series. Here I validate the internal error estimates for the individual $latex F_{\ell,\reg}$ using a large set of numerical self-force computations of widely-varying accuracies. I present numerical evidence that the actual numerical errors in $latex F_{\ell,\reg}$ for different~$latex \ell$ are at most weakly correlated, so the usual statistical error estimates are valid for computing the self-force. I show that the tail fit is numerically ill-conditioned, but this can be mostly alleviated by renormalizing the basis functions to have similar magnitudes.
Using AMR, fixed mesh refinement, and extended-precision floating-point arithmetic, I obtain the (contravariant) radial component of the self-force for a particle in a circular geodesic orbit of areal radius $latex r = 10M$ to within $latex 1$~ppm relative error.