REVIEW
Second order gauge invariant gravitational perturbations of a Kerr black hole
General Relativity and Quantum Cosmology · Submitted 5 Nov 1998 · 45 references
Pith's verdict: A real construction — the tetrad-invariant piece checks out and the second-order Kerr wave equation is genuinely new — but the gauge-vector integrability is assumed, not shown, and the paper says so.
T0 · AI review
Why Pith says this
Highlight any text to comment on it.
Referee report
The paper
The paper claims
The paper builds a second-order gravitational waveform around a rotating black hole that does not change when coordinates or the null frame are changed, and shows it obeys the same Teukolsky wave equation as the first-order field with a new
The reading
What carries the argument
The load-bearing object is the corrected waveform of Eq. (41), $psi_I = psi_4^{(2)} + Q - Q_{AF} + psi_4^{(1)},_{mu}(xi^mu_{ORG} - xi^mu)$, where $Q$ is the quadratic tetrad-correcting term and $Q_{AF}$ is its counterpart evaluated in the outgoing radiation gauge. The tetrad part of the correction works because $psi_3^{(1)}$ shifts by $3a,psi_2^{(0)}$ under a null rotation of type I, which lets the rotation parameter be solved for and replaced; the coordinate part works because the first-order gauge vector can be recovered from variations of the tetrad-invariant $psi_2^{(1)}$ (giving $xi^r$ and $xi^vartheta$) and from variations of the metric and extrinsic curvature (giving the derivatives of $xi^t$ and $xi^phi$). The payoff is that the same Teukolsky operator $hat{T}$, written explicitly in Boyer-Lindquist coordinates and the Kinnersley tetrad, acts on the second-order invariant field, so an existing first-order evolution code can be reused with an added source term.
Extended reading notes
Core claim
The central claim is an existence proof by construction: for second-order perturbations about a Kerr black hole there exists a waveform $psi_I$ that is invariant under first-order coordinate (gauge) transformations and under tetrad rotations, and that coincides with the second-order $psi_4$ in an asymptotically flat gauge at the radiation zone. The construction proceeds in two stages. First, $psi_4^{(2)}$ is corrected by a quadratic term involving the first-order tetrad components and the first-order Weyl scalars, $2psi_4^{(1)}(l^r_{(1)}/l^r_{(0)} + m^vartheta_{(1)}/m^vartheta_{(0)}) - (2/3)(psi_3^{(1)})^2/psi_2^{(0)}$, which cancels the tetrad dependence, with a further subtraction in the outgoing radiation gauge so that the added piece does not survive at large $r$. Second, the residual coordinate dependence, which is linear in the first-order gauge vector $xi^mu$, is cancelled by adding $psi_4^{(1)}$ times the difference between the gauge vector in the outgoing radiation gauge and the actual one. The resulting $psi_I$ obeys $hat{T}psi_I = S_I$, where $hat{T}$ is exactly the operator governing first-order perturbations and $S_I$ is a source assembled from first-order perturbations and their derivatives. The authors also lay out how to impose initial data for $psi_I$, how to extract radiated energy, linear momentum and angular momentum from it, and how to fix the two residual integration constants by an a posteriori gauge choice.
Load-bearing premise
The construction assumes that the time and azimuthal components of the first-order gauge vector can be recovered by integrating their derivatives in Eqs. (45)-(46) with consistent integrability conditions, and that the outgoing-radiation-gauge reference vectors exist so that the difference between the two is well defined — the paper states in Sec. IV.C that these components are 'assumed a priori' and tied to the existence of the outgoing radiation gauge proved in Ref. [32].
Editorial extensions
If this is right
- A second-order radiation field can be computed without committing to a coordinate system, so its waveform and energy flux can be compared directly with the covariant $psi_4$ produced by full numerical integrations of Einstein's equations.
- Because the wave operator is unchanged from first order, the (2+1)-dimensional time-domain Teukolsky evolution code already used for first-order perturbations can be extended to second order by adding a source term built from the first-order solution.
- The ratio of second-order to first-order radiated energy supplies the otherwise missing error estimate for the close-limit approximation, giving a systematic measure of how small the perturbative parameter must be for the approximation to hold.
- The formalism reduces to the nonrotating case and is intended to reproduce the known second-order Schwarzschild close-limit results, while also giving access to higher multipoles such as $ell = 4$.
- Radiated energy, linear momentum and angular momentum at future null infinity can all be computed from $psi_I$ using the standard Bondi and Landau-Lifshitz expressions, since $psi_I$ equals $psi_4^{(2)}$ in an asymptotically flat gauge.
Reading between the lines
- The construction amounts to a machine that converts any first-order gauge into an asymptotically flat one, which suggests that the outgoing radiation gauge may be the most economical setting for the actual numerical evolution, rather than evolving in an arbitrary gauge and correcting afterwards.
- The method is presented for Kerr but the authors note it should carry over to any Petrov type D (or even type II) background, so the same template might be reused for second-order perturbations of nonrotating neutron stars or of cosmological backgrounds.
- The non-uniqueness of the quadratic correction — one can always add products of first-order invariants — means that different implementations will produce waveforms differing at high order in $1/r$; a useful cross-check would be to verify that two such choices give the same radiated energy at infinity.
- Because the correction requires integrating $xi^t$ and $xi^phi$ from their derivatives, the resulting waveform is nonlocal in the first-order data even though the evolution equation is local; this asymmetry between the equation and the observable is worth watching in any implementation.
Original abstract
We investigate higher than the first order gravitational perturbations in the Newman-Penrose formalism. Equations for the Weyl scalar $\psi_4,$ representing outgoing gravitational radiation, can be uncoupled into a single wave equation to any perturbative order. For second order perturbations about a Kerr black hole, we prove the existence of a first and second order gauge (coordinates) and tetrad invariant waveform, $\psi_I$, by explicit construction. This waveform is formed by the second order piece of $\psi_4$ plus a term, quadratic in first order perturbations, chosen to make $\psi_I$ totally invariant and to have the appropriate behavior in an asymptotically flat gauge. $\psi_I$ fulfills a single wave equation of the form ${\cal T}\psi_I=S,$ where ${\cal T}$ is the same wave operator as for first order perturbations and $S$ is a source term build up out of (known to this level) first order perturbations. We discuss the issues of imposition of initial data to this equation, computation of the energy and momentum radiated and wave extraction for direct comparison with full numerical approaches to solve Einstein equations.
The compiled record
The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.
The abstract, sentence by sentence
We investigate higher than the first order gravitational perturbations in the Newman-Penrose formalism. Equations for the Weyl scalar $\psi_4,$ representing outgoing gravitational radiation, can be uncoupled into a single wave equation to any perturbative order. For second order perturbations about a Kerr black hole, we prove the existence of a first and second order gauge (coordinates) and tetrad invariant waveform, $\psi_I$, by explicit construction. This waveform is formed by the second order piece of $\psi_4$ plus a term, quadratic in first order perturbations, chosen to make $\psi_I$ totally invariant and to have the appropriate behavior in an asymptotically flat gauge. $\psi_I$ fulfills a single wave equation of the form ${\cal T}\psi_I=S,$ where ${\cal T}$ is the same wave operator as for first order perturbations and $S$ is a source term build up out of (known to this level) first order perturbations. We discuss the issues of imposition of initial data to this equation, computation of the energy and momentum radiated and wave extraction for direct comparison with full numerical approaches to solve Einstein equations.
background contribution method result interpretation main claim
| # | sentence | role | the authors' own work | states the main conclusion | strength |
|---|---|---|---|---|---|
| s0 | We investigate higher than the first order gravitational perturbations in the Newman-Penrose formalism. | contribution 0.82 | 0.90 | 0.03 | asserted 0.88 |
| s1 | Equations for the Weyl scalar $\psi_4,$ representing outgoing gravitational radiation, can be uncoupled into a single wave equation to any perturbative order. | background 0.80 | 0.38 | 0.09 | asserted 0.91 |
| s2 | For second order perturbations about a Kerr black hole, we prove the existence of a first and second order gauge (coordinates) and tetrad invariant waveform, $\psi_I$, by explicit construction. | result 0.60 | 0.96 | 0.73 | asserted 1.00 |
| s3 | This waveform is formed by the second order piece of $\psi_4$ plus a term, quadratic in first order perturbations, chosen to make $\psi_I$ totally invariant and to have the appropriate behavior in an asymptotically flat gauge. | method 0.67 | 0.91 | 0.24 | asserted 1.00 |
| s4 | $\psi_I$ fulfills a single wave equation of the form ${\cal T}\psi_I=S,$ where ${\cal T}$ is the same wave operator as for first order perturbations and $S$ is a source term build up out of (known to this level) first order perturbations. | result 0.83 | 0.88 | 0.32 | asserted 1.00 |
| s5 | We discuss the issues of imposition of initial data to this equation, computation of the energy and momentum radiated and wave extraction for direct comparison with full numerical approaches to solve Einstein equations. | contribution 0.63 | 0.59 | 0.06 | unclear 0.52 |
What was examined, and what was not
- arXiv listing abstract v2
- full text
- figures and tables
- references and citing sentences
- relations to other papers
- P abstract
- S sections
- T body sentences
- C citations
- R cross-paper relations
References
How later papers used it
88 later papers cite it
1 later paper by other authors is listed below: 1 used its method or tool. 3 more papers by the same authors are listed: 3 used its method or tool.
Uses its method or tool 4 of 43 papers with this label · 3 by the same authors
-
Defined the derivatives following this paper.
Second Order Perturbations of Kerr Black Holes: Reconstruction of the Metric (2020) · 0.92 · paraphrased: the licence does not allow quoting
By the same authors
-
Performed the extraction of gravitational radiation from numerical relativity simulations using formulas (22), (23) and (24) from this paper for the energy and linear and angular momentum radiated in terms of the extracted Weyl scalar Ψ4 at the observer location sphere.
The maximum radiated energy and final spin of high speed collision of two black holes (2024) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Performed the extraction of gravitational radiation from numerical relativity simulations using formulas (22) and (23) from this paper for the energy and linear momentum radiated, and formulas from another work for angular momentum radiated, in terms of the extracted Weyl scalar Ψ4 at the observer location sphere.
Parameter Estimation with Targeted Eccentric Numerical-Relativity Simulations for GW200208_22 and GW190620 (2025) (same authors) · 1.00 · paraphrased: the licence does not allow quoting
-
Derived extrapolation formulas for the radiated energy and momenta calculated from the Weyl scalar ψ4 as in this paper.
Perturbative extraction of gravitational waveforms generated with Numerical Relativity (2015) (same authors) · 0.81 · paraphrased: the licence does not allow quoting
How the literature uses it
83 arXiv papers cite this in 177 sentences
Papers citing it · 86 sentences without a confident label
Contests it 1 paper
-
Our result is consistent with [ref] up to the overall minus sign due to different signatures, which corrects some errors in [ref] .
Perturbations of spinning black holes in dynamical Chern-Simons gravity I. Slow rotation equations · Perturbed NP quantities · contest 0.82
Builds on it 1 paper
-
We further show that our formalism can be generalized to higher order in both $ζ$ and $ϵ$, which is thus a bGR extension of the higher-order Teukolsky formalism in GR developed in [ref] .
Perturbations of spinning black holes beyond General Relativity: Modified Teukolsky equation · Executive summary · 0.99
Uses its method or tool 42 papers
-
From the waveforms we compute the radiated energy and radiated linear and angular momentum using the formulas given in [ref] .
Evolutions of unequal mass, highly spinning black hole binaries · Results · 1.00
-
Thus, [equation] where $S_i$ and $S_f$ are the initial and final angular momentum before and after scattering, respectively. $J_GW$ is the angular momentum radiated away by s, which we compute as [ref] , [equation] Here, $h^+$ and $h^×$ are polarizations of the strain, and a dot denotes derivatives with respect to time.
Spin-up and mass-gain in hyperbolic encounters of spinning black holes · Extraction of observables · 0.97
-
The extraction of gravitational radiation from the numerical relativity simulations is performed using the formulas (22),(23),(24) from [ref] for the energy and linear and angular momentum radiated, respectively, in terms of the extracted Weyl scalar $Ψ_4$ at the observer location sphere $R_obs=113M$.
The maximum radiated energy and final spin of high speed collision of two black holes · Techniques · 1.00
Compares with it 6 papers
-
Here we provide a complete listing of these relations (compare also to Eq. (A4) of [ref] ):
Second Order Perturbations of Kerr Black Holes: Reconstruction of the Metric · Linearized NP spin coefficients in terms of the linearized metric · 0.86
-
Similarly to <ref>, we expand $Ψ_4$ as follows: [equation] This is in contrast to other works on QQNMs [ref] , which define their perturbative expansion similarly to Campanelli and Lousto [ref] , as [equation] The background and first-order perturbations are identical in both expansions.
Quadratic quasinormal modes at null infinity on a Schwarzschild spacetime · Perturbative expansion · 0.93
-
In [ref] , the Teukolsky equation was successfully extended to second- and higher-order, so let us show now that these higher-order equations are very similar to what we obtained in this paper.
Perturbations of spinning black holes beyond General Relativity: Modified Teukolsky equation · Connection to the second-order Teukolsky formalism in GR · 0.84
Cites it as support 2 papers
-
Even if we do not specialize the perturbed tetrad or gauge, $Ψ^(2)_4L$ and $Ψ^(2)_4$ contain identical GW content [ref] .
Quadratic quasinormal modes at null infinity on a Schwarzschild spacetime · General case · 0.99
-
On the one hand, the inability of the linear-order Teukolsky equation to estimate the errors due to the use of a perturbative expansion makes the study of higher-order Teukolsky equations necessary [ref] .
Perturbations of spinning black holes beyond General Relativity: Modified Teukolsky equation · Connection to the second-order Teukolsky formalism in GR · 0.97
Cites it as background 18 papers
-
The equation for first order perturbations was originally derived by Teukolsky [ref] , and was later generalized to $n$-th order perturbations by Campanelli & Lousto [ref] .
Second Order Perturbations of Kerr Black Holes: Reconstruction of the Metric · Master equations for perturbations of a Petrov Type D spacetime · 0.87
-
This involves the same Teukolsky operator as in Eq. [ref] acting on the second-order curvature correction, and a complicated source $𝒮$ that depends quadratically on the linear perturbations [ref] .
Nonlinearities in Black Hole Ringdowns · Introduction · 0.89
-
There have been a number of investigations on this topic, starting with analyses on Schwarzschild black holes [ref] , which characterized the QQNM frequency spectrum and the sources that drive these quadratic modes, followed by generalizations to Kerr black holes [ref] .
Generation and propagation of nonlinear quasi-normal modes of a Schwarzschild black hole · Introduction · 0.92
Discussion
Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.