The Synergy between Numerical and Perturbative Approaches to Black Holes

Edward Seidel

DOI: https://doi.org/10.48550/arXiv.gr-qc/9806087

1998-06-23

Abstract:I describe approaches to the study of black hole spacetimes via numerical relativity. After a brief review of the basic formalisms and techniques used in numerical black hole simulations, I discuss a series of calculations from axisymmetry to full 3D that can be seen as stepping stones to simulations of the full 3D coalescence of two black holes. In particular, I emphasize the interplay between perturbation theory and numerical simulation that build both confidence in present results and tools to aid and to interpret results of future simulations of black hole coalescence.

General Relativity and Quantum Cosmology

What problem does this paper attempt to address?

The problem that this paper attempts to solve is: How to study the dynamical evolution of black - hole spacetime by combining numerical relativity and perturbation theory, especially during the process of two black holes colliding. Specifically, the author explores the step - by - step method from axisymmetric to full three - dimensional simulations and emphasizes the synergy between perturbation theory and numerical simulations to enhance confidence in current results and provide tools for future binary - black - hole merger simulations.

### Main problems
1. **Challenges in numerical simulations**:
- Black - hole numerical simulations are very difficult because one must simultaneously handle the highly nonlinear region inside the singularity, near the event horizon, and the linear region in the radiation zone.
- Using the traditional singularity - avoiding time - slicing method in 3D will lead to numerical instability, making the simulation difficult to carry out.

2. **Application of perturbation theory**:
- Perturbation theory plays an important role in verifying and interpreting numerical simulation results, especially when extracting waveforms and dealing with black - hole initial data sets.
- By applying perturbation theory to numerical relativity, one can better understand the dynamics of distorted and colliding black holes.

### Solutions
- **Combination of perturbation theory and numerical simulation**:
- Use perturbation theory to test and interpret fully nonlinear numerical simulation results.
- Transition gradually from the simple axisymmetric case to the complex 3D simulation, forming a "confidence ladder".

- **Specific applications**:
- **Waveform extraction**: Extract waveforms by decomposing the numerically generated metric into a spherically symmetric part and a perturbed part and using spherical - harmonic - function projection.
- **Axisymmetric distorted black hole**: Study the evolution of a single axisymmetric black hole under the influence of gravitational waves and verify the consistency of linear and nonlinear results.
- **Axisymmetric black - hole collision**: Use perturbation theory to simulate black - hole collisions in two limit cases (approximate limit and long - distance limit).
- **3D test platform**: Extend the existing axisymmetric initial data set to 3D simulations and verify the consistency of results of different codes.

### Conclusion
By combining numerical relativity and perturbation theory, the author shows the effectiveness of this method in studying black - hole dynamics, especially for highly distorted and colliding black - hole systems. This not only enhances confidence in existing results but also provides a powerful tool for future complex simulations.

### Formula summary
- **Metric decomposition**:
\[
g_{\alpha\beta}=\bar{g}_{\alpha\beta}+h_{\alpha\beta}
\]
where \(\bar{g}_{\alpha\beta}\) is the Schwarzschild background metric and \(h_{\alpha\beta}\) is a small perturbation.

- **Zerilli equation**:
\[
\frac{\partial^{2}\psi^{(\ell m)}}{\partial r_{*}^{2}}-\frac{\partial^{2}\psi^{(\ell m)}}{\partial t^{2}}+V^{(\ell)}(r)\psi^{(\ell m)} = 0
\]
where
\[
V^{(\ell)}(r)=\left(1-\frac{2M}{r}\right)\left[\frac{1}{\Lambda^{2}}\left(\frac{72M^{3}}{r^{5}}-\frac{12M}{r^{3}}(\ell - 1)(\ell + 2)\left(1-\frac{3M}{r}\right)\right)+\frac{\ell(\ell - 1)(\ell + 2)(\ell + 1)}{r^{2}\Lambda}\right]
\]
\(r_{*}\) is the so - called tortoise coordinate, defined as \(r_{*}=r + 2M\ln(r/2M - 1)\), and \(\Lambda\) is a function depending on \(r\).

These formulas show the core content of perturbation theory and are used to describe the gravitational perturbation of Schwarzschild black holes.