Linear Stability of Schwarzschild Spacetime: The Cauchy Problem of Metric Coefficients
In this paper, we study the theory of linearized gravity and prove the linear stability of Schwarzschild black holes as solutions of the vacuum Einstein equations. In particular, we prove that solutions to the linearized vacuum Einstein equations centered at a Schwarzschild metric, with suitably regular initial data, remain uniformly bounded and decay to a linearized Kerr metric on the exterior region. Our method employs Hodge decomposition to split the solution into closed and co-closed portions, respectively identified with even-parity and odd-parity solutions in the physics literature. For the co-closed portion, we extend previous results of the first two authors, deriving Regge-Wheeler type equations for two gauge-invariant quantities without the earlier paper’s need of axisymmetry. For the closed portion, we build upon earlier work of Zerilli and Moncrief, wherein the authors derive a master equation for a gauge-invariant quantity in a spherical harmonic decomposition. We work with gauge-invariant quantities at the level of perturbed connection coefficients, with the initial value problem formulated on Cauchy data sets. With the choice of an appropriate gauge in each of the two portions, decay estimates on these decoupled quantities are used to establish decay of the metric coefficients of the solution, completing the proof of linear stability. Our approach complements Dafermos-Holzegel-Rodnianski’s in , where the linear stability of Schwarzschild is established for characteristic initial data sets.
The Schwarzschild solution of the vacuum Einstein equation in general relativity is the unique static solution that represents an isolated gravitating system of a single black hole. Studies, both theoretically and experimentally, of such a system are modeled on the Schwarzschild solution and its perturbation. The stability of the Schwarzschild solution is thus of utmost importance. More than two decades after the celebrated work of nonlinear stability of Minkowski space by Christodoulou and Klainerman , the nonlinear stability of Schwarzschild remains open. This paper addresses the linear stability of the Schwarzschild solution, which has a long history and rich literature involving the works of both physicists and mathematicians, and culminating in the recent breakthrough of Dafermos-Holzegel-Rodnianski . This paper provides a different and simpler proof that reveals the underlying geometric structure of the vacuum Einstein equation at a more elementary level.
The question of linear stability is formulated in the following way. Consider the vacuum Einstein equation , of which the Schwarzschild metric is solution. Let be a solution of the linearization of the vacuum Einstein equation at Schwarzschild:
is a smooth symmetric tensor on the background Schwarzschild spacetime. Since the Einstein equation is invariant under diffeomorphisms of the spacetime, infinitesimally any smooth co-vector field generates another solution of (1). Here is the deformation tensor of ; see Lemma 9. On the other hand, the Schwarzschild solution lies in the larger family of Kerr solutions and there are solution of (1) that correspond to Kerr perturbations. To prove the linear stability amounts to:
Firstly, formulate (1) as the initial value problem of a linear hyperbolic system of PDE’s of components of under a suitable gauge condition. Secondly, prove that the linearized metric coefficients of , after being normalized by the gauge condition, decay through a suitable foliation to a Kerr perturbation under appropriate initial conditions.
There are two main approaches to the perturbation problem:
1. Perturbation of metric coefficients. This was initiated by the seminal work of Regge-Wheeler . It was shown that (1) under a suitable gauge decouples into even-parity and odd-parity perturbations, which correspond to axial and polar perturbations in Chandrasekhar . Originally, Regge-Wheeler  discovered the master Regge-Wheeler equation in the axial or odd-parity case; more than a decade later, Zerilli  derived the eponymous equation in the polar or even-parity case. Later work by Moncrief  phrased the decoupling in terms of a gauge-invariant, connection-level quantity, which we refer to as the Zerilli-Moncrief function. See also [18, 26].
2. Perturbation in Newman-Penrose (N-P) formalism. In particular, it is known that the extreme linearized Weyl curvature components satisfy the Teukolsky equation , which can be further solved by separation of variables.
On the Schwarzschild background, there is a transformation theory, initiated by Wald  with further refinements by Aksteiner et al. , that relates the Regge-Wheeler equation and the Teukolsky equation. The Regge-Wheeler equation has a favorable potential which allows for a direct analysis, while this is not clear for the Teukolsky equation.
In the Schwarzschild and Kerr settings, the estimates and techniques from the study of the scalar wave equation, regarded as a “poor man’s” linearization of the vacuum Einstein equations, are expected to prove an essential ingredient in further progress on linear stability, with developments in linear stability playing a similar role in non-linear stability. A complete theory, including uniform boundedness and decay estimates, is now in place for scalar waves, first appearing in the seminal works of  in the Schwarzschild setting and  in sub-extremal Kerr, with contributions and refinements also appearing in [14, 16, 3, 10, 2, 19, 24, 17].
To get a complete theory of linear stability, one needs to prove such decay estimates for each component of a solution under a suitable gauge, and modulo the Kerr perturbations. Finster-Smoller  prove the decay estimates of Teukolsky equation on the Schwarzschild background (see also the Kerr case ). However, the decay estimates of the perturbed metric coefficients do not follow from this. Note the extreme linearized Weyl curvature components are gauge invariant quantities, while the decay of metric coefficients holds true only after a gauge condition. It seems that the reconstruction of metric coefficients from the extreme components of Weyl curvatures in the N-P formalism remains unsolved, with partial results appearing in . See also .
The authors of  make use of the aforementioned transformation theory and shows that a certain second derivative of extreme linearized Weyl curvature satisfies a Regge-Wheeler type equation with a favorable potential. A double null gauge is imposed to derive that all perturbed metric coefficients decay modulo the Kerr perturbation. The initial value problem is formulated on characteristic data sets.
The current paper gives a complete theory of linear stability at the level of the metric perturbation. In the space of linearized solutions of symmetric tensors , we identify the Kerr perturbation in the subspace of lower angular modes. For the higher angular modes, we work with gauge-invariant quantities at the level of perturbed connection coefficients which satisfy Regge-Wheeler type equations with favorable potentials. The Regge-Wheeler gauge and an interpolated Chandrasekhar gauge are adopted to prove the decay of all perturbed metric coefficients and complete the proof. In contrast to , the initial value problem is formulated on Cauchy data sets.
In order to identify the Kerr perturbation, we first decompose any smooth, symmetric tensor according to angular modes
see Proposition 4.
Summarizing our results, we have the following linear stability theorem on the Schwarzschild spacetime, stated in Theorem 1 and Theorem 2:
Let be a smooth, symmetric tensor on the Schwarzschild spacetime, satisfying the linearized vacuum Einstein equations (1).
For the component of , there exists a unique smooth co-vector (modulo Killing fields) on the Schwarzschild spacetime and constants , such that
where are smooth symmetric tensors that correspond to linearized Kerr solutions specified in Definition 8.
Under the same assumption for as in Theorem 1. Assuming moreover is compactly supported away from the bifurcation sphere on the time-slice , there exists a smooth co-vector such that
with the gauge-normalized solution decaying pointwise through a suitable foliation.
Another decomposition (the Hodge type decomposition) of the space of symmetric tensors is adopted to study . Any is decomposed into the closed and co-closed portions, which generalize the even-odd or axial-polar decompositions in physics literature, without any symmetry or mode assumptions. Note that belongs to the closed part while belong to the co-closed part.
For each portion of , we decouple gauge-invariant quantities satisfying Regge-Wheeler type equations (71, 74, 97); analysis of these Regge-Wheeler type equations shows that each quantity decays to zero through a suitable foliation. The rest of the proof consists of the reconstruction of components of from these quantities under suitable gauge choice of and the deduction of decay of all components of . The choice of the gauge turned out to be quite subtle. The well-known Regge-Wheeler gauge  and Chandrasekhar gauge  both use the gauge freedom to remove four components of . However, neither seem sufficient to imply the decay of the other six components under the boundary condition. At the end, the decay of all components of is achieved by imposing the Regge-Wheeler gauge and an interpolated Chandrasekhar gauge.
The paper is organized as follows. In Section 2, we present the Schwarzschild spacetimes as a family of static, spherically symmetric spacetimes satisfying the vacuum Einstein equations. In Section 3, we discuss linearized gravity about such spherically symmetric spacetimes. In particular, we discuss Hodge decomposition on the spheres of symmetry and decomposition into spherical harmonics. In Section 4 we present the well-known linearized Kerr family of solutions, along with the pure gauge solutions. Using such solutions, we treat the analysis of and the proof of Theorem 1 in Section 5. Subsequent sections deal with the analysis of the closed and co-closed portions of the remainder . In Section 6, we prove decay of the co-closed portion in the Regge-Wheeler gauge, extending results from the previous . In Section 7, we present the well-known Zerilli-Moncrief function as a gauge-invariant quantity, satisfying the Zerilli equation, the analysis of which is the subject of Section 8. In Section 9, we introduce the Chandrasekhar gauge, and prove decay of the closed solution under a suitable modification of the gauge. We summarize our results on in Section 10, wherein we prove Theorem 2.
2. The Schwarzschild Spacetime
The Schwarzschild spacetimes comprise a family of static, spherically symmetric spacetimes, parametrized by mass . Each such spacetime is vacuum; i.e., each metric satisfies the vacuum Einstein equations .
The staticity and spherical symmetry of the Schwarzschild family are encoded in a number of Killing fields. In particular, we have the static Killing field, denoted , and the rotational Killing fields, denoted , with . For convenience in what follows, we collect the rotation Killing fields in the set . Moreover, we denote by the full set of Killing fields.
Our results concern the Schwarzschild exterior region, up to and including the future event horizon. In the course of our analysis, various coordinate systems will prove useful; we enumerate them below.
The Schwarzschild exterior, not including the event horizon, is covered by a coordinate patch with . In these standard Schwarzschild coordinates, the Schwarzschild metric has the form
is the round metric on the unit sphere. Often we use the shorthand
Alternatively, we can cover this region with the Regge-Wheeler coordinates , with tortoise coordinate normalized as
such that the metric
is defined on , with on the photon sphere .
A variant of the above takes , with
now defined for coordinates satisfying In contrast with the previous two, this coordinate system covers both the exterior region and the black hole region.
Finally, we shall refer to the double-null coordinate system , with null coordinates and related to the Regge-Wheeler coordinates by
With this relation, the Schwarzschild metric takes the form
Note that each of the above coordinate systems covers only a portion of the maximally extended Schwarzschild spacetime, globally parametrized by the Kruskal coordinates . As we work only on the exterior region and event horizon, we do not employ this cumbersome coordinate system in what follows.
3. Linearized Gravity in a Spherically Symmetric Background
3.1. Spherically symmetric background
The analysis in this section applies to a spherically symmetric spacetime such that the group acts by isometry.
Let be a two-dimensional Lorentzian manifold with local coordinates . Let be the unit two-sphere with the standard Riemannian metric in local coordinates . We adopt the convention that any repeated index is summed. Each point on represents an orbit sphere, with a positive function which represents the areal radius of each orbit sphere. We consider a general spherically symmetric spacetime in local coordinates :
The following indices notations are adopted throughout the paper: for quotient indices, for spherical indices, and for spacetime indices.
For example, the metric on the Schwarzschild spacetime with coordinates has
The Christoffel symbols of a spherically symmetric spacetime are
where and are the Christoffel symbols of and , respectively.
For tensors defined on , we consider two types of differential operators, and . When applied to functions, and are just differentiation with respect to coordinate variables and , respectively. We then define
The operators can then be extended to associated tensor bundles, particularly the bundle of symmetric two-tensors on . We use the notation and for the quotient d’Alembertian and the spherical Laplacian operators. Furthermore, we denote the volume forms for the quotient space and the unit sphere by and , respectively. The operator can be used to characterize functions that are supported at lower angular modes.
A smooth function on a spherically symmetric spacetime is said to be supported at if
A symmetric tensor on a general spherically symmetric spacetime is of the form
in our calculation. Each component of depends on all spacetime variables.
Any symmetric tensor on a spherically symmetric spacetime can be decomposed into in which the components of
are characterized by the vanishing of the integrals
with respect to any function that is supported at (Definition 3).
For more details on the proposition above, refer to the later subsection on spherical harmonic decomposition. The above discussion applies to any spherically symmetric spacetime. In the following, we collect several formulae specific to the Schwarzschild spacetime that will be used in further calculations. The Christoffel symbols of the quotient metric (16) are
and in particular .
On the other hand the radial function satisfies
Th Gauss curvature of the quotient metric (16) is .
3.2. The Linearized Vacuum Einstein Equations
For an Einstein manifold with , we consider an infinitesimal solution of the linearization of . Suppose a symmetric (0,2)-tensor is a linear perturbation of . We recall that the perturbation of the Ricci curvature satisfies
On a spherically symmetric spacetime, we further decompose the last component of (19)
into trace and traceless parts. The trace is regarded as a function on , while the traceless part is a symmetric traceless two-tensor with respect to .
Expressed in this way, the linear perturbation takes the form
with each of depending upon all spacetime variables.
Perturbing about a spherically symmetric spacetime, with radial function , we compute the linearized Ricci tensor:
Specializing to the Schwarzschild spacetime, we rewrite and in alternative forms which will be used later. For , the formulas from the preceding section yield
The last line in (24) is simplified by using
and observe that
holds when the quotient spacetime is two-dimensional. We also use on Schwarzschild and the commutation formula
3.3. Hodge Decomposition
Recall that a general symmetric tensor on a spherically symmetric spacetime is of the form
in our calculation. Each component of depends on all spacetime variables. Such a decomposition is at the pointwise level. In the following, we consider a Hodge type decomposition on the spheres of symmetry.
The following Hodge type decomposition of components of of the form (28) follows from the appendix:
Let be the second component of a symmetric 2-tensor of the form (28). For each , there exist functions and on such that
Let be the last component a symmetric 2-tensor of the form (28) on that is traceless in the sense that . There exist functions , on such that
Any symmetric (0, 2) tensor of the form (28) can be decomposed as where
We note that the total number of components of remains ten, given by , , , , , and .
We refer to the portions and into which was decomposed above as the closed and co-closed portions, respectively. In the context of linearized gravity, each portion is itself a solution of (1) above. As we work within a linear theory, there is no trouble in studying, separately, the closed and co-closed pieces of a solution satisfying the linearized vacuum Einstein equations above; we simply add the two together to recover the original.
The closed and co-closed solutions generalize the even-parity (polar) and the odd-parity (axial) solutions in the physics literature, respectively. In the course of the paper, we will simply refer to them as closed and co-closed solutions.
3.4. Spherical Harmonics
The metric perturbation is comprised of objects which can be regarded as scalar, co-vector, and symmetric two-tensor quantities on the spheres of symmetry. In the course of our work, we find it necessary to decompose these objects into spherical harmonics. Owing to the Hodge decomposition of the previous subsection, wherein each of these quantities is decomposed in constituent scalar potentials, the familiar scalar spherical harmonics are seen to underly the harmonics of each object.
The scalar spherical harmonics , indexed by integers and , are eigenfunctions of the spherical Laplacian, with eigenvalue . That is, the satisfy
The eigenfunctions form a complete, orthonormal basis of . Note that we are normalizing
For the co-vectors on the sphere, we have the closed harmonics
and the co-closed harmonics
Regarding symmetric two-tensors, closed harmonics have the form
while, for the co-closed harmonics,
We note that, for the co-vectors
with support on . For the symmetric-two tensors, we have
with support on .
Subsequently, we shall find it necessary to project the linearized solution to particular harmonics. For the closed solution , we use the notation
and for the co-closed solution ,
with , and all objects on the quotient space . The spherical decomposition above is done in on the spheres of symmetry, so that, for example,
holds for each pair in the sense of .
In the decomposition of Proposition 4
is the portion of supported on the spherical harmonics , with defined similarly.
4. Linearized Kerr Solutions and Pure Gauge Solutions
4.1. Linearized Kerr Solutions
Considering the Boyer-Lindquist coordinates as an extension of the standard Schwarzschild coordinates, we write the Kerr metric in the suggestive form (28)
We treat separately the linearized change in mass and change in angular velocity below.
Linearized Change in Mass
In the expression above, linearized mass solutions have the form
giving infinitesimal change in mass within the Schwarzschild family. We can verify directly that the symmetric two-tensor satisfies the linearized vacuum Einstein equations (1). Note that linearized Schwarzschild solutions are closed solutions, supported at the lowest harmonic .
Linearized Change in Angular Velocity
Infinitesimal change in angular velocity appears in the linearized Kerr solution as
for . Again, direct computation shows that is a solution of the linearized vacuum Einstein equations (1). Such linearized Kerr solutions are co-closed solutions, supported at the harmonic .
Taken together, the above form the four dimensional family of linearized Kerr solutions.
The Kerr solutions of the linearized vacuum Einstein equation on Schwarzschild are given by:
4.2. Pure Gauge Solutions
First we have the following calculation lemma to identity pure gauge solutions.
Suppose is a co-vector on Schwarzschild with in which
Then the deformation tensor of decomposes into closed and co-closed parts,