Published December 31, 2024 | Published
Journal Article Open

Recipes for computing radiation from a Kerr black hole using a generalized Sasaki-Nakamura formalism: Homogeneous solutions

  • 1. ROR icon University of Copenhagen
  • 2. ROR icon California Institute of Technology

Abstract

Central to black hole perturbation theory calculations is the Teukolsky equation that governs the propagation and the generation of radiation emitted by Kerr black holes. However, it is plagued by a long-ranged potential associated with the perturbation equation and hence a direct numerical integration of the equation is challenging. Sasaki and Nakamura devised a formulation that transforms the equation into a new equation that is free from the issue for the case of outgoing gravitational radiation. The formulation was later generalized by Hughes to work for any type of radiation. In this work, we revamp the generalized Sasaki-Nakamura (GSN) formalism and explicitly show the transformations that convert solutions between the Teukolsky and the GSN formalism for both in- and outgoing radiation of scalar, electromagnetic, and gravitational types. We derive all necessary ingredients for the GSN formalism to be used in numerical computations. In particular, we present a new numerical implementation of the formalism, GeneralizedSasakiNakamura.jl, that computes homogeneous solutions to both perturbation equations in the Teukolsky and the GSN formalism. The code works well at low frequencies and is even better at high frequencies by leveraging the fact that black holes are highly permeable to waves at high frequencies. This work lays the foundation for an efficient scheme to compute gravitational radiation from Kerr black holes and an alternative way to compute quasinormal modes of Kerr black holes.

Copyright and License

© 2024 American Physical Society

Acknowledgement

The author would like to thank Yanbei Chen, Manu Srivastava, Shuo Xin, Emanuele Berti, Scott Hughes, Aaron Johnson, Jonathan Thompson, and Alan Weinstein for the valuable discussions and insights when preparing this work. The author would like to especially thank Manu Srivastava for the read of an early draft of this manuscript and Shuo Xin for performing the scaled Wronskian calculations using the fortran code in Refs. [47, 48]. R. K. L. L. acknowledges support from the National Science Foundation Awards No. PHY-1912594 and No. PHY-2207758. R. K. L. L. also acknowledges support from the research Grants No. VIL37766 and No. VIL53101 by the Villum Fonden, the DNRF Chair program Grant No. DNRF162 by the Danish National Research Foundation, the European Union’s Horizon 2020 research and the innovation programme under the Marie Sklodowska-Curie Grant Agreement No. 101131233.

Funding

R. K. L. L. acknowledges support from the National Science Foundation Awards No. PHY-1912594 and No. PHY-2207758. R. K. L. L. also acknowledges support from the research Grants No. VIL37766 and No. VIL53101 by the Villum Fonden, the DNRF Chair program Grant No. DNRF162 by the Danish National Research Foundation, the European Union’s Horizon 2020 research and the innovation programme under the Marie Sklodowska-Curie Grant Agreement No. 101131233.

Data Availability

Mathematica notebooks deriving and storing all the expressions shown here are publicly available from the Zenodo repository [66].

Files

PhysRevD.110.124070.pdf
Files (2.4 MB)
Name Size Download all
md5:1cd5b2380a162d6e5b7bee2448e71d6f
2.4 MB Preview Download

Additional details

Created:
January 2, 2025
Modified:
January 2, 2025