The radial component of a perturbation outside the event horizon satisfies the following wave equation,
For “axial” perturbations, which is a single potential barrier with a peak around , which is the location of the unstable photon orbit. The form (23) is true even if we consider scalar or electromagnetic test fields as perturbations. The parameter takes the values 1 for scalar perturbations, 0 for electromagnetic perturbations, and –3 for gravitational perturbations and can be expressed as , where is the spin of the perturbing field.
For “polar” perturbations the effective potential was derived by Zerilli  and has the form has shown that one can transform the Equation (21) for “axial” modes to the corresponding one for “polar” modes via a transformation involving differential operations. It can also be shown that both forms are connected to the Bardeen–Press  perturbation equation derived via the Newman–Penrose formalism. The potential decays exponentially near the horizon, , and as for .
From the form of Equation (21) it is evident that the study of black hole perturbations will follow the footsteps of the theory outlined in Section 2.
Kay and Wald  have shown that solutions with data of compact support are bounded. Hence we know that the time independent Green function is analytic for . The essential difficulty is now to obtain the solutions (cf. Equation (10)) of the equation
As the horizon () is a regular singular point of (26), a representation of as a converging series exists. For it reads:. (The analytic extension of is investigated in .) The result is that has an extension to the complex plane with poles only at negative real integers. The representation of is more complicated: Because infinity is a singular point no power series expansion like (28) exists. A representation coming from the iteration of the defining integral equation is given by Jensen and Candelas , see also . It turns out that the continuation of has a branch cut due to the decay for large .
The most extensive mathematical investigation of quasi-normal modes of the Schwarzschild solution is contained in the paper by Bachelot and Motet-Bachelot . Here the existence of an infinite number of quasi-normal modes is demonstrated. Truncating the potential (23) to make it of compact support leads to the estimate (16).
The decay of solutions in time is not exponential because of the weak decay of the potential for large . At late times, the quasi-normal oscillations are swamped by the radiative tail [166, 167]. This tail radiation is of interest in its own right since it originates on the background spacetime. The first authoritative study of nearly spherical collapse, exhibiting radiative tails, was performed by Price [166, 167].
Studying the behavior of a massless scalar field propagating on a fixed Schwarzschild background, he showed that the field dies off with the power-law tail,[66, 67, 68]. Today it is apparent in any simulation involving evolutions of various fields on a black hole background including Schwarzschild, Reissner–Nordström , and Kerr [132, 133]. It has also been observed in simulations of axial oscillations of neutron stars , and should also be present for polar oscillations. Leaver  has studied in detail these tails and associated this power low tail with the branch-cut integral along the negative imaginary axis in the complex plane. His suggestion that there will be radiative tails observable at and has been verified by Gundlach, Price, and Pullin . Similar results were arrived at recently by Ching et al.  in a more extensive study of the late time behavior. In a nonlinear study Gundlach, Price, and Pullin  have shown that tails develop even when the collapsing field fails to produce a black hole. Finally, for a study of tails in the presence of a cosmological constant refer to , while for a recent study, using analytic methods, of the late-time tails of linear scalar fields outside Schwarzschild and Kerr black holes refer to [36, 37].
Using the properties of the waves at the horizon and infinity given in Equation (27) one can search for the quasi-normal mode frequencies since practically the whole problem has been reduced to a boundary value problem with being the complex eigenvalue. The procedure and techniques used to solve the problem will be discussed later in Section 6, but it is worth mentioning here a simple approach to calculate the QNM frequencies proposed by Schutz and Will . The approach is based on the standard WKB treatment of wave scattering on the peak of the potential barrier, and it can be easily shown that the complex frequency can be estimated from the relation
|0||0.37367||–0.08896 i||0.59944||–0.09270 i||0.80918||–0.09416 i|
|1||0.34671||–0.27391 i||0.58264||–0.28130 i||0.79663||–0.28443 i|
|2||0.30105||–0.47828 i||0.55168||–0.47909 i||0.77271||–0.47991 i|
|3||0.25150||–0.70514 i||0.51196||–0.69034 i||0.73984||–0.68392 i|
Figure 2 shows some of the modes of the Schwarzschild black hole. The number of modes for each harmonic index is infinite, as was mathematically proven by Bachelot and Motet-Bachelot . This was also implied in an earlier work by Ferrari and Mashhoon , and it has been seen in the numerical calculations in [25, 157]. It can be also seen that the imaginary part of the frequency grows very quickly. This means that the higher modes do not contribute significantly in the emitted gravitational wave signal, and this is also true for the higher modes (octapole etc.).
As is apparent in Figure 2 that there is a special purely imaginary QNM frequency. The existence of “algebraically special” solutions for perturbations of Schwarzschild, Reissner–Nordström and Kerr black holes were first pointed out by Chandrasekhar . It is still questionable whether these frequencies should be considered as QNMs  and there is a suggestion that the potential might become transparent for these frequencies . For a more detailed discussion refer to .
As a final comment we should mention that as the order of the modes increases the real part of the frequency remains constant, while the imaginary part increases proportionally to the order of the mode. Nollert  derived the following approximate formula for the asymptotic behavior of QNMs of a Schwarzschild black hole, and .
For large values of the distribution of QNMs is given by [164, 86, 85, 113].
The perturbations of Reissner–Nordström black holes, due to the spherical symmetry of the solution, follow the footsteps of the analysis that we have presented in this section. Most of the work was done during the seventies by Zerilli , Moncrief [153, 154] and later by Chandrasekhar and Xanthopoulos [55, 209]. For an extensive discussion refer to . We have again wave equations of the form (21), one for each parity with potentials which are like (23) and (24) plus extra terms which relate to the charge of the black hole. An interesting feature of the charged black holes is that any perturbation of the gravitational (electromagnetic) field will also induce electromagnetic (gravitational) perturbations. In other words, any perturbation of the Reissner–Nordström spacetime will produce both electromagnetic and gravitational radiation. Again it has been shown that the solutions for the odd parity oscillations can be deduced from the solutions for even parity oscillations and vice versa . The QNM frequencies of the Reissner–Nordström black hole have been calculated by Gunter , Kokkotas, and Schutz , Leaver , Andersson , and lately for the nearly extreme case by Andersson and Onozawa .
© Max Planck Society and the author(s)