In this paper, a direct rigorous mathematical proof of the Gregory–Laflamme instability for the five-dimensional Schwarzschild black string is presented. Under a choice of ansatz for the perturbation and a gauge choice, the linearized vacuum Einstein equation reduces to an ordinary differential equation (ODE) problem for a single function. In this work, a suitable rescaling and change of variables is applied, which casts the ODE into a Schrödinger eigenvalue equation to which an energy functional is assigned. It is then shown by direct variational methods that the lowest eigenfunction gives rise to an exponentially growing mode solution, which has admissible behavior at the future event horizon and spacelike infinity. After the addition of a pure gauge solution, this gives rise to a regular exponentially growing mode solution of the linearized vacuum Einstein equation in harmonic/transverse-traceless gauge.
I. INTRODUCTION
The main topic of this paper is the study of the stability problem for the Schwarzschild black string solution to the Einstein vacuum equation in five dimensions. In 1993, the work of Gregory–Laflamme1 gave strong numerical evidence for the presence of an exponentially growing mode instability. This phenomenon has since been known as the Gregory–Laflamme instability. This work has been widely invoked in the physics community to infer instability of many higher dimensional spacetimes, for example, black rings, ultraspinning Myers–Perry black holes, and black Saturns. The interested reader should consult Refs. 2 and 3 and references therein, as well as Ref. 4 and Refs. 5 and 6, which give a general approach to stability problems. The purpose of the present paper is to provide a direct, self-contained, and elementary mathematical proof of the Gregory–Laflamme instability of the 5D Schwarzschild black string.
A. Schwarzschild black holes, black strings, and black branes
The most basic solution to the vacuum Einstein equation
giving rise to the black hole phenomena is the Schwarzschild black hole solution (Schn, gs). It arises dynamically as the maximal Cauchy development of the following initial data: an initial hypersurface , a first fundamental form (in isotropic coordinates)
and the second fundamental form K = 0, where is the metric on the unit (n − 2)-sphere . This spacetime is asymptotically flat and spherically symmetric. The Penrose diagram in Fig. 1 represents the causal structure of (Schn, gs) arising from this initial data, restricted to the future of Σ0. The metric on the exterior (see Fig. 1) of the n-dimensional Schwarzschild black hole in traditional Schwarzschild coordinates (t, r, φ1, …, φn−2) takes the form
where t ∈ [0, ∞), , and is the metric on the unit (n − 2)-sphere.
The Penrose diagram of the Schwarzschild spacetime (Schn, gs). Here, is future null infinity, and are future timelike infinity and spacelike infinity, respectively, is the distinguished exterior region, is another exterior region, is the black hole region, is the future event horizon, and is the bifurcation sphere. The wavy line denotes a singular boundary, which is not part of the spacetime (Schn, gs) but towards which the Kretchmann curvature invariant diverges. It is in this sense that (Schn, gs) is singular. Note that every point in this diagram is, in fact, an (n − 2)-sphere.
The Penrose diagram of the Schwarzschild spacetime (Schn, gs). Here, is future null infinity, and are future timelike infinity and spacelike infinity, respectively, is the distinguished exterior region, is another exterior region, is the black hole region, is the future event horizon, and is the bifurcation sphere. The wavy line denotes a singular boundary, which is not part of the spacetime (Schn, gs) but towards which the Kretchmann curvature invariant diverges. It is in this sense that (Schn, gs) is singular. Note that every point in this diagram is, in fact, an (n − 2)-sphere.
The Lorentzian manifold that is the main topic of this paper is the Schwarzschild black string spacetime in five dimensions, which is constructed from the 4D Schwarzschild solution (Sch4, gs). Before focusing on this spacetime explicitly, it is of interest to discuss more general spacetimes constructed from the n-dimensional Schwarzschild black hole solution (Schn, gs). Let denote the circle of radius R, and let with its associated p-dimensional Euclidean metric δp. If one has the n-dimensional Schwarzschild black hole spacetime (Schn, gs) and takes its Cartesian product with Fp, then one realizes the (n + p)-dimensional Schwarzschild black brane (Schn ×Fp, gs ⊕ δp). This means that the (n + p)-dimensional Schwarzschild black brane (Schn ×Fp, gs ⊕ δp) is a product manifold made from Ricci-flat manifolds, which is again Ricci-flat and hence satisfies the vacuum Einstein equation (1.1). Note that in contrast to (Schn, gs), the spacetimes (Schn ×Fp, gs ⊕ δp) are not asymptotically flat but are called “asymptotically Kaluza–Klein.”
The Schwarzschild black brane spacetimes (Schn ×Fp, gs ⊕ δp) arise dynamically as the maximal Cauchy development of suitably extended Schwarzschild initial data, i.e., (Σ0 ×Fp, hs ⊕ δp, K = 0). Hence, the above Penrose diagram in Fig. 1 can be reinterpreted as the Penrose diagram for the Schwarzschild black brane, but instead of each point representing a (n − 2)-sphere, it represents a . In particular, the notation will be used henceforth to denote the distinguished exterior region of (Schn ×Fp, gs ⊕ δp).
Taking p = 1 gives rise to the (n + 1)-dimensional Schwarzschild black string spacetime or alternatively . The topic of the present paper is the 5D Schwarzschild black string spacetime or alternatively . The metric on the exterior in standard Schwarzschild coordinates is
where t ∈ [0, ∞), r ∈ (2M, ∞), and or .
Finally, to analyze the subsequent problem of linear stability on the exterior region up to the future event horizon , one requires a chart with coordinate functions that are regular up to this hypersurface , where now denotes the bifurcation surface. A good choice is ingoing Eddington–Finkelstein coordinates defined by
The (n + p)-dimensional Schwarzschild black brane metric becomes
B. Previous works
For a good introduction to the Gregory–Laflamme instability and the numerical result of Ref. 1, see Ref. 7. A detailed survey of the key work8 related to the present paper is undertaken in Sec. III. A brief history of the problem is presented here:
In 1988, Gregory–Laflamme examined the Schwarzschild black string spacetime and stated that it is stable.9 However, an issue in the analysis arose from working in Schwarzschild coordinates, which lead to incorrect regularity assumptions for the asymptotic solutions.
In 1993, Gregory–Laflamme used numerics to give strong evidence for the existence of a low-frequency instability of the Schwarzschild black string and branes in harmonic gauge.1
In 1994, Gregory–Laflamme generalized their numerical analysis to show instability of “magnetically-charged dilatonic” black branes10 (see Refs. 10 and 11 for a discussion of these solutions).
In 2000, Gubser–Mitra discussed the Gregory–Laflamme instability for general black branes. They conjectured that a necessary and sufficient condition for stability of the black brane spacetimes is thermodynamic stability of the corresponding black hole.12,13
In 2000, Reall,14 with the aim of addressing the Gubser–Mitra conjecture, explored further the relation between the stability of black branes arising from static, spherically symmetric black holes and thermodynamic stability of those black holes. In particular, the work of Reall argues that there is a direct relation between the “negative mode” of the Euclidean Schwarzschild instanton solution (this mode was initially identified in a paper by Gross, Perry, and Yaffe15) and the threshold of the Gregory–Laflamme instability. This idea was further explored in a work of Reall et al.,16 which extended the idea that “negative modes” of the Euclidean extension of a Myers–Perry black hole (the generalization of the Kerr spacetime to higher dimensions, see Refs. 17 and 2 for details) correspond to the threshold for the onset of a Gregory–Laflamme instability.
In 2006, Hovdebo and Myers8 used a different gauge (which was introduced in Ref. 18) to reproduce the numerics from the original work of Gregory and Laflamme. This gauge choice will be called spherical gauge and will be adopted in the present work. This work discusses the presence of the Gregory–Laflamme instability for the “boosted” Schwarzschild black string and the Emparan–Reall black ring (for a discussion of this solution, see Refs. 2, 19, and 20).
In 2010, Lehner and Pretorius numerically simulated the non-linear evolution of the Gregory–Laflamme instability; see the review21 and references therein.
In 2011, Figueras, Murata, and Reall4 put forward the idea that a local Penrose inequality gives a stability criterion. Furthermore, Ref. 4 showed numerically that this local Penrose inequality was violated for the Schwarzschild black string for a range of frequency parameters, which closely match those found in the original work of Gregory–Laflamme.1
In 2012, Hollands and Wald5 and, later in 2015, Prabu and Wald6 developed a general method applicable to many linear stability problems, which encompasses the problem of linear stability of the Schwarzschild black string exterior . References 5 and 6 are explored in detail in Sec. I E.
C. Statement of the main theorem
The purpose of this paper is to give a direct, self-contained, elementary proof of the Gregory–Laflamme instability for the 5D Schwarzschild black string.
For the statement of the main theorem, one should have in mind the Penrose diagram in Fig. 2 for the 5D Schwarzschild black string spacetime.
The Penrose diagram for the 5D Schwarzschild black string illustrating the set up for the linear instability problem. Indicated is a spacelike asymptotically flat hypersurface , which extends from spacelike infinity to intersect the future event horizon to the future of the bifurcation surface . Furthermore, or , is the black hole region, is the exterior region, is future null infinity, and is future timelike infinity. The hypersurface Σ can be expressed as Σ = {(t, r*, θ, φ, z) : t = f(r*)} such that f ∼ 1 for r* → ∞. An explicit example would be a hypersurface of constant t*, where t* = t + 2M log(r − 2M).
The Penrose diagram for the 5D Schwarzschild black string illustrating the set up for the linear instability problem. Indicated is a spacelike asymptotically flat hypersurface , which extends from spacelike infinity to intersect the future event horizon to the future of the bifurcation surface . Furthermore, or , is the black hole region, is the exterior region, is future null infinity, and is future timelike infinity. The hypersurface Σ can be expressed as Σ = {(t, r*, θ, φ, z) : t = f(r*)} such that f ∼ 1 for r* → ∞. An explicit example would be a hypersurface of constant t*, where t* = t + 2M log(r − 2M).
A way of establishing the linear instability of an asymptotically flat black hole is exhibiting a mode solution of the linearized Einstein equation (1.7), which is smooth up to and including the future event horizon and decays toward spacelike infinity and such that μ > 0.
Suppose R > 4M, then one can choose k such that there exists an integer , and therefore, h induces a smooth solution on the exterior of the Schwarzschild black string . Moreover, the initial data for such a mode solution on the exterior of have finite energy.
Hence, the exterior of the Schwarzschild black string or for R > 4M is linearly unstable as a solution of the vacuum Einstein equation (1.7), and the instability can be realized as a mode instability in harmonic/transverse-traceless gauge (1.10), which is not pure gauge.
One can construct a gauge invariant quantity, the tztz-component of the linearized Weyl tensor , which is non-vanishing for a non-trivial mode solution h with k ≠ 0 and μ ≠ 0 and exhibits exponential growth in t when μ > 0. This allows one to show that the mode solution constructed in Theorem 1.1 is not pure gauge. Hence, one expects that the above mode solution persists in any “good” gauge, not just (1.10).
The reader should note that the lower bound on the frequency parameter k should not be interpreted as ruling out the existence of unstable modes with arbitrarily long wavelengths. The lower bound on k in Theorem 1.1 results from the use of a test function in the variational argument (see Proposition 4.5 in Sec. IV C). The numerics of Gregory–Laflamme and Hovdebo–Myers1,8 both provide evidence that there are unstable modes for k arbitrarily small.
D. Difficulties and main ideas of the proof
It may seem natural to directly consider the problem in harmonic gauge since the equation of study (1.7) reduces to a tensorial wave equation
The above equation (1.11) results from the linearization of the gauge reduced non-linear vacuum Einstein equation (1.1), which is strongly hyperbolic and therefore well-posed. Equation (1.11) reduces to a system of ordinary differential equations (ODEs) under the mode solution ansatz (1.8) with (1.9). This system can be reduced to a single ODE of the form
where u = Htt, Htr, Hrr, or Hθθ, and Pμ,k(r) and Qμ,k(r) depend on μ, k, and r. However, if one insists on this decoupling, one introduces a regular singular point in the range r ∈ (0, ∞). For certain ranges of μ and k, this value occurs on the exterior , i.e., the regular singular point occurs in r ∈ (2M, ∞). In particular, this regular singularity occurs on the exterior for the numerical values of k and μ for which Gregory–Laflamme identified instability. In the original works of Gregory and Laflamme, the decoupled ODE for Htr was studied (see Refs. 1, 7, and 9)
It turns out that in looking for an instability, one can make a different gauge choice called spherical gauge. As shown in Sec. III, the linearized vacuum Einstein equation (1.7) for a mode solution (1.8) in spherical gauge can be reduced to a second-order ODE of the form (1.12), where, in contrast to harmonic/transverse-traceless gauge, Pμ,k(r) = Pk(r) and Qμ,k(r) = Qk(r) depend only on k and r. Hence, the existence of solution to ODE (1.12) becomes a simple eigenvalue problem for μ. Spherical gauge was originally introduced in Ref. 18 and has another advantage over harmonic/transverse-traceless gauge, which is that all r ∈ (2M, ∞) are ordinary points of ODE (1.12). Hence, the spherical gauge choice also avoids the issues of a regular singularity at some r ∈ (2M, ∞). However, in contrast to harmonic gauge, for this gauge choice, well-posedness is unclear. If one were trying to prove stability, then exhibiting a well-posed gauge would be key since well-posedness of the equations is essential for understanding general solutions. For instability, it turns out that it is sufficient to exhibit a mode solution of the non-gauge reduced Eq. (1.7), which is not pure gauge. One expects then that such a mode solution will persist in all “good” gauges, of which harmonic gauge is an example. The discussion of pure gauge mode solutions in spherical gauge in Sec. III C provides a proof that if k ≠ 0 and μ ≠ 0, then a mode solution in spherical gauge is not pure gauge. This can be shown directly or from the computation of a gauge invariant quantity, namely, the tztz-component of the linearized Weyl tensor, . Furthermore, it is shown that if a non-trivial mode solution in spherical gauge grows exponentially in t, then is non-zero and grows exponentially t.
An issue with spherical gauge is that mode solutions in the spherical gauge do not, in general, extend smoothly to the future event horizon , even when they represent physically admissible solutions. However, as shown in Sec. III D, one can detect what are the admissible boundary conditions at the future event horizon in spherical gauge by adding a pure gauge perturbation to the metric perturbation to try and construct a solution that indeed extends smoothly to . In fact, the pure gauge perturbation found is precisely one that transforms the metric perturbation to harmonic/transverse-traceless gauge (1.10). Hence, after also identifying the admissible boundary conditions at spacelike infinity in Sec. III D, proving the existence of an unstable mode solution to the linearized vacuum Einstein equation (1.7) that is not pure gauge is reduced to showing the existence of a solution to ODE (1.12) with μ > 0 and k ≠ 0, which satisfies the admissible boundary conditions that are identified in this work.
In this paper, ODE problem (1.12) is approached from a direct variational point of view in Sec. IV. To run a direct variational argument, the solution u of ODE (1.12) is rescaled and change of coordinates is applied. It is shown in Sec. IV A that Eq. (1.12) can be cast into a Schrödinger form
with Vk independent of μ. ODE (1.13) can be interpreted as an eigenvalue problem for −μ2; finding an eigenfunction, in a suitable space, with a negative eigenvalue will correspond to an instability. As shown in Sec. IV B, this involves assigning the following energy functional to the Schrödinger operator on the left-hand side of (1.13):
Using a suitably chosen test function, one can show that the infimum over functions in of this functional is negative for a range of k. One then needs to argue that this infimum is attained as an eigenvalue by showing that this functional is lower semicontinuous and that the minimizer is non-trivial. The corresponding eigenfunction is then a weak solution in to ODE (1.13) with μ > 0 for a range of . Elementary one-dimensional elliptic regularity implies that the solution is indeed smooth away from the future event horizon, , and therefore corresponds to a classical solution of the problem (1.13). Finally, the solution can be shown to satisfy the admissible boundary conditions by the condition that the solution lies in .
This paper is organized in the following manner. The remainder of the present section contains additional background on the Gregory–Laflamme instability. In Sec. II, linear perturbation theory is reviewed and the linearized Einstein equation (1.7) is derived. In Sec. III, the analysis in spherical gauge is presented. The decoupled ODE (1.12) resulting from the linearized Einstein equation (1.7) is derived, and it is established that the problem can be reduced to the existence of a solution to the decoupled ODE with μ > 0 and k ≠ 0 satisfying admissible boundary conditions. In Sec. IV, the proof of the existence of such a solution is presented via the direct variational method.
Appendix A contains a list of the Riemann tensor components and the Christoffel symbols for the Schwarzschild black string spacetime or . Appendix B collects results on singularities in second order ODE relevant for the discussion of the boundary conditions for the decoupled ODE (1.12). Appendix C provides a method of transforming a second order ODE into a Schrödinger equation. Appendix D collects some useful results from analysis that are needed in the Proof of Theorem 1.1. Appendix E compliments Theorem 1.1 with some stability results.
E. The canonical energy method
The reader should note that there are two papers5,6 concerning a very general class of spacetimes, which are of relevence to the stability problem for the Schwarzschild black string. In particular, it follows from Refs. 5 and 6 that there exists a linear perturbation of the Schwarzschild black string spacetime, which is not pure gauge and grows exponentially in the Schwarzschild t-coordinate. The following describes the results of these works.
In 2012, a paper of Hollands and Wald5 gave a criterion for linear stability of stationary, axisymmetric, vacuum black holes and black branes in D ≥ 4 spacetime dimensions under axisymmetric perturbations. They define a quantity called the “canonical energy” of the perturbation, which is an integral over an initial Cauchy surface of an expression quadratic in the perturbation. It can be related to thermodynamic quantities by
where M and JB are the ADM mass and ADM angular momenta in the Bth plane and A is the cross-sectional area of the horizon. Note that the right-hand side of (1.15) refers to the second variation of thermodynamic quantities. It is remarkable that the combination of these second variations is, in fact, determined by linear perturbations.
Reference 5 considers initial data for a perturbation of either a stationary, axisymmetric black hole or black brane with the following properties: (i) the linearized Hamiltonian and momentum constraints are satisfied, (ii) that δM = 0 = δJA and that the ADM momentum vanishes, and (iii) specific gauge conditions and finiteness/regularity conditions at the future horizon and infinity are satisfied. In what follows, initial data satisfying (i)–(iii) will be referred to as admissible. Hollands and Wald showed that if for all admissible initial data, then one has mode stability. The work also establishes that if there exist admissible initial data such that , then there exist admissible initial data for a perturbation, which cannot approach a stationary perturbation at late times, i.e., one has failure of asymptotic stability.
For the Schwarzschild black hole, one can take initial data, which corresponds simply to a change of the mass parameter M ↦ M + α, and therefore, by Eq. (1.15) and since the cross-sectional area of the horizon is given by A = 16π(M + α)2, it follows that . This is the “thermodynamic instability” of the Schwarzschild black hole. However, the initial data for a change of mass perturbation is manifestly not admissible (the family of Schwarzschild black holes is, after all, dynamically stable).
The work of Hollands and Wald5 also shows an additional result relevant specifically to the problem of stability of black branes. Suppose that there exist initial data for a perturbation of the ADM parameters of a black hole such that . Reference 5 shows that starting from such a perturbation of the black hole, one can infer the existence of admissible initial data, which depend on a parameter l, for a perturbation (which is not pure gauge) of the associated black brane such that again . One should note that this argument does not give an explicit bound on l. This criterion formalized a conjecture by Gubser–Mitra that a necessary and sufficient condition for stability of the black brane spacetimes is thermodynamic stability of the corresponding black hole.12,13 Since the change of mass perturbation of the Schwarzschild black hole produces , this argument implies that the Schwarzschild black string fails to be asymptotically stable.
The reader should note that the Hollands and Wald paper5 also showed that a necessary and sufficient condition for stability, with respect to axisymmetric perturbations, is that a “local Penrose inequality” is satisfied. The idea that a local Penrose inequality gives a stability criterion was originally discussed in the work of Figueras, Murata, and Reall,4 which gave strong evidence in favor of sufficiency of this condition for stability. Furthermore, Ref. 4 showed numerically that this local Penrose inequality was violated for the Schwarzschild black string for a range of frequency parameters that closely match those found in the original work of Gregory–Laflamme.1
The failure of asymptotic stability does not in itself imply that perturbations grow. However, the results of Ref. 5 were strengthened in 2015 by Prabhu and Wald.6 They showed, using some spectral theory, that if there exist admissible initial data for a perturbation such that for a black brane, then there exists initially well-behaved perturbations that are not pure gauge and that grow exponentially in time. Having established that there exist admissible initial data for a perturbation such that for the Schwarzschild black string in Ref. 5, the existence of a linear perturbation, which is not pure gauge and has exponential growth, follows.
The present work differs from the above as it gives a direct, self-contained, elementary proof of the Gregory–Laflamme instability following the original formulation of Refs. 1 and 7–9, which is completely explicit. In particular, it gives an exponentially growing mode solution with an explicit growth rate of the form defined by equations (1.8) and (1.9) in harmonic/transverse-traceless gauge, which is not pure gauge.
It would also be of interest to see if Theorem 1.1 in the form stated could be inferred from the canonical energy method of Hollands, Wald, and Prabu5,6 in an explicit way bypassing some of the functional calculus applied there. In particular, it would be interesting to explore the possible relation between the variational theory applied to and that applied here (see Sec. IV B).
F. Outlook
This paper brings together what is known about the Gregory–Laflamme instability as well as providing a direct elementary mathematically rigorous proof of its existence without the use of numerics and with an explicit bound on μ and k. Note that while only the 5D Schwarzschild black string was considered here, the result of instability readily extends to higher dimensions with the replacement of kz in the exponential factor with ∑ikizi.
Further directions of work could be to study the non-linear problem, the extension to or , the extension to charged black branes of the work,10 and the extension to black rings or ultraspinning Myers–Perry black holes.
G. Contextual remarks
1. Motivation for the study of higher dimensions
The study of higher dimensions merits a few words of motivation since, from a physical standpoint, only 3 + 1 are perceived classically. First, from a purely mathematical perspective, it is of interest to see how general relativity differs in higher dimensions from the 4D case. This throws light on how general Lorentzian manifolds obeying the vacuum Einstein equation (1.1) behave. Second, the physics community is very interested in higher-dimensional gravity from the point of view of string theory. Understanding how general relativity behaves in higher dimensions is therefore of relevance to the low energy limit of string theory.2
2. Some differences in higher dimensions
In higher dimensions, many results from 4D general relativity no longer hold. As shown by Hawking, in 4D, the cross sections of the event horizon of an asymptotically flat stationary black hole spacetime must be topologically (under the dominant energy condition).22 In higher dimensions, it is possible to construct explicit examples of black hole spacetimes with non-spherical cross-sectional horizon topology. For example, the black ring solution with horizon topology is .19 In higher dimensions, there also exists a generalized Kerr solution known as the Myers–Perry black hole,17 which has cross-sectional horizon topology . Hawking’s theorem has been generalized to higher dimensions,23 which shows that the horizon topology must be of positive scalar curvature. In 5D, under the assumptions of stationarity, asymptotic flatness, two commuting axisymmetries and “rod structure” black holes are unique, and further the horizon topology is either , , or lens space.24
In 4D, it is conjectured that maximal developments of “generic” asymptotically flat initial data sets can asymptotically be described by a finite number of Kerr black holes. This “final state conjecture” cannot generalize immediately since there exist at least two distinct families of black hole solutions that can have the same mass and angular momentum: the Myers–Perry black hole and the black ring. Moreover, there exist distinct black ring solutions with the same mass and angular momentum.2,20 The final state conjecture may need to be modified to include the property of stability.
3. Related works
A few other works are of relevance to this discussion. The review paper2 and book chapter20 discuss the black ring solution19 in great detail. This relates to the work presented here since the Gregory–Laflamme instability is often heuristically invoked when discussing higher-dimensional black hole solutions. In particular, if the black ring of study has a large radius and is sufficiently thin, then it “looks like” a Schwarzschild black string and therefore would be susceptible to the Gregory–Laflamme instability. There have been heuristic and numerical results to give evidence to this claim.8,25 Finally, note that in 2018, Ref. 26 produced the first mathematically rigorous result on the stability problem for the black ring spacetime.
II. LINEAR PERTURBATION THEORY
A. Linearized vacuum Einstein equation
Consider a Lorentzian manifold (M, g) with metric satisfying the vacuum Einstein equation
In this section, a “perturbation” of the spacetime metric will be discussed. This will be represented by a new metric of the form g + ϵh with ϵ > 0. h here is a symmetric bilinear form on the fibers of TM. In the following, a series of results on how various quantities change to O(ϵ) (the linear level) are derived. This will result in an expression for the Ricci tensor under such a perturbation to linear order.
An important point to note is that indices are raised and lowered here with respect to g.
Direct computation.
□
If one assumes that g satisfies the vacuum Einstein equation (1.1) and g + ϵh satisfies the vacuum Einstein equation (1.1) to O(ϵ), then it follows from Proposition 2.1 that h must satisfy
to O(ϵ). In what follows, Eq. (2.4) will be called the linearized vacuum Einstein equation. This will be the main equation of interest, with g being the Schwarzschild black string metric,
B. Pure gauge solutions in linearized theory
The vacuum Einstein equation (1.1) is a system of second order quasilinear partial differential equations of the pair (M, g), which are invariant under the diffeomorphisms of M. This means that for given initial data, the vacuum Einstein equation (1.1) only determines a spacetime unique up to diffeomorphism, i.e., if there exists a diffeomorphism Φ : M → M, then (M, g) and (M, Φ*(g)) are equivalent solutions of the vacuum Einstein equation (1.1). For constructing spacetimes, one often imposes conditions on local coordinates called a gauge choice. For linearized theory, this can be formulated as follows.
Consider a Lorentzian manifold with ϵ > 0. Let {Φτ} be a one-parameter family of diffeomorphisms generated by a vector field X and define ξ ≔ τX ∈ TM. Then, from the definition of the Lie derivative, one has
if one treats . Hence, in the context of linearized theory, one considers two solutions to the linearized vacuum Einstein equation (2.4), h1 and h2, as equivalent if
for some vector field ξ ∈ TM.
Showing that a solution h to the linearized vacuum Einstein equation (2.4) is not pure gauge is tantamount to showing that h is not equivalent to the trivial solution. It is thus essential that the solution constructed in this paper not be pure gauge. The following propositions establish that the tztz-component of the linearized Weyl tensor is invariant under gauge transformation. This means that if is non-zero for a solution h to the linearized vacuum Einstein equation (2.4), then h cannot be pure gauge.
Direct computation.
□
For the 5D Schwarszchild black string, evaluated on a pure gauge solution vanishes.
□
III. ANALYSIS IN SPHERICAL GAUGE
In this section, a mode solution, h, of the linearized vacuum Einstein equation (2.4) on the exterior of the Schwarzschild black string spacetime or is considered. One makes the additional assumption that this mode solution preserves the spherical symmetry of Sch4. Hence, in particular, the solution can be expressed in (t, r, θ, φ, z) coordinates as
where α, β ∈{t, r, θ, φ, z}. Moreover, in search of instability, the most interesting case for the present work is μ > 0.
This section contains the analysis of the ODEs resulting from the linearized Einstein vacuum equation (2.4) for a mode solution of the form (3.1) when it is expressed in spherical gauge.
The terminology “spherical gauge” is motivated by the fact that a mode solution of this form preserves the area of the spheres of the original spacetime.
First, it is shown in Sec. III A that one can impose the gauge consistently at the level of modes, i.e., if there is a mode solution of the form (3.1), with μ ≠ 0 and either k ≠ 0 or , then there is a mode solution of the form (3.2) differing from the original one by a pure gauge solution. In the case where Htz = 0, Hrz = 0, and Hzz = 0, this consistency condition is already implicit in Refs. 18 and 8. In Sec. III B, the original decoupling of the ODEs resulting from the linearized vacuum Einstein equation (2.4) and the spherical gauge ansatz (3.2) is reproduced from Ref. 8. This decoupling results in a single ODE for the component Hz(r) in Eq. (3.2). It is then shown, in Sec. III C, that if k ≠ 0 and μ ≠ 0, then mode solutions in spherical gauge (3.2) are not pure gauge. This is proved by examining the tztz-component of the linearized Weyl tensor associated with a mode solution in spherical gauge, which is gauge invariant by Proposition 2.3. In this section, it is also proved that if a non-trivial mode solution in spherical gauge has μ > 0 (i.e., it grows exponentially in t) and k ≠ 0, then is non-zero and also grows exponentially. By the gauge invariance of , this behavior will persist in all gauges. Next, in Sec. III D, the admissible boundary conditions for the solution at the future event horizon and finiteness conditions at spacelike infinity are identified. Note this issue is subtle since, in general, both “basis” elements for a mode solution h of the form (3.2) are, in fact, singular at the future event horizon in this gauge. By adding a pure gauge perturbation, the admissible boundary conditions for the solution h in the form (3.2) can be identified. Moreover, this pure gauge solution can be chosen such that, after adding it, the harmonic/transverse-traceless gauge (1.10) conditions are satisfied. Finally, in Sec. III E, the problem of constructing a linear mode instability of the form (3.1) is reduced to show that there exists a solution to the decoupled ODE for Hz(r), with μ > 0 and k ≠ 0, that satisfies the admissible boundary conditions at the future event horizon and spacelike infinity (see Proposition 3.8).
A. Consistency
In Ref. 8, it is stated that any mode solution of the form in Eq. (3.1) with Htz = 0, Hrz = 0, and Hzz = 0 can be brought to the spherical gauge form (3.2) by the addition of a pure gauge solution. Slightly more generally, one, in fact, has the following proposition:
(Consistency of the Spherical Gauge). Consider a mode solution h to the linearized Einstein vacuum equation (2.4) on the exterior of the Schwarzschild black string spacetime or of the form (3.1) with μ ≠ 0. Further suppose that either k ≠ 0 or . Then, there exists a pure gauge solution hpg such that h + hpg is of form (3.2). It is in this sense that the spherical gauge (3.2) can be consistently imposed on the exterior of the Schwarzschild black string or .
From Sec. II B, a pure gauge solution is given by for a vector field ξ. Hence, is the new mode solution. Consider a diffeomorphism generating vector field of the form ξ = eμt+ikz(ζt(r), ζr(r), 0, 0, ζz(r)).
□
B. Reduction to ODE
Under a spherical gauge ansatz (3.2) with μ ≠ 0 and k ≠ 0, the linearized vacuum Einstein equation (2.4) reduces to a system of coupled ODEs for the components Ht, Hv, Hr, and Hz. This system of ODEs can be decoupled to the single ODE for ,
with
The following proposition establishes this decoupling of the linearized vacuum Einstein equation (2.4) to ODE (3.5) and the construction of a mode solution h in spherical gauge (3.2) from a solution to ODE (3.5).
Given a mode solution h in spherical gauge (3.2) with μ ≠ 0 and k ≠ 0 on the exterior of the Schwarzschild black string or , ODE (3.5) is satisfied by hzz. Conversely, given a C2((2M, ∞)) solution to ODE (3.5) with k ≠ 0 and μ ≠ 0, one can construct a mode solution h in spherical gauge (3.2) to the linearized vacuuum Einstein equation (2.4) on the exterior of the Schwarzschild black string . If , then h induces a mode solution on .
□
A mode solution h in spherical gauge with Ht(r) ≡ 0 satisfying the linearized vacuum Einstein equation (2.4) on the exterior of the Schwarzschild black string is then again equivalent to the system of ODE (3.8)–(3.14) [with k = 0 and Ht ≡ 0] being satisfied. Equations (3.8) and (3.10) are automatically satisfied by k = 0. Eq. (3.12) automatically gives the decoupled equation (3.5) for Hz. Then, Eq. (3.9) can be solved for Hr in terms of Hz and . This gives the relation in Eq. (3.15) for Hr with k = 0. Equation (3.11) can be used to solve for Hv in terms of Hz and . At this point, the equations (3.13) and (3.14) are automatically satisfied. Therefore, again a solution to ODE (3.5) induces a mode solution in spherical gauge with Ht = 0.
C. Excluding pure gauge perturbations
This section contains two proofs that if k ≠ 0 and μ ≠ 0, then a non-trivial mode solution h of the form (3.2) cannot be a pure gauge solution. One can prove this directly via the following proposition:
□
Perhaps more satisfactorily, one can establish that if h is a non-trival mode solution in spherical gauge (3.2) with k ≠ 0 and μ ≠ 0, then the tztz-component of the linearized Weyl tensor is non-vanishing. Moreover, if h has μ > 0, then grows exponentially. Since is gauge invariant, this behavior persists in all gauges. More precisely, one has the following proposition:
Suppose k ≠ 0, μ ≠ 0, and h is a non-trivial mode solution in spherical gauge (3.2) of the linearized vacuum Einstein equation (2.4) on the exterior of the Schwarzschild black string or . Then, is non-vanishing and h is not pure gauge. Moreover, if μ > 0 then also grows exponentially.
□
D. Admissible boundary conditions
One can construct two sets of distinguished solutions to ODE (3.5) associated with the “end points” of the interval (2M, ∞). Note that by Definition B.1 from Appendix B, r = 2M is a regular singularity, as 2M is not an ordinary point and
are analytic near r = 2M. By Definition B.3, ODE (3.5) has an irregular singularity at infinity since there exist convergent series expansions
in a neighborhood of infinity with p0 = 0, p1 = −4, q0 = −(k2 + μ2), and q1 = −2M(k2 + 2μ2). The asymptotic analysis of the ODEs around these points is examined in Secs. III D 1 and III D 2. This analysis of ODE (3.5) near r = 2M and r = ∞ will lead to the identification of the admissible boundary conditions for a mode solution h in spherical gauge (3.2) of the linearized Einstein vacuum equation (2.4).
1. The future event horizon
The goal of this section is to identify the admissible boundary conditions for a solution to ODE (3.5) near r = 2M. This requires one to understand the behavior near r = 2M of the mode solution h in spherical gauge (3.2) of the linearized vacuum Einstein equation (2.4), which results (through the construction in Proposition 3.2) from .
Associated with the future event horizon , there exists a basis for solutions to ODE (3.5). From , one can examine the behavior near r = 2M of any mode solution h in spherical gauge (3.2) with μ ≠ 0 and k ≠ 0 through Proposition 3.2. A mode solution h in spherical gauge (3.2) with μ > 0 and k ≠ 0 constructed from never smoothly extends to the future event horizon. A mode solution h in spherical gauge (3.2) with μ > 0 and k ≠ 0 constructed from also does not smoothly extend to the future event horizon unless μ satisfies particular conditions. However, if h is a mode solution in spherical gauge (3.2) with μ > 0 and k ≠ 0 constructed from , then after the addition of a pure gauge solution hpg, it turns that out one can smoothly extend h + hpg to the future event horizon. Moreover, it will be shown that h + hpg satisfies the harmonic/transverse-traceless gauge (1.10) conditions. This will be the content of Proposition 3.5.
First, some preliminaries: The coefficients of ODE (3.5) extend meromorphically to r = 2M and behave asymptotically as
Hence, one may write ODE (3.5) as
From Appendix B, the indicial equation associated with the ODE (3.38) is
which has roots
If , then one can deduce from Theorem B.1 the asymptotic basis for solutions near r = 2M. If , then the relevant result for the asymptotic basis of solutions is Theorem B.2. Combining the results of Theorems B.1 and B.2, one has the following basis for solutions for μ > 0:
where the coefficients , , and the anomalous term CN can be calculated recursively (see Theorems B.1 and B.2). A general solution to ODE (3.5) will be of the form
with .
If 4Mμ is not an integer or 4Mμ is an integer and CN = 0, then the asymptotic basis for solutions for μ > 0 reduces to
The main result of this section is the following:
Suppose μ > 0, k ≠ 0, and let be a solution to ODE (3.5). Let h be the mode solution on the exterior of the Schwarzschild black string constructed from in Proposition 3.2. Then, there exists a pure gauge solution hpg such that h + hpg extends to a smooth solution of the linearized vacuum Einstein equation (2.4) at the future event horizon if k2 = 0, where k2 is defined in Eq. (3.43). Moreover, h + hpg can be chosen to satisfy the harmonic/transverse-traceless gauge (1.10) conditions.
Before proving the statement of Proposition 3.5, it is useful to prove the following lemma:
□
□
First, if 4Mμ is a positive integer and the coefficient CN does not vanish, then by Eq. (3.89), the basis element has an essential logarithmic divergence and is therefore always singular at the future event horizon .
Note that the full mode solution h constructed from Proposition 3.2 involves a factor of eμt = eμve−μr(r − 2M)−2Mμ, so after multiplication by this exponential factor, one can see that the basis elements are always singular, i.e., a solution with k2 ≠ 0 is always singular at the future event horizon. The components and are unconditionally smooth. However, in general, the components and remain singular at the future event horizon unless 4Mμ = 1 or or −2 + 2Mμ > 2. [In Appendix E, it is shown that for existence of a solution with μ > 0, which has k2 = 0 and is finite at infinity (see Sec. III D 2), then .] Hence, neither basis perturbation h± in spherical gauge (3.2) extends, in general, smoothly across the future event horizon .
2. Spacelike infinity
The goal of this section is to identify the admissible boundary conditions for a solution to ODE (3.5) as r → ∞. This requires one to understand the behavior as r → ∞ of the mode solution h in spherical gauge (3.2) of the linearized vacuum Einstein equation (2.4), which results (through the construction in Proposition 3.2) from .
In this section, a basis for solution associated with r → ∞ is constructed. This basis captures the asymptotic behavior of any solution to ODE (3.5) as r → ∞. In particular, as r → ∞, grows exponentially and decays exponentially. It will be shown that after the addition of the pure gauge solution hpg defined in equations (3.64) and (3.65), h + hpg is a mode solution in harmonic/transverse-traceless gauge (1.10) to the linearized Einstein vacuum equation, which is a linear combination of solutions that grow or decay exponentially as r → ∞. The admissible boundary condition will be that the solution should decay exponentially, from which it will follows that .
One should note that the functions Pk(r) and admit convergent series expansions in a neighborhood of r = ∞,
with p0 = 0, p1 = −4, q0 = −(k2 + μ2), and q1 = −2M(k2 + 2μ2). Therefore, r = ∞ is an irregular singular point of ODE (3.5) according to the discussion of Appendix B. Eqs. (B18) and (B19) from Appendix B give
From Theorem B.3, there exists a unique basis for solutions to ODE (3.5) satisfying
Therefore, a general solution will be of the form
with .
Let be a solution to ODE (3.5). Let h be the mode solution to the linearized vacuum Einstein equation (2.4) in spherical gauge (3.2) associated with the solution , and let hpg be the pure gauge solution defined by Eqs. (3.64) and (3.65) such that h + hpg satisfies the harmonic/transverse-traceless gauge (1.10) conditions. Then, the solution h + hpg to ODE (3.5) decays exponentially towards spacelike infinity if c1 = 0, where c1 is defined by Eq. (3.100).
□
E. Reduction of the proof of theorem 1.1
This section summarizes Propositions 3.2–3.5 and 3.7 to give a full description of the permissible asymptotic behavior of a mode solution h in spherical gauge (3.2), which is not pure gauge. This provides a reduction of Theorem 1.1 to proving that there exists a solution to ODE (3.5), which has μ > 0, k ≠ 0, and obeys the admissible boundary conditions: k2 = 0 and c1 = 0.
Under the additional assumption that , the mode solution h defined above can be interpreted as a mode solution to the linearized vacuum Einstein equation (2.4) on the exterior of the Schwarzschild black string . Hence, if , the above statement applies to the exterior of .
Section IV (see Proposition 4.1) will prove the existence of a solution to ODE (3.5) satisfying the properties of Proposition 3.8. In particular, for all , a solution to ODE (3.5) with , k2 = 0, and c1 = 0 is constructed. If R > 4M, then there exists an integer . Hence, one can choose k such that the constructed gives rise to a mode solution on . Moreover, on , h will manifestly have finite energy in the sense that and are finite. (Note that on , h will not have finite energy due to the periodic behavior in z on .) Thus, Theorem 1.1 follows from Propositions 3.8 and 4.1.
IV. THE VARIATIONAL ARGUMENT
By Proposition 3.8, the Proof of Theorem 1.1 has now been reduced to exhibiting a solution to (3.5) with μ > 0, k ≠ 0, k2 = 0, and c1 = 0. This section establishes the required proposition thus completing the proof.
For all , there exists a C∞((2M, ∞)) solution to ODE (3.5) with μ > 0, and in the language of Proposition 3.8, k2 = 0 and c1 = 0.
In order to exhibit such a solution to ODE (3.5), it is convenient to rescale the solution and change coordinates in ODE (3.5) so as to recast as a Schrödinger equation for a function u. This transformation is given in Sec. IV A. In Sec. IV B, an energy functional is assigned to the resulting Schrödinger operator. With the use of a test function (constructed in Sec. IV C), a direct variational argument can be run to establish that for , there exists a weak solution with such that μ > 0. The Proof of Proposition 4.1 concludes by showing that the solution u is indeed smooth for r ∈ (2M, ∞) and satisfies the conditions of Proposition 3.8, i.e., k2 = 0 and c1 = 0.
A. Schrödinger reformulation
To reduce the number of parameters in ODE (3.5), one can eliminate the mass parameter with , , and to find
with
Following Proposition C.1 from Appendix C, one can now transform Eq. (4.1) into the regularized Schrödinger form by introducing a weight function and changing coordinates to . This will produce a Schrödinger operator with a potential, which decays to zero at the future event horizon and tends to the constant at spatial infinity. From Proposition C.1, the weight function must satisfy the ODE,
The desired solution for the weight function is
ODE (4.1) becomes
where can be found from Eq. (C10) to be
where x is understood as an implicit function of x*.
As a trivial consequence of Proposition 3.8 in Sec. III D on asymptotics of the solution to ODE (3.5), one has the following proposition for the asymptotics of the Schrödinger equation (4.6).
B. Direct variational argument
This section establishes a variational argument which, will be used to infer the existence of a negative eigenvalue to the Schrödinger operator in Eq. (4.6).
□
□
C. The test function and existence of a minimizer
ODE (4.6) is now in a form where a direct variational argument can be used to prove that there exists an eigenfunction of the Schrödinger operator associated with the left-hand side of ODE (4.6) with a negative eigenvalue, i.e., . The following proposition constructs a suitable test function such that it is in the correct function space, , and, for all , implies that the infimum of the energy functional in Eq. (4.26) is negative. (As will be apparent, the negativity is inferred via complicated but purely algebraic calculations.)
Define , where x is an implicit function of x*, n is a finite non-zero natural number, and define E and E0 as in Eq. (4.26) of Corollary 4.4. Then, and for n = 100 and , .
□
D. Proof of proposition 4.1
To prove Proposition 4.1, one can clearly reformulate as follows:
For all , there exists a solution to the Schrödinger equation (4.6) with , and in the language of Proposition 4.2, and .
From the regularity Theorem D.3, any , which weakly solves the Schödinger equation (4.6), is, in fact, smooth. Therefore, for all , there exists a solution to the Schrödinger equation (4.6) with .
Therefore, taking and gives a solution to the Schrödinger equation (4.6) with , , and .
□
ACKNOWLEDGMENTS
First and foremost, I would like to express my gratitude to my supervisor Professor Mihalis Dafermos for introducing me to this project, for his guidance, and his comments on this manuscript. In particular, I would also like to thank Harvey Reall for his comments on this manuscript and Claude Warnick, Pierre Raphael, Yakov Shlapentokh-Rothman, Christoph Kehle, Rita Teixeira da Costa, and Renato Velozo-Ruiz for many useful discussions. I would like to thank an anonymous reviewer for their useful comments on this manuscript. In particular, I would like to thank them for suggesting that the tztz-component of the linearized Weyl tensor could be invariant under gauge transformation. Finally, I would like to thank Kasia for her support. I am supported by the UK Engineering and Physical Sciences Research Council (EPSRC) Grant No. EP/L016516/1 for the University of Cambridge Centre for Doctoral Training, the Cambridge Centre for Analysis.
DATA AVAILABILITY
Data sharing is not applicable to this article as no new data were created or analyzed in this study.
APPENDIX A: CHRISTOFFEL AND RIEMANN TENSOR COMPONENTS FOR THE 5D SCHWARZSCHILD BLACK STRING
To compute □ghab, one requires the Christoffel symbols and the Riemann tensor components; the non-zero Christoffel symboles are listed as follows:
The others are obtained from symmetry of lower indices. Note, . Hence, the Riemann tensor components that are relevant are the ones with spacetime indices μ ∈ {0, …, 3}, which are just the usual Schwarzschild Riemann tensor components; the non-zero ones are listed below for completeness,
Any others can be found from the symmetry.
APPENDIX B: SINGULARITIES IN SECOND ORDER ODE
This section is heavily based on the book of Olver.31 In particular, see Chap. 5 Secs. IV and V and Chap. 7 Sec. II.
The singular behavior of z = ∞ is determined by making the change of variables in ODE (B1). This case will be considered explicitly in Appendix B 2.
In the following, general results for ODE are presented.
1. Regular singularities
In this paper, solutions of a second order ODE in a neighborhood |z − z0| < r of a regular singular point are required. The classical method is to search for a convergent series solution in such a neighborhood.
The following two theorems deal with the asymptotic behavior of solutions in the neighborhood of a regular singularity.
If the roots of the indicial equation do not differ by an integer, then Theorem B.1 gives a basis of solutions for the ODE in a neighborhood of the singular point. Equation (B10) determines the coefficients of the series expansion recursively from an arbitrarily assigned a0 ≠ 0, which can be taken to be 1. This process runs into difficulty if, and only if, the two roots differ by a positive integer. To see this, let α+ be the root of the indicial equation with the largest real part, the other root is then α+ − s for some . Then, since I((α+ − s) + s) = 0, one cannot determine as via Eq. (B10) for this power series. In this case, one solution can be found with the above method by taking the root of the indicial equation with the largest real part.
The following theorem investigates the case where the roots differ by an integer. Let α+ be the root of the indicial equation with the largest real part, and the other root is then α+ − s for some .
2. Irregular singularities
This section summarizes the key result for constructing a basis of solutions to ODE (3.5) associated with r → ∞. [The results presented can, in fact, be applied to any irregular singular point of an ODE (B1) since without loss of generality, the irregular singularity can be assumed to be at infinity after a change of coordinates.] The following definition makes precise the notion of an irregular singularity at infinity.
The main Theorem B.3 of this section can be motivated by the following discussion. Consider a formal power series
Substituting the expansions into the ODE and equating coefficients yield
and
Now, Eq. (B15) has two roots,
These give rise to
The two values of a0, can be, without loss of generality, set to 1, and the higher order coefficients determined iteratively from Eq. (B17) unless one is in the exceptional case where (for further information on this case, see Sec. I C of Chap. 7 in Ref. 31). The issue that arises is that in most cases, the formal series solution (B14) does not converge. However, the following theorem characterizes when (B14) provides an asymptotic expansion for the solution for sufficiently large |z|.
APPENDIX C: TRANFORMATION TO SCHRÖDINGER FORM
□
APPENDIX D: USEFUL RESULTS FROM ANALYSIS
1. Sobolev embedding
2. The multiplication operator is compact
Let with lim|x|→∞q(x) = 0 and s > 0. Then, Mq : u → qu is a compact operator from to .
□
3. A regularity result
(Regularity for the Schrödinger Equation). Let be a weak solution of the equation (−Δ + V)u = λu, where V is a measurable function and . Then, if V ∈ C∞(Ω) with open, not necessarily bounded, then u ∈ C∞(Ω) also.
(Ref. 33, Vol. II, p. 55). Note one can argue this from standard elliptic regularity results and Sobolev embeddings. In this paper, only the one-dimensional case of this is applied, which is completely elementary.
□
APPENDIX E: A RESULT ON STABILITY IN SPHERICAL GAUGE
This section contains a few technical results on where the instability may lie in frequency space. This helped guide the search for a suitable test function and the subsequent instability.
(regions of stability in frequency space). Let μ > 0 and k ≠ 0. There does not exist a solution of ODE (3.5) with c1 = 0, k2 = 0, and or .
□
By an almost identical argument, one can make the bound for even sharper and show that and .