Producing self-force data on Kerr backgrounds requires solving the Teukolsky and related perturbation equations at the limit of available numerical techniques. The required accuracy is severe — phase precision of one part in a million is typical — and the parameter space is large, with each waveform model depending on hundreds of thousands of individual field-equation solutions.
We develop solvers spanning a range of complementary techniques. Frequency-domain methods, including MST series and spectral expansions, give high accuracy for bound orbits. Hyperboloidal slicing enables long-time evolutions free of artificial boundary conditions. Our solvers are integrated into the open-source Black Hole Perturbation Toolkit and are used by groups worldwide.
Current focus:
- GPU-accelerated frequency-domain Teukolsky solvers
- Hyperboloidal frequency-domain solutions at second order
- Spectral methods for puncture-regularized sources
- Continued maintenance and release of solvers through BHPToolkit