Higher Order Linear Stability and Instability of Reissner-Nordstr\"om's Cauchy Horizon
Jo\~ao L. Costa, Pedro M. Gir\~ao

TL;DR
This paper analyzes the stability of wave solutions near the Cauchy horizon of Reissner-Nordström black holes, establishing criteria for regularity and blow-up based on surface gravities, with implications for black hole interior stability.
Contribution
It provides new criteria relating surface gravities to wave regularity and blow-up at the Cauchy horizon in Reissner-Nordström spacetimes.
Findings
Criteria for wave regularity up to the Cauchy horizon.
Conditions for wave blow-up in $C^1$ and $H^1$.
Applicability to various Reissner-Nordström black hole models.
Abstract
We consider smooth solutions of the wave equation, on a fixed black hole region of a subextremal Reissner-Nordstr\"om (asymptotically flat, de Sitter or anti-de Sitter) spacetime, whose restrictions to the event horizon have compact support. We provide criteria, in terms of surface gravities, for the waves to remain in , , up to and including the Cauchy horizon. We also provide sufficient conditions for the blow up of solutions in and .
Peer Reviews
No public reviews on file for this paper yet. If you reviewed it on a platform where reviews are public (OpenReview, ICLR, NeurIPS, ICML), you can paste yours below so the community can read it here.
Videos
No videos yet. Explain this paper in a talk, walkthrough, or lecture? Add one.
**Higher Order Linear Stability and Instability of Reissner-Nordström’s Cauchy Horizon
**
João L. Costa*†‡* and Pedro M. Girão‡
†ISCTE - Instituto Universitário de Lisboa, Lisboa, Portugal.
‡Center for Mathematical Analysis, Geometry and Dynamical Systems,
Instituto Superior Técnico, Universidade de Lisboa,
Av. Rovisco Pais, 1049-001 Lisbon, Portugal.
Abstract.
We consider smooth solutions of the wave equation, on a fixed black hole region of a subextremal Reissner-Nordström (asymptotically flat, de Sitter or anti-de Sitter) spacetime, whose restrictions to the event horizon have compact support. We provide criteria, in terms of surface gravities, for the waves to remain in , , up to and including the Cauchy horizon. We also provide sufficient conditions for the blow up of solutions in and .
Key words and phrases:
Wave equation, black holes, positive cosmological constant
2010 Mathematics Subject Classification:
Primary: 35L05; Secondary: 35R01, 58J45, 83C57
1. Introduction
Cauchy horizons are the spacetime boundary of the maximal Cauchy development of initial value problems for the Einstein field equations. Whenever non-empty, their existence and stability puts into question global uniqueness, and consequently challenges the deterministic character of General Relativity. To understand how perturbations of a static charged black hole behave at the Cauchy horizon that lies in its interior, we will study solutions of the wave equation on the black hole region of fixed subextremal Reissner-Nordström (asymptotically flat, de Sitter or anti-de Sitter) spacetimes. In this framework, it is natural to consider that Cauchy horizons that allow solutions with higher regularity are more stable than the ones that do not.
The stability of Cauchy horizons is a classical problem in General Relativity and, in recent years, considerable progress has been made in its understanding through the mathematical analysis of wave equations. Stability results can be found in [12, 24, 3, 13, 14, 19, 16, 17] and instability results in [20, 23, 9, 10], and the references therein. For developments concerning the analysis of the full Einstein equations we refer to [4, 5, 6, 7, 21, 22, 8, 25].
Most of the literature about the linear problem focuses on stability-regularity at the and levels, in line with the modern formulations of the Strong Cosmic Censorship Conjecture. There are however some notable exceptions. In [13], Gajic provides criteria for the and extendibility of spherically symmetric waves on (asymptotically flat) extremal black holes. In the subextremal de Sitter setting, Hintz and Vasy [17] have shown that solutions of the wave equation arising from smooth Cauchy data have regularity up to the Cauchy horizon, with the degree of regularity being dictated by , the spectral gap of the operator (which also controls the decay rate of solutions along the event horizon), and , the Cauchy horizon’s surface gravity. However, recent numerical computations of the spectral gap [2] suggest that the regularity never exceeds .
Here we present criteria for higher order linear stability of the Cauchy horizon, meaning with , in a subextremal Reissner-Nordström spacetime, as well as criteria for linear instability, in both and . We will achieve this by considering waves, without symmetry assumptions, whose restrictions to the event horizon have compact support. Although, in view of the results in [11, 1, 18], this behavior on the event horizon cannot arise from generic Cauchy data, it provides a class of bona fide characteristic initial value problems for the wave equation. We will show that an arbitrarily high regularity at the Cauchy horizon can be obtained by increasing the order to which the wave vanishes in a direction transverse to the event horizon. Moreover, for this initial value problem, the role of the surface gravities in determining the degree of stability of the Cauchy horizon becomes particularly transparent. For instance, we will prove that, as a consequence of a well known relation between surface gravities, if the wave only vanishes to zeroth order at the event horizon then, in spite of having compact support on the event horizon, it cannot be extended in to any neighborhood of any point on the Cauchy horizon. In particular, this shows that we cannot expect to obtain arbitrarily high regularity for waves up to and including the Cauchy horizon by simply increasing their decay rate along the event horizon.
1.1. Statement of the main results
Let us set some basic terminology and notation. Let be a connected component of the black hole region of a subextremal Reissner-Nordström (asymptotically flat, de Sitter or anti-de Sitter) spacetime. Denote by and the surface gravities of the future event horizon and the future Cauchy horizon , respectively, and let and denote the “right side” components of these horizons (see Figure 1). Let be a future increasing affine parameter of the generators of , constant on each symmetry sphere, and let denote an ingoing null hypersurface that intersects , at . Letting be a smooth vector field which is tangent to and transverse to , we will say that vanishes to order at if
[TABLE]
We are interested in properties of functions that belong to the space
[TABLE]
for a fixed and some .
We may now state our four main theorems. In all of them belongs to .
Theorem 1.1**.**
If and , then belongs to . Moreover, the second mixed null derivatives of belong to , the restriction of to symmetry spheres is , and satisfies the wave equation on the Cauchy horizon.
Theorem 1.2**.**
Let . If and , then belongs to .
Theorem 1.3**.**
If the spherical mean of (see (35)) belongs to and , then does not belong to , for any open set .
Since the inequality is valid in the entire subextremal range of Reissner-Nordström we conclude that, if , then it cannot be extended in to any neighborhood of any point on the Cauchy horizon.
It is an easy consequence of [24] that if with , then belongs to . We prove that this result is essentially sharp.
Theorem 1.4**.**
If the spherical mean of belongs to with , then does not belong to , for any open set .
The organization of this paper is as follows. In Section 2 we explain the basic setup of our problem. In Section 3 we recall three energy estimates due to Sbierski. In Section 4 we upgrade the previous to pointwise estimates. In Section 5 we prove Theorem 1.1 which establishes the existence of a classical solution up to and including the Cauchy horizon. In Section 6 we prove Theorem 1.2 concerning solutions with higher regularity. Finally, in Section 7 we prove Theorems 1.3 and 1.4 about blow up in and in .
2. Setup
2.1. Some useful coordinate systems.
We will study solutions of the wave equation on a fixed background consisting of the black hole region of a subextremal Reissner-Nordström (asymptotically flat, de Sitter or anti-de Sitter) spacetime. This spacetime has a metric given in a local coordinate system by
[TABLE]
where is the round metric on the -sphere, and
[TABLE]
Here is the mass, is the charge parameter and is the cosmological constant. We will assume that the function has at least two positive roots, the smallest of which are
[TABLE]
The values and correspond to the values of at the Cauchy horizon and at the event horizon , respectively. The Penrose diagram of this spacetime for positive is given in Figure 1.
The surface gravities of the Cauchy and event horizons, defined by
[TABLE]
are of fundamental importance to us here. Throughout we will assume that the surface gravities do not vanish, which restricts the scope of our analysis to the subextremal setting.
For , we have
[TABLE]
Moreover, any tortoise coordinate
[TABLE]
satisfies, for ,
[TABLE]
The black hole region corresponds to
[TABLE]
a region where the function is negative, and where varies in .
We will often rely on the double null coordinates given in terms of and by
[TABLE]
In these coordinates the metric takes the form
[TABLE]
Clearly, we have
[TABLE]
Note that the event horizon corresponds to and the Cauchy horizon corresponds to . Since our double null coordinates are singular at these horizons, at the event horizon we change from coordinates to coordinates using
[TABLE]
In these coordinates the metric becomes
[TABLE]
At the Cauchy horizon we change from coordinates to coordinates using
[TABLE]
In these coordinates the metric is written as
[TABLE]
Note that to change from coordinates to coordinates we can use
[TABLE]
By abuse of notation, we will write .
It is important to note that the vector field is Killing. We also denote by , for , the generators of spherical symmetry, and just by any one of the three. The vector fields are also Killing.
2.2. The wave equation
Define
[TABLE]
[TABLE]
with denoting the spherical laplacian of ,
[TABLE]
The wave equation,
[TABLE]
is equivalent to both
[TABLE]
and
[TABLE]
2.3. The energy-momentum tensor
Recall that to a scalar function we may associate the energy-momentum tensor
[TABLE]
whose relevance for the study of solutions of the wave equation stems from the fact that its divergence satisfies
[TABLE]
Our energy estimates for will be obtained by applying the Divergence Theorem to certain currents, which are contractions of the energy-momentum tensor with appropriate vector fields. It will be useful to have the expression of the energy-momentum tensor in coordinates. One readily checks that
[TABLE]
Again, denotes the spherical gradient of ,
[TABLE]
and
[TABLE]
2.4. Energy identities and the Divergence Theorem
We will apply the Divergence Theorem in regions bounded by hypersurfaces , where is constant equal to , hypersurfaces , where is constant equal to , and hypersurfaces , where the geometric variable is constant equal to . Denoting by , and the corresponding normals, with unit and all three future directed, and denoting by , and the corresponding volume elements, we have
[TABLE]
where is the volume form associated to . Note that along the null hypersurfaces there is no natural choice of normal or volume form, so one can just choose a convenient normal and then let the Divergence Theorem determine the volume form.
Our currents will be vector fields of the form
[TABLE]
with timelike and future pointing, so that if is a solution of the wave equation, then
[TABLE]
Our choices of will be such that is nonnegative. We denote by
[TABLE]
Applying the Divergence Theorem to the current in the region
[TABLE]
(see Figure 2) we get the energy identity
[TABLE]
For a hypersurface , the integral
[TABLE]
controls first order derivatives of . Let us give an example by defining, near the Cauchy horizon, . This choice leads to
[TABLE]
[TABLE]
[TABLE]
Note that the expressions inside the square parentheses above are nonnegative as required by the fact that energy-momentum tensor satisfies the Dominant Energy Condition.
For we have
[TABLE]
and
[TABLE]
so that for sufficiently close to we have
[TABLE]
[TABLE]
[TABLE]
The notation means that there exist positive constants and such that . The blue-shift vector field, which will play a relevant role below, satisfies
[TABLE]
and so expressions (11), (12) and (13) are also valid if we replace by .
3. Basic energy estimates
We denote by
[TABLE]
We will now recall some basic energy estimates. The first one applies to the red-shift region. According to [24, p. 113, (4.5.5)] we have
Lemma 3.1**.**
For every , there exists a future directed timelike time invariant vector field in , a and a constant such that we have
[TABLE]
for all belonging to .
At the event horizon and satisfies for .
The second energy estimate applies to the no-shift region. According to the proof of [24, Lemma 4.5.6] we have
Lemma 3.2**.**
Given and a future directed timelike time invariant vector field in , there exists a constant such that
[TABLE]
for all satisfying .
Note that, in the previous case, might be negative but we can apply the Divergence Theorem with vector field and sufficiently large so that is nonnegative.
The third energy estimate applies to the blue-shift region. According to the proof of [24, Proposition 4.5.8] we have
Theorem 3.3**.**
Assume . For every sufficiently small , , there exists a future directed timelike time invariant vector field in , a and a constant such that, for all , we have (see [24, p. 114, last line])
[TABLE]
and (see [24, p. 117, (4.5.10)] together with the previous inequality)
[TABLE]
for all belonging to .
Moreover, the function belongs to .
At the Cauchy horizon and satisfies for .
4. Pointwise estimates for , and
Let . Throughout this section we will assume that . We will fix satisfying . The objective of the next three subsections is to prove that the three estimates
[TABLE]
hold for and . We choose a satisfying .
4.1. Estimates for
In a parallel manner to (13), there exists such that
[TABLE]
holds for .
To obtain a uniform bound on we use the following five ingredients:
(i) From (14) and the fact that (or from [24, p. 113, (4.5.5)]), for , we have that
[TABLE]
for a red-shift vector field that satisfies . Since and the vector fields are Killing, and implies that , for every and every multi-index (see Remark 6.1), we see that estimate (21) holds with replaced by .
(ii) We now apply Sobolev’s inequality in symmetry spheres and (20) to obtain
[TABLE]
and then we use (21) to conclude that
[TABLE]
(iii) We recall [3, Lemma 4.5].
Lemma 4.1**.**
Let and assume that for some , and for all ,
[TABLE]
Then, for all ,
[TABLE]
Below we will take and .
(iv) If we take squares of both sides of
[TABLE]
and then apply Hölder’s inequality we get
[TABLE]
We can then use (22) and Lemma 4.1 to conclude that
[TABLE]
(v) From [3, Lemma 4.2], we know that
Lemma 4.2**.**
Let , and . Let be continuous and such that
[TABLE]
for all . Then, for \alpha<\min\bigl{\{}\frac{b}{M},\Delta\bigr{\}}, we have
[TABLE]
for all .
If we consider the function and add a large multiple of (which satisfies (22)) to both sides of (4.1), we can apply the previous lemma to obtain the pointwise estimate
[TABLE]
for and .
Note that the constant , in the last estimate, is uniform in because is a bounded function of . For we also have (17) since and .
Since estimate (21) also holds with replaced by , for any multi-index , we can repeat the proceeding of the previous paragraphs to obtain (18) and (19), at , with constants , once again, uniform in .
4.2. Estimates for
Let be a timelike future directed and time independent vector field. Similarly to (12), we have
[TABLE]
Since is Killing, and implies that , and , for , by (15) we have
[TABLE]
This together with (17), applied with , implies that
[TABLE]
Thus, we obtain (17) for .
Applying (15) to , , and , and using (18) and (19) for , we obtain (18) and (19) for . The constants do not depend on or .
4.3. Estimates for
In this region, according to (16), we have
[TABLE]
Applying (24) to , and implies that extends continuously to along segments of constant . Moreover, is the uniform limit of as (for , where ). Arguing as in [6, Proposition 5.2, Step 2], is continuous in . The same reasoning can be used to show that , and extend continuously to . A simple argument implies that the derivative of the continuous extension of with respect to exists and coincides with the continuous extension of . Analogous statements apply to and . Using (24), and reasoning as we did in the region , we see that (17), (18) and (19) hold for .
5. Existence of a classical solution up to the Cauchy horizon
Henceforth, by “up to the Cauchy horizon” we mean up to and including the Cauchy horizon. In this section, we will use the energy estimates of the previous sections, together with Lemma 4.2, to obtain a pointwise bound for , for a fixed . This together with the previously established pointwise bounds for other derivatives of , which are valid up to the Cauchy horizon, can then be used to integrate (10) and obtain a pointwise bound , up to the Cauchy horizon. Finally, the control of this quantity in will allow us to extend as a classical solution of the wave equation, all the way up to the Cauchy horizon.
Proof of Theorem 1.1.
We proceed in four steps.
(i) Bounding for . Assume . We now fix satisfying and, as before, choose satisfying .
We will start by showing that, for a fixed , we have
[TABLE]
Indeed, this follows by the procedure developed in Section 4.1: we start by realizing that for any , with and a mutli-index, we have, in view of (16),
[TABLE]
which implies that
[TABLE]
and allows one to estimate
[TABLE]
We can now apply Lemma 4.2 to to finish the proof of (25).
(ii) Bounding up to the Cauchy horizon. Let
[TABLE]
so that . Integrating the wave equation, in its form (10), along a segment with fixed , from to , we get
[TABLE]
Choose and such that , for . Using (25), and (17) and (19) to estimate , yields
[TABLE]
for and .
(iii) Continuity of up to the Cauchy horizon. We define by
[TABLE]
Let . To prove the uniform convergence of to , as , for the first variable belonging to , we write
[TABLE]
Let . Again using estimates (17) and (19) to control , we can fix and sufficiently big so that , for . Now, using (25), fix sufficiently large so that , for and . Finally, invoking the uniform convergence, in , of to and of to , as , we are allowed to fix such that , for . For and , we have
[TABLE]
This proves the stated uniform convergence. Again, arguing as in [6, Proposition 5.2, Step 2], is continuous in . A simple argument implies that the derivative of the continuous extension of with respect to exists and coincides with the continuous extension of .
(iv) The wave equation is satisfied on the Cauchy horizon. To justify that the wave equation (10) is satisfied on the Cauchy horizon we just have to differentiate the right-hand side of (27) with respect to . Note that we are not claiming that is up to the Cauchy horizon but merely that exists, is continuous and satisfies (10). We can also guarantee that exists and is continuous. Indeed, define
[TABLE]
The wave equation can also be written as
[TABLE]
This can be integrated to
[TABLE]
Note that is continuous up to the Cauchy horizon because it is equal to . Therefore, equation (29) holds with . Another application of the Fundamental Theorem of Calculus guarantees that exists and is continuous up to the Cauchy horizon and that (28) is satisfied also on the Cauchy horizon. Of course, the fact that exists and is continuous up to the Cauchy horizon is enough to guarantee that (28) is satisfied on the Cauchy horizon. ∎
6. Solutions with higher regularity
This section is devoted to the proof of Theorem 1.2.
6.1. Wave equations for
We start by deriving the inhomogeneous wave equations satisfied by higher order -derivatives of . The commutators of and , and of and , are
[TABLE]
and
[TABLE]
where
[TABLE]
and
[TABLE]
Consequently, the function satisfies the inhomogeneous wave equation
[TABLE]
the function satisfies the inhomogeneous wave equation
[TABLE]
and, in general, satisfies the inhomogeneous wave equation
[TABLE]
with
[TABLE]
We also define
[TABLE]
Because () is a differential operator of order one in (), () involves a sum of derivatives of whose order with respect to () is at most .
Taking into account that the sequence defined by (7) satisfies
[TABLE]
we can write (31) as
[TABLE]
Similarly, we have that
[TABLE]
6.2. Derivatives of in
Recall the definition of in (2).
Remark 6.1**.**
Let . Then .
Proof.
For the result follows immediately from the fact that the vector fields are Killing and tangent to the sphere . For the remaining cases it suffices to note that:
- (i)
Obviously, . 2. (ii)
. Recalling that , the wave equation (9) shows that if vanishes on the event horizon and , then . Differentiating both sides of (9) once with respect to we then conclude that . If we keep differentiating both sides of (9) with respect to , we see that . 3. (iii)
. Before handling the general case, lets us also go over the case in detail. Suppose now that vanishes on the event horizon, and . Using (30) together with the previous paragraph we conclude that . One can also argue that , for all . 4. (iv)
. More generally, using (31), if vanishes on the event horizon and its first derivatives with respect to vanish at , then , for all .
∎
6.3. Bounding higher derivatives of
The wave equation (10) can be used to bound when . If we integrate the inhomogeneous wave equation for and use the bounds for derivatives of , whose order with respect to is at most equal to one, then we can bound up to the Cauchy horizon if . The comes from the fact that in (32), with , the function in front of is and . Moreover, when we will be able to generalize the previous procedure and establish boundedness of .
Proof of Theorem 1.2.
We proceed in three steps.
(i) . Since implies , using (26) applied to , when we have
[TABLE]
for . Then the wave equation (10) shows that on the hypersurface we have that
[TABLE]
Integrating the wave equation (32) with we obtain
[TABLE]
The derivatives of appearing inside the integral have order at most one with respect to , namely they are , , , and . Moreover,
[TABLE]
for . Therefore, under the assumption that we obtain that
[TABLE]
for . As, in addition, when , , and are continuous on we have that when .
(ii) . Let us consider another specific case, , before analyzing the general situation: repeating the previous argument, using (33) applied to , when we have
[TABLE]
for . The wave equation (30) shows that on the hypersurface we have that
[TABLE]
Since the above mentioned derivatives, and are controlled when , integrating the wave equation (32) with , when we obtain that
[TABLE]
for . All other third order derivatives of are continuous when and we are able to conclude that when .
(iii) . The general case now clearly follows by induction: in fact, when we obtain that
[TABLE]
for . In conclusion, provided . ∎
7. Blow up in and in
7.1. Blow up in
This subsection is devoted to the proof of Theorem 1.3. Let us start by sketching the main ideas for this proof. Suppose that is spherically symmetric. If for all large we have that is positive, then it turns out that and are positive for . This fact can be used to propagate a lower bound for , at , all the way up to the Cauchy horizon. We can then obtain a lower bound for and a negative upper bound for , which can be used to obtain the desired blow up result.
Proof of Theorem 1.3.
We proceed in six steps.
(i) Initial data for . Assume first that is spherically symmetric and that the restriction of to the ingoing null hypersurface , through the event horizon, vanishes to order and does not vanish to order , on the event horizon. Then there exist constants such that (eventually replacing by )
[TABLE]
for . As , using (4) we have
[TABLE]
for . According to (5) and (6), for , there exist constants such that
[TABLE]
Thus, we get
[TABLE]
for .
(ii) in . Since is spherically symmetric the wave equation reduces to
[TABLE]
According to [6, Lemma B.1], and are positive for . So is an increasing function. This implies that
[TABLE]
for .
(iii) for . Now we integrate the following (spherically symmetric) version of the wave equation
[TABLE]
between and . Taking into account that in this region
[TABLE]
we obtain
[TABLE]
On the hypersurface we have , and so
[TABLE]
As is positive, it follows that
[TABLE]
for .
(iv) for . In the region , according to (5) and (6), there exist constants such that
[TABLE]
Therefore,
[TABLE]
(v) Blow up. For the right-hand side of (34) goes to as goes to . In this case does not extend to a function up to the Cauchy horizon.
Suppose now that is not spherically symmetric. Its spherically mean
[TABLE]
is also a solution of the wave equation. According to our hypotheses, is not up to the Cauchy horizon. Therefore cannot be up to the Cauchy horizon.
(vi) Uniform blow up. The previous analysis only provides blow up of the norm of along null rays , with . To extend the result to all outgoing null rays intersecting , assume that is bounded along some . Then we would be able to solve the spherically symmetric wave equation sideways, with characteristic initial data provided by the spherical mean along , for . But since the spacetimes region is compact, local well posedness (see for instance [4, Theorem 4.5]) and the regularity of the initial data would imply boundedness of along . This is a contradiction. ∎
7.2. Blow up in
To prove that does not belong to it is enough to prove that its spherically symmetric part does not belong to . Moreover, the negative upper bound (34) applied to the spherical mean can be used to obtain a lower bound for the norm of .
Proof of Theorem 1.4.
We proceed in three steps.
(i) The norm of . To define a norm on , we define a Riemannian metric on in the usual way. After choosing a unit timelike vector field , we let
[TABLE]
Our choice of is
[TABLE]
This leads to
[TABLE]
The square of the norm of the gradient of is
[TABLE]
The volume element on is
[TABLE]
Now we may define the norm of to be
[TABLE]
(ii) Decomposition of the norm. Again, let be the spherical mean of . We remark that
[TABLE]
Indeed, this follows from
[TABLE]
For the second equality we used (35) and the fact that belongs to .
(iii) Blow up. From (36), to prove that does not belong to it is enough to prove that does not belong to . Using (34), the norm of is bounded below by
[TABLE]
where and . This integral is infinite provided . The extension of the blow up to any neighborhood of follows once again by local well posedness, as in the end of Subsection 7.1. ∎
Acknowledgements
We thank J. Natário and J.D. Silva for useful comments on a preliminary version of this paper. We also thank Anne Franzen for sharing and allowing us to use Figure 1. This work was partially supported by FCT/Portugal through UID/MAT/04459/2013 and grant (GPSEinstein) PTDC/MAT-ANA/1275/2014.
The reference list from the paper itself. Each links out to its DOI / PubMed record.
- 1[1] Angelopoulos, Y., Aretakis, S. and Gajic, D., Late-time asymptotics for the wave equation on spherically symmetric stationary spacetimes , Advances in Mathematics, Vol. 323 (2018) 529-621.
- 2[2] Cardoso, V., Costa, J.L., Kyriakos, D., Hintz, P. and Jansen, A., Quasinormal modes and Strong Cosmic Censorship , Phys. Rev. Lett. 120 (2018) 031103.
- 3[3] Costa, J.L. and Franzen, A., Bounded energy waves on the black hole interior of Reissner-Nordström-de Sitter , Ann. Henri Poincaré (2017) 18: 3371.
- 4[4] Costa, J.L., Girão, P.M., Natário, J. and Silva, J.D., On the global uniqueness for the Einstein-Maxwell-scalar field system with a cosmological constant. Part 1: Well posedness and breakdown criterion , Class. Quant. Gravity 32 (2015) 015017.
- 5[5] Costa, J.L., Girão, P.M., Natário, J. and Silva, J.D., On the global uniqueness for the Einstein-Maxwell-scalar field system with a cosmological constant. Part 2: Structure of the solutions and stability of the Cauchy horizon , Comm. in Math. Phys. 339 (2015), 903–947.
- 6[6] Costa, J.L., Girão, P.M., Natário, J. and Silva, J.D. On the global uniqueness for the Einstein-Maxwell-scalar field system with a cosmological constant. Part 3: Mass inflation and extendibility of the solutions , Ann. PDE (2017) 3: 8.
- 7[7] Costa, J.L., Girão, P.M., Natário, J. and Silva, J.D., On the occurrence of mass inflation for the Einstein-Maxwell-scalar field system with a cosmological constant and an exponential Price law , Comm. Math. Phys. (2018) 361:289.
- 8[8] Dafermos, M. and Luk, J., The interior of dynamical vacuum black holes I: The C 0 superscript 𝐶 0 C^{0} -stability of the Kerr Cauchy horizon , ar Xiv:1710.01722 v 1.
