[Paper Review] Second-Order Black Hole Perturbations: A Computer Algebra Approach, I - The Schwarzschild Spacetime

George Davies|ArXiv.org|Oct 18, 1998

Relativity and Gravitational Theory1 references3 citations

TL;DR

This paper presents a computer algebra-based derivation of second-order perturbations in Schwarzschild spacetime, focusing on the quadrapole-quadrapole coupling. Using GRTensorII on Maple, the authors compute the effective source term for the second-order Zerilli equation, confirming earlier results by Gleiser with a direct method that bypasses standard function transformations, enabling scalable higher-order analysis with publicly available tools.

ABSTRACT

This article outlines our derivation of the second order perturbations to a Schwarzschild black hole, highlighting our use of, and necessary reliance on, computer algebra. The particular perturbation scenario that is presented here is the case of the linear quadrapole seeding the second order quadrapole. This problem amounts to finding the second order Zerilli wave equation, and in particular the effective source term due to the linear quadrapole. With one minor exception, our calculations confirm the earlier findings of Gleiser, et.al. On route to these results we also illustrate that, with the aid of computer algebra, the linear Schwarzschild problem can be solved in a very direct manner (i.e., without resorting to the usual function transformations), and it is this ``direct method'' that drives the higher order perturbation analysis. The calculations were performed using the GRTensorII computer algebra package, running on the Maple V platform, along with several new Maple routines that we have written specifically for these types of problems. Although we have chosen to consider only the ``quadrapole-quadrapole'' calculation in this article, the GRTensor environment, with the inclusion of these new routines, would allow this analysis to be repeated for a far more general problem. These routines, along with Maple worksheets that reproduce our calculations, are publicly available at the GRTensor website: www.astro.queensu.ca/~grtensor . The interested reader is invited to download and use them to reproduce our results and experiment.

Motivation & Objective

• To develop a systematic, computationally efficient method for calculating second-order black hole perturbations in the Schwarzschild geometry.
• To demonstrate the feasibility and advantages of using computer algebra systems—specifically GRTensorII on Maple—for handling complex tensorial calculations in general relativity.
• To provide a direct, transformation-free approach to solving linear perturbations, which facilitates extension to higher-order problems.
• To derive the effective source term for the second-order Zerilli wave equation in the case of linear quadrapole seeding a second-order quadrapole response.
• To make the computational tools and Maple workbooks publicly available for reproducibility and broader application.

Proposed method

• Employing the GRTensorII computer algebra package on the Maple V platform to perform symbolic tensor computations in curved spacetime.
• Implementing custom Maple routines specifically designed for second-order perturbation theory in black hole spacetimes.
• Using a direct method to solve the linear Schwarzschild perturbation problem without relying on function transformations such as those involving the Teukolsky or Zerilli functions.
• Deriving the second-order Zerilli equation by computing the nonlinear source term arising from the product of linear perturbations.
• Validating results against earlier work by Gleiser through consistent symbolic computation and cross-checking.
• Publishing all computational code and worksheets to ensure reproducibility and enable extension to other multipole modes.

Experimental results

Research questions

• RQ1 How can second-order perturbations in Schwarzschild spacetime be systematically computed using symbolic computation?
• RQ2 What advantages does a direct, transformation-free method offer over conventional approaches in solving linear and second-order perturbations?
• RQ3 Can the effective source term for the second-order Zerilli equation be accurately derived using computer algebra for the quadrapole-quadrapole case?
• RQ4 To what extent can the computational framework be generalized to other multipole modes and perturbation scenarios?
• RQ5 How can publicly available computational tools enhance reproducibility and scalability in higher-order black hole perturbation theory?

Key findings

• The second-order Zerilli wave equation is successfully derived for the quadrapole-quadrapole case using a computer algebra approach.
• The effective source term for the second-order perturbation is computed and found to be in agreement with earlier results by Gleiser, confirming their findings.
• The direct method of solving linear perturbations—without function transformations—proves effective and scalable for higher-order analysis.
• The use of GRTensorII with custom Maple routines enables efficient, accurate, and reproducible computation of complex tensorial expressions in second-order perturbation theory.
• All computational code, including Maple worksheets, is made publicly available at the GRTensor website for reuse and extension by the research community.
• The framework is extensible and can be applied to a broader class of perturbation problems beyond the quadrapole-quadrapole case.