Naked Singularity Censoring With Anisotropic Apparent Horizon
Naked Singularity Censoring With Anisotropic Apparent Horizon
Naked Singularity Censoring With Anisotropic Apparent Horizon
XINLIANG AN
Abstract. Employing the Einstein-scalar field system, we demonstrate an approach for prov-
ing high co-dimensional nonlinear instability of naked-singularity solutions as constructed by
Christodoulou in [18]. We further investigate the censorship of Christodoulou’s naked singu-
larity and show that a tiny anisotropic perturbation arising from the outgoing characteristic
initial data would lead to the emergence of an anisotropic apparent horizon, which covers and
censors the naked singularity. Our approach advances the hyperbolic short-pulse method by
not requiring the aid of additional large parameters, by permitting the use of initial pertur-
bations for the shear tensor and the derivative of scalar field to be with finite BV and C 0
arXiv:2401.02003v1 [gr-qc] 3 Jan 2024
norms, and by allowing the initial perturbation to be arbitrarily small in scale-critical norms.
New elliptic arguments based on non-perturbative methods are also developed.
1. Introduction
In 1994 Christodoulou [18] constructed his famous naked-singularity solution in 3 + 1 dimen-
sions for the Einstein-scalar field system
1
Rµν − Rgµν = Tµν ,
2 (1.1)
1 λ
Tµν =Dµ φDν φ − gµν D φDλ φ.
2
3+1
Here φ : R → R is a real-valued scalar function and (M, g) is the 3+1 dimensional spacetime.
The stability or instability of these naked singularities is closely related to the weak cosmic
censorship [20, 34]:
Conjecture 1. (Weak cosmic censorship conjecture). For generic asymptotically flat initial
data, the maximal development of Einstein’s field equations possesses a complete future null
infinity I + and hides the possibly formed singularities in a (black hole) region causally discon-
nected from I + .
The above is one of the greatest open problems in classic general relativity. When Penrose
formulated the original weak cosmic censorship conjecture in [34], the caveat generic was not
included. Christodoulou’s example in [18] surprisingly shows that naked singularity could form
in the evolution of Einstein’s field equations. Within spherical symmetry, for the Einstein-scalar
field system (1.1), by a celebrated series of works [16,17,19] Christodoulou proved that the initial
datum leading to the naked-singularity formation is rare, and for generic initial data the possibly
formed singularities are indeed hidden in the trapped (black hole) region. This motivates him
to add the term generic for Conjecture 1 in [20].
Christodoulou’s proofs in [16,17,19] reveal novel mathematical and physical properties of the
spherically symmetric gravitational collapse. For the initial datum leading to naked singularity,
in [19] he identified 2 instability mechanisms and deduced co-dimensional 2 instability. These
instability theorems show that each rare initial datum, which leads to naked singularity forma-
tion, is always associated with 2-dimensional perturbations and each of the perturbations would
give rise to a trapped surface and hence a black hole formation. Here a trapped surface is a
2-sphere in M with both null expansions negative.
In above arguments, the spherically symmetric perturbation required in [19] is a bit restric-
tive, since it imposes a global condition in all angles. Can we relax this requirement by only
placing a small anisotropic-in-angle perturbation? If it is possible, this may also significantly
raise the co-dimensions of the instability theorems. In this article, we give such an anisotropic
result. In particular, we prove that for Christodoulou’s constructed naked singularity in [18], by
imposing a tiny perturbation from any angle, we can trigger a trapped region and an apparent
horizon to censor the naked singularity. Compared to Christodoulou’s instability argument that
the singularity is covered by later-formed trapped surfaces, here we further show that the naked
singularity is censored by an achronal apparent horizon emerging from it. Subjected to a general
anisotropic perturbation, the naked singularity in [18] is not naked even locally. It is completely
censored for observers outside the apparent horizon.
1.1. Main Theorem. Here we use the double-null foliation for the spacetime (M, g). Let u
and u be optical functions increasing towards the future and satisfy the eikonal equations
g µν ∂µ u∂ν u = 0, g µν ∂µ u∂ν u = 0.
Set Hu to be the level sets of u and H u to be the level sets of u. And we require that Hu and H u
correspond to the outgoing and incoming null cones, respectively. The intersections of Hu and
H u are topologically 2-spheres and are denoted by Su,u . We further fix a frame e1 , e2 tangent
to Su,u and denote ∇ to be the angular derivative on it. In this paper, we employ the double
null coordinates (u, u, θ1 , θ2 ) for (M, g), where ∂/∂θ1 and ∂/∂θ2 are tangential to Su,u .1 For
A, B = 1, 2, the spacetime metric g takes the form
g = −2Ω2 (du ⊗ du + du ⊗ du) + g/AB (dθA − dA du) ⊗ (dθB − dB du).
Here Ω(u, u, θ1 , θ2 ) is the lapse function, g/AB is the induced metric on Su,u , and A, B = 1, 2.
For notational simplicity, in this article we also use ω ∈ S2 to denote (θ1 , θ2 ). With the help of
Ω, we choose the null pair e3 , e4 in this article as
∂ ∂ ∂
e3 = Ω−1 , e4 = Ω−1 + dA A
∂u ∂u ∂θ
with dA obeying dA = 0 on H−1 . And it holds
g(e3 , e3 ) = 0, g(e4 , e4 ) = 0, g(e3 , e4 ) = −2.
Define
χAB := g(DeA e4 , eB )
and set χ̂ to be its traceless part, i.e., the shear tensor. In this article, we first prove
Theorem 1.1. Consider the characteristic initial value problem for the Einstein-scalar field
system (1.1). Assigning Christodoulou’s naked-singularity initial data in [18] along H 0 with
−1 ≤ u ≤ 0 and prescribing perturbed initial data along H−1 satisfying
X
uj k(∂u )j ∇i χ̂kL∞ 2
u L (S−1,u )
+ uj k(∂u )j ∇i ∇4 φkL∞ 2
u L (S−1,u )
≤ B, (1.2)
i≤5,j≤3
then for each B there exist a sufficient small δ = δ(B) and the Einstein-scalar field system
admits a unique regular solution in the spacetime region (u, u, ω) satisfying 0 ≤ u ≤ δ and
u ≤ |u|Ω2−δ̃ (u, 0) ≤ 1 with 0 < δ̃ ≪ 1 a fixed constant.
Moreover, if we further assume
Z u
|χ̂|2 (−1, u′ , ω) + [∇4 φ(−1, u′ , ω) − ∇4 φ(−1, 0, ω)]2 du′ = f (ω, u)u, (1.3)
0
with f (ω, u) obeying
0 ≤ f ≤ 1 on S2 × (0, δ] and f (·, u) ≥ m on Bp (ǫ) (1.4)
2
for a point p ∈ S and constants m ∈ (0, 1) and ǫ ∈ (0, π/2), then within the solved hyperbolic
region, there exists a unique MOTS (marginally outer trapped surface) Mu on each H u with
0 < u ≤ δ and the collection of MOTS {Mu } emerges and censors the central singularity and
forms an achronal apparent horizon.
1See Chapter 1 and Chapter 2 in [21] for a detailed discussion of the double null coordinates.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 3
For the shear tensor and the derivative of scalar field, the above theorem prescribes initial
[0,δ]
perturbation along H−1 with a finite BV norm satisfying (1.3). Our next theorem further
[0,δ]
allows initial perturbation along H−1 to be with finite C 0 norm.
Theorem 1.2. The above statements in Theorem 1.1 still hold, if we replace the requirements
for f (ω, u) in (1.3) and (1.4) by
1 1
f (ω, u) = g(u)f˜(ω, u) with g(u) ≈ [ln(ln )]− 2 obeying lim+ g(u) = 0,
u u→0
and f˜(ω, u) satisfying 0 ≤ f˜ ≤ 1 on S2 × (0, δ], f˜(·, u) ≥ m on Bp (ǫ) for a point p ∈ S2 with
constants m ∈ (0, 1) and ǫ ∈ (0, π/2).
,ũ
,δ
µ∗
µ∗
1−
1−
a)
a)
(ũ
(δ
−
−
S
S
O
Hδ
(u
H ũ
=
δ)
(u
=
H0
ũ)
1)
(u
−
=
=
0)
1(
u
−
H
Figure 1
Remark 2. Compared with the previous short-pulse results in the scale-critical regime, here
we don’t need the additional large parameter a and we can view that the parameter a is set to
be a = a(B) = 1. This is an important advance for the short-pulse method invented in [21] by
1
Christodoulou with a = δ − 2 ≫ 1 and generalized in [2] and [4] with a = a(B) ≫ 1 by the author
and Luk, and by the author, respectively. With the aid of this parameter a, we could consider the
1
strength of the short pulse to be of order a 2 . The larger the short-pulse perturbation, the easier
to form a trapped surface. A novel insight of this paper is that, via employing the instability
4 XINLIANG AN
nature of the interior naked-singularity solution, to form an anisotropic trapped surface and the
corresponding apparent horizon, we can set the upper bound of the short-pulse perturbation to
1
be of order a(B) 2 = 1, and according to (1.4) we could even set the size of short pulse to be
a universal small number m. As pointed out in Remark 1, we can further choose its size to be
g(u) with g(u) → 0 as u → 0+ . The perturbation allowed can be arbitrarily small measured at
the scale-critical level. This is markedly different from [2] by the author and Luk, and [6] by the
author and Han, where the constant a = a(B) has to be a large parameter and the fact 1/a ≪ 1
is utilized crucially throughout [2, 6] to control the nonlinear terms. This paper contains other
hyperbolic and elliptic advancements. To capture the key features of the anisotropic instability
mechanism, we also develop quite a few new ingredients. These are summarized in Section 1.4.
Remark 3. In this paper, all the estimates still hold if we restore sharp a-weights and set
a = a(B) to be a parameter greater than or equal to 1. In particular, we allow a = a(B) ≫ 1,
but we do not need a ≫ 1 to close the arguments. To be consistent with [2,6] and to contain more
analytic information for further use, in below sections we work under the following assumptions
for initial perturbation:
X 1 1
uj a− 2 k(∂u )j ∇i χ̂kL∞ 2
u L (S−1,u )
+ uj a− 2 k(∂u )j ∇i ∇4 φkL∞ 2
u L (S−1,u )
≤ B, (1.5)
i≤5,j≤3
and
Z u
|χ̂|2 (−1, u′ , ω) + [∇4 φ(−1, u′ , ω) − ∇4 φ(−1, 0, ω)]2 du′ = f (ω, u)ua. (1.6)
0
O
Anisotropic Apparent Horizon
Anisotropic Perturbation
along u = −1
S−1,0
Figure 2
Theorem 1.3. Consider the space of prescribed initial data for χ̂(−1, u, ω)q
and ∂u φ(−1, u, ω)
√ 1−µ∗
along u = −1 with u ≥ 0. With |χ̂|(−1, 0, ω) = 0, ∂u φ(−1, 0, ω) = 2 · 1−µ
∗
µ∗ · 4 cor-
responding to their values of Christodoulou’s naked-singularity solution in [18], then for any
k ∈ Z+ , we can find functions f1 , f2 , ..., fk and g1 , g2 , ..., gk , which are continuous with respect
to the u variable and smooth with respect to the ω variable, such that the below prescribed
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 5
would lead to the trapped surface and the apparent horizon formation when
k
X k
X 2
λ2i + λ′i 6= 0 with λk , λ′k ∈ R.
1 1
Moreover, if it holds
0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ′1 g1 + λ′2 g2 +···+ λ′k gk
= 0 + λ̃1 f1 + λ̃2 f2 + · · · + λ̃k fk , ∂u φ(−1, 0, ω) + λ̃′1 g1 + λ̃′2 g2 +···+ λ̃′k gk ,
then it renders
λi = λ̃i and λ′i = λ̃i with 1 ≤ i ≤ k.
+
Hence, we may say that, for any k ∈ Z , Christodoulou’s naked-singularity solution has at least
co-dimension 2k nonlinear instability subject to outgoing continuous characteristic perturbations
at (−1, 0, ω).
Remark 5. In Section 13, we obtain two instability theorems. With respect to the shear tensor
and the derivative of the scalar field, Theorem 13.1 proves instability subject to perturbations
with finite BV norms and Theorem 1.3 proves instability subject to perturbations with finite
C 0 norms. In [2, 3, 6], to trigger the short pulse, BV perturbations for ∂u g/|S−1,0 are prescribed.
In [2, 3, 6], when u ≤ 0, the shear tensor χ̂ vanishes, and when u > 0, the shear tensor χ̂ obeys
a large lower bound related to the large parameter a. For this article, in Section 13.1 we first
prove Theorem 13.1, which is consistent with [2, 3, 6] by imposing the BV perturbations. In
Section 13.2, we further allow the perturbations to be C 0 . Additional gains in both hyperbolic
and elliptic parts enable us to do so. Detailed discussions are provided in Section 13.
Remark 6. Our above theorems prove arbitrarily large co-dimensional nonlinear instability for
a naked-singularity solution. This is not reported before and we spot the importance of the
anisotropic arguments. Together with Remark 1, we can see that this paper contains a desired
scale-critical nonlinear instability result.
Remark 7. The proofs of above two theorems and the approach we demonstrate here have
the potential to be generalized to other Einstein field equations. For the hyperbolic part, we
employ Ω(u, 0) → 0 as u → 0. For the elliptic part, we utilize Ωω(u, 0) = −∂u log Ω(u, 0)/2 >
0, which also holds for Rodnianski-Shlapentokh-Rothman’s naked-singularity in [35] for the
Einstein vacuum equations.
Our proofs of above two theorems are rooted in recent progress on trapped surface formation
and apparent horizon emergence. In below, we review these developments.
1.2. Formation of Trapped Surface. With a 589-page proof, Christodoulou [21] proved the
first trapped-surface-formation result for Einstein’s equations with no symmetry. As depicted in
Figure 3 below, Christodoulou considered a characteristic initial-value problem for the Einstein
vacuum equations in the shadowed region. Minkowskian initial data are prescribed along the
incoming null cone H 0 and short-pulse initial data are prescribed along the outgoing null cone
H−1 . Along H−1 , denoting χ to be its associated outgoing null second fundamental form and
setting χ̂ to be its traceless part, in 2008 Christodoulou proved the following2:
2In [21] Christodoulou’s short-pulse initial data are prescribed at past null infinity. Here we focus on its
correspondence in a finite region.
6 XINLIANG AN
Hδ
∗
u
H
(u
=
δ)
H0
1)
(u
−
=
=
0)
1(
u
−
H
Figure 3
Theorem (Christodoulou [21]). Set H 0 to coincide with a backward light cone in the Minkowskian
spacetime for −1 ≤ u ≤ 0. For every B > 0 and −1 < u∗ < 0, there exists a δ = δ(B, u∗ ) > 0
sufficiently small such that the following holds: for 0 ≤ u ≤ δ, if the prescribed χ̂ on H−1
satisfies X 1
δ 2 +j k(∂u )j ∇i χ̂kL∞ 2
u L (S−1,u )
≤B (1.7)
i≤5, j≤3
with ∇ being the angular derivative on S−1,u , then the solution to Einstein vacuum equations
remains regular in the spacetime region −1 ≤ u ≤ u∗ , 0 ≤ u ≤ δ. Furthermore, if the below
lower-bound assumption also satisfies
Z δ
inf2 |χ̂|2 (−1, u′ , ω) du′ ≥ M∗ > 2|u∗ |, (1.8)
ω∈S 0
then the Einstein vacuum equations admit a unique regular solution in the spacetime region
1
−1 ≤ u ≤ −bδa 2 and 0 ≤ u ≤ δ. Moreover, if we further require
Z δ
inf2 |χ̂|2 (−1, u′ , ω)du′ ≥ 4δa,
ω∈S 0
then S−δa,δ is the formed 2-dimensional trapped surface.
Later, via employing only the signature for decay rates and a spacetime rescaling, a reproof
1
and an extension of Theorem 1.2 to the past null infinity was obtained in [4] for the case b = a 2
by the author.
Remark 8. The scale-critical bounds in [2] inspire the hyperbolic estimates of this article.
However, there are at least three major differences:
(1) In [2] the formed trapped surface is S−δa,δ and it is in the approximately self-similar
regime. The designed short-pulse hierarchy and the hyperbolic estimates are all con-
sistent with the considerations of scaling criticality. This is the underlying reason that
in [4] the author can give an alternative proof in the far field regime and very few bor-
derline terms appear. In this article, the trapped surface is located around S−(δa)1−µ∗ ,δ
with 0 < µ∗ < 1 and it is not in the approximately self-similar regime. Hence, designing
a suitable hierarchy and closing the bootstrap arguments are more delicate. We need
to divide the spacetime into two regions and deal with many new borderline terms.
Moreover, finding a proper order to conduct the estimates for recovering the bootstrap
constants is crucial. The newly designed hierarchy is given in Section 2.4 and the order
of estimates is explained in Section 2.6.
(2) To carry out the hyperbolic estimates in [2] and [4], in the spacetime region considered,
we apply the below inequality ubiquitously
1
ua 2 1
≤ 1 ≪ 1.
|u| a2
The above inequality hinges on the largeness of the additional parameter a. While in
this article our conclusions do not depend on the largeness of a and we can set a = 1,
which advances the short-pulse method.
(3) In this paper, we can even set the short-pulse to be small in scale-critical norms. To
fulfill the requirement (13.15) in Theorem 13.3, we can choose
1 1
|χ̂(−1, u, ω)| + |∇4 φ(−1, u, ω) − ∇4 φ(−1, 0, ω)| ≈ [ln(ln )]− 4 .
u
1
Note that the function [ln(ln 1/u)]− 4 → 0 as u → 0+ . Hence, we can let |χ̂(−1, u, ω)| and
∇4 φ(−1, u, ω) − ∇4 φ(−1, 0, ω) to be absolutely continuous functions. Their BV norms
in u variable (and other scale-critical norms) can be set small. While in all previous
trapped surface formation works outside of spherical symmetry, the BV perturbation of
|χ̂(−1, u, ω)| in u variable requires to be large.
1.3. Emergence of Apparent Horizon. Besides trapped surface formation, the rise of the
apparent horizon in the gravitational collapse of Einstein’s field equations was proved in [3] by
the author and in [6] by the author and Han. In [3], the following result was proved
Theorem (An [3]). Under the same initial-data assumption as in Theorem 1.2, the Einstein
1
vacuum equations admit a unique regular solution in the region 0 < u ≤ δ and −1 ≤ u ≤ −bua 2 .
In addition, assuming
Z u
|χ̂|2 (−1, u′ , ω)du′ = f (ω, u)ua for each 0 < u ≤ δ, (1.10)
0
with f (ω, u) being a smooth function in S2 × (0, δ] satisfying 20/21 ≤ f (ω, u) ≤ 22/21. Then
there exists a unique MOTS Mu along each H u (0 < u ≤ δ). Furthermore, requiring
1
a− 2 kuj ∇je4 ∇i χ̂kL∞ 2
u L (S−1,u )
≤ B for and i, j ∈ Z+ ,
then the MOTSs {Mu } together form a smooth apparent horizon.
8 XINLIANG AN
,δ
ũ
δa
,
ũa
−
−
S
S
O
Hδ
(u
H ũ
=
δ)
(u
=
H0
ũ)
1)
(u
−
=
=
0)
1(
u
−
H
Figure 4
The next work by An-Han generalized the above result to the fully anisotropic scenario.
Theorem (An-Han [6]). With the same prescribed initial data as in Theorem 1.2, for each
1
0 < u ≤ δ, if we require a− 2 kuj ∇je4 ∇i χ̂kL∞ 2
u L (S−1,u )
≤ B for any i, j ∈ Z+ and
Z u
|χ̂|2 (−1, u′ , ω)du′ = f (ω, u)ua, (1.11)
0
with f (ω, u) satisfying
0≤f ≤1 on S2 × (0, δ] and f (·, u) ≥ m on Bp (ǫ)
2
for a point p ∈ S and constants m ∈ (0, 1) and ǫ ∈ (0, π/2). Then in a large region of each
H u with 0 < u ≤ δ, there exists a unique MOTS Mu and the collection of MOTS {Mu } emerge
from the central singularity and form a smooth apparent horizon. Furthermore, if we require
1
|u∂u f (ω, u)| ≤ a− 3 , then the apparent horizon is a spacelike hypersurface, hence a dynamical
horizon.
In [6], this theorem was further extended to the multi-valley scenario and concrete examples
of dynamically formed apparent horizon with single and multi valleys are provided. In Figure
5, one can see the shape of MOTS along H u with two valleys, which are bounded by the upper
and lower barriers.
Hu
B
C
upper barrier
A
MOTS Mu (ω)
lower barrier
p2 p1
S2
Figure 5
- deriving apriori estimates via working through details of Moser’s iteration and Schauder’s
estimates in the order of
1 1
C 0 → C 0,α → C 1, 3 → C 2, 3 with 0 < α ≪ 1;
- conducting the method of continuity based on the explicit structure of the linearized
operator.
Compared with [6], in this article the geometry along the initial incoming cone is completely
different and the main terms in the elliptic equation for R(u, ω) (with u = R(u, ω) being the
MOTS along H u ) differentiate, which raises extra challenges. In [6] we employ the ansatz
e e In this
R = uae−φ to translate the elliptic equation for R(u, ω) into an elliptic equation for φ.
article, we design and utilize two forms of ansatz:
∗ ∗ e
(1) R = (ua)1−µ e−(1−µ )φ . We apply this form of R(u, ω) to perform the apriori estimates
via Moser’s iteration and Schauder’s estimates.
1−µ∗
1−µ∗
(2) R = (ua) φ̌ µ∗
. We employ this ansatz to eliminate the additional nonlinear bor-
derline term (rising from naked-singularity initial data) in the elliptic equation of R(u, ω)
and to conduct the method of continuity for the corresponding new linearized equation.
1.4. Difficulties and New Ingredients in the Proof. To prove Theorem 1.1, we divide the
considered spacetime region into two parts:
1
- the exterior {(u, u, ω)| − 1 ≤ u ≤ u1 and 0 ≤ ua 2 ≤ |u1 |Ω2−δ̃ (u1 , 0)},
1
- the interior {(u, u, ω)| u1 ≤ u ≤ 0 and 0 ≤ ua 2 ≤ |u|Ω2−δ̃ (u, 0)}.
Here |u1 | is a small constant to be fixed later. For notational simplicity, in the rest of this article,
we set constant δ̃ to obey 0 < δ̃ ≪ 1. In reality, if needed we can choose δ̃ to be any number
satisfying 0 < δ̃ < 1/2 and the below proof still works.
1.4.1. Hyperbolic Part. Compared with [2] by An-Luk, here more challenges arise. We list some
of them and state several strategies:
- For the position of MOTS along each H u , in contrast to the results in [2, 4, 6], where
∗
the MOTS lies around u = −ua, in this paper the MOTS is around u = −(ua)1−µ
with 0 < µ∗ < 1. The hyperbolic estimates in [2,4,6] are consistent with asymptotically
self-similar bounds and crucially rely on
1
ua 2 1
. 1 ≪1
|u| a2
with a being a large parameter. While here the MOTS is not located in the asymp-
totically self-similar regime, and to obtain a desired instability result we would allow
general (tiny) anisotropic perturbations and would set the parameter a to be of size 1.
These bring many new difficulties, especially on how to design and prove the short-pulse
hierarchy when |u| is small.
- In this paper, we generalize our constructed ansatz in [2,4] by inserting extra |u| weights
to the hierarchy. Employing the fact
1 µ∗
Ω(u, 0) ≤ Ω(u, u, θ1 , θ2 ) ≤ 2Ω(u, 0) and Ω2 (u, 0) = |u| 1−µ∗ ,
2
the suitable |u| weights are added via imposing the corresponding Ω weights. Espe-
cially, the Ω-weight top-order energy estimates and top-order elliptic estimates are
critical for our proof. In the energy estimates, because of the additional Ω weights,
basic integration-by-part and pairing arguments are revised. In comparison to [2],
the top-order elliptic estimatesPare more subtle. As shown in Section 4 and Sec-
tion 7, enhanced estimates for 1≤i≤4 kui ∇i (Ωtrχ)kL2 (Su,u ) , ku6 ∇5 (Ωtrχ)kL2u L2 (Su,u ) ,
precise estimates involving borderline terms for ku5 ∇5 ηkL2u L2 (Su,u ) , ku6 ∇5 ηkL2u L2 (Su,u ) ,
10 XINLIANG AN
ku5 ∇5 (trχ, χ̂)kL2u L2 (Su,u ) are derived and are employed in an essential way. Subtle struc-
tures of the Einstein-scalar field system are explored.
- When bounding the nonlinear terms in the interior region where u1 ≤ u ≤ 0, we utilize
the inequality
1
ua 2 1 µ∗ 1 µ∗ 1
3 ≤ Ω 2 −δ̃ (u, 0) = |u| 1−µ∗ · 2 −δ̃ ≤ |u1 | 1−µ∗ · 2 −δ̃ ≪ 1, (1.12)
|u|Ω 2
where 0 < δ̃ ≪ 1 and |u1 | ≪ 1. Notice that the largeness of a is no longer required, but
the power of Ω has to be discreetly tracked.
- To conduct energy estimates for the scalar field, to derive energy estimates for curvature
components and to obtain elliptic estimates, we encounter terms, which only obey the
upper bound
1 1 1 1
u2 a4 ua 2 2
1 = .
|u| 2 Ω |u|Ω2
By contract to (1.12), this quotient is no longer . 1. The terms satisfying only this
bound serve as the new borderline terms to be estimated. We then add the weights
1 1
1 + u 1a
2 4
to the bootstrap assumptions. Finding the correct order to retrieve the
|u| 2 Ω
constants (which do not depend on other bootstrap assumptions) in these borderline es-
1 1
|u| 2 Ω
strap assumptions for Õ5,2 (η), Õ5,2 (Ωtrχ, trχ), Õ5,2 (Ωχ̂, χ̂), S(∇A φ), S(∇3 φ), R(K, σ̌)
and R(β) listed in Section 2.4. To improve the bootstrap assumption for each of these
terms, the correctly designed order of arguments and enhanced estimates are vitally
employed.
- Compared with An-Luk [2], there are extra borderline terms arising from Ωω and ∂u φ.
Along u = 0, they satisfy
r
µ∗ 1 √ µ∗ 1
Ωω(u, 0) = · and ∂u φ(u, 0) = 2 .
4(1 − µ∗ ) |u| 1 − µ∗ |u|
They behave similarly to Ωtrχ and cause multiple new borderline terms.
- To carry out energy estimates, besides σ̌ defined in [2], we also use renormalized quan-
tities βA − 21 ∇4 φ∇A φ (equivalent to βA − 12 R4A ) and β A + 21 ∇3 φ∇A φ (equivalent to
β A + 12 R3A ) to combine contributions of ∇3 βA , ∇3 R4A or ∇4 β A , ∇4 R3A in null Bianchi
equations (2.13) together.3 Furthermore, we also introduce and employ the renormal-
ized Gaussian curvature K − |u|1 2 − 41 ∇A φ∇A φ. This new renormalization comes from
the below identity for the Weyl curvature
1 1 1
∇A WA434 = − ∇4 R43 + ∇3 R44 − ∇A R(gA3 g44 − gA4 g43 ).
2 2 12
It gives
1 1 1 1
∇4 (ρ − R43 − R) = ∇A WA434 − ∇3 R44 + non-top-order terms
4 12 2 4
with ρ = W3434 /4. Further calculations via applying the constraint equation imply that
1 1 1 1 1
ρ − R43 − R = −(K − ∇A φ∇A φ) + χ̂ · χ̂ − trχtrχ.
4 12 4 2 4
For using vanishing initial data along H 0 , we then utilize the renormalized Gaussian
curvature K − |u|1 2 − 41 ∇A φ∇A φ to carry out the energy estimates. In our estimates,
3A similar idea was used in [5] by An-Athanasiou for the Einstein-Maxwell system.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 11
1.4.2. Elliptic Part. Because of the naked-singularity initial data, here the key structure of the
main elliptic equation for MOTS is quite different from its counterpart in [6]. In this article, we
design and utilize two forms of ansatz:
∗ ∗ e e
- R = (ua)1−µ e−(1−µ )φ . This is a generalization of [6], where R = uae−φ . We use this
form to construct the anisotropic trapped surfaces and to conduct Moser’s iteration and
Schauder’s estimates as in [6].
∗ 1−µ∗
- R = (ua)1−µ φ̌ µ∗ . The ansatz here is new and is indispensable for this paper. We
employ this form to carry out the key continuity method for the new linearized equation
to solve for the MOTS.
In below, we will give a more comprehensive explanation on the utilization of these two forms.
To construct the trapped and untrapped barriers and to derive the regularity estimates, we
∗ ∗ e
employ the ansatz R = (ua)1−µ e−(1−µ )φ . We observe that
- Owning to Christodoulou’s initial data, the constructions of the multi- and single-valley
anisotropic trapped surfaces as in [6] and [28] still hold. One can see the details in
Section 9.
- For the apriori regularity estimates, we find that the approach developed in [6] by An-
Han is robust enough and can be adapted to this new setting. Consider the main elliptic
equation (10.1):
1 e e2
0 =∆γ φe + 1 − f (u, ω)eφ − 4(1 − µ∗ )RΩω|∇γ φ|
2
R e 2 + 2γ ij ηi ∂j φe
− (1 + Ωtrg χ)(1 − µ∗ )|∇γ φ|
2
1 1 e
+[ ∗
RΩ−1 trg χ − 1 + f (u, ω)eφ ]
2(1 − µ ) 2
∗ e e
− 2(1 − µ )RΩχ̂ · ∇γ φ · ∇γ φ.
e
On the right side, besides the main terms ∆γ φe + 1 − 21 f (u, ω)eφ we denote
∗ 1−µ∗
The new ansatz R = (ua)1−µ φ̌ µ∗ plays a vital role in linearization arguments.
12 XINLIANG AN
∗ ∗ e
- If we proceed with the ansatz R = (ua)1−µ e−(1−µ )φ . The main elliptic equation for
MOTS is then reduced to
1 e
0 = ∆γ φe + 1 − f (ω, u)eφ − 4(1 − µ∗ )RΩω|∇γ φ|
e 2 + small terms.
2
In [6], since RΩω is small, the main terms in above equation are the first three terms
on the right and the continuity argument via linearization can be carried out accord-
e 2 term). While in this paper, due to the naked-
ingly (without troubles from the |∇γ φ|
singularity initial data, we have
µ∗ 1
−4(1 − µ∗ )RΩω(u, u, ω) ≈ −4(1 − µ∗ )R · ∗
· = −4µ∗
4(1 − µ ) |u|
being of size 1. It is no longer negligible. Here, by introducing and employing the new
∗ 1−µ∗
ansatz R = (ua)1−µ φ̌ µ∗ , we eliminate this quadratic gradient borderline term and
the equation of MOTS becomes
µ∗ 1
0 = ∆γ φ̌ − µ∗ φ̌ + f (ω, u)φ̌1− µ∗ + small terms.
2
This enables us to complete the method of continuity and to construct the unique MOTS
along each H u in the considered spacetime region.
1.4.3. Instability Results. Our proven anisotropic trapped surface and apparent horizon forma-
tion criteria enable us to obtain instability theorems of high co-dimensions.
• Allowing perturbations with finite BV norms along u = −1 with u > 0, in Theorem
13.2 we consider the space of initial data for χ̂(−1, u, ω) and ∂u φ(−1, u, ω), that can
be prescribed along u = −1 with u > 0. Given Christodoulou’s naked-singularity data
in [18] along u = 0, then for any k ∈ Z+ , we can find smooth functions f1 , f2 , ..., fk and
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 13
g1 , g2 , ..., gk with variables u > 0 and ω, such that the below prescribed initial data for
|χ̂| and ∂u φ with u > 0
|χ̂|(−1, u, ω), ∂u φ(−1, u, ω)
′ ′ ′
= 0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ1 g1 + λ2 g2 + · · · + λk gk
would lead to the trapped surface and the apparent horizon formation, when
k
X k
X 2
λ2i + λ′i 6= 0 with λk , λ′k ∈ R.
1 1
Moreover, if it holds
0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ′1 g1 + λ′2 g2 +···+ λ′k gk
= 0 + λ̃1 f1 + λ̃2 f2 + · · · + λ̃k fk , ∂u φ(−1, 0, ω) + λ̃′1 g1 + λ̃′2 g2 +···+ λ̃′k gk ,
then it renders
λi = λ̃i and λ′i = λ̃i with 1 ≤ i ≤ k.
We may therefore say that, for any k ∈ Z+ , Christodoulou’s naked-singularity solution
in [18] has at least co-dimensional 2k nonlinear instability subject to outgoing BV char-
acteristic perturbations at (−1, 0, ω).
• In Theorem 13.2, with mi and m′i being positive constants, the constructed fi and gi
obey the requirement: for u > 0
1 1
fi (ω, u) ≥ mi a 2 on Bpi (ǫi ) ⊂ Di , gi (ω, u) ≥ m′i a 2 on Bp′i (ǫ′i ) ⊂ Di′ ,
and
fi (ω, u) = 0, gi (ω, u) = 0 for u ≤ 0.
Hence, functions fi (ω, u) and gi (ω, u) are not continuous in u, when u → 0+ . In Section
13.2, by exploiting the extra Ω-weight gained in the hyperbolic estimates, we extend
this instability theorem by allowing to choose continuous functions fi (ω, u) and gi (ω, u)
for the u variable with u ≥ 0. In Theorem 1.3, we impose the below requirements for fi
and gi :
1 1 1
With 1 ≤ i ≤ k, C1 > 2, Λ > 0 and g(u) = (C1 Λ + C1 ) 2 · [ln(ln )]− 2 ,
u
we request the prescribed fi and gi to obey
1 1 1 1
0 ≤ fi ≤ g(u) 2 a 2 on S2 × (0, δ], fi (·, u) ≥ mi g(u) 2 a 2 on Bpi (ǫi ) ⊂ Di ,
1 1 1 1
0 ≤ gi ≤ g(u) 2 a 2 on S2 × (0, δ], gi (·, u) ≥ m′i g(u) 2 a 2 on Bp′i (ǫ′i ) ⊂ Di′
with points pi ∈ Di ⊂ S2 , p′i ∈ Di′ ⊂ S2 and constants mi , m′i ∈ (0, 1), ǫi , ǫ′i ∈ (0, π/2).
Note that we can set a = 1. And as u → 0+ , it holds
1
lim g(u) 2 = 0.
u→0+
1.5.1. Trapped Surfaces, MOTS and Short-Pulse Method. After Christodoulou’s monumental
work [21], a systematical approach by introducing signature for the short pulse was developed
by Klainerman-Rodnianski in [29] to simplify and to extend [21] in a finite region. Yu in [38]
extended Klainerman-Rodnianski’s approach to the Einstein-Maxwell system. The author in [1]
introduced a new signature, called the signature for decay rates, and extended Klainerman-
Rodnianski’s results to past null infinity. Based on the hyperbolic estimates in [29], Klainerman-
Luk-Rodnianski in [28] found a fully anisotropic mechanism for trapped surface formation.
In [2] the author and Luk designed a new set of short-pulse hierarchies and proved the first
scale-critical trapped surface formation criterion for the Einstein vacuum equations. With only
the signature for decay rates, in [4] the author gave a short proof of the scale-critical trapped
surface formation in the far-field regime and related the short-pulse hierarchies to the peeling
properties of the gravitational waves. In [5], the author and Athanasiou further extended the
new approach in [4] to the Einstein-Maxwell system. With spherically symmetric singular initial
data prescribed along the incoming null hypersurface, Li and Liu in [32] studied the Einstein-
scalar field system and proved an almost scale-critical trapped surface formation criterion. Their
result showed that the spherically symmetric singularities considered in [32] are unstable (as a
result of trapped surface formation) subject to close-to-isotropic gravitational perturbations
and the instability in [32] is in the sense of the first category, slightly weaker than co-dimension
1. Very recently, the methods in [4] by the author and in [5] by the author and Athanasiou
are also extended to the Einstein-Yang-Mills system and the Einstein-scalar field system, by
Athanasiou-Mondal-Yau [15] and Zhao-Hilditch-Kroon [39], respectively.
Based on the hyperbolic estimates in [2], in the short-pulse regime, the MOTS along each
incoming null hypersurface and the corresponding apparent horizon are constructed in [3] by the
author and in [6] by the author and Han for the general anisotropic scenarios. For solving the
MOTS and the apparent horizon based on spacelike foliations, interested readers are referred
to [11–14] by Andersson, Eichmair, Mars, Metzger, Simon and references therein.
1.5.2. Weak Cosmic Censorship within Spherical Symmetry and Naked Singularity Formation.
In a series of works [16]- [19], Christodoulou achieved the proof of weak cosmic censorship for the
spherically symmetric Einstein-(real) scalar field system. After constructing the first example of
the naked-singularity solution under spherical symmetry, he ingenuously obtained 4-integrated
theorems:
(1) BV extension principle in triangle and in rectangle regions;
(2) a sharp trapped surface formation criterion;
(3) the first instability theorem;
(4) the second instability theorem.
Christodoulou’s arguments in [16]- [19] rely on the special structures of the spherically symmet-
ric Einstein-(real) scalar field system and the associated monotonicity properties. His 4-step
argument is unified. Each step is exactly matched with the other steps. In recent years, each of
Christodoulou’s steps has been separately extended. In [8] the author and Lim extended [16] to
the spherically symmetric Einstein-Maxwell-charged (complex) scalar field system. See also the
extension in [23] by Costa to the spherically symmetric Einstein-(real) scalar field system with
positive cosmological constant. For the instability argument of the Einstein-(real) scalar field
system, Li-Liu in [31] gave a proof with apriori estimates. Their instability theorem is slightly
different from Christodoulou’s.
Extending the proof of the weak cosmic censorship to the other spherically symmetric Einstein-
matter field is still very challenging. The Einstein-Maxwell-charged (complex) scalar field system
is the next to be considered. However, the corresponding weak cosmic censorship has remained
open since the late 1990s. Several difficulties are easily seen:
• Due to the presence of non-constant charge Q, almost all previously employed mono-
tonicity formulas fail to hold. These monotonicity formulas are critically used in Christodoulou’s
BV extension/trapped surface formation/instability arguments when Q = 0. A new
strategy to incorporate non-constant charge Q remains to be developed.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 15
• For the trapped surface formation criterion, with η being the initial mass input and
δ being the small deformation parameter, Christodoulou needs η & δ log δ1 . Following
Christodoulou’s method, for the charged case, in [8] the author and Lim obtained a
ω
result requiring η & δ 1− 2 (> δ log δ1 ) with 0 < ω ≪ 1. But Christodoulou’s “blue-shift”
ω
γ in the instability argument is around log δ1 rather than δ − 2 . How to reconcile?
• Christodoulou’s way of proving his instability theorems is very delicate. Can one obtain
the same sharp arguments as Christodoulou did?
In a new preprint [9], the author and Tan address these questions and we show that the weak
cosmic censorship still holds (in the sense of Christodoulou) for the gravitational collapse of the
spherically symmetric Einstein-Maxwell-charged (complex) scalar field system.
Recently, there is also the exciting development of constructing the naked-singularity solutions
for Einstein’s field equations. Under spherical symmetry and self-similar ansatz, in [10] the
author and Zhang showed naked-singularity formation for Einstein vacuum equations in higher
dimensions, and in [25] Guo-Hadzic-Jang extended [18] to the Einstein-Euler system. Outside
spherical symmetry, important progress has been made by Rodnianski-Shlapentokh-Rothman in
[35] and by Shlapentokh-Rothman in [36]. There they prove the first naked-singularity formation
result for the Einstein vacuum equations in 3 + 1 dimensions.
1.6. Organizations of the Paper. The paper is organized as follows. In Section 2 we present
the settings. In Section 3, we derive the preliminary estimates. In Section 4, we provide the
L2 (Su,u ) estimates for Ricci coefficients. In Section 5, we derive the L2 (Su,u ) estimates for the
scalar field. In Section 6, we prove energy estimates for the scalar field. In Section 7, we establish
the elliptic estimates. In Section 8, we conduct the energy estimates for curvature components.
In Section 9, we construct the sub and super barriers for the MOTS. In Section 10, we carry
out Moser’s iteration and the Schauder estimates. In Section 11, we demonstrate the crucial
linearization and continuity arguments. In Section 12, we obtain the existence and properties
for the apparent horizon. In Section 13, we prove the instability theorems.
1.7. Acknowledgement. XA would like to thank Yu Deng and Willie Wong for valuable cor-
respondences. XA is supported by MOE Tier 1 grants A-0004287-00-00, A-0008492-00-00 and
MOE Tier 2 grant A-8000977-00-00.
2. Setting
2.1. Double Null Coordinates. In this paper, we adopt the double null coordinates (u, u, θ1 , θ2 )
for (M, g). For A, B = 1, 2, the spacetime metric g takes the form
2.2. Equations. The equations studied here are the Einstein-scalar field system:
1
Rµν − Rgµν = Tµν ,
2
1
Tµν = Dµ φDν φ − gµν Dλ φDλ φ.
2
16 XINLIANG AN
Utilizing the null frame e3 , e4 and an frame e1 , e2 tangent to the 2-spheres Su,u , with the indices
A, B ∈ {1, 2}, we define the Ricci coefficients:
1 1
ω̂ = − ∇4 (log Ω), ω = − ∇3 (log Ω),
2 2 (2.3)
ηA = ζA + ∇A (log Ω), η A = −ζA + ∇A (log Ω).
1 1
Wiklm := Riklm + (Rim gkl − Ril gkm + Rkl gim − Rkm gil ) + R(gil gkm − gim gkl )
2 6
and define the following null curvature components
where /ǫ is the volume form associated to the metric g/. For totally symmetric tensors, we also
employ the below div and curl operators
Let χ̂ and χ̂ be the traceless parts of χ and χ respectively. Employing above notations,
rewriting Rµν − 12 Rgµν = Tµν according to the null frame {e3 , e4 , eA , eB } with A, B = 1, 2, we
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 17
And the Ricci coefficients ηA , η A , ω, ω̂ obey the remaining null structure equations:
1
∇4 ηA = −χAB · (η − η)B − βA − DA φD4 φ,
2
1
∇3 η A = −χAB · (η − η)B + β A − DA φD3 φ,
2 (2.6)
1 1 1 1
∇4 ω = 2ω̂ω − ηA · η A + |η|2 + ρ + D3 φD4 φ + DA φDA φ,
2 2 6 12
1 1 1 1
∇3 ω̂ = 2ω̂ω − ηA · η A + |η|2 + ρ + D3 φD4 φ + DA φDA φ.
2 2 6 12
On Su,u , with the null frame, the Gauss-Codazzi equations boil down to the constraint equations:
1 1 1 1
div χ̂A = ∇A trχ − (ηB − η B ) · (χ̂BA − trχgBA ) − βA + DA φD4 φ,
2 2 2 2
1 1 1 1
div χ̂A = ∇A trχ + (ηB − η B ) · (χ̂BA − trχgBA ) + β A + DA φD3 φ,
2 2 2 2
1
curl η = − curl η = σ + χ̂ ∧ χ̂, (2.7)
2
1 1 1 1
K = − (ρ + R) − trχtrχ + χ̂ · χ̂ + g AB DA φDB φ
6 4 2 2
1 1 1 1
= − ρ − trχtrχ + χ̂ · χ̂ + D3 φD4 φ + DA φDA φ.
4 2 6 3
Here K is the Gauss curvature of the spheres Su,u .
For the system (1.1), the scalar field φ obeys ✷g φ = 0. And we rewrite the wave equation as
a transport equation
0 = ✷g φ =g λµ Dλ Dµ φ = g λµ Dλ (Dµ φ) − g λµ DDλ eµ φ
=g 34 D3 (D4 φ) − g 34 DD3 e4 φ + g 43 D4 (D3 φ) − g 43 DD4 e3 φ + g AB DA (DB φ) − g AB DDA eB φ
1 1
= − D3 D4 φ − D4 D3 φ + g AB DA (DB φ) − g AB D∇A eB + 12 χAB e3 + 21 χ e4 φ
2 2 AB
1 1
= − D3 D4 φ + ∆g φ − trχe3 φ − trχe4 φ.
2 2
It hence holds
1 1
D3 D4 φ + trχD3 φ + trχD4 φ = g AB ∇A ∇B φ. (2.8)
2 2
18 XINLIANG AN
For convenience, we also employ the following two equivalent forms (see Section 6)
1 1
(Ωe3 )(Ωe4 φ) + ΩtrχΩe4 φ = Ω2 ∆g φ − ΩtrχΩe3 φ+2Ω2 η A eA φ, (2.9)
2 2
1 1
(Ωe4 )(Ωe3 φ) + ΩtrχΩe4 φ = Ω2 ∆g φ − ΩtrχΩe3 φ+2Ω2 η A eA φ. (2.10)
2 2
Using the property of covariant derivatives, the scalar field φ also satisfies
1 1 1
∇i Wiklm = − ∇m Rlk + ∇l Rmk − (∇l Rgmk − ∇m Rglk ). (2.12)
2 2 12
Employing the null frame {e3 , e4 , eA , eB } with A, B = 1, 2, we rewrite the above (2.12) as the
below null Bianchi equations:
1 1 1
∇3 α + trχα = ∇⊗β b + 4ωα − 3(χ̂ρ +∗ χ̂σ) + (ζ + 4η)⊗β
b + (D3 R44 − D4 R43 )gAB + ∇4 RgAB ,
2 4 6
1
∇4 β + 2trχβ = div α − 2ω̂β + η · α − (DA R44 − D4 R4A ),
2
1 1 1
∇3 β + trχβ = ∇ρ + 2ωβ +∗ ∇σ + 2χ̂ · β + 3(ηρ +∗ ησ) + ∇A (g CD RCD ) − ∇B RBA + ∇A R + ∇3 R4A ,
2 6 2
3 1 1
∇4 σ + trχσ = −div ∗ β + χ̂ ·∗ α − ζ ·∗ β − 2η ·∗ β − (Dµ R4ν − Dν R4µ )ǫµν 34 ,
2 2 4
3 1 1
∇3 σ + trχσ = −div ∗ β + χ̂ ·∗ α − ζ ·∗ β − 2η ·∗ β + (Dµ R3ν − Dν R3µ )ǫµν 34 ,
2 2 4
3 1 1 1
∇4 ρ + trχρ = div β − χ̂ · α + ζ · β + 2η · β − (D3 R44 − D4 R43 ) + ∇4 R,
2 2 4 12
3 1 1 1
∇3 ρ + trχρ = −div β − χ̂ · α + ζ · β − 2η · β + (D3 R34 − D4 R33 ) + ∇3 R,
2 2 4 12
1 1 1
∇4 β + trχβ = −∇ρ +∗ ∇σ + 2ω̂β + 2χ̂ · β − 3(ηρ −∗ ησ)− ∇A (g CD RCD ) + ∇B RBA − ∇A R − ∇4 R3A ,
2 6 2
1
∇3 β + 2trχ β = −div α − 2ωβ + η · α + (DA R33 − D3 R3A ),
2
1 1 1
∇4 α + trχα = −∇⊗β b + 4ω̂α − 3(χ̂ρ −∗ χ̂σ) + (ζ − 4η)⊗β
b + (D4 R33 − D3 R34 )gAB + ∇3 RgAB .
2 4 6
(2.13)
4For more general incoming initial data, the renormalized quantity for K is K − K(u, 0) − 1 ∇A φ∇ φ.
4 A
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 19
1
∇3 (β· − ∇4 φ∇· φ)A + trχβA = − ∇K +∗ ∇σ̌ + 2ωβ + 2χ̂ · β − 3(ηK − ∗ ησ̌)
2
1 3
+ (∇(χ̂ · χ̂) −∗ ∇(χ̂ ∧ χ̂)) − ηtrχtrχ
2 4
3 ∗ 1 (2.14)
+ (η χ̂ · χ̂ − η χ̂ ∧ χ̂) − (∇trχtrχ + trχ∇trχ)
2 4
+ ∇A (g CD RCD ) − ∇B RBA
1
+ ηA ∇3 φ∇4 φ + ηA ∇C φ∇C φ.
2
3 1 1 1
∇4 σ̌ + trχσ̌ = − div ∗ (β − ∇4 φ∇φ) − ζ ·∗ β − 2η ·∗ β + χ̂ ·∗ (∇⊗η)
b + χ̂ ·∗ (η ⊗η)
b
2 2 2 2 (2.15)
1 1 C AB
− ∇1 R42 + ∇2 R41 + (∇A φ∇B φ − gAB ∇C φ∇ φ) ∧ χ̂ ,
2 2
1 1
∇4 (K − ∇A φ∇A φ) + trχ(K − ∇A φ∇A φ)
4 4
1 1
b + χ̂ · (η ⊗η)
= − div β − ζ · β − 2η · β + χ̂ · ∇⊗η b
2 2 (2.16)
1 1 1 1
− trχdiv η − trχ|η|2 + η A ∇A φ∇4 φ + ∆g φ∇4 φ
2 2 2 2
1 1
+ (∇A φ∇B φ − gAB ∇C φ∇C φ) · χ̂AB ,
2 2
3 1 1 1
∇3 σ̌ + trχσ̌ = − div ∗ (β + ∇3 φ∇φ) − ζ ·∗ β − 2η ·∗ β + χ̂ ·∗ (∇⊗η)
b + χ̂ ·∗ (η ⊗η)
b
2 2 2 2 (2.17)
1 1
+ ∇1 R32 − ∇2 R31 + χ̂AB ∧ (∇A φ∇B φ − gAB ∇C φ∇C φ),
2 2
1 1 3 1 1
∇3 (K − − ∇A φ∇A φ) + trχ(K − 2 − ∇A φ∇A φ)
|u|2 4 2 |u| 4
1 1 A
=div (β + ∇3 φ∇φ) − ∇ ∇3 φ∇A φ − ζ · β + 2η · β
2 2
1 1 1 1
+ χ̂ · ∇⊗η b − trχ|η|2 + η A ∇A φ∇3 φ
b + χ̂ · (η ⊗η) (2.18)
2 2 2 2
1 1 C AB 1
+ (∇A φ∇B φ − gAB ∇C φ∇ φ) · χ̂ − trχ∇A φ∇A φ
2 2 8
1 1 1 A Ω−1 2
+ trχ(−div η + K − 2 − ∇ φ∇A φ) − (Ωtrχ + ),
2 |u| 4 |u|2 |u|
1 1 1
∇4 (β · + ∇3 φ∇· φ)A + trχβ A =∇(K − 2 − ∇A φ∇A φ) +∗ ∇σ̌ + 2ωβ + 2χ̂ · β + 3(−ηK +∗ ησ̌)
2 |u| 4
1 A 1 1 1
+ ∇(− ∇ φ∇A φ + trχtrχ − χ̂ · χ̂) − ∗ ∇(χ̂ ∧ χ̂)
4 4 2 2
+ ∇A (g CD RCD ) + ∇B RBA
1
− η A ∇3 φ∇4 φ − η A ∇C φ∇C φ.
2
(2.19)
Remark 10. In below, if no further clarification, we employ the capital Latin letters A ∈ {1, 2}
as frame indices on the spheres Su,u , and use Greek letters µ ∈ {1, 2, 3, 4} as frame indices in
the whole 3+1 dimensional spacetime.
20 XINLIANG AN
u
O =
Naked Singularity δ
u
=
0
r=0
u
Anisotropic Perturbation
=
0
r=0
1
−
=
u
1
−
=
u
Christodoulou’s solution solves the Einstein-scalar field system 1.1. In [18] a spacetime sin-
gularity forms in gravitational collapse and the whole spacetime is free of trapped surfaces and
the apparent horizon.
In this article, we consider a characteristic initial value problem with data prescribed along
u = 0 with −1 ≤ u ≤ 0 and u = −1 with 0 ≤ u ≤ δ. For u ≤ 0, we adopt the spacetime
(Mint , gint ) constructed in [18] with Mint standing for the interior spacetime region u ≤ 0.
1
Setting v = 4u, we let gint = − Ω2 (u, v)(du ⊗ dv + dv ⊗ du) + r2 (u, v)[(dθ1 )2 + sin2 θ1 (dθ2 )2 ].
2
Here r(u, v) is the radius function.
We will add general (anisotropic) perturbation along u = −1 starting from u = 0 and for
u ≥ 0 the spacetime metric g is of the form (2.1). In particular, along u = 0, we prescribe
Christodoulou’s spherically-symmetric naked-singularity initial data: via Christodoulou’s con-
struction in [18], for Hawking mass m = 2r (1 + 4Ω−2 ∂u r∂v r), it holds that
2m
(u, 0) = µ∗ with µ∗ being a positive constant less than 1.
r
We also set
1 − µ∗
r(−1, 0) = 1, ∂v r(−1, 0) = and ∂u r(u, 0) = −1.
4
2m
Since along u = 0 it holds µ∗ = r = 1 + 4Ω−2 ∂u r∂v r, we get
1 − µ∗
Ω−2 ∂v r(u, 0) = along u = 0.
4
Employing the equation5
2m
∂u r r
∂u [ln(∂v r)] = · ,
r 1 − 2m
r
with r(−1, 0) = 1 and ∂v r(−1, 0) = (1 − µ∗ )/4, we obtain
1 − µ∗ 1−µ
µ∗ µ∗
∂v r(u, 0) = r ∗ (u, 0) and Ω2 (u, 0) = r 1−µ∗ (u, 0).
4
Using v = 4u, these further imply
µ∗
∂u r(u, 0) = (1 − µ∗ ) · r 1−µ∗ (u, 0) and Ω−2 ∂u r(u, 0) = 1 − µ∗ .
5It is equivalent to r∂ ∂ r = −∂ r∂ r − 1 Ω2 .
u v u v 4
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 21
Recall trχ(u, 0) = g AB g(DA e4 , eB )(u, 0), trχ(u, 0) = g AB g(DA e3 , eB )(u, 0) and along u = 0 it
holds e3 = Ω−1 ∂u , e4 = Ω−1 ∂u . We hence derive
2 −2 2(1 − µ∗ )
Ω−1 trχ(u, 0) = Ω−1 g AB · ΓD
A4 · gBD = Ω ∂u r(u, 0) = ,
r r(u, 0)
2 2
Ωtrχ(u, 0) = Ωg AB · ΓD
A3 · gBD = ∂u r(u, 0) = − .
r r(u, 0)
Utilizing ω = − 12 ∇3 log Ω, we also have
1 1 µ∗ µ∗ 1
Ωω(u, 0) = − ∂u log Ω(u, 0) = − ∂u log[r 2(1−µ∗ ) (u, 0)] = · . (2.20)
2 2 4(1 − µ∗ ) r(u, 0)
We proceed to calculate the initial value of K along u = 0. Applying the below null structure
equation and the Gauss equation
1 1 1
∇3 trχ + trχtrχ =2ωtrχ + 2ρ − χ̂ · χ̂ + 2div η + 2|η|2 − D3 φD4 φ + DA φDA φ,
2 3 3
1 1 1 1 A
K = − ρ − trχtrχ + χ̂ · χ̂ + D3 φD4 φ + DA φD φ,
4 2 6 3
1 −1
further employing ω = − 2 Ω ∇3 Ω and restricting the equations to u = 0, we obtain
2Ω2 K(u, 0) + ∂u (Ω2 · Ω−1 trχ)(u, 0) + Ω2 · Ω−1 trχ · Ωtrχ(u, 0) = 0.
Together with the fact
2(1 − µ∗ ) 2 µ∗
Ω−1 trχ(u, 0) = , Ωtrχ(u, 0) = − , Ω2 (u, 0) = r 1−µ∗ (u, 0), ∂u r(u, 0) = −1,
r(u, 0) r(u, 0)
we deduce
µ∗ µ∗ 4(1 − µ∗ )
2K(u, 0) + r− 1−µ∗ · ∂u (r 1−µ∗ −1 ) · 2(1 − µ∗ )(u, 0) − = 0,
r(u, 0)2
which infers
1 1
K(u, 0) = −(2 − 4µ∗ − 4 + 4µ∗ ) · = .
2r(u, 0)2 r(u, 0)2
Quoting the results in [18], for the scalar field along u = 0 it holds
r
√ µ∗ 1
Ω∂3 φ(u, 0) =∂u φ(u, 0) = 2 ,
1 − µ∗ r(u, 0)
r (2.21)
√ 1 − µ∗ ∂u r(u, 0)
Ω∂4 φ(u, 0) =∂u φ(u, 0) = 2 · .
µ∗ r(u, 0)
Remark 11. Along u = 0, employing ∂u r(u, 0) = −1 and ∂u (Ω−2 ∂u r) = − 21 rΩ−2 (∂u φ)2 , via
(2.21) we also obtain
1 |r∂u φ|2 µ∗ 1
Ωω(u, 0) = − ∂u log Ω(u, 0) = (u, 0) = ∗
· ,
2 8r 4(1 − µ ) r(u, 0)
which is consistent with (2.20).
√
Remark 12. The 2 on the right of (2.21) is owing to in [18] Christodoulou’s use of
1
Ricµν − Rgµν = 2Tµν ,
2
1
Tµν = ∂µ φc ∂ν φc − gµν ∂ σ φc ∂σ φc .
2
q
µ∗
Here we denote φc to be the scalar field function employed in [18]. Denoting k := 1−µ∗ , on
page 644 of [18], we have that along s = s∗ (i.e. u = 0) it holds
∂u φc ∂u φc 1
r (u, 0) = −k and r (u, 0) = . (2.22)
∂u r ∂u r k
22 XINLIANG AN
Remark 13. At the end of this subsection, we point out that, along u = 0 by employing
equations r∂u ∂v φ = −∂u r∂v φ − ∂v r∂u φ and r∂u ∂v r = −∂u r∂v r − 14 Ω2 , it holds
∂ ∂v φ Ω2 ∂v φ ∂u φ
− · + = 0. (2.23)
∂u ∂v r 4r∂v r ∂v r r
Treat ∂v φ(u, 0) as it is determined by other geometric quantities. Applying the fact
r
Ω2 −1 1 ∂u φ √ µ∗ 1
− (u, 0) = ∗
· and (u, 0) = 2 · ,
4r∂v r 1 − µ r(u, 0) r 1 − µ r(u, 0)2
∗
we can solve an ODE and we see that the below equality as stated in (2.21)
r
∂v φ(u, 0) ∂u φ(u, 0) √ 1 − µ∗ 1 1
= = 2 · ≈
∂v r(u, 0) ∂u r(u, 0) µ∗ r(u, 0) r(u, 0)
is debt to a special initial choice of ∂v φ/∂v r at S−1,0 . If we do not use this choice (by adding a
general perturbation), back to (2.23), we would have that ∂v φ(u, 0)/∂v r(u, 0) behaviors like
−1 −µ∗ µ∗
r(u, 0) 1−µ∗ = r(u, 0)−1 · r(u, 0) 1−µ∗ . With ∂v r(u, 0) ≈ r 1−µ∗ (u, 0), it yields ∂v φ(u, 0) ≈
r(u, 0)−1 . In below sections, we demonstrate an approach, which can deal with these general
scenarios.
2.4. Norms. For geometric quantities defined on Su,u , we now set norms that we will work
with. We use ψ and ψ to denote the below geometric quantities
ψ ∈ {η, η, Ωω − Ωω(u, 0), Ωχ̂, Ωtrχ − Ωtrχ(u, 0), ∇e1 φ, ∇e2 φ, Ω∂3 φ − Ω∂3 φ(u, 0)}. (2.25)
For conventional simplicity, we also write
Ωω(u, 0) =(Ωω)0 , Ω∂3 φ(u, 0) = (Ω∂3 φ)0 ,
−1 (2.26)
Ωtrχ(u, 0) = (Ωtrχ)0 , Ω trχ(u, 0) = (Ω−1 trχ)0 , Ω(u, 0) = Ω0 .
1 1 i+3 i
Ri,∞ (u, u) := 1 kui+2 ∇i (Ωβ)kL∞ (Su,u ) +ku ∇ ρkL∞ (Su,u )
a 2 ua
1 i+3 i 1 i+3 i 1
+ 1 ku ∇ σ̌kL∞ (Su,u ) + 1 ku ∇ (K − 2 )kL∞ (Su,u ) (2.28)
ua 2 ua 2 |u|
1 i+3 i
+ 1 ku ∇ (Ωβ)kL∞ (Su,u ) with 0 ≤ i ≤ 1.
ua 2
1 1
Oi,2 (u, u) := 1 kui ∇i ψ kL2 (Su,u ) + 1 kui+1 ∇i ψ kL2 (Su,u ) with 0 ≤ i ≤ 4. (2.29)
a 2 ua 2
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 23
We also set
1 5 5
Õ5,2 (u, u) := 1 1 ku ∇ (Ωtrχ, Ωω̂, Ωχ̂)kL2 2
u L (Su,u )
u a2
2
1 1 −1
u2 a4 1 5 5 1
+ 1 1 1 ku ∇ ηkL2 2
u L (Su,u )
+ 3 3 ku6 ∇5 ηkL2u L2 (Su,u )
|u| 2 Ω u2 a4 u2 a4
(2.30)
Ω 1
+ 1 1 ku5 ∇5 ηkL2u L2 (Su,u ) + 1 1 kΩu5 ∇5 ηkL2u L2 (Su,u )
u2 a2 u2 a2
1 1 −1
u2 a4 1 5 5
+ 1+ 1 1 1 ku ∇ (Ωtrχ, Ωχ̂, trχ, χ̂)kL2 2
u L (Su,u )
,
|u| 2 Ω u2 a4
′ 1
Õ5,2 (u, u) := 3 ku6 ∇5 (Ωtrχ)kL2u L2 (Su,u ) . (2.31)
u2 a
For 1 ≤ i ≤ 4, we further let
1 1
Ri,1 (u, u) := 1 1 kui+1 ∇i (ΩβA − Ω∇4 φ∇A φ)kL2u L2 (Su,u )
u a
2 2 2
(2.32)
1 1 1
+ 1 1 kΩui+1 ∇i (K − 2 − ∇A φ∇A φ, σ̌)kL2u L2 (Su,u ) ,
u2 a2 |u| 4
1 1 −1
u2 a4 1 i+2 i 1 1
Ri,2 (u, u) := 1 + 1 3 3 ku ∇ (K − 2 − ∇A φ∇A φ, σ̌)kL2u L2 (Su,u )
|u| 2 Ω u2 a4 |u| 4
1 1 −1
(2.33)
u2 a4 1 i+2 i 1
+ 1+ 1 3 3 ku ∇ (β A + ∇3 φ∇A φ)kL2u L2 (Su,u ) .
|u| 2 Ω u2 a4 2
For 1 ≤ i ≤ 5, we then set
1 1
Si,1 (u, u) := 1 1 kui ∇i (Ωe4 φ)kL2u L2 (Su,u ) + 1 1 kΩui ∇i (eA φ)kL2u L2 (Su,u ) , (2.34)
u a
2 2 u a2
2
1 1 −1
u2 a4 1 i i
Si,2 (u, u) := 1 + 1 1 1 ku ∇ (eA φ)kL2 2
u L (Su,u )
|u| 2 Ω u2 a4
1 1 −1
(2.35)
u2 a4 1 −1 i i
+ 1+ 1 1 1 kΩ u ∇ (Ωe3 φ)kL2u L2 (Su,u ) .
|u| 2 Ω u2 a4
And we further define O(u, u), R(u, u), S(u, u) as
X X X
O(u, u) := Oi,∞ (u, u) + Oi,2 (u, u), R(u, u) := Ri,1 (u, u) + Ri,2 (u, u),
i≤2 i≤4 1≤i≤4
X
S(u, u) := Si,1 (u, u) + Si,2 (u, u).
1≤i≤4
′
Next we denote Oi,∞ , Oi,2 , Õ5,2 , Õ5,2 , Ri,∞ , Ri,1 , Ri,2 , Si,1 , Si,2 to be the supremum over u, u
′
in the considered spacetime region of Oi,∞ (u, u), Oi,2 (u, u), Õ5,2 (u, u), Õ5,2 (u, u), Ri,∞ (u, u), Ri,1 (u, u), Ri,2 (u, u), Si
respectively. Finally, we set O, R, S to be:
X X X X
O := Oi,∞ + Oi,2 , R := Ri,1 + Ri,2 , S := Si,1 + Si,2 .
i≤2 i≤4 1≤i≤4 1≤i≤4
with
X X 1
I (0) (u) := 1 + 1 k(∇)m (∇4 )k (χ̂, ∇4 φ)kL2 (S−1,u ) .
0≤m≤6, 0≤k≤5
a2
24 XINLIANG AN
Remark 14. Following the arguments carried out in detail in Chapter of [21], we have
X X X
Oi,∞ (−1, u) + Oi,2 (−1, u) + Ri,∞ (−1, u) . I (0) .
i≤2 i≤4 i≤2
In Section 7, we will employ the initial values of Ωκ, µ, Ω−1 trχ along H−1 , with
1 1 1 1
κ := ∇ω̂ + ∗ ∇ω † − (β − ∇4 φ∇φ), µ := −div η + K − 2
− ∇A φ∇A φ.
2 2 |u| 4
Note that we set ω † to be 0 and Ω to be 1 along H−1 . And it holds
1 1 1
ku4 ∇4 (Ωκ)kL2 (S−1,u ) . a 2 , ku5 ∇4 µkL2 (S−1,u ) . ua 2 , ku5 ∇5 (Ω−1 trχ)kL2 (S−1,u ) . ua 2 .
2.5. Exterior and Interior Regions. To establish the hyperbolic existence result, we first
choose a u1 < 0 satisfying 0 < |u1 | ≪ 1 and then divide the spacetime M into two regions:
1
· the exterior: −1 ≤ u ≤ u1 and 0 ≤ ua 2 ≤ |u1 |Ω2−δ̃ (u1 , 0),
1
· the interior: u1 ≤ u ≤ 0 and 0 ≤ ua 2 ≤ |u|Ω2−δ̃ (u, 0).
Remark 15. For notational simplicity, in the rest of this paper, we set constant δ̃ obeying
0 < δ̃ ≪ 1. In reality, if needed we can set δ̃ to be any number satisfying 0 < δ̃ < 1/2 and all
the below proof stays the same.
In this article, we will use the below inequality frequently
1 1 1
ua 2 ua 2 ua 2
≤ ≤ 3 ≪ 1. (2.36)
|u| |u|Ω |u|Ω 2
This inequality is correct because (with a = 1)
1 3
ua 2 |u1 |Ω2−δ̃ (u1 ,0) |u1 |Ω 2 (u1 ,0) 1
· in the exterior we have 3 ≤ 3 = 3 · Ω 2 −δ̃ (u1 , 0) ≪ 1,
|u|Ω 2 |u|Ω 2 |u|Ω 2
1
ua 2 1 1
· in the interior we have 3 ≤ Ω 2 −δ̃ (u, 0) . Ω 2 −δ̃ (u1 , 0) ≪ 1.
|u|Ω 2
With above discussion, we can further choose a constant B, such that in both regions
1 1 1
ua 2 ua 2 ua 2 1
≤ ≤ 3 ≤ ≪ 1. (2.37)
|u| |u|Ω |u|Ω 2 B
2.6. An Approach of Bootstrap. In this article, we carry out a bootstrap argument to prove
uniform upper bounds for O, R and S. For initial data along u = 0 and u = −1, we have
O(0) + R(0) + S (0) ≤ C0 I (0) ,
where C0 depends only on initial data. For notational simplicity, we rewrite the above line as
O(0) + R(0) + S (0) . I (0) .
For the rest of this article, we keep the same notation: if P ≤ C0 Q with constant C0 depending
only on the initial data, we could write it as
P . Q.
(0)
Under this notation, we have I . 1.
Our goal in below sections is to show that in the considered spacetime region, it holds
O + R + S ≤ I (0) + (I (0) )2 + 1.
To achieve this, we make bootstrap assumptions
O(u, u) ≤ O, (2.38)
R(u, u) ≤ R, (2.39)
S(u, u) ≤ S, (2.40)
′
Õ5,2 (u, u) ≤ Õ, Õ5,2 (u, u) ≤ Õ′ . (2.41)
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 25
Recall (2.37). For constants in bootstrap assumptions (2.38)-(2.41), for the rest of this paper,
we assume that they satisfy
1 1
1 1 ua 2 1 ua 2 − 32
1
1 ≤ O, R, S, Õ ≪ (Õ′ ) 8 ≪ B 32 ≤ ( 3 )− 32 ≤ ( ) . (2.42)
|u|Ω 2 |u|Ω
In Section 5 we obtain
P P P
i≤4 Oi,2 (∇4 φ) . 1, i≤4 Oi,2 (∇A φ) . S, i≤4 Oi,2 (∇3 φ) . S.
In Section 6 we derive
S5,1 . 1, S5,2 . 1.
With the bounds obtained in Section 5 and Section 6, we revisit the estimates in Section 4.
Hence we now get
P P
i≤4 Oi,2 (ψ ) . 1, i≤4 Oi,2 (ψ ) . 1 + Õ5,2 (ω̂) + R(β).
In Section 7 we show
Õ5,2 (ω̂) . 1 + R(β), Õ5,2 (trχ, χ̂) . 1 + R(β),
′
Õ5,2 (trχ) . 1,
Õ5,2 (η) . 1, Õ5,2 (η) . 1 + R(K, σ̌) + R(β) + R(K, σ̌),
Õ5,2 (Ωtrχ, trχ) . 1, Õ5,2 (Ωχ̂) . 1, Õ5,2 (χ̂) . R(β).
In Section 8 we derive
R(β) + R(K, σ̌) . 1, R(K, σ̌) + R(β) . 1.
Remark 16. The proof strategy and the main theorems of this paper also hold for perturbed
Christodoulou’s initial data prescribed along u = 0. Denote
∂u r ∂u r
Ω̊(u, 0) = Ω(u, 0) and Γ̊(u, 0) ∈ { (u, 0), (u, 0), ∂u φ(u, 0), ∂u φ(u, 0)}.
r r
We further set Ω̊c (u, 0) and Γ̊c (u, 0) to be the corresponding values of Christodoulou’s naked-
singularity solution along u = 0. Our below proof strategy also extends to spherically-symmetric
initial data along u = 0 satisfying
with constants 0 < δ̃ ≪ 1 and a ≥ 1.7 For simplicity of the proof and for more precise
statements of the main theorems, in the main body of this paper we stick to Christodoulou’s
naked-singularity initial data with Footnotes 4, 9, 13, 14, 15, 16, 17 and Remark 20 added to
explain how to deal with the perturbed scenario.8
6A more detailed explanation of the bootstrap method can be found in Section 3.1 of [4].
7The reader is also referred to a spherically-symmetric extension of Christodoulou’s naked-singularity solution
by Singh in [37].
8The methods developed in this paper also allow non-spherically-symmetric initial data prescribed along
u = 0, the author will provide the details in a separate paper.
26 XINLIANG AN
3. Preliminary Estimates
3.1. Estimates for metric components. We now start the hyperbolic estimates. We first
obtain a bound for Ω under the bootstrap assumptions:
Proposition 3.1. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), we have
Ω 1
k − 1kL∞ (Su,u ) ≤ B − 2 .
Ω(u, 0)
Proof. We employ the equation
Ω 1 ∂ ∂
Ωω̂ = − ∇4 log Ω = − ( + dA A ) log Ω. (3.1)
2 2 ∂u ∂θ
Integrating equation (3.1) along e4 direction, we obtain
Z u 1
Ω ua 2 O
|| log ||L∞ (Su,u ) . ||Ωω̂||L∞ (Su,u′ ) du′ . .
Ω(u, 0) 0 |u|
Here we use the bootstrap assumption (2.38). This implies that
1
Ω ua 2 O 1
k − 1kL∞ (Su,u ) . ≤ B− 2 .
Ω(u, 0) |u|
We then control g/ under the bootstrap assumptions. Following the same argument as for
Proposition 5.2 in [2], we have
Proposition 3.2. Fix a point (u, ω) on the initial hypersurface H 0 . Along the outgoing charac-
teristic emanating from (u, ω), we define Λ(u) and λ(u) to be the larger and smaller eigenvalue
of g/−1 (u, u = 0, ω)g/(u, u, ω). Under the assumptions of Theorem 1.1 and the bootstrap assump-
tions (2.38), (2.39), (2.40) and (2.41), then it holds
1
|Λ(u) − 1| + |λ(u) − 1| ≤ B − 2
for every u ∈ [0, δ]. Consequently, for every u ∈ [0, δ], this also implies
1 dvolg/u
|ξ(u) − 1| ≤ B − 2 with ξ(u) = .
dvolg/0
As corollaries, we also obtain the below area estimate and the Sobolev embedding.
Proposition 3.3. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds that
1
sup |Area(Su,u ) − Area(Su,0 )| ≤ |u|2 B − 2 .
u,u
Proposition 3.4. With the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40) and (2.41), it holds that
X
kφkL∞ (Su,u ) . kui−1 ∇i φkL2 (Su,u ) . (3.2)
i≤2
Note that the detailed proof of the Sobolev embedding is provided in Propositions 5.6-5.9 in [2].
3.2. Estimates for Transport Equations. In this article, we will employ the below two
propositions frequently:
Proposition 3.5. With the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40) ,(2.41), for an Su,u -tangent tensor field φ of arbitrary rank it holds
Z u
||φ||L2 (Su,u ) . ||φ||L2 (Su,u′ ) + ||Ω∇4 φ||L2 (Su,u′′ ) du′′ .
u′
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 27
Proposition 3.6. Set φ and F to be Su,u -tangent tensor fields of rank k obeying the below
transport equation:
Ω∇3 φA1 ...Ak + λ0 ΩtrχφA1 ...Ak = ΩFA1 ...Ak
and denote λ1 = 2(λ0 − 21 ). Under the assumptions of Theorem 1.1 and the bootstrap assump-
tions (2.38), (2.39), (2.40), (2.41), then it holds that
Z u
|u|λ1 kφkL2 (Su,u ) . kφkL2 (S−1,u ) + |u′ |λ1 kΩF kL2 (Su′ ,u ) du′ .
−1
The proof of the above two propositions are the same as in Proposition 5.4 and Proposition 5.5
of [2] and hence are omitted.9
3.3. Commutation Formula. We quote the following formula from [27]:
Proposition 3.7. The commutator [∇4 , ∇] acting on an (0, r) S-tensor satisfies
[∇4 , ∇B ]φA1 ...Ar
=[D4 , DB ]φA1 ...Ar + (∇B log Ω)∇4 φA1 ...Ar − (g/−1 )CD χBD ∇C φA1 ...Ar
r
X r
X
− (g/−1 )CD χBD η A φA1 ...Âi C...Ar + (g/−1 )CD χAi B η D φA1 ...Âi C...Ar .
i
i=1 i=1
Remark 17. For the rest of this paper, when there is no danger of confusion, the constants on
the left of the equations are kept precise. The uniform constants in front of coefficients on the
right are usually omitted.
3.4. Elliptic Estimates for Hodge Systems. We first recall the definition of the divergence
and curl of a symmetric covariant tensor with an arbitrary rank:
(div φ)A1 ...Ar := ∇B φBA1 ...Ar ,
(curl φ)A1 ...Ar := /ǫ BC ∇B φCA1 ...Ar .
Here /ǫ is the volume form associated to the metric g/ on Su,u . Recall that the trace is defined
via
(tr φ)A1 ...Ar−1 = (g/−1 )BC φBCA1 ...Ar−1 .
In this article, we will use the following L2 (Su,u ) elliptic estimates proved in Chapter 7 of [21].
Proposition 3.11. Under the assumptions of Theorem 1.1 and the bootstrap assumptions
(2.38), (2.39), (2.40), (2.41), denoting φ to be a totally symmetric r + 1 covariant tensorfield on
a 2-sphere (Su,u , g/) satisfying
div φ = f˜, curl φ = g̃, trφ = h̃,
then, for 1 ≤ i ≤ 5, it holds that
i−1
X
kui ∇i φ||L2 (Su,u ) . (||uj+1 ∇j (f˜, g̃)||L2 (Su,u ) + kuj ∇j h̃kL2 (Su,u ) + ||uj ∇j φ||L2 (Su,u ) ).
j=0
As a special case with φ being a symmetric traceless 2-tensor, it holds that curl φ =∗ div φ.
Using the above proposition, we then get
Proposition 3.12. Assume that φ is a symmetric traceless 2-tensor obeying
div φ = f˜.
Then, under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38), (2.39),
(2.40), (2.41), for 1 ≤ i ≤ 5, it holds
i−1
X
i i
||u ∇ φ|| L2 (S u,u )
. (||uj+1 ∇j f˜||L2 (Su,u ) + ||uj ∇j φ||L2 (Su,u ) ).
j=0
Lemma 4.1. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds that
X X
kui1 +i2 +1 ∇i1 ψ i2 +1 kL2 (Su,u ) . kui1 +1 ∇i1 ψ kL2 (Su,u )
i1 +i2 ≤4 i1 ≤4
and
X X
kui1 +i2 +2 ∇i1 ψ i2 +1 kL∞ (Su,u ) . kui1 +2 ∇i1 ψ kL∞ (Su,u ) .
i1 +i2 ≤2 i1 ≤2
Proposition 4.2. With the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
X 1
kui ∇i (Ωχ̂)kL2 (Su,u ) . a 2 .
i≤4
i+1
Ω∇3 ∇i (Ωχ̂) + Ωtrχ∇i (Ωχ̂)
2
X X X 1 i1 i2 +1 i3
=Ω2 ∇i+1 η + ∇i1 ψ i2 +2 + ∇i1 ψ i2 +1 ∇i3 ψ + ∇ ψ ∇ ψ.
i1 +i2 =i i1 +i2 +i3 =i i1 +i2 +i3 =i−1
|u|
Note that the fourth term only presents when i ≥ 1. To bound kui ∇i (Ωχ̂)kL∞ ∞ 2
u Lu L (Su,u )
, we
then apply Proposition 3.6 with λ0 = i+1
2 and estimate each term on the right.
For the ∇i+1 η term, when i = 4 we have
1 1
Ω u2 a2 1
− 14
ku4 Ω2 ∇5 ηkL1u L2 (Su,u ) ≤ k kL2u kΩu5 ∇5 ηkL2u L2 (Su,u ) . 1 Õ5,2 ≤ a B
2 .
|u| |u| 2
The rest terms are lower order. For example, we control the last term as
X X uaO2 1 1
k ui−1 ∇i1 ψ i2 +1 ∇i3 ψ kL1u L2 (Su,u ) . ≤ a 2 B− 4 . (4.1)
i1 +i2 +i3 =i−1
|u|
i≤4
1 1
All the remaining terms on the
P right obey the same upper bound a 2 B − 4 .
1
i
For initial data, we have i≤4 k∇ (Ωχ̂)kL∞ 2
u L (S−1,u )
≤ a 2 . Gathering all above estimates,
we get
X 1
kui ∇i (Ωχ̂)kL∞ 2
u L (Su,u )
. a2 .
i≤4
We then derive the estimates for Ωtrχ and its derivatives. We will also point out that
∇i (Ωtrχ) obeys a better bound than the standard estimate for ∇i ψ .
30 XINLIANG AN
Proposition 4.3. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
X 1
kui+1 ∇i (Ωtrχ−Ω2 (Ω−1 trχ)0 )kL2 (Su,u ) . uaO2 ≪ |u|a 2 ,
i≤4
X uaO2 1
kui ∇i (Ωtrχ)kL2 (Su,u ) . ≪ a2 .
|u|
i≤4
We then employ Proposition 3.5 and control the right hand side in the kui+1 · kL1u L2 (Su,u ) norm
to bound kui+1 ∇i (Ω−1 trχ − (Ω−1 trχ)0 )kL∞ 2
u L (Su,u )
.
Remark 18. In below estimates and for the rest of this article, when there is no danger of
confusion, we use Ω−1 to stand for kΩ−1 kL∞ (Su,u ) .
Employing the bootstrap assumption (2.38) together with Sobolev embedding in Proposition
3.4, for the case i2 = 0, we first control
X X
k ui+1 ∇i3 (χ̂, ∂4 φ)∇i4 (χ̂, ∂4 φ)kL1u L2 (Su,u )
i≤4 i3 +i4 =i
X X
.u kui3 +1 ∇i3 (χ̂, ∂4 φ)kL∞
u Lu L (Su,u )
∞ ∞ kui4 ∇i4 (χ̂, ∂4 φ)kL∞ ∞ 2
u Lu L (Su,u )
. Ω−2 uaO2 .
i3 ≤2 i4 ≤4
P
Similarly, contribution from i1 +i2 +i3 +i4 =i ∇i1 ψ i2 ∇i3 ψ ∇i4 trχ with i2 = 0 obeys
X X
k Ω−1 ui+1 ∇i3 ψ ∇i4 trχkL1u L2 (Su,u )
i≤4 i3 +i4 =i
ua O X i+1 i −1
1
2
. ku ∇ (Ω trχ − (Ω−1 trχ)0 )kL∞ ∞ 2
u Lu L (Su,u )
+ Ω−2 uaO2 .
|u|
i≤4
P 1 i1 i2 i3
When i2 = 0, we can also control the third term i1 +i2 +i3 =i |u| ∇ ψ ∇ ψ by
X X 1 Ω2 O 1
kui ∇i ψ kL1u L2 (Su,u ) . u kui ∇i ψ kL∞ ∞ 2
u Lu L (Su,u )
.ua 2 O = Ω−2 ua · 1 ≤ Ω−2 uaB − 2 .
i≤4 i≤4
a 2
For the case i2 ≥ 1, using the bootstrap assumption (2.38) together with Proposition 4.1, for
the first two terms we also have
X X 3
u2 a 2 O 3 1
i+1 i1 i2 +1 i3 i4
k u ∇ ψ ∇ ψ ∇ ψ kL1u L2 (Su,u ) . ≤ Ω−2 uaB − 4 .
i +i +i +i =i−1
|u|
i≤4 1 2 3 4
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 31
ua 2 O X i+1 i −1
1
We proceed to prove the estimates for ω̂ and its derivatives. And ω̂ also obeys a slightly
better bound than a general component ψ .
Proposition 4.4. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
X 1
X
kui ∇i (Ωω̂)kL2 (Su,u ) . a 2 + kui ∇i (Ω∇4 φ)kL2 (Su,u ) .
i≤4 i≤4
For the terms containing K, if needed, we also employ K . (K − |u|1 2 )+ |u|1 2 . Applying bootstrap
assumptions (2.40),(2.38) and Proposition 4.1, their contributions obey
1 1 1 1 1
u 2 a 2 (R + O) ua 2 O′ 1 1 a 2 B− 4
. 3 + 5 · u 2 a 2 (R + O) ≤ .
|u| 2 |u| 2 |u|
We proceed to control the contribution from |u|1 2 . Note that for this term, the case i2 = 0
happens only when i = 0. For this scenario, the term obeys
1 1
ku−1 2 kL1u L2 (Su,u ) . . (4.2)
|u| |u|
For i2 ≥ 1, it holds
X X 1
1 1
ua 2 O a 2 B − 4
1
i−1 i1 i2 +1
k u ∇ ψ k L 1 L2 (S ) . ≤ .
i1 +i2 =i−1
|u|2 u u,u
|u|2 |u|
i≤4
Here in the second inequality, we employ the bootstrap assumption (2.38). Collecting all the
above estimates, we prove
X X 1 1
1
a 2 B− 4
1
2 i−1 i1 i2 +1 i3
k Ω u ∇ ψ ∇ (K − 2 )kL1u L2 (Su,u ) . + . (4.3)
i +i +i =i−1
|u| |u| |u|
i≤4 1 2 3
For the remaining term, using the bound for ∇i (Ωtrχ) in Proposition 4.3, we have
X ua 1 1 X i i
kui−2 ∇i (Ωtrχ, Ω∇4 φ)kL1u L2 (Su,u ) . 2 + + ku ∇ (Ω∇4 φ)kL2 (Su,u )
|u| |u| |u|
i≤4 i≤4
(4.4)
1 X i i
1 3
u a 2 4
. 3 + ku ∇ (Ω∇4 φ)kL2 (Su,u ) .
|u| 2 |u|
i≤4
Remark 19. Note that in Proposition 5.1, we will prove
X 1
kui ∇i (Ω∇4 φ)kL2 (Su,u ) . a 2 .
i≤4
In particular, it holds
X 1 1
i+1 i
1 ku ∇ ψ kL2 (Su,u ) . 1+O(ω̂) + 1 1 ku5 ∇5 (Ωω̂)kL∞ 2 2
u Lu L (Su,u )
+ R(β).
i≤4
ua 2 u a2
2
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 33
Proof. For ψ ∈ {η, ∇ log Ω, Ωtrχ − (Ωtrχ)0 , Ωχ̂, Ωω − (Ωω)0 }, the below schematic transport
equation holds
1 11
Ω∇4 ψ = Ωβ+∇(Ωω̂) + Ω2 K + Ω2 ∇η + ψ ψ + ψ ψ + ψ. (4.6)
|u|
By commuting the above equation with i angular derivatives, we have
X
Ω∇4 ∇i ψ =∇i (Ωβ) + ∇i+1 (Ωω̂, Ω2 η) + Ω2 ∇i1 ψ i2 ∇i3 K
i1 +i2 +i3 =i
X X 1 i1 i2 i3
+ ∇i1 ψ i2 +1 ∇i3 (ψ , ψ ) + ∇ ψ ∇ ψ.
i1 +i2 +i3 =i i1 +i2 +i3 =i
|u|
Notice ∇(Ωω)0 = 0. Due to the naked-singularity initial data along u = 0, all the quantities ψ
are initially 0. By Proposition 3.5, hence to estimate the quantity kui ∇i ψ kL∞ 2
u L (Su,u )
, we only
need to bound the kui · kL1u L2 (Su,u ) norm of the right hand side. We proceed to control each of
the terms in the equation. We bound the Ωβ term first
X u 2 X i+1 i
1 1
ua 2
i i
ku ∇ (Ωβ)k L1u L2 (Su,u ) ≤ ku ∇ (Ωβ)kLu L (Su,u ) .
2 2 R(β). (4.7)
|u| |u|
i≤4 i≤4
u X i+1 i+1
1
u2 5 5
. ku ∇ (Ωω̂)kL∞ 2
u L (Su,u )
+ ku ∇ (Ωω̂)kL2u L2 (Su,u ) (4.8)
|u| |u|
i≤3
1
ua 2
. (Oi,2 (ω̂) + Õ5,2 (ω̂)).
|u|
For the ∇i+1 η term, using the bootstrap assumptions (2.38) and (2.41), it holds
X
kΩ2 ui ∇i+1 ηkL1u L2 (Su,u )
i≤4
u X i+2 i+1
1
u2
. 2 ku ∇ ψ kL∞ 2
u L (Su,u )
+ kΩ2 u5 ∇5 ηkL2u L2 (Su,u ) (4.9)
|u| |u|
i≤3
1 1 1 1 3 3 3 1
u2 a 2 O u 2 1 1 u2 a4 u2 a 4 O u 2 a 4 ua 2
. 2
+ · Ω2 · u 2 a 4 · 1 · Õ ≤ 2
+ 3 ≤ .
|u| |u| |u| 2 Ω |u| |u| 2 |u|
Employing the bootstrap assumptions (2.38), (2.40), we then bound the term containing the
Gauss curvature K and we write
1 1 1 1
K = (K − 2 − ∇A φ∇A φ) + 2 + ∇A φ∇A φ
|u| 4 |u| 4
to deduce the estimate. It holds
X X ua 2
1
Applying (2.38) to control ψ and using Proposition 4.1 to estimate the product of ψ , as well
as using Sobolev embedding in Proposition 3.4, for i ≤ 4 we bound the fourth term and the
11On the right of the equation (4.6), we have that Ω2 K obeys at least the same and even better estimates
compared with K. When there is no danger of confusion, we sometimes also write K instead of Ω2 K for the
notational convenience and then the Ω2 could be viewed as a coefficient like 1 and is usually omitted in the
schematic writing.
34 XINLIANG AN
X 1
ua 2 u2 aO2
k ui−1 ∇i1 ψ i2 ∇i3 ψ kL1u L2 (Su,u ) . (1 + Oi,2 (ω̂)) + . (4.12)
i1 +i2 +i3 =i
|u| |u|2
1
Notice that here kui ∇i ψ kL∞ 2
u L (Su,u )
≤ a 2 (1 + Oi,2 (ω̂)) was derived before.
Collecting all above estimates, we hence obtain
X 1
ua 2 X
kui ∇i ψ kL∞ 2
u L (Su,u )
. [1+ Oi,2 (ω̂) + Õ5,2 (ω̂) + R(β)].
|u|
i≤4 i≤4
For future use, we also prove an improved bound for ∇i (Ωtrχ − (Ωtrχ)0 ) for i ≥ 1:
Proposition 4.6. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds that
X 3
u2 a4
3 1
ua 2 B − 2
1
i i
ku ∇ (Ωtrχ − (Ωtrχ)0 )kL2 (Su,u ) . 3 ≤ .
|u| 2 |u|
1≤i≤4
Here besides using Proposition 4.1, we also employ the improved bound for ∇i (Ωtrχ) derived in
Propositions 4.3. Collecting this with (4.9) and (4.11), we then obtain the desired conclusion.
We also notice that the ∇4 equation for K − |u|1 2 − 41 ∇A φ∇A φ contains exactly the same type
of terms as the ∇4 equation for ∇ψ . Hence, for ∇i (K − |u|1 2 − 41 ∇A φ∇A φ) we have the same
estimates as for ∇i+1 ψ for i ≤ 3.
Proposition 4.7. With the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40) and (2.41), it holds
X 1 1 1
X
kui+2 ∇i (K − 2 − ∇A φ∇A φ)kL∞ L ∞ L2 (S
u,u ) . ua 2 [1 + Oi,2 (ω̂) + Õ5,2 (ω̂) + R].
|u| 4 u u
i≤3 i≤4
In particular,
X 1 1 A
1
ua 2 X
i+1 i
ku ∇ (K − 2 − ∇ φ∇A φ)kL∞ L ∞ L2 (S
u,u ) . [1 + Oi,2 (ω̂) + Õ5,2 (ω̂) + R] ≪ 1.
|u| 4 u u
|u|
i≤3 i≤4
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 35
Note that the last term only appears when i ≥ 1.Using Proposition 3.6 with λ0 = i+1 2 , we have
i i i 1 2
that ku ∇ (Ωe4 φ)kL2 (Su,u ) is controlled by the ku · kLu L (Su,u ) norm of the right. We first
estimate
X
kui ∇i1 (η + η)i2 ∇i3 (Ω2 ∆g φ)kL1u L2 (Su,u )
i1 +i2 +i3 =i
X
≤kΩ u ∇i+1 (eA φ)kL1u L2 (Su,u ) +kΩ2 ui
2 i
∇i1 (η + η)i2 +1 ∇i3 ∆g φkL1u L2 (Su,u )
i1 +i2 +i3 =i−1
1 1
Ω ua 2 O ua 2 S 2
≤kΩui+1 ∇i+1 (eA φ)kL2u L2 (Su,u ) k kL2u L∞ (Su,u ) + · |u|
|u| |u| |u|3
1 1
u2 a2 S u2 aOS 1 1
≤ 1 + ≤ a 2 B− 2 .
|u| 2 |u|2
Here for the last line we use bootstrap assumption
1 1
ku5 ∇5 (Ωe4 φ)kL2u L2 (Su,u ) + kΩu5 ∇5 (eA φ)kL2u L2 (Su,u ) . u 2 a 2 S.
The terms including − 21 Ωtrχ · Ωe3 φ are lower order and they obey
X 1 ua 1
kui ∇i1 (η + η)i2 ∇i3 (− Ωtrχ · Ωe3 φ)kL1u L2 (Su,u ) ≤ ≤ a 2 B −1 ,
i +i +i =i
2 |u|
1 2 3
where to control ku ∇i (Ωtrχ) · Ωe3 φkL1u L2 (Su,u ) we employ improved estimates in Proposition
i
(4.3).
Using the schematic Codazzi equation, we can rewrite Ωβ = ∇ψ + ψ (Ωtrχ + ψ ). All the rest
1
terms are already estimated in Proposition 4.2 and they are bounded by a 2 . Hence we prove
1
kui ∇i (Ωe4 φ)kL∞ 2
u L (Su,u )
≤ a2 . (5.2)
To derive estimates for eA φ, we employ the equation
∇Ωe4 ∇A φ = eA (Ωe4 φ) − ΩχA B eB φ.
36 XINLIANG AN
Via Proposition 3.5, we then bound kui ∇i (eA φ)kL2 (Su,u ) by the kui ·kL1u L2 (Su,u ) norm of the right
hand side. With the aid of bootstrap assumption 2.40, the term with the top-order derivative
can be controlled by
1
kui ∇i+1 (Ωe4 φ)kL1u L2 (Su,u ) ≤kui+1 ∇i+1 (Ωe4 φ)kL2u L2 (Su,u ) k kL2u L∞ (Su,u )
|u|
1 1
1 1 u2 ua 2
≤u 2 a 2 S · = S.
|u| |u|
All the rest terms can be estimated as in Proposition 4.5. There is neither ∇i (Ωβ) nor ∇i+1 (Ωω̂).
1
And they obey the upper bound ua 2 /|u|. Thus, we arrive at
4
X 1
ua 2 S
kui ∇i (eA φ)kL∞ 2
u L (Su,u )
. .
i=0
|u|
To estimate (Ωe3 φ) − (Ωe3 φ)0 , using (2.9) and commutation formula, we get
(Ω∇4 )∇i (Ωe3 φ − (Ωe3 φ)0 )
X 1 1
= ∇i1 (η + η)i2 ∇i3 (Ω2 ∆g φ − ΩtrχΩe3 φ − ΩtrχΩe4 φ + 2Ω2 η A eA φ)
i1 +i2 +i3 =i
2 2
X X
+ ∇i1 (η + η)i2 ∇i3 (Ωχ)∇i4 (Ωe3 φ) + ∇i1 (η + η)i2 ∇i3 (Ωβ)∇i4 (Ωe3 φ).
i1 +i2 +i3 +i4 =i i1 +i2 +i3 +i4 =i−1
i i
With Proposition 3.5, we control ku ∇ [(Ωe3 φ) − (Ωe3 φ)0 ]kL∞ 2
u L (Su,u )
by the kui · kL1u L2 (Su,u )
norm of the right. The term with the top-order derivative is bounded by
Ω2
kΩ2 ui ∇i+1 (eA φ)kL1u L2 (Su,u ) ≤kui+1 ∇i+1 (eA φ)kL2u L2 (Su,u ) k kL2 L∞ (Su,u )
|u| u
1 1 1 1
1 1 u2 a4 u2 2 ua 2 S
≤u 2 a 4 1 + 1 S · Ω ≤ .
|u| 2 Ω |u| |u|
Here for the second inequality, we use bootstrap assumption 2.40
1 1
u2 a4 1 1
ku5 ∇5 (eA φ)kL2u L2 (Su,u ) + kΩ−1 u5 ∇5 (Ωe3 φ)kL2u L2 (Su,u ) . 1 + 1 u 2 a 4 S.
|u| 2 Ω
Two borderline terms on the right are (Ωe3 φ)∇i (Ωχ) and (Ωtrχ)∇i (Ωe4 φ) with i = 4. And
they obey
1
i i ua 2
ku (Ωe3 φ)∇ (Ωχ)kL1u L2 (Su,u ) ≤Oi,2 (Ωχ) ,
|u|
1
i i ua 2
ku (Ωtrχ)∇ (Ωe4 φ)kL1u L2 (Su,u ) ≤Oi,2 (Ωe4 φ) .
|u|
Note that by Proposition 4.2, Proposition 4.3 and inequality (5.2), we have that Oi,2 (Ωχ),
Oi,2 (Ωe4 φ) depend only on initial data.
All the rest terms can be estimated as in Proposition 4.5. And they are
1 1 1 1
ua 2 ua 2 ua 2 B − 2
≤ · · (O + S) ≤ .
|u| |u| |u|
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 37
We hence arrive at
4
X 1
i i ua 2 S
ku ∇ (Ωe3 φ − (Ωe3 φ)0 )kL2 (Su,u ) . .
i=0
|u|
1 1
DA eB = ∇A eB + χAB e3 + χAB e4 , (6.1)
2 2
D3 eA = ∇3 eA + ηA e3 , D4 eA = ∇4 eA + η A e4 , (6.2)
♯B ♯B
D A e 3 = χA eB + ζA e3 , DA e4 = χA eB − ζA e4 , (6.3)
♯A ♯A
D3 e4 = 2η eA + 2 ω e4 , D4 e3 = 2η eA + 2 ω̂ e3 , (6.4)
D3 e3 = −2ωe3 , D4 e4 = −2ω̂e4 . (6.5)
Next we start to use equation (2.8)
1 1
D3 D4 φ = D4 D3 φ = ∆g φ − trχe3 φ − trχe4 φ.
2 2
With the fact −2ω̂ = ∇4 log Ω = Ω−1 ∇4 Ω and the properties of covariant derivative
D4 D3 φ = D4 (D3 φ) − DD4 e3 φ = D4 (D3 φ) − 2η A eA (φ) + Ω−1 D4 Ω e3 φ,
we can write (2.8) as
1 1
D4 (ΩD3 φ) − 2Ωη A eA φ = ΩD4 D3 φ = ΩD3 D4 φ = Ω∆g φ − Ωtrχe3 φ − Ωtrχe4 φ.
2 2
This implies
1 1
ΩD4 (ΩD3 φ) = Ω2 ∇A ∇A φ − Ωtrχ · Ωe3 φ − Ωtrχ · Ωe4 φ + 2Ω2 η A eA φ. (6.6)
2 2
Similarly, via employing the fact −2ω = ∇3 log Ω = Ω−1 ∇3 Ω and using
D3 D4 φ =D3 (D4 φ) − DD3 e4 φ = D3 (D4 φ) − D2ηA eA +2ωe4 φ
=D3 (D4 φ) − 2η A eA φ − 2ωe4 φ = D3 (D4 φ) − 2η A eA φ + Ω−1 D3 Ωe4 φ,
we also obtain
1 1
ΩD3 (ΩD4 φ) = Ω2 ∇A ∇A φ − Ωtrχ · Ωe3 φ − Ωtrχ · Ωe4 φ + 2Ω2 η A eA φ. (6.7)
2 2
We then proceed to derive the desired form for the transport equation of ∇A φ. For f being a
/f to be the 1-form restricted to each surface Su,u of df (the differential
scalar function, we define d
of f ). Denote L, X being vector fields, we first have
/f )(X) = L((df
(LL d /f )[L, X] = L(Xf ) − [L, X]f = X(Lf ).
/ )(X)) − (d
∂
Letting f = φ, L = ∂u and X = eA , we hence obtain
∂
(L ∂ /dφ)A = d
/A ( φ). (6.8)
∂u ∂u
For (L ∂ /
dφ)A , we also have
∂u
∂ ∂ ∂ ∂
(L ∂ /dφ)A = (eA φ)−[ , eA ]φ = (eA φ)−(D ∂ eA ) φ+(DeA ) φ = D ∂ DeA φ+DeA (Ωe3 )φ.
∂u ∂u ∂u ∂u ∂u ∂u ∂u
(6.9)
For the second pair, multiplying ∇A φ for the first equation and multiplying Ω−2 (Ωe3 φ) for the
second equation, it holds
1
∇A φ∇Ωe3 ∇A φ + Ωtrχ∇A φ∇A φ = ∇A φeA (Ωe3 φ) − Ωχ̂A B ∇B φ∇A φ,
2
Ω−2 (Ωe3 φ)Ωe4 (Ωe3 φ) =Ω−2 (Ωe3 φ)Ω2 ∇A ∇A φ
1 1
+ Ω−2 (Ωe3 φ) · − Ωtrχ · Ωe4 φ − Ωtrχ · Ωe3 φ + 2Ω2 η A · eA φ .
2 2
(6.22)
If we add the two equations in (6.21) and integrate on Su,u , we will encounter the potential
top-order-derivative terms
Z
Ω2 (∇A φ)eA (Ωe4 φ) + (Ωe4 φ)Ω2 ∇A ∇A φ. (6.23)
Su,u
Conducting integration by parts with respect to ∇A for both (6.23) and (6.24), the potential
top-order-derivative terms will be canceled.
Proposition 6.1. Suppose (1) φ and (2) φ are r tensorfields, then it holds
Z Z
Ω−1 (1) φ∇4 (2) φ + Ω−1 (2) φ∇4 (1) φ
Du,u Du,u
Z Z Z
= Ω−1 (1) φ(2) φ − Ω−1 (1) φ(2) φ − Ω−1 trχ(1) φ(2) φ.
H u (−1,u) H 0 (−1,u) Du,u
(1)
Set φ = (2) φ = φ in the above proposition. Together with (6.25) we then obtain
Proposition 6.2. Suppose φ is a r tensorfield, then it holds
Z uZ Z uZ Z uZ u Z Z u Z u Z
2 2
|φ| = |φ| + 2φ(Ω∇4 φ) + Ωtrχφ2 .
−1 Su′ ,u −1 Su′ ,0 0 −1 Su′ ,u′ 0 −1 Su′ ,u′
Proposition 6.3. For an r tensorfield (1) φ and an r − 1 tensorfield (2) φ, they satisfy
Z Z
(1) A1 A2 ...Ar
φ ∇Ar (2) φA1 ...Ar−1 + ∇Ar (1) φA1 A2 ...Ar (2) φA1 ...Ar−1
Du,u Du,u
Z
=− (η + η)(1) φ(2) φ.
Du,u
2λ1 (−Ωe3 u) 2
+ (1 − 2λ0 )Ωtrχ = (1 − 2λ0 )(Ωtrχ + ).
|u| |u|
By integrating with respect to du du and applying the fundamental theorem of calculus in u, we
then prove this proposition.
We now proceed to prove
Proposition 6.5. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
1 1
ku5 ∇5 (Ωe4 φ)kL2u L2 (Su,u ) + kΩu5 ∇5 (eA φ)kL2u L2 (Su,u ) . u 2 a 2 . (6.26)
Proof. Taking equation (6.18)
1 1
Ωe3 (Ωe4 φ) = Ω2 ∇A ∇A φ − Ωtrχ · Ωe4 φ − Ωtrχ · Ωe3 φ + 2Ω2 η A · eA φ,
2 2
and commuting it with i angular derivatives, we obtain
i+1
Ω∇3 ∇i (Ωe4 φ) + Ωtrχ∇i (Ωe4 φ)
2
X 1
= ∇i1 (η + η)i2 ∇i3 Ω2 ∆g φ − Ωtrχ · Ωe3 φ + 2Ω2 η A eA φ
i1 +i2 +i3 =i
2
X
+ ∇i1 (η + η)i2 ∇i3 (Ωtrχ)∇i4 (Ωe4 φ)
i1 +i2 +i3 +i4 =i, (6.27)
i4 <i
X 2
+ ∇i1 (η + η)i2 ∇i3 (Ωχ̂, Ωtrχ + )∇i4 (Ωe4 φ)
i1 +i2 +i3 +i4 =i
|u|
X
+ ∇i1 (η + η)i2 ∇i3 (Ωβ)∇i4 (Ωe4 φ).
i1 +i2 +i3 +i4 =i−1
Note that the term with the top-order derivative on the right is Ω2 ∇i ∆g φ.
Similarly, commuting (6.17) with i angular derivatives, we obtain
X
Ω∇4 ∇i (Ω∇A )φ = ∇i1 (η + η)i2 ∇i3 ΩeA (Ωe4 φ) − Ω2 χA B eB φ + Ωe4 (Ω)∇A φ
i1 +i2 +i3 =i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωχ)∇i4 (Ω∇A )φ
i1 +i2 +i3 +i4 =i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωβ)∇i4 (Ω∇A )φ.
i1 +i2 +i3 +i4 =i−1
(6.28)
i
Notice that the term with the top-order derivative on the right is ∇ [ΩeA (Ωe4 φ)].
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 41
We now consider
Z uZ uZ
i+1
u2i ∇A1 ···Ai (Ωe4 φ)[Ω∇3 ∇A1 ···Ai (Ωe4 φ) + Ωtrχ∇A1 ···Ai (Ωe4 φ)]
0 u0 Su,u 2
Z uZ uZ (6.29)
+ u2i ∇A1 ···Ai (Ω∇A )φ[Ω∇4 ∇A1 ···Ai (Ω∇A )φ].
0 u0 Su,u
Employing integration by parts with respect to ∇A , we then cancel the top-order derivatives
∇i+1 (∇φ), ∇i+1 (Ωe4 φ). Applying Proposition 6.2 for (6.29), together with (6.27) and (6.28) we
arrive at
It is a straight-forward check that the non-top-order-derivative terms in (6.30) obey the same
(and even smaller) upper bounds, compared with the estimates above. There is no borderline
term and their upper bounds are much smaller than ua. Putting all the estimates together, we
then finish the proof of this current proposition.
Next, we move to establish the energy estimate for the other pair and we have
Proposition 6.6. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
1 1
5 5 −1 5 5 u2 a4 1 1
ku ∇ (eA φ)kL2u L2 (Su,u ) + kΩ u ∇ (Ωe3 φ)kL2u L2 (Su,u ) . 1+ 1 u2 a4 .
|u| 2 Ω
In addition, we also have the improved estimate
1 1 1
u2 a4 1 1 ua 2 1 1 1
ku5 ∇5 (eA φ)kL2u L2 (Su,u ) + kΩ−1 u5 ∇5 (Ωe3 φ)kL2u L2 (Su,u ) . 1 · u2 a4 + ( 3 )4 · u2 a4 .
|u| 2 Ω |u|Ω 2
Proof. Employing the paired equations (6.19), (6.20)
1
∇Ωe3 ∇A φ + Ωtrχ∇A φ =eA (Ωe3 φ) − Ωχ̂A B ∇B φ,
2
1 1
∇Ωe4 [(Ωe3 φ) − (Ωe3 φ)0 ] =Ω2 ∆φ − Ωtrχ · Ωe4 φ − Ωtrχ · Ωe3 φ + 2Ω2 η A · eA φ.
2 2
and commutating them with i angular derivatives, for i ≥ 1 we arrive at
i+1
(Ω∇3 )∇i (eA φ) + Ωtrχ∇i (eA φ)
2
X
= ∇i1 (η + η)i2 ∇i3 ∇A (Ω∇3 φ) − Ωχ̂A B ∇B φ
i1 +i2 +i3 =i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωtrχ)∇i4 (eA φ)
i1 +i2 +i3 +i4 =i,
i4 <i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωχ̂, Ωtrχ − (Ωtrχ)0 )∇i4 (eA φ)
i1 +i2 +i3 +i4 =i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωβ)∇i4 (eA φ),
i1 +i2 +i3 +i4 =i−1
i
(Ω∇4 )∇ (Ωe3 φ)
=(Ω∇4 )∇i [(Ωe3 φ) − (Ωe3 φ)0 ]
X
i1 i2 i3 2 1 1 2 A
= ∇ (η + η) ∇ Ω ∆g φ − Ωtrχ · Ωe4 φ − Ωtrχ · Ωe3 φ + 2Ω η eA φ
i1 +i2 +i3 =i
2 2
X
+ ∇i1 (η + η)i2 ∇i3 (Ωχ)∇i4 [(Ωe3 φ) − (Ωe3 φ)0 ]
i1 +i2 +i3 +i4 =i
X
+ ∇i1 (η + η)i2 ∇i3 (Ωβ)∇i4 [(Ωe3 φ) − (Ωe3 φ)0 ].
i1 +i2 +i3 +i4 =i−1
With i ≥ 1, we then consider
Z uZ u Z
i+1
u2i ∇A1 ···Ai (eA φ)[Ω∇3 ∇A1 ···Ai (eA φ) + Ωtrχ∇A1 ···Ai (eA φ)]
0 −1 Su,u 2
Z uZ u Z
+ u2i ∇A1 ···Ai [(Ωe3 φ) − (Ωe3 φ)0 ]{Ω−2 · Ω∇4 ∇A1 ···Ai [(Ωe3 φ) − (Ωe3 φ)0 ]}.
0 −1 Su,u
Using (6.19) and (6.20) into the above expression, we see that the terms involving the top-order
derivatives are of the form
Z
u2i ∇A1 ···Ai (eA φ)∇A1 ···Ai [eA (Ωe3 φ)] + u2i Ω−2 Ω2 ∇A1 ···Ai (Ωe3 φ)∇A1 ···Ai (∆g φ).
Su,u
44 XINLIANG AN
ku10 (Ωχ̂, Ωtrχ − (Ωtrχ)0 )∇5 (eA φ)∇5 (eA φ)kL1u L1u L1 (Su,u )
≤ku5 ∇5 (eA φ)kL∞ 2 2
u Lu L (Su,u )
ku5 ∇5 (eA φ)kL∞ 2 2
u Lu L (Su,u )
kΩχ̂, Ωtrχ − (Ωtrχ)0 kL1u L∞
u L (Su,u )
∞
1 1 1 1 1 1
1 1 u2 a4 1 1 u2 a4 ua 2 1 ua 2 1
≤u 2 a 4 · 1 + 1 S · u2 a4 · 1 + 1 S· O ≤ ua 2 · ( 3 ) ,
2
|u| 2 Ω |u| 2 Ω |u| |u|Ω 2
1 1 1 1 1 1
1 1 u2 a4 1 1 u2 a4 a2 1 ua 2 1
≤u 2 a 4 1 + 1 S · u 2 a4 1 +
1 S · uO ≤ ua 2 · ( )2,
|u| 2 Ω |u| 2 Ω |u| |u|Ω
1 1 1 1 1 1
1 1 u2 a4 1 1 u2 a4 a2 1 ua 2 1
≤u 2 a 4 1 + 1 Õ · u 2 a4 1 +
1 S · uO ≤ ua 2 · (
3 ) .
2
|u| 2 Ω |u| 2 Ω |u|Ω |u|Ω 2
For the top-order-derivative terms containing ∇5 (Ωe3 φ) and η, we have
ku10 η A ∇5 (eA φ)∇5 (Ωe3 φ)kL1u L1u L1 (Su,u )
≤ku5 ∇5 (eA φ)kL∞ 2 2
u Lu L (Su,u )
kΩ−1 u5 ∇5 (Ωe3 φ)kL∞ 2 2
u Lu L (Su,u )
kΩη A kL2u L2u L∞ (Su,u )
1 1 1 1 1 1
1 1 u2 a4 1 1 u2 a4 ua 2 1 1 1 ua 2
≤u 2 a 4 1 + 1 S · u2 a4 1 + 1 S· 2
|u| 2 u 2 O ≤ ua 2 · ,
|u| 2 Ω |u| 2 Ω |u| |u|Ω
The top-order-derivative terms involving ∇5 (Ωe3 φ) and Ωe3 φ − (Ωe3 φ)0 obey
kΩ−2 u10 [(Ωe3 φ) − (Ωe3 φ)0 ]∇5 (Ωtrχ)∇5 (Ωe3 φ)kL1u L1u L1 (Su,u )
≤ku5 ∇5 (Ωtrχ)kL∞ 2 2
u Lu L (Su,u )
kΩ−1 u5 ∇5 (Ωe3 φ)kL∞ 2 2
u Lu L (Su,u )
kΩ−2 u10 [(Ωe3 φ) − (Ωe3 φ)0 ]∇5 (Ωχ̂)∇5 (Ωe3 φ)kL1u L1u L1 (Su,u )
≤ku5 ∇5 (Ωχ̂)kL∞ 2 2
u Lu L (Su,u )
kΩ−1 u5 ∇5 (Ωe3 φ)kL∞ 2 2
u Lu L (Su,u )
kΩ−1 [(Ωe3 φ) − (Ωe3 φ)0 ]kL2u L2u L∞ (Su,u )
1 1 1 1
1 1 1 1 u2 a4 ua 2 1 1 1 ua 2
≤u a Õ · u a
2 2 2 4 1+ 1 S · 2 |u| u O ≤ ua ·
2 2 2 .
|u| 2 Ω |u| Ω |u|Ω
And the top-order-derivative term with ∇5 (Ωe3 φ) and ∇4 (Ωβ) obey
kΩ−2 u10 [(Ωe3 φ) − (Ωe3 φ)0 ]∇4 (Ωβ)∇5 (Ωe3 φ)kL1u L1u L1 (Su,u )
≤ku5 ∇4 (Ωβ)kL∞ 2 2
u Lu L (Su,u )
kΩ−1 u5 ∇5 (Ωe3 φ)kL∞ 2 2
u Lu L (Su,u )
Note that in above inequality Ωω † obeys the same bounds as for ψ . Hence we regard Ωω † as
one of the ψ terms and we employ the notation ψ ∈ {Ωtrχ, Ωχ̂, Ωω̂, Ωω † }. Using this notation,
we write the transport equation for Ωκ as
1
Ω∇3 (Ωκ) + Ωtrχ(Ωκ)
2
X 2 X
= ψ i1 ∇i2 (Ωtrχ + , Ωχ̂) · ∇i3 ψ + Ω2 ψ i1 +1 ∇i2 (η, η, eA φ)
i +i +i =1
|u| i +i =1
1 2 3 1 2
1 Ω 1
+ ψψ + β + ψ K + ∇(Ωe3 φ)ψ + ∇(Ωe4 φ).
|u| |u| |u|
Commuting
P with angular derivatives for 4 times and employing the schematic Codazzi equation
Ωβ = i1 +i2 =1 ψ i1 ∇i2 ψ , we obtain
5
Ω∇3 ∇4 (Ωκ) + Ωtrχ∇4 (Ωκ)
2
X 2 X
i1 i2 i3
= ∇ ψ ∇ (Ωtrχ + , Ωχ̂) · ∇i4 ψ + Ω2 ∇i1 ψ i2 +1 ∇i3 (η, η, eA φ)
i +i +i +i =5
|u| i
1 2 3 4 1 +i2 +i3 =5
i1 ≤4
1 4 1 X X
+ ∇ (Ωβ) + ∇i1 ψ i2 +1 ∇i3 ψ + ∇i1 +1 ψ i2 +1 ∇i3 ψ
|u| |u| i
1 +i2 +i3 =4 i1 +i2 +i3 =4
i1 ≤3
X 1 5
+ ∇i1 ψ i2 ∇i3 ψ ∇i4 K+ψ ∇5 (Ωe3 φ) + ∇ (Ωe4 φ).
i1 +i2 +i3 +i4 =4
|u|
X 3
u 2 aO2
1
u 2 a 2 B− 2
1 1
.ua(OÕ + O2 ).
This gives
3 3
ku6 ∇5 (Ωtrχ)kL2u L2 (Su,u ) . u 2 a(OÕ + O2 ) ≪ u 2 aÕ′ . (7.3)
It also implies
1 1
ku5 ∇5 (Ωtrχ)kL2u L2 (Su,u ) . u 2 a 2 . (7.4)
We further take the L2 norm in u. With the bound (7.3) for ∇5 trχ, we deduce
X
ku5 ∇5 (Ωχ̂)kL2u L2 (Su,u ) . kui ∇i (Ωtrχ)kL2u L2 (Su,u )
i≤5
X Ω X
+ kui+1 ∇i (Ωβ− ∇4 φ∇A φ)kL2u L2 (Su,u ) + kui+1 ∇i (Ω∇4 φ∇A φ)kL2u L2 (Su,u )
2
i≤4 i≤4
X X X
i+1 i1 i2
+ ku ∇ ψ ∇ ψ kL2u L2 (Su,u ) + kui ∇i (Ωχ̂)kL2u L2 (Su,u )
i≤4 i1 +i2 =i i≤4
3
2
1 1 1 1 u aO2 1 1
.u 2 a 2 + u 2 a 2 R(β) + . u 2 a 2 [1 + R(β)].
|u|
With the above estimate for u5 ∇5 (Ωχ̂), together with Proposition 4.3, Proposition 5.1, Propo-
sition 6.5, we then revisit (7.2) and obtain
1
ku6 ∇5 (Ωtrχ)kL2 (Su,u ) ≤ u 2 k(Ωtrχ, Ωω̂)kL∞
u L (Su,u )
∞ ku6 ∇5 (Ωtrχ)kL2u L2 (Su,u )
1
+ u 2 kuΩtrχkL∞
u L (Su,u )
∞ ku5 ∇5 (Ωω̂)kL2u L2 (Su,u )
1
+ u 2 kuΩχ̂kL∞
u L (Su,u )
∞ ku5 ∇5 (Ωχ̂)kL2u L2 (Su,u )
1
+u 2 kuΩ∇4 φkL∞
u L (Su,u )
∞ ku5 ∇5 (Ω∇4 φ)kL2u L2 (Su,u )
X
+ ku6 ∇i1 ψ i2 ∇i3 ψ ∇i4 ψ kL∞ 1 2
u Lu L (Su,u )
i1 +i2 +i3 +i4 =5
i3 ,i4 ≤4
1 1
u2 a2 3 1 1 1 1 1 1 1 1 1 1 1
≤ · u 2 aÕ + u 2 Ω2 · u 2 a 4 [1 + R(β)] + u 2 a 2 · u 2 a 2 [1 + R(β)]+u 2 a 2 · u 2 a 2
|u|
≤ua[1 + R(β)].
3
ku6 ∇5 (Ωtrχ)kL2u L2 (Su,u ) . u 2 a[1 + R(β)]. (7.5)
Proposition 7.3. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
1 1
1 1 u2 a4
ku5 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
. u2 a4 · 1 .
|u| Ω2
1 1
div η = −µ + K − 2
− ∇A φ∇A φ, curl η = −σ̌.
|u| 4
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 51
1 1 3 1 1
∇3 (K − − ∇A φ∇A φ) + trχ(K − 2 − ∇A φ∇A φ)
|u|2 4 2 |u| 4
1 1
=div (β + ∇3 φ∇φ) − ∇A ∇3 φ∇A φ − ζ · β + 2η · β
2 2
1 1 1 1
b − trχ|η|2 + η A ∇A φ∇3 φ
b + χ̂ · (η ⊗η)
+ χ̂ · ∇⊗η
2 2 2 2
1 1 1
+ (∇A φ∇B φ − gAB ∇C φ∇C φ) · χ̂AB − trχ∇A φ∇A φ
2 2 8
1 1 1 Ω−1 2
+ trχ(−div η + K − 2 − ∇A φ∇A φ) − 2
(Ωtrχ + ).
2 |u| 4 |u| |u|
Ω
Ω∇3 µ + Ωtrχ µ =ψ ∇(η, η) + ψ ψ ψ + ψ ∇(Ωχ̂, Ωtrχ) + Ωtrχdiv η + ΩtrχK +
|u|3
+Ωtrχ ψ ψ + ∆g φ(Ωe3 φ) + eA φeA (Ωe3 φ) + ηeA φ(Ωe3 φ).
Ω∇3 ∇4 µ + 3Ωtrχ∇4 µ
Ω 1 X
= ∇5 η + Ωe3 φ · ∇5 (eA φ) + eA φ · ∇5 (Ωe3 φ) + ψ ∇5 (η, η, Ωtrχ, Ωχ̂) + ∇i1 ψ i2 +1 ∇i3 ψ
|u| |u| i +i +i =4
1 2 3
1 X X X
+ ∇i1 ψ i2 ∇i3 K + ∇i1 ψ i2 +2 ∇i3 ψ + ∇i1 ψ i2 +1 ∇i3 K
|u| i +i +i =4 i1 +i2 +i3 =4 i1 +i2 +i3 =4
1 2 3
X X
i1 i2 i3 +1 i4
+ ∇ ψ ∇ eA φ∇ (Ωe3 φ) + i1 i2 i3
∇ ψ ∇ eA φ∇i4 +1 (Ωe3 φ).
i1 +i2 +i3 =4 i1 +i2 +i3 =4
i3 ≤3 i4 ≤3
1 1
1 1 u2 a4
The remaining lower-order terms are all bounded by u 2 a 4 · 1 . Hence, combining the
|u| 2 Ω
above estimates, we arrive at
1 1
1 1 u2 a4
ku5 ∇4 µkL2u L∞ 2
u L (Su,u )
≤ u2 a4 · 1 .
|u| 2 Ω
By using the div-curl system
1 1
div η = −µ + K − − ∇A φ∇A φ, curl η = −σ̌
|u|2 4
and employing elliptic estimates from Proposition 3.11, we further deduce
1 1
ku5 ∇5 ηkL2u L2 (Su,u ) .ku5 ∇4 µkL2u L2 (Su,u ) + ku5 ∇4 (K − 2 − ∇A φ∇A φ, σ̌)kL2u L2 (Su,u )
|u| 4
X 1 1
+ kui+1 ∇i (−µ + K − 2 − ∇A φ∇A φ, σ̌)kL2u L2 (Su,u )
|u| 4
i≤3
1 1 3 3 1 1 1 1
1 1 u2 a4
u2 a4 u2 a4 1 1 u2 a4
≤u a ·2
1
4 + 1+ 1 R ≤ u2 a4 · 1 .
|u| Ω
2 |u| |u| Ω
2 |u| 2 Ω
We proceed to obtain the highest order derivative estimates for η. In the proof we will derive
an improved bound for ∇i µ with 1 ≤ i ≤ 4.
Proposition 7.4. With the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41) and the proved O5,2 (η) . 1, we have
1 1
3 3 u2 a4 3 3 1
ku6 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
. u2 a4 1 + 1 [1 + R(K, σ̌)] + u 2 a 4 · a− 4 R(β),
|u| 2 Ω
1 1
ku5 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
. Ω−1 u 2 a 2 [1 + R(K, σ̌)],
1 1
ku5 Ω∇5 ηkL∞ 2 2
u Lu L (Su,u )
. u 2 a 2 [1 + R(K, σ̌)].
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 53
Employing Proposition 4.1 together with the bootstrap assumption (2.39), we then bound the
1 1
3 u2 a4
other lower order terms and they obey . ua 4 · 1 .
|u| 2 Ω
Combining all the estimates above, we obtain
X 3
1
u2 a4
1
kui+2 ∇i µkL∞ 2
u L (Su,u )
.ua 4 · 1 . (7.7)
1≤i≤4
|u| Ω2
3
1 1
u2 a4 X 1 1 1
.ua 4 1 + 1 · R(K) + kui+2 ∇i (K − 2
− ∇A φ∇A φ, σ̌)kL2 (Su,u ) + ua 2 R(β),
|u| Ω2 |u| 4
i≤4
We then take the L2 norm in u and u respectively. Via using (7.7), we further obtain
1 1
3 3 u2 a4 3 3 1
ku6 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
.u 2 a 4 1 + 1 · R(K, σ̌) + u 2 a 4 · a− 4 R(β),
|u| Ω 2
1 1 1 1 1 1
kΩu5 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
. u 2 a 2 R(K, σ̌) + u 2 a 2 . u 2 a 2 [1 + R(K, σ̌)]. (7.8)
This also infers
1 1
ku5 ∇5 ηkL∞ 2 2
u Lu L (Su,u )
≤ Ω−1 u 2 a 2 · [1 + R(K, σ̌)].
We hence finish the proof of this proposition.
Finally, we move to bound the highest order estimates for trχ and χ̂.
Proposition 7.5. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
1 1
5 5 1 u2 a4
1
ku ∇ (Ωtrχ, Ωχ̂)kL∞ 2 2
u Lu L (Su,u )
.u a 1 +2 4
1 ,
|u| 2 Ω
1 1
1 1 u2 a4
ku5 ∇5 trχkL∞ 2 2
u Lu L (Su,u )
.u a 1 +
2 4
1 ,
|u| 2 Ω
1 1
1 1 u2 a4
ku5 ∇5 χ̂kL∞ 2 2
u Lu L (Su,u )
.u 2 a 4 1 + 1 · R(β).
|u| 2 Ω
Proof. We utilize the following equation for trχ:
1
Ω∇3 (Ω−1 trχ) + Ωtrχ(Ω−1 trχ) = −|χ̂|2 − ∇3 φ∇3 φ.
2
Via commuting the above equation with 5 angular derivatives, we obtain
Ω∇3 ∇5 (Ω−1 trχ) + 3Ωtrχ∇5 (Ω−1 trχ)
X
=ψ ∇5 (Ω−1 trχ, Ω−1 χ̂) + Ω−1 e3 φ∇5 (Ωe3 φ) + Ω−2 ∇i1 ψ i2 +1 ∇i3 ψ
i1 +i2 +i3 =5
i1 ,i3 ≤4 (7.9)
X Ω −2 X Ω −2
+ ∇i1 ψ i2 +1 ∇i3 ψ + 2
∇i1 ψ i2 +1 .
i1 +i2 +i3 =4
|u| i
|u|
1 +i2 =4
X 1
u2 a4
1 1
u2 a4
1
i1 ,i3 ≤4
X Ω−2 i1 i2 +1 i3 u2 a4 u2 a4
1 1 1 1
X Ω−2 i1 i2 +1 u2 a4 u2 a4 1
1 1 1 1
1 1
1 1 u a
2 4
ku5 ∇5 (Ωtrχ)kL∞ 2 2
u Lu L (Su,u )
≤ u 2 a 4 (1 + 1 ),
|u| 2 Ω
1 1 1 1 1 1 1 1 1 1 1
(7.10)
5 5 u2 a4 u2 a4 ua 2 1 u2 a4 u2 a4 u2 a4 1
ku ∇ trχkL2 (Su,u ) ≤ 1 · 1 +( 3 ) · 4
1 + 1 ( 1 1 ) , 2
|u| 2 Ω |u| 2 |u|Ω 2 |u| 2 |u| 2 |u| 2 Ω 2
1 1
1 1 u2 a4
ku5 ∇5 trχkL∞ 2 2
u Lu L (Su,u )
≤ u 2 a 4 (1 + 1 ).
|u| 2 Ω
To bound the fifth angular derivatives of χ̂, we now employ the Codazzi equation
1 1 1
div χ̂ = β + ∇trχ − (η − η) · (χ̂ − trχ).
2 2 2
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 57
X 1
ua 2 O(η, η)
5 5 i+1 i+1 i+1 i
ku ∇ χ̂kL2 (Su,u ) . ku ∇ trχkL2 (Su,u ) + ku ∇ βkL2 (Su,u ) + .
|u|Ω
i≤4
2
We then take the L norm in u and obtain
ku5 ∇5 (Ωχ̂)kL∞ 2 2
u Lu L (Su,u )
1 1
1 1
u2 a4 X 1
ua 2 O
.u 2 a 4 1 + 1 + kui+2 ∇i βkL∞ 2 2
u Lu L (Su,u )
ku−1 ΩkL∞
u Lu L (Su,u )
∞ ∞ + 1
|u| 2 Ω i≤4
|u| 2
1 1 1 1 1 1 1
1 1 u2 a4 3 3 u2 a4 1 ua 2 O 1 1 u2 a4
.u 2 a 4 (1 + 1 ) + u2 a4 1 + 1 R· + 1 ≤ u 2 a 4 (1 + 1 ),
|u| Ω 2 |u| 2 Ω |u| |u| 2 |u| 2 Ω
ku5 ∇5 χ̂kL∞ 2 2
u Lu L (Su,u )
X 1
i+1 i+1 i+1 i
ua 2 O(η, η)
. ku ∇ trχkL∞ 2 2
u Lu L (Su,u )
+ ku ∇ βkL∞ 2 2
u Lu L (Su,u )
+ 1
i≤4
|u| 2 Ω
X X 1
ua 2 R(β)
i+1 i+1 i+2 i −1
. ku ∇ trχkL∞ 2 2
u Lu L (Su,u )
+ ku ∇ βkL∞ 2 2
u Lu L (Su,u )
· ku kL∞
u Lu L (Su,u )
∞ ∞ + 1
i≤4 i≤4
|u| 2 Ω
1 1 1 1
1 1 u2 a4 1 1 u2 a4
.u 2 a 4 (1 + 1 ) + u2 a4 · 1 · R(β),
|u| 2 Ω |u| 2 Ω
where we use Proposition 4.5 and Proposition 7.1.
Collecting all above estimates, we hence prove this proposition.
We conclude this section by summarizing all the proved estimates in this section with the aid
of the Õ5,2 norm:
Proposition 7.6. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
Õ5,2 . 1 + R.
3 1
∇4 (Ωσ̌) + trχΩσ̌ = − Ωdiv ∗ (β − ∇4 φ∇φ) + Ω(−∇1 R42 + ∇2 R41 ) − 2Ωω̂σ̌
2 2
1 1
+ Ω[−ζ · β − 2η ·∗ β + χ̂ ·∗ (∇⊗η)
∗ b + χ̂ ·∗ (η ⊗η)]
b (8.2)
2 2
1 1
+ Ω(∇A φ∇B φ − gAB ∇C φ∇C φ) ∧ χ̂AB ,
2 2
1 1 1 1
∇4 [Ω(K − − ∇A φ∇A φ)] + trχΩ(K − 2 − ∇A φ∇A φ)
|u|2 4 |u| 4
1 1 1 1 1
= − Ωdiv (β − ∇4 ∇φ) − Ω∇A ∇4 φ∇A φ − 2Ωω̂(K − 2 − ∇A φ∇A φ) − Ωtrχ · 2
2 2 |u| 4 |u| (8.3)
1 1 1 1
b − trχdiv η − trχ|η|2 ]
b + χ̂ · (η ⊗η)
+ Ω[−ζ · β − 2η · β + χ̂ · ∇⊗η
2 2 2 2
1 A 1 1 C AB
+ Ω[ η ∇A φ∇4 φ + (∇A φ∇B φ − gAB ∇C φ∇ φ) · χ̂ ].
2 2 2
The second pair are
3 1 1 1
∇3 σ̌ + trχσ̌ = − div ∗ (β + ∇3 φ∇φ) − ζ ·∗ β − 2η ·∗ β + χ̂ ·∗ (∇⊗η)
b + χ̂ ·∗ (η ⊗η)
b
2 2 2 2 (8.4)
1 1
+ ∇1 R32 − ∇2 R31 + χ̂AB ∧ (∇A φ∇B φ − gAB ∇C φ∇C φ),
2 2
1 1 3 1 1
∇3 (K − 2
− ∇A φ∇A φ) + trχ(K − 2 − ∇A φ∇A φ)
|u| 4 2 |u| 4
1 1 A
=div (β + ∇3 φ∇φ) − ∇ ∇3 φ∇A φ − ζ · β + 2η · β
2 2
1 1 1 1
+ χ̂ · ∇⊗η b − trχ|η|2 + η A ∇A φ∇3 φ
b + χ̂ · (η ⊗η) (8.5)
2 2 2 2
1 1 1
+ (∇A φ∇B φ − gAB ∇C φ∇C φ) · χ̂AB − trχ∇A φ∇A φ
2 2 8
1 1 1 Ω−1 2
+ trχ(−div η + K − 2 − ∇A φ∇A φ) − 2
(Ωtrχ + ),
2 |u| 4 |u| |u|
1 1
∇4 (β · + ∇3 φ∇· φ)A + trχ(β A + ∇3 φ∇A φ)
2 2
1 1
=∇(K − 2 − ∇A φ∇A φ) +∗ ∇σ̌ + 2ω̂β + 2χ̂ · β + 3(−ηK +∗ ησ̌)
|u| 4
(8.6)
1 1 1 1 1
+ ∇(− ∇A φ∇A φ + trχtrχ − χ̂ · χ̂) − ∗ ∇(χ̂ ∧ χ̂)+ trχ∇3 φ∇A φ
4 4 2 2 2
CD B 1 C
+ ∇A (g RCD ) + ∇ RBA − η A ∇3 φ∇4 φ − η A ∇C φ∇ φ.
2
Via applying Propositions 6.3-6.4, for the first pair we have
Proposition 8.1. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
X 1 1 1 1 1
(kui+1 ∇i (Ωβ − Ω∇4 φ∇φ)kL∞ 2 2
u Lu L (Su,u )
+ kΩui+1 ∇i (K − 2 − ∇A φ∇A φ, σ̌)kL∞ 2 2
u Lu L (Su,u )
) . u2 a2 .
2 |u| 4
i≤4
X Ω−1 X X
+ ∇i1 ψ i2 +1 ∇i3 ψ + ∇i1 ψ i2 +1 ∇i3 ψ + ∇i1 ψ i2 +2
i1 +i2 +i3 =i+1
|u| i1 +i2 +i3 =i i1 +i2 =i+1
i1 ,i3 ≤i i1 ≤i
1 X
+ ∇i1 ψ i2 +1 +Ω∇i+2 φ∇φ,
|u|2 i1 +i2 =i
1
∇4 ∇i (Ωσ̌) + Ωdiv ∗ ∇i (β − ∇4 φ∇φ) = G2,i , (8.8)
2
where G2,i is
X
G2,i =ψ ∇i+1 ψ + ψ ∇i+1 η + ∇i1 ψ i2 ∇i3 ψ ∇i4 σ̌
i1 +i2 +i3 +i4 =i
X
i1 i2 +1
+ ∇ ψ ∇i3 ψ +∇i+1 (Ωe4 φ)∇φ,
i1 +i2 +i3 =i+1
i1 ,i3 ≤i
1 1 1
∇4 ∇i [Ω(K − − ∇A φ∇A φ)] + Ωdiv ∇i (β − ∇4 φ∇φ) = G3,i (8.9)
|u|2 4 2
with G3,i satisfying
X 1
G3,i =ψ ∇i+1 ψ + ψ ∇i+1 η + ∇i1 ψ i2 ∇i3 ψ ∇i4 (K − , σ̌)
i1 +i2 +i3 +i4 =i
|u|2
X 1 X
+ ∇i1 ψ i2 +1 ∇i3 ψ + ∇i1 ψ i2 ∇i3 ψ + ∇i+1 (Ωe4 φ)∇φ.
i1 +i2 +i3 =i+1
|u|2 i1 +i2 +i3 =i
i1 ,i3 ≤i
Note that by the structure of the renormalized null Bianchi equation, we will have −∇1 R42 +
∇2 R41 , which enables us to avoid the troublesome term ∇i+2 φe4 φ in above two equations.
We now employ Proposition 6.2 for φ = ui+1 ∇i [Ω(K− |u|1 2 − 41 ∇A φ∇A φ)] and φ = ui+1 ∇i [Ωσ̌],
and apply Proposition 6.4 for φ = ∇i [Ω(β − 12 ∇4 φ∇φ)] with λ1 = i + 1. We then add these
identities together and obtain
X 1
kui+1 ∇i (Ωβ)kL2u L2 (Su,u ) + kΩui+1 ∇i (K − 2 )kL2 (H) + kΩui+1 ∇i σ̌kL2u L2 (Su,u )
|u|
i≤4
X X
. kui+1 ∇i (Ωβ)kL2u L2 (S0,u ) + kΩui+1 G1,i kL1u L2u L2 (Su,u ) + kΩui+1 G2,i kL1u L2u L2 (Su,u )
i≤4 i≤4
X
i+1
+ kΩu G3,i kL1u L2u L2 (Su,u ) .
i≤4
(8.10)
60 XINLIANG AN
Note that the highest derivatives of the scalar field in in G1,i satisfies
kΩui+1 ∇i+1 ∇φ∇φkL1u L2u L2 (Su,u ) ≤ kui+1 ∇i+1 ∇φkL∞ 2 2
u Lu L (Su,u )
kΩ∇φkL1u L∞
u L (Su,u )
∞
1 1 1
1 u2 a4
1 ua 2 O 1 1
≤u a · 1 +
2 2
1 S· ≤ u2 a2 .
|u| Ω
2 |u|
The corresponding terms in G2,i and G3,i obey
kΩui+1 ∇i+1 (Ωe4 φ)∇φkL1u L2u L2 (Su,u ) ≤ kui+1 ∇i+1 (Ωe4 φ)kL∞ 2 2
u Lu L (Su,u )
kΩ∇φkL2u L2u L∞ (Su,u )
1 3 1
1 ua 2 1 1 1 1 1 u2 a2 1 1
≤u a S ·
2 |u| 2 u 2 O = u 2 a 2
2
3 SO ≤ u a .
2 2
|u|2 |u| 2
We then proceed to bound other terms. Among the terms in Ωui+1 G1,i , we first control the
contributions from the fifth derivatives the Ricci coefficients. These terms are
1 5
ψ ∇5 (Ωtrχ, Ωχ̂), ∇ (Ωtrχ), ψ ∇5 (Ωχ̂, Ωtrχ) (8.11)
|u|
and in ku5 · kL1u L2u L2 (Su,u ) norm they satisfy
1
1 1 ua 2 O 1 1
.u 2 a 2 Õ · ≤ u2 a2 ,
|u|
3
4 5 6 5 1 u 2 aÕ 1 1
ku ∇ (Ωtrχ)kL1u L2u L2 (Su,u ) .ku ∇ (Ωtrχ)kL∞ 2 2
u Lu L (Su,u )
k 2 kL1u . ≤ u2 a2 .
|u| |u|
Here we use the improved bound for u6 ∇5 (Ωtrχ) in Proposition 7.2, which is better than other
∇5 (Ωψ ) components. For the terms containing (Ωχ̂, Ωtrχ), employing the bootstrap assumption
(2.40), we get
ku5 ψ ∇5 (Ωχ̂, Ωtrχ)kL1u L2u L2 (Su,u ) . ku5 ∇5 (Ωχ̂, Ωtrχ)kL∞ 2 2
u Lu L (Su,u )
kψ kL2u L2u L∞ (Su,u )
1 1 1 1 1 1 1 1
1 1 u2 a4 u2 a2 O 1 1 u2 a4 u2 a4 1 1
.u 2 a 4 1 + 1 Õ · 1 ≤ u 2 a2
1 1 + 1 ÕO ≤ u 2 a 2 .
|u| 2 Ω |u| 2 |u| 2 |u| 2 Ω
For kΩui+1 G1,i kL1u L2u L2 (Su,u ) , the other terms obey
X 1 ua 2 O 12 12 1 1
1
1 1
kΩui+1 ∇i1 ψ i2 +1 ∇i3 (K − 2
, σ̌)kL1u L2u L2 (Su,u ) ≤ |u| u · u 2 a 2 ≤ u 2 a 2 ,
i1 +i2 +i3 =i
|u| |u|2
X 1
ua 2 O a 2 O 1
1
1 1
kΩui+1 ∇i1 ψ i2 +1 ∇i3 ψ kL1u L2u L2 (Su,u ) ≤ 2
· · |u|2 u 2 ≤ u 2 a 2 ,
i1 +i2 +i3 =i+1
|u| |u|
i1 ,i3 ≤i
X 1
ua 2 O a 2 O 1 1
1
1
i i1 i2 +1 i3
ku ∇ ψ ∇ ψ kL1u L2u L2 (Su,u ) ≤ · · |u|2 u 2 ≤ u 2 a 2 ,
i1 +i2 +i3 =i
|u|2 |u|
X 1
ua 2 · ua 2 O2 1
1
1 1
i+1 i1 i2 +2
kΩu ∇ ψ kL1u L2u L2 (Su,u ) ≤ · |u|2 u 2 ≤ u 2 a 2 ,
i1 +i2 =i+1
|u|4
i1 ≤i
X ua 2 O
1
1 1 1
i−1 i1 i2 +1
kΩu ∇ ψ k L1u L2u L2 (Su,u ) ≤ · |u|2 u 2 ≤ u 2 a 2 .
i1 +i2 =i
|u|3
We then proceed to derive the estimates for G2,i and G3,i . Since all the terms of G2,i are
contained in the expression for G3,i , hence we only need to bound the terms in G3,i . The terms
containing the 5 derivatives of the Ricci coefficients are
Ωu5 ψ ∇5 ψ , Ωu5 ψ ∇5 η,
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 61
The rest terms in kΩui+1 (G2,i , G3,i )kL1u L2u L2 (Su,u ) are lower order and they satisfy
X 1
1
ua 2 O 1 1 1 1
kΩui+1 ∇i1 ψ i2 ∇i3 ψ ∇i4 (K − 2
, σ̌)kL1u L2u L2 (Su,u ) ≤ · u2 a2 ≤ u2 a2 ,
i1 +i2 +i3 +i4 =i
|u| |u|
X 1
ua 2 O a 2 O
1
3 1 1
kΩui+1 ∇i1 ψ i2 +1 ∇i3 ψ kL1u L2u L2 (Su,u ) ≤ 2
· · |u| 2 u ≤ u 2 a 2 ,
i1 +i2 +i3 =i+1
|u| |u|
i1 ,i3 ≤i
X a2
1
3 1 1
kΩui−1 ∇i1 ψ i2 ∇i3 ψ kL1u L2u L2 (Su,u ) ≤ |u| 2 uO ≤ u 2 a 2 .
i1 +i2 +i3 =i
|u|2
Proposition 8.2. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40), (2.41), it holds
X 1 1 A 1
i+2 i i+2 i
ku ∇ (K − 2 − ∇ φ∇A φ, σ̌)kL∞ 2 2
u Lu L (Su,u )
+ ku ∇ (β + ∇3 φ∇φ)kL∞ 2 2
u Lu L (Su,u )
|u| 4 2
1≤i≤4
1 1
3 3 u2 a4
.u 2 a 4 1 + 1 .
|u| 2 Ω
By further employing the L2 (Su,u ) estimates proved in Section 5.1 for φ, it also holds
X 1
3 3
1
u2 a4
1
i+2 i i+2 i
ku ∇ (K − 2 , σ̌)kL∞ 2 2
u Lu L (Su,u )
+ ku ∇ βkL∞ 2 2
u Lu L (Su,u )
.u a 1+
2 4
1 .
|u| |u| 2 Ω
1≤i≤4
Proof. We first commute the renormalized null Bianchi equations with i angular derivatives and
get
3+i 1
∇3 ∇i σ̌ + trχ∇i σ̌ + div ∗ ∇i (β + ∇3 φ∇φ) = F1,i , (8.12)
2 2
where F1,i is
X Ω−1 X
F1,i =∇i+1 (e3 φ)∇φ + Ω−1 ∇i1 ψ i2 +1 ∇i3 ψ + ∇i1 ψ i2 +1 ∇i3 σ̌,
i1 +i2 +i3 =i+1
|u| i1 +i2 +i3 =i−1
1 1 1 3+i 1 1
∇3 ∇i (K − − ∇A φ∇A φ)−div ∇i (β + ∇3 φ∇φ)+ trχ∇i (K − 2 − ∇A φ∇A φ) = F2,i
|u|2 4 2 2 |u| 4
(8.13)
62 XINLIANG AN
1 1 1
∇4 ∇i (β + ∇3 φ∇φ) − div ∇i (K − 2 − ∇A φ∇A φ) −∗ ∇∇i σ̌ = F3,i , (8.14)
2 |u| 4
where F3,i satisfies
F3,i =Ω−2 (Ωtrχ)∇i+1 (Ωtrχ) + Ω−2 · Ωχ̂∇i+1 (Ωχ̂, Ωtrχ) + Ω−2 · Ωχ̂∇i+1 (Ωχ̂)
+ Ω−2 (Ωtrχ)∇i+1 (Ωtrχ) + Ω−2 · Ωω̂∇i+1 (Ωχ̂, Ωtrχ − (Ωtrχ)0 ) + ∇i+1 (∇φ)∇φ
X
+ Ω−2 ∇i1 ψ i2 ∇i3 (Ωχ̂, Ωtrχ − (Ωtrχ)0 )∇i4 (Ωtrχ, Ωχ̂, Ωω̂)
i1 +i2 +i3 +i4 =i+1
i3 ,i4 ≤i (8.15)
X Ω−2 X
+ ∇i1 ψ i2 +1 ∇i3 (K, σ̌) + ∇i1 ψ i2 +1 ∇i3 ψ
i1 +i2 +i3 =i
|u| i1 +i2 +i3 =i
X X
i1 +1 i2 +1 −2
+ ∇ ψ∇ ψ +Ω ∇i1 ψ ∇i2 ψ ∇i3 ψ .
i1 +i2 =i−1 i1 +i2 +i3 =i
We then applying Proposition 6.2 for φ = ui+2 ∇i (β + 21 ∇3 φ∇φ) and Proposition 6.4 for
φ = ∇i (K − |u|1 2 − 41 ∇A φ∇A φ) and φ = ∇i σ̌ with λ1 = i + 2 and add these identities together.
We then obtain
1 1 1
kui+2 ∇i (K − − ∇A φ∇A φ, σ̌)k2L2u L2 (Su,u ) + kui+2 ∇i (β+ ∇3 φ∇φ)k2L2u L2 (Su,u )
|u|2 4 2
1 1
.k∇i (K − 2 − ∇A φ∇A φ, σ̌)k2L2u L2 (S0,u ) + ku2i+4 ∇i σ̌·ΩF1,i k2L1u L1u L1 (Su,u )
|u| 4
1 1 1
+ ku2i+4 ∇i (K − 2 − ∇A φ∇A φ)·ΩF2,i kL1u L1u L1 (Su,u ) + ku2i+4 ∇i (β+ ∇3 φ∇φ)·ΩF3,i kL1u L1u L1 (Su,u ) .
|u| 4 2
(8.16)
We proceed to control terms on the right. For the term involving F3,i we have
1
ku2i+4 ∇i (β+ ∇3 φ∇φ)·ΩF3,i kL1u L1u L1 (Su,u )
2
i+2 i 1
≤ku ∇ (β+ ∇3 φ∇φ)kL∞ 2 2
u Lu L (Su,u )
kΩui+2 F3,i kL1u L2u L2 (Su,u ) .
2
In kΩui+2 F3,i kL1u L2u L2 (Su,u ) , we note that the borderline terms are from the first two lines of
(8.15). These borderline terms obey
ku6 ∇5 (Ωtrχ)trχkL1u L2u L2 (Su,u ) .kΩ−1 u5 ∇5 (Ωtrχ)kL1u L2u L2 (Su,u )
1
≤ku6 ∇5 (Ωtrχ)kL∞ 2 2
u Lu L (Su,u )
·k kL2 L2 L2 (Su,u )
|u|Ω u u
1 1 1 1 1
3 u2 3 3 u2 a4 3 3 u2 a4
≤u aÕ(trχ) ·
2
1 ≤u a ·
2 4
1 R(β) . u a ·
2 4
1 ,
|u| 2 Ω |u| 2 Ω |u| 2 Ω
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 63
ku6 ∇5 (trχ, χ̂)(Ωtrχ, Ωχ̂, Ωω̂)kL1u L2u L2 (Su,u ) + ku6 ∇5 (Ωχ̂, Ωtrχ)χ̂kL1u L2u L2 (Su,u )
≤ku5 ∇5 (trχ, χ̂)kL∞ 2 2
u Lu L (Su,u )
· kuψ kL1u L∞
u L (Su,u )
∞
In the similar manner, we also control terms in F2,i and F1,i . In kΩui+2 F2,i kL1u L2u L2 (Su,u ) and
kΩui+2 F1,i kL1u L2u L2 (Su,u ) , the top-order derivatives appear in the next four terms. The first term
satisfies
kui+2 ∇i+1 (Ωe3 φ)∇φkL1u L2u L2 (Su,u )
≤kΩ−1 ui+1 ∇i+1 (Ωe3 φ)kL∞ 2 2
u Lu L (Su,u )
kuΩ∇φkL2u L2u L∞ (Su,u )
1 1 1 1 1 1
1 1 u2 a4 ua 2 1 1 3 3 u2 a4 u2 3 3
≤u 2 a 4 1 + 1 S· |u| 2 u 2 ≤ u 2 a 4 1 + 1 1 S ≤ u a .
2 4
For the rest terms in kΩui+2 F2,i kL1u L2u L2 (Su,u ) and kΩui+2 F1,i kL1u L2u L2 (Su,u ) , we have
X 1 1
ua 2 ua 2 O2 2 12 3 3 ua 4 O2 3 3
1
X 1 1
kΩ−1 ui+1 ∇i1 ψ i2 ∇i3 (K − − ∇A φ∇A φ, σ̌)kL1u L2u L2 (Su,u )
i1 +i2 +i3 =i−1
|u|2 4
1 1 1
ua 2 ua 2 O2 2 1 3 3 ua 4 3 3
≤Ω−1 |u|2 · 2
· 3
|u| u 2 ≤ u 2 a 4 · ≤ u2 a4 ,
|u| |u| |u|Ω
X
kui ∇i1 ψ i2 ∇i3 (Ωtrχ − (Ωtrχ)0 )kL1u L2u L2 (Su,u )
i1 +i2 +i3 =i
X 1
u2
i+ 23 i1 i2 i3 − 32 3 3 3 3
≤ku ∇ ψ ∇ (Ωtrχ − (Ωtrχ)0 )k L∞ ∞ 2
u Lu L (Su,u )
ku k L1u L2u L∞ (Su,u ) ≤u a ·
2 4
1 ≤ u2 a4 ,
i1 +i2 +i3 =i |u| 2
Gathering the above estimates, together with Proposition (8.1) we then obtain
X 1
kui+2 ∇i (K − 2 , σ̌)kL2u L2 (Su,u ) + kui+2 ∇i βkL2u L2 (Su,u )
|u|
i≤4
1 1 1 1
3 3 u2 a4 3 3 u2 a4
.u 2 a 4 1 + 1 (1 + R[β]) . u 2 a 4 1 + 1 .
|u| 2 Ω |u| 2 Ω
Combining Propositions 8.2, 8.1, 8.3, we thus conclude
R . 1. (8.18)
Substituting this bound (8.18) into Proposition 7.6, we furthermore acquire
Õ5,2 . 1. (8.19)
This proves the hyperbolic part of Theorem 1.1.
1
We further derive an additional estimate for K − |u|2 − 14 ∇A φ∇A φ.
Proposition 8.3. Under the assumptions of Theorem 1.1 and the bootstrap assumptions (2.38),
(2.39), (2.40),(2.41), it holds
1 1 3 3
ku2 (K − 2 − ∇A φ∇A φ)kL∞ 2 2
u Lu L (S)
.u 2 a 4 .
|u| 4
Proof. Applying Proposition 4.7 and the bootstrap assumption (2.40), we get
1 1 1 1 3 3
ku2 (K − 2 − ∇A φ∇A φ)kL∞ 2 2
u Lu L (S)
. u 2 · ua 2 = u 2 a 4 .
|u| 4
The above sections finish the hyperbolic estimates. Together with a standard local existence
and extension argument, we have obtained a desired spacetime existence region. We now move
to the elliptic part and construct the apparent horizon within this obtained existence regime.
Remark 20. For perturbed Christodoulou’s initial data prescribed along u = 0, we can define
the corresponding φ via letting
φ
R · Ω(−R, 0)2 = uae− 1−µ∗ .
66 XINLIANG AN
Hu
u = −R(u, ω)
Mu (ω)
ω ∈ S2 u = −1
Figure 6
Remark 21. Compared with [6] by An-Han, there µ∗ = 0 and a ≫ 1. In this paper, we have
0 < µ∗ < 1 and we can set a = 1. Two major observations in below elliptic arguments are that
i) owing to Christodoulou’s initial data, a potentially dangersome term is with a favored sign
and hence the approach to construct the multi- and single-valley anisotropic trapped surface as
in [6] and [28] still works. ii) even though the borderline quadratic nonlinear term with gradient
is no longer negligible as 1/a fails to be the small parameter, the method of deriving apriori
estimates via Moser’s iteration developed in [6] by An-Han is robust enough to be extended to
this new setting. Besides these observations, a novel ingredient in the step of linearization is a
new ansatz for the to-be-solved solution and the use of it plays a vital role.
In below, we provide the details. Note that on Su,u it holds
∆g φ = R−2 ∆γ φ.
And equation G(R, u) = 0 can be rewritten as
0 = ∆γ φ − |∇γ φ|2 − 4ΩωR|∇γ φ|2 + 2γ ij ηi ∂j φ + 2R−1 Ωχ̂kl γ ik γ jl ∇i φ∇j φ
R R (9.2)
− Ωtrg χ|∇γ φ|2 + Ω−1 trg χ.
2 2
We point out that the main contribution in Ω−1 trg χ is
2 uaf (u, ω)
· (1 − µ∗ ) − 2 2 · (1 − µ∗ )
R R Ω0 (R)
µ∗
and Ω20 (R) := Ω2 (−R, 0) = R 1−µ∗ . We then write 12 RΩ−1 trg χ as
1 uaf (u, ω) · (1 − µ∗ ) 1 uaf (u, ω) · (1 − µ∗ )
RΩ−1 trg χ = 1 − µ∗ − 2 + [ RΩ−1 trg χ − (1 − µ∗ ) + ].
2 2RΩ0 (R) 2 2RΩ20 (R)
Together with the fact
µ∗ 1 ∗ 1 φ
RΩ20 (R) = R · R 1−µ∗ = R 1−µ∗ = [(ua)1−µ e−φ ] 1−µ∗ = uae− 1−µ∗ , (9.3)
we can rewrite G(R, u) = 0 as
1 φ
0 = ∆γ φ + (1 − µ∗ ) − f (u, ω)(1 − µ∗ )e 1−µ∗ − 4ΩωR|∇γ φ|2
2
R
− (1 + Ωtrg χ)|∇γ φ|2 + 2γ ij ηi ∂j φ + 2R−1 Ωχ̂kl γ ik γ jl ∇i φ∇j φ
2
1 1 φ
+ [ RΩ−1 trg χ − (1 − µ∗ ) + f (u, ω)(1 − µ∗ )e 1−µ∗ ].
2 2
Multiplying 1/(1 − µ∗ ) on both sides of the equation and setting
φ
φe := ,
1 − µ∗
we then transfer G(R, u) = 0 into
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 67
∗ ∗ e ∗
Hence, the 2-sphere Su,u with u = −(ua)1−µ e−(1−µ )φ
= ( 23 ua)1−µ is an untrapped surface.
¯
To construct the supersolution (trapped surface), we need S(φ̃, u) < 0. Here we apply a key
property from the hyperbolic parts:
1 µ∗ 1
Ωω(u, u) ≥ Ωω(u, 0) = ∗
· > 0.13
2 8(1 − µ ) r(u, 0)
¯ it suffices to require
Thus, back to (9.4) to find the supersolution φ̃,
¯ 1 ¯ R ¯
∆γ φ̃ + 1 − f (u, ω)eφ̃ − (1 + Ωtrg χ)(1 − µ∗ )|∇γ φ̃|2
2 2
¯ ¯ ¯ R 1 ¯
+ 2γ ij ηi ∂j φ̃ + 2(1 − µ∗ )R−1 Ωχ̂kl γ ik γ jl ∇i φ̃∇j φ̃+[ Ω−1 trg χ − 1 + f (u, ω)eφ̃ ] < 0.
2(1 − µ∗ ) 2
Notice that the one-valley supersolution φ̃¯ given in Lemma 3.7 of [6] and the multi-valley superso-
¯ ¯
lution φ̃ constructed in Theorem 7.4 of [6] fulfill this requirement. Since −4(1−µ∗ )ΩωR|∇γ φ̃|2 ≤
¯ ¯
0, with these φ̃ we also have S(φ̃, u) < 0. As listed in Lemma 3.7 of [6], for any τ > 2.2, the
¯
constructed φ̃ obeys
In below, for equation (9.4) we seek the solution φe bounded between φ̃ = ln( 23 ) and φ̃¯
constructed above as in [6]. We require that φe satisfies
3
0 < ln( ) ≤ φe ≤ κ (9.13)
2
¯ Here κ = κ(m, ǫ) depends on constants m
with κ being the maximum of the supersolution φ̃.
and ǫ stated in Theorem 1.1. We further define constant K via
ē C 2 (S2 ) }
K := max{C(κ), |φ| (9.14)
with
4Λκeκ
C(κ) = (Ceκ )Ce (9.15)
being explicitly given in Theorem 10.4.
Throughout this paper, we will use parameters a, function Ω and constant κ, which satisfy
4Λκ ¯ 1 1
K = max{(Ceκ )Ce , kφ̃kC 2 } ≤ a 4 · Ω−δ̃/2 ≪ a 3 · Ω−δ̃ . (9.16)
Here, C is a uniform positive constant defined in Theorem 10.4, and Λ is introduced in (9.6).
Note that for the interior region, we have a = 1 and 0 < Ω . Ω(u, 0) ≤ Ω(u1 , 0) ≪ 1. When
κ = κ(m, ǫ) is fixed, in both regions we choose either parameter a being sufficiently large or
0 < Ω ≤ Ω(u, 0) being sufficiently small, such that (9.16) holds.
Remark 23. We can set a = 1 and the above arguments hold. Moreover, the hypersurface
∗
u = −(ua · mǫτ +2 )1−µ locates within the hyperbolic existence region.
13Note that the key property Ωω(u, 0) > 0 also holds for perturbed initial data.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 69
In this paper, F1 is of size 1 and serves as the main term. A key observation of this paper
(also pointed out one line below (4.3) in [6]) is that the regularity estimates carried out in [6]
also work for the scenario
e 2 + 1) with c = 1.
|F1 | + |F2 | ≤ c(|Dφ| (10.6)
In below, we start to conduct these elliptic estimates. We then write (10.5) in divergence
form as Z Z p h i
1 e e Dφ)
e ψ̃ dω,
aij (ω, φ̃)Di φ̃Dj ψ̃ dω = |γ| 1 − f eφ + F (ω, φ, (10.7)
2
with ψ̃ being any function in H01 (B1 ) ∩ L∞ (B1 ). We also assume φ̃ = φ̃(ω) ∈ H 1 (B1 ) and
require
0 ≤ φ̃ ≤ κ on B1 (10.8)
with κ being some positive constant, which is related to the anisotropicity via the construction
of upper and lower barriers.
Our first goal is to bound the Hölder norm of φ̃ via Moser’s iteration. For the estimates, we
pay attention to their explicit dependence on κ. Following the steps in the proof of Theorem
4.2 in [6], we obtain
Proposition 10.1. Let φ̃ ∈ H 1 (B1 ) be a positive solution to (10.7) with (10.8) satisfied. Then,
it holds
[φ̃]C α (B1/2 ) ≤ Ceκ (10.9)
with C being a positive constant depending only on Λ. And α is given by
α = ǫ0 e−4Λκ ∈ (0, 1). (10.10)
Here ǫ0 > 0 is a small constant depending on Λ.
Proof. For s ∈ (0, 1), we set
M (s) = max φ̃, m(s) = min φ̃.
Bs Bs
Writing m = m(s) and M = M (s) for brevity, proceeding the same as in the reasoning for
Theorem 4.2 in [6], we get
s h s i
M − m(s) ≤ Ce4Λκ m − m(s) + Ce4Λκ+κ s2 , (10.11)
2 2
s h s i
M (s) − m ≤ Ce4Λκ M (s) − M + Ce4Λκ+κ s2 . (10.12)
2 2
For any s ∈ (0, 1), define
h(s) := M (s) − m(s). (10.13)
Adding (10.11) and (10.12), we then obtain
s h s i
h(s) + h ≤ Ce4Λκ h(s) − h + Ce4Λκ+κ s2 .
2 2
This implies
s Ce4Λκ − 1
h ≤ Cγ h(s) + eκ s2 with Cγ = < 1.
2 Ce4Λκ + 1
Using this Cγ , we pick up µ ∈ (0, 1) satisfying
α := (1 − µ) log Cγ / log(1/2) < 2µ.
Further applying Lemma 8.23 in [24] or Lemma 4.19 in [26], for any s ∈ (0, 1/2], we hence prove
h(s) ≤ Csα h(1) + eκ s2µ ≤ Ceκ sα , (10.14)
where we employ h(1) ≤ M (1) ≤ κ ≤ eκ for the second inequality. This concludes (10.9).
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 71
Next, we derive the gradient estimate for φ̃. By the hyperbolic part, we have
|Dω akl (ω, φ̃)| + |∂φ̃ akl (ω, φ̃)| ≤ Λ1 (10.15)
with Λ1 being some positive constant. The constant in the below theorem depends on Λ1 . And
we have
Proposition 10.2. Set φ̃ ∈ H 1 (B1 ) to be a positive solution to (10.7) with (10.8) satisfied.
Then, it holds
4Λκ
|φ̃|C 1,1/3 (B1/2 ) ≤ (Ceκ )Ce (10.16)
with C being a positive constant depending only on Λ and Λ1 .
Proof. The detailed steps for proving Theorem 4.3 in [6] can also be carried out here. For any
function ψ̃ belonging to H01 (B1 ) ∩ L∞ (B1 ) with ψ̃ ≥ 0, by (10.7) and (10.6), it holds
Z Z p h i Z
1
aij Di φ̃Dj ψ̃ dω ≤ |γ| 1 − f eφ̃ + |Dφ̃|2 + 1 ψ̃ dω ≤ [4 + 2|Dφ̃|2 ]ψ̃ dω.
2
Taking ξ(τ ) = 2−1 Λ−1 (e2Λτ − 1), we have
Z Z
aij Di (ξ(φ̃))Dj ψ̃ dω ≤ 4ξ ′ (φ̃)ψ̃ dω.
This is equivalent to
Z Z
aij e2Λφ̃ Di φ̃Dj ψ̃ dω ≤ 4e2Λφ̃ ψ̃ dω.
Taking ψ̃ = φ̃µ̃2 for some µ̃ ∈ C01 (B1 ) with µ̃ = 1 on B3/4 , via using the Cauchy-Schwarz
inequality, we get Z
|Dφ̃|2 dω ≤ Cκ2 e2Λκ (10.17)
B3/4
with C being a positive constant depending only on Λ.
To derive a bound for the Hölder norm of Dφ, e we need to employ the regularity of aij (ω0 , φ̃(ω0 ))−
aij (ω, φ̃(ω)), which was already obtained in Proposition 10.1. Following the procedures of prov-
ing Theorem 4.3 in [6], for any 0 < ρ ≤ s we then obtain
Z h ρ 2 iZ
|Dφ̃|2 dω ≤ C + h(s) + [h(s)]2 + s2 |Dφ̃|2 dω + Ce2κ s4 ,
Bρ (ω0 ) s Bs (ω0 )
where
e
h(s) = h(ω0 ; s) := sup φe − inf φ.
Bs (ω0 ) Bs (ω)
Take any δ ∈ (0, 1). For a small ǫδ > 0, we choose sδ small such that Bsδ (ω0 ) ⊂ B3/4 and
h(sδ ) + [h(sδ )]2 + s2δ < ǫδ . (10.18)
κ −C2 e4Λκ
By employing (10.14) and κ ≫ 1, if we set sδ = (C1 e ) , it satisfies the above requirement.
Applying Lemma 3.4 in [26], for any s ≤ sδ we then get
Z h Z i
|Dφ̃|2 dω ≤ Cs2δ s−2δ
δ |Dφ̃|2 dω + e2κ .
Bs (ω0 ) Bsδ (ω0 )
Proof. We first fix an arbitrary ball B on S2 . In local coordinates within B, we have (10.5) and
p
aij (ω, φ̃) = |γ(u, ω)|γ ij (u, ω).
∗ ∗
Here u = −R = −(ua)1−µ e−(1−µ )φ̃ . Recall that F = F1 + F2 with F1 = −4(1 − µ∗ )RΩω|∇γ φ̃|2
and F2 being a linear combination of ∇i φ̃∇j φ̃, ∇i φ̃, 1 with the below coefficients
Thus, with the aid of (9.1), the equation of MOTS 0 = G(R, u) is changed into
1 1
0 =∆g R + Ωtrg χ|∇g R|2 + 4Ωω|∇g R|2 − Ω−1 trg χ + 2g(η, ∇g R)+2Ωχ̂(∇g R, ∇g R)
2 2
µ∗ 1
=[−λ + λ2 ]Rφ̌−2 |∇g φ̌|2 + λRφ̌−1 ∆g φ̌ − Ω−1 trg χ + 2λRφ̌−1 g(η, ∇g φ̌)
1 − µ∗ 2
1 1 µ∗ 1
+( Ωtrg χ + )λ2 R2 φ̌−2 |∇g φ̌|2 + [4Ωω − · ]λ2 R2 φ̌−2 |∇g φ̌|2 + 2Ωχ̂(∇g R, ∇g R).
2 R 1 − µ∗ R
We now pick
1 − µ∗ µ∗
λ= ∗
, which satisfies − λ + λ2 = 0.
µ 1 − µ∗
This renders the above equation to be
1 1 1
0 =λRφ̌−1 ∆g φ̌ − Ω−1 trg χ + 2λRφ̌−1 g(η, ∇g φ̌) + ( Ωtrg χ + )λ2 R2 φ̌−2 |∇g φ̌|2
2 2 R
µ∗ 1 2 2 −2 2
+[4Ωω − · ]λ R φ̌ |∇g φ̌| + 2Ωχ̂(∇g R, ∇g R).
1 − µ∗ R
74 XINLIANG AN
µ∗ 1
With Christodoulou’s naked-singularity initial data, we note that 4Ωω(−R, 0) = 1−µ∗ ·R . Hence
we further define α1 as below and recall α2 given in (10.26)
1 1
α1 := −( Ωtrg χ + )−[4Ωω − (4Ωω)0 ],14
2 R
1 −1 1 − µ∗ uaf (ω, u)
α2 := Ω trg χ − + (1 − µ∗ ).
2 R 2Ω20 R2
∗ µ∗
Together with R = (ua)1−µ φ̌λ , γ = R−2 g, ∆g φ̌ = R−2 ∆γ φ̌ and Ω20 (R) = Ω2 (u, 0) = R 1−µ∗ ,
the above equation infers
f (ω, u) − 1−µ
λ
λφ̌−1 ∆γ φ̌ − (1 − µ∗ ) + (1 − µ∗ ) φ̌ ∗
+2λR2 φ̌−1 g(η, ∇g φ̌)
2
− λ2 R3 φ̌−2 |∇g φ̌|2 α1 − Rα2 +2RΩχ̂(∇g R, ∇g R) = 0.
With λ = (1 − µ∗ )/µ∗ , we now define
µ∗ 1 2µ∗
S(φ̌, ω) :=∆γ φ̌ − µ∗ φ̌ + f (ω, u)φ̌1− µ∗ + λR2 g(η, ∇g φ̌)
2 1 − µ∗
(11.1)
µ∗ µ∗ 2µ∗
− ∗
λ2 R3 φ̌−1 |∇g φ̌|2 α1 − ∗
φ̌Rα2 + φ̌RΩχ̂(∇g R, ∇g R).
1−µ 1−µ 1 − µ∗
And for 2-sphere S 1−µ∗ , we have that, if
−(ua)1−µ∗ φ̌ µ∗ ,u
< 0, it corresponds to an untrapped surface,
S(φ̌, ω) = 0, it corresponds to a MOTS,
> 0, it corresponds to a trapped surface.
For the operator S(φ̌, ω) defined in (11.1), it further holds
Proposition 11.1. Suppose φ̌ is a given function on S2 . The linearized operator ∂φ̌ S(φ̌) can
be expressed as
1 − µ∗ 1
∂φ̌ S(φ̌)[w] = ∆γ w − µ∗ w − f (ω)φ̌− µ∗ w + ci ∇i w + cw. (11.2)
2
Here the connection is with respect to the metric γ, and functions c = c(ω, φ̌, u) and ci =
ci (ω, φ̌, u) satisfy
1
|c| + |ci | ≤ Ca− 3 Ωδ̃ (|φ̌−1 ∇2γ φ̌| + |φ̌−1 ∇γ φ̌|2 + 1).
Proof. We first write
S(φ̌) = I1 + I2 + I3 + I4 + I5 , (11.3)
with
µ∗ 1− µ1∗
I1 = ∆γ φ̌ − µ∗ φ̌ + f φ̌ , (11.4)
2
2µ∗ 2 −1 −1 ik jl 2µ∗
I2 = λ φ̌ R Ωχ̂ γ γ ∇i φ̌∇j φ̌, I3 = λγ ij ηi ∇j φ̌,
1 − µ∗ kl 1 − µ∗
µ∗ µ∗
I4 = − ∗
λ2 φ̌−1 Rα1 |∇γ φ̌|2γ , I5 = − φ̌Rα2 . (11.5)
1−µ 1 − µ∗
∗ 1−µ∗
Then we take ∂φ̌ with respect to each Ii with i = 1, 2, 3, 4, 5. Recall R = (ua)1−µ φ̌ µ∗ and we
set
∗ 1−µ∗
R(ǫ) = (ua)1−µ (φ̌ + ǫw) µ∗ .
14For perturbed initial data, in α there is an additional term, i.e., 4Ωω(u, 0) − µ∗ · 1 . Its absolute value
1 1−µ∗ |u|
1 Ω(u,0)δ̃
is bounded by |u|
· 1 and in later arguments, this additional term obeys a desired control.
a3
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 75
d
Denoting ˙ = dǫ ǫ=0 , it then holds
1 − µ∗ 1 − µ∗
Ṙ = Rφ̌−1 w, u̇ = −Rφ̌−1 w, φ̌˙ = w.
µ∗ µ∗
We also have
1 − µ∗ −1
γ̇ ij = −γ ik γ jl γ̇kl and γ̇kl = u̇∂u γkl = − φ̌ w[(RΩtrg χ + 2)γkl + 2R−1 Ωχ̂kl ].
µ∗
µ∗ 1− µ1∗
For I1 = ∆γ φ̌ − µ∗ φ̌ + 2 f φ̌ , we further apply Lemma 9.1 to γ. Together with φ̇ = w, we
obtain
1
(∆γ φ̌)· =∆γ w − γ ik γ jl γ̇kl ∇2ij φ̌ − γ ij γ kl ∇i γ̇jl ∇k φ̌ + γ ij γ kl ∇l γ̇ij ∇k φ̌ (11.6)
2
and
1 − µ∗ 1
I˙1 = ∆γ w − µ∗ w − f (ω)φ̌− µ∗ w.
2
2µ∗ ∗
For I2 = 1−µ ∗ λ2 −1 −1
φ̌ R Ωχ̂ kl
γ ik jl
γ ∇i φ̌∇j φ̌, with λ = 1−µ
µ∗ it holds
I˙2 = − 2λφ̌−2 wR−1 Ωχ̂kl γ ik γ jl ∇i φ̌∇j φ̌ − 2λφ̌−1 R−2 Rφ̌−1 λwΩχ̂kl γ ik γ jl ∇i φ̌∇j φ̌
+ 2λφ̌−1 R−1 ∂u (Ωχ̂kl ) · (−Rφ̌−1 λw) · γ ik γ jl ∇i φ̌∇j φ̌ + 2λφ̌−1 R−1 Ωχ̂kl (γ̇ ik γ jl + γ ik γ̇ jl )∇i φ̌∇j φ̌
+ 2λφ̌−1 R−1 Ωχ̂kl γ ik γ jl (∇i w∇j φ̌ + ∇i φ̌∇j w),
where the coefficient c for w and the coefficient ci for ∇i w satisfying
Ωδ̃ −1 Ωδ̃
|c| ≤ 1 |φ̌ ∇j φ̌|2 , |ci | ≤ 1 |φ̌−1 ∇φ̌|.
a3 a3
∗
2µ ij
For I3 = − 1−µ ∗γ ηi ∇j φ̌, we get
˙ 2µ∗ ik jl 1 − µ∗ −1 −1
I3 = γ γ − φ̌ w[(RΩtrg χ + 2)γkl + 2R Ωχ̂kl ] ηi ∇i φ̌
1 − µ∗ µ∗
2µ∗ ij 1 − µ∗ 2µ∗ ij
− ∗
γ ∂u ηi · (−Rφ̌−1 ∗
w)∇j φ̌ − γ ηi ∇j w.
1−µ µ 1 − µ∗
And on the right the coefficient c for w and the coefficient ci for ∇i w satisfy
Ωδ̃ Ωδ̃
|c| ≤ 1 |φ̌−1 ∇j φ̌|, |ci | ≤ 1 .
a 3 a3
∗
µ 2 −1
For I4 = − 1−µ∗ λ φ̌ Rα1 |∇γ φ̌|2 , we have the expression
15For perturbed initial data, after the cancellation of the leading terms, these worrisome items are of size
1
Ω(−R, 0)δ̃ R−2 a− 3 , which are sufficiently small.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 77
1 Ωδ̃
|∂u α2 | ≤ .
|u|2 a 13
Back to I˙5 , we conclude the coefficient c for w and the coefficient ci for ∇i w obey
Ωδ̃
|c| + |ci | ≤ 1 .
a3
This implies Z Z
w2 dω ≤ C |L0 w|2 dω.
S2 S2
We further let
Z
∗ 1 − µ∗ − µ1∗
< w1 , w2 >:= ∇γ w1 · ∇γ w2 + µ w1 w2 + f (ω, u)φ̌ w1 w2 dω.
S2 2
Since 0 < µ∗ < 1 and f ≥ 0, we can use < w1 , w2 > as an inner-product on H 1 (S2 ), equivalent
to the standard H 1 (S2 ) inner product. The first eigenvalue of −L0 can be expressed as
Z Z
1 − µ∗ 1
µ1 = inf{ [|∇γ w|2 + µ∗ w2 + f (ω, u)φ̌− µ∗ w2 ]dω; w2 dω = 1}.
S2 2 S2
Immediately, we see that µ1 ≥ µ∗ . And via using Theorem 8.38 in [24], we also have
Lemma 11.2. Assume f and φ̌ to be given nonnegative functions on S2 with (1.4) satisfied.
Then the equation L0 w = 0 only admits the trivial solution. Moreover, the first eigenvalue µ1
of −L0 is positive and simple and it has a positive eigenfunction.
With this property, we obtain the crucial invertibility of the operator ∂φ̌ S(φ̌):
Lemma 11.3. Let β ∈ (0, 1) be a fixed constant. Set f and φ̌ to be given nonnegative functions
on S2 with (1.4) satisfied and |φ̌|C 2 (S2 ) ≤ K with K obeying (9.14). Then within the hyperbolic
region, it holds that ∂φ̌ S(φ̌) : C 2,β (S2 ) → C 0,β (S2 ) is invertible.
Proof. Applying Lemma 11.2, we first take a positive eigenfunction ψ1 associated with the first
eigenvalue µ1 > 0, i.e.,
1 − µ∗ − µ1∗
∆γ ψ1 − µ∗ ψ1 − f φ̌ ψ1 = −µ1 ψ1 .
2
Using Proposition 11.1, we have
1 − µ∗ − µ1∗
∂φ̌ S(φ̌)[ψ1 ] = ∆γ ψ1 − µ∗ ψ1 − f φ̌ ψ1 + ci ∇i ψ1 + cψ1
2
= −µ1 ψ1 + ci ∇i ψ1 + cψ1 .
From derived apriori estimates and constructed lower barrier for φ̌, we have |φ̌−1 ∇2γ φ̌|+|φ̌−1 ∇γ φ̌|2 +
1 is independent of u, u and it satisfies
1
|c| + |ci | ≤ Ca− 3 Ωδ̃ (|φ̌−1 ∇2γ φ̌| + |φ̌−1 ∇γ φ̌|2 + 1).
With Ω being sufficiently small depending on u, u in the interior region, we have |ci | + |c| ≪ 1.
This implies ∂φ̌ S(φ̌)[ψ1 ] < 0 on S2 . Employing the calculation in the proof of Theorem 2.11
in [26], for any w ∈ C 2 (S2 ), we consider an elliptic equation for ψw1 . For notional simplicity, let
us rewrite the elliptic operator ∂φ̌ S(φ̌) as
∂φ̌ S(φ̌)[w] = aij Dij w + bi Di w + cw
w
with aij , bi , c being functions. We then calculate and directly check that ψ1 satisfies
Since ∂φ̌ S(φ̌)[ψ1 ] < 0 and ψ1 > 0, the above elliptic equation is invertible. For fixed ψ1 > 0,
giving ∂φ̌ S(φ̌)[w]/ψ1 , we can solve for the unique solution ψw1 . This further implies that ∂φ̌ S(φ̌) :
C 2,β (S2 ) → C 0,β (S2 ) is invertible and it holds
w ∂φ̌ S(φ̌)[w]
| |C 2,β (S2 ) ≤ C| |C 0,β (S2 ) and |w|C 2,β (S2 ) ≤ C|∂φ̌ S(φ̌)[w]|C 0,β (S2 ) (11.10)
ψ1 ψ1
with some positive constant C independent of w and depending on function φ̌ and its associated
eigenfunction ψ1 .
In the next lemma, we eliminate such dependence on φ̌.
Lemma 11.4. Assume f and φ̌ to be given nonnegative functions on S2 with (1.4) satisfied and
|φ̌|C 2 (S2 ) ≤ K with K obeying (9.14). Then, within the hyperbolic region, for any w ∈ C 2 (S2 ),
we have that
max
2
|w| ≤ C max
2
|∂φ̌ S(φ̌)[w]|. (11.11)
S S
Here the positive constant C depends only on m and ǫ in (1.4) and K, independent of φ̌ and
u ∈ (0, δ].
Proof. We take a smooth nonnegative function f0 = f0 (ω) on S2 satisfying f0 ≥ m/2 on Bp (ǫ)
and f (·, u) ≥ f0 on S2 for all u ∈ (0, δ]. With the standard spherical metric γ0 on S2 , we first
consider
1 − µ∗
L∗ w = ∆γ0 w−µ∗ w − f0 w.
2
Employing Lemma 11.2 to L∗ , we can find a positive eigenfunction ψ1 associated with its first
eigenvalue µ1 > 0, i.e.,
1 − µ∗
∆γ0 ψ1 −µ∗ ψ1 − f ψ1 = −µ1 ψ1 .
2
Here µ1 and ψ1 depend only on m, ǫ and are independent of u ∈ (0, δ]. With the fact 0 < φ̌ ≤ 1,
f ≥ f0 , applying Lemma 9.2, and Proposition 11.1, we then get
1 − µ∗ − µ1∗
∂φ̌ S(φ̌)[ψ1 ] = ∆γ0 ψ1 − µ∗ ψ1 − f φ̌ ψ1 + ∆γ ψ1 − ∆γ0 ψ1 + ci ∇i ψ1 + cψ1
2
1 − µ∗
≤ ∆γ0 ψ1 − µ∗ ψ1 − f0 ψ1 + ∆γ ψ1 − ∆γ0 ψ1 + ci ∇i ψ1 + cψ1
2
= −µ1 ψ1 + ∆γ ψ1 − ∆γ0 ψ1 + ci ∇i ψ1 + cψ1 ≤ −c(m, ǫ, K) < 0.
For the last inequality, we choose a being suitably large in the exterior and Ω being sufficiently
small in the interior. We then repeat the proof of Lemma 11.3 with this ψ1 , which is independent
of φ̌. And we prove the current lemma.
In below we further prove that the equation S(φ̌) = 0 admits a unique solution φ̌ satisfying
|φ̌|C 2 (S2 ) ≤ K.
Proposition 11.5. Set f and φ̌1 , φ̌2 to be given nonnegative functions on S2 with (1.4) satisfied.
Require |φ̌1 |C 2 (S2 ) ≤ K and |φ̌2 |C 2 (S2 ) ≤ K with K obeying (9.14). If S(φ̌1 ) = S(φ̌2 ), then it
holds that φ̌1 ≡ φ̌2 for all ω ∈ S2 .
Proof. For the given φ̌1 and φ̌2 , we define the below L and view it as a linear operator of φ̌1 − φ̌2 .
Z 1
L(φ̌1 − φ̌2 ) := S(φ̌1 ) − S(φ̌2 ) = ∂φ̌ S(tφ̌1 + (1 − t)φ̌2 )[φ̌1 − φ̌2 ]dt.
0
The same reasoning as for Lemma 5.5 of [6] yields the desired result.
R1
Here operator 0 ∂φ̌ S(tφ̌1 +(1−t)φ̌2 )dt possesses a similar structure to the linearized operator
∂φ̌ S(φ̌) in Proposition 11.1. The corresponding Lemmas 11.2-11.4 can be obtained in the same
manner. By the assumption, we first use S(φ̌1 ) ≥ S(φ̌2 ) and it gives L(φ̌1 − φ̌2 ) ≥ 0. With
the same proof as in Lemma 11.3, for the operator L, applying Theorem 2.11 in [26], the strong
maximum principle holds here. Hence it holds either φ̌1 < φ̌2 or φ̌1 ≡ φ̌2 on S2 . We then
80 XINLIANG AN
use S(φ̌2 ) ≥ S(φ̌1 ). Repeating the argument, we now have either φ̌2 < φ̌1 or φ̌1 ≡ φ̌2 on S2 .
Combining these together, we have proved the desired uniqueness result.
As a corollary of the above lemma, we also prove that
Theorem 11.6. Take u to be a fixed constant in (0, δ]. Let f be a given smooth function on
S2 , satisfying (1.3) and (1.4), and φ̌ and φ̌ be constructed as in Section 9. Then, there exists a
unique smooth solution φ̌ = φ̌(ω, u) of S(φ̌, u) = 0, satisfying φ̌ < φ̌(·, u) < φ̌ on S2 .
Proof. Here, for any constant λ ∈ [0, 1] and any C 2 -function φ̌ on S2 , we set
µ∗ 1
Sλ (φ̌) = ∆γ φ̌ − µ∗ φ̌ + λf + (1 − λ)]φ̌1− µ∗ + λ(I2 + I3 + I4 + I5 ),
2
where Ii , with 2 ≤ i ≤ 5, are defined in (11.5). With the same line of reasoning as for the detailed
proof of Theorem 6.1 in [6], by applying the method of continuity for parameter λ ∈ [0, 1],
together with the utilization of the invertibility of ∂φ̌ Sλ (φ̌), we conclude the existence. The
uniqueness follows from Proposition 11.5.
The desired estimates for ∂u α1 and ∂u α2 were already obtained in Proposition 11.1. To bound
∂u α1 and ∂u α2 , we utilize the explicit formulas for ∂u [ 12 Ωtrg χ + R1 ] and ∂u (Ωω), the estimates
for ∂u α1 can be obtained straightforwardly. And for ∂u α2 , we have
1 1
∂u α2 (u, u, ω) = − × ( trχtrχ + |χ̂|2g + ∂4 φ∂4 φ)(u, u, ω)
2 2
1 Ω−2 (u, u, ω) 2 2
+ × |Ω χ̂| g (−1, u, ω) + [Ωe 4 φ(−1, u, ω) − Ωe 4 φ(−1, 0, ω)] .
2 |u|2
The desired estimate for above expression can be got in the same fashion as for (11.8).
Proceeding similarly to the argument for Corollary 6.4 in [6], we then obtain
Corollary 12.3. For each u ∈ (0, δ], let φ̌ = φ̌(ω, u) be the solution to S(φ̌, u) = 0. Under the
assumptions in Theorem 11.6, it holds
1 1
max
2
|∂u φ̌| ≤ C u−1 a− 3 Ωδ̃ max
2
(|φ̌| + |∇γ φ̌| + |∇2ij φ̌| + |φ̌−1 ∇γ φ̌∇γ φ̌|) + max
2
|∂u f φ̌1− µ∗ | ,
S S S
(12.3)
where C is a positive constant independent of a, u and φ̌.
This further implies
Proposition 12.4. For each u ∈ (0, δ], set φ̌ = φ̌(ω, u) to be the solution of S(φ̌, u) = 0. Under
the assumptions in Theorem 11.6, if we further require
1
|u∂u f (ω, u)| ≤ a− 3 Ωδ̃ (u1 ), (12.4)
∗
the apparent horizon {(ω, u, u)| u = −(ua)1−µ φ̌λ (ω, u)} is spacelike.
∗ 1−µ∗
Proof. Using R = (ua)1−µ φ̌λ (ω, u) and λ = µ∗ , we have ∂u R = u−1 R[(1 − µ∗ ) + λuφ̌−1 ∂u φ̌].
∗
− 31 δ̃
Applying Corollary 12.3, if |u∂u f (ω, u)| ≤ a Ω , then we get ∂u R ≥ 1−µ
2 u
−1
R. Assume g ′ is
2
the induced metric to the apparent horizon u = −R(ω, u). Rewrite ω ∈ S as ω = (θ1 , θ2 ). As
shown in [3], then it holds
gθ′ i θj = gθi θj , gθ′ i u = 2∂θi R, gu′ u = 4∂u R.
Suppose X = λ1 ∂θ1 + λ2 ∂θ2 + λ3 ∂u to be an arbitrary nonzero tangent vector along the
apparent horizon with constants λ1 , λ2 , λ3 . It holds
g ′ (X, X) = λi λj gij + 4λi λ3 ∂i R + 4λ23 ∂u R,
where i, j = 1, 2 are summed over. The case λ3 = 0 implies g ′ (X, X) = λi λj gij > 0. For the
case λ3 6= 0, using gij = R2 γij and ∂i R = λRφ̌−1 ∇i φ̌, we get
g ′ (X, X) =λi λj gij + 4λi λ3 ∂i R + 4λ23 ∂u R
1 − µ∗ −1
≥R2 λi λj γij + 4λi λ3 λRφ̌−1 ∇i φ̌ + 4λ23 · u R
2 (12.5)
≥2λ23 (1 − µ∗ )u−1 R − 4λ23 λ2 (φ̌−1 ∇i φ̌)2
∗ ∗
>2λ23 (1 − µ∗ )u−µ (amǫτ +2 )1−µ − 4λ3 λ2 (φ̌−1 ∇i φ̌)2 .
For the last inequality, we use φēmax = − log(mǫτ +2 ) + C(τ ) in (9.12) and hence
∗
R ≥ (ua · mǫτ +2 )1−µ with τ > 2.2 fixed.
∗
At the same time, by (9.11) we have u = −( 23 ua)1−µ is an untrapped surface. For the con-
∗ 1−µ∗ 1−µ∗ ∗
structed MOTS locating at u = −R = −(ua)1−µ φ̌ , it thus holds φ̌ µ∗ ≥ ( 23 )1−µ . That is
µ∗
2 µ∗
φ̌ ≥ ( 3 ) . Therefore, back to (12.5), when u is sufficiently small, we have g ′ (X, X) > 0. This
renders our constructed apparent horizon spacelike and hence a dynamical horizon.
82 XINLIANG AN
13.1. Instability Subject to Perturbations with Finite BV Norms. We first recall equa-
tion (5.1)
1 1
(Ωe3 )(Ωe4 φ) + ΩtrχΩe4 φ = Ω2 ∆g φ − ΩtrχΩe3 φ+2Ω2 η A eA φ.
2 2
With e3 = Ω−1 ∂/∂u, this equation is equivalent to
∂ 1 2 1
(uΩe4 φ) + (Ωtrχ + )uΩe4 φ = uΩ2 ∆g φ − uΩtrχ · Ωe3 φ + 2uΩ2 η A eA φ.
∂u 2 |u| 2
Via integration, we hence obtain
Z u
1
uΩe4 φ(u, u, ω) =uΩe4 φ(−1, u, ω) + − u′ Ωtrχ · Ωe3 φ(u′ , u, ω)du′
−1 2
Z u (13.1)
1 2
+ [− (Ωtrχ + ) · uΩe4 φ + uΩ2 ∆g φ + 2uΩ2 η A eA φ](u′ , u, ω)du′ .
−1 2 |u|
For the first two terms on the right of (13.1), we can rewrite them as
Z u
1
uΩe4 φ(−1, u, ω) + − u′ Ωtrχ · Ωe3 φ(u′ , u, ω)du′
−1 2
Z
1 u ′
=uΩe4 φ(−1, u, ω) − uΩe4 φ(−1, 0, ω) + uΩe4 φ(−1, 0, ω) − u Ωtrχ · [Ωe3 φ − (Ωe3 φ)0 ](u′ , u, ω)du′
2 −1
Z Z
1 u ′ 1 u ′
− u [Ωtrχ − (Ωtrχ)0 ] · (Ωe3 φ)0 (u′ , u, ω)du′ − u (Ωtrχ)0 · (Ωe3 φ)0 (u′ , u, ω)du′ .
2 −1 2 −1
Employing (13.1), with Christodoulou’s initial data along u = 0 we have
Z u
1
uΩe4 φ(u, 0, ω) = uΩe4 φ(−1, 0, ω) + − u′ (Ωtrχ)0 · (Ωe3 φ)0 (u′ , 0, ω)du′ .
−1 2
Hence, the underlined two terms can be replaced by uΩe4 φ(u, 0, ω). Back to (13.1) again, we
now obtain
uΩe4 φ(u, u, ω) − [uΩe4 φ(−1, u, ω) − uΩe4 φ(−1, 0, ω)]
Z
1 u ′
=uΩe4 φ(u, 0, ω) − u Ωtrχ[Ωe3 φ − (Ωe3 φ)0 ](u′ , u, ω)du′
2 −1
Z
1 u ′
− u [Ωtrχ − (Ωtrχ)0 ] · (Ωe3 φ)0 (u′ , u, ω)du′
2 −1
Z u
1 2
+ [− (Ωtrχ + ) · uΩe4 φ + uΩ2 ∆g φ + 2uΩ2 η A eA φ](u′ , u, ω)du′ .
−1 2 |u|
Utilizing the fact
r
√ 1 − µ∗ 1 − µ∗ µ∗
uΩe4 φ(u, 0, ω) = − 2 · · |u| 1−µ∗
µ∗ 4
and the hyperbolic estimates derived in previous sections, we get
ua µ∗
|uΩe4 φ(u, u, ω) − [Ωe4 φ(−1, u, ω) − Ωe4 φ(−1, 0, ω)]| . + |u| 1−µ∗ . (13.2)
|u|
Now we consider the key equation that drives the trapped surface formation
1
Ω∇4 (Ω−1 trχ − (Ω−1 trχ)0 ) = −|χ̂|2 − (∇4 φ)2 − (trχ)2 .
2
17For perturbed initial data, analogous instability theorems can be proved correspondingly.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 83
then for each B there exist a sufficient small δ = δ(B) and constant a = a(B)18, and the
Einstein-scalar field system admits a unique regular solution in the spacetime region (u, u, ω),
where 0 ≤ u ≤ δ and u ≤ |u|Ω2−δ̃ (u, 0) ≤ 1 with 0 < δ̃ ≪ 1.
Note that for u < 0 we require our spacetime to coincide with Christodoulou’s naked-
singularity solution in [18]. Along u = −1 we consider the possible extensions of χ̂(−1, u, ω), ∂u φ(−1, u, ω)
for u > 0. Via the next theorem, we can see that Christodoulou’s naked-singularity solution is
associated with nonlinear instability of high co-dimensions:
Theorem 13.2. Consider the space of prescribed initial data for χ̂(−1, u, ω) and ∂u φ(−1, u, ω)
√ q 1−µ∗
along u = −1 with u > 0. With |χ̂|(−1, 0, ω) = 0, ∂u φ(−1, 0, ω) = 2 · 1−µ
∗
µ∗ · 4 cor-
responding to their values of Christodoulou’s naked-singularity solution in [18], then for any
k ∈ Z+ , we can find smooth functions f1 , f2 , ..., fk and g1 , g2 , ..., gk with variables u and ω, such
that the below prescribed initial data for |χ̂| and ∂u φ with u > 0
|χ̂|(−1, u, ω), ∂u φ(−1, u, ω)
= 0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ′1 g1 + λ′2 g2 + · · · + λ′k gk
would lead to the trapped surface and the apparent horizon formation, when
k
X k
X 2
λ2i + λ′i 6= 0 with λk , λ′k ∈ R.
1 1
Moreover, if it holds
0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ′1 g1 + λ′2 g2 + · · · + λ′k gk
= 0 + λ̃1 f1 + λ̃2 f2 + · · · + λ̃k fk , ∂u φ(−1, 0, ω) + λ̃′1 g1 + λ̃′2 g2 + · · · + λ̃′k gk ,
then it renders
λi = λ̃i and λ′i = λ̃i with 1 ≤ i ≤ k.
+
We may therefore say that, for any k ∈ Z , Christodoulou’s naked-singularity solution in
[18] has at least co-dimensional 2k nonlinear instability subject to outgoing BV characteristic
perturbations at (−1, 0, ω).
Proof. For any k ∈ Z+ and 0 ≤ u ≪ 1, on the 2-sphere S−1,u we first fix 2k disjoint discs:
D1 , D2 , ..., Dk and D1′ , D2′ , ..., Dk′ . We prescribe smooth functions fi being compactly supported
in Di and smooth functions gi being compactly supported in Di′ . For 1 ≤ i ≤ k, we further
require
1 1
0 ≤ fi ≤ a 2 on S2 × (0, δ], fi (·, u) ≥ mi a 2 on Bpi (ǫi ) ⊂ Di ,
1 1
0 ≤ gi ≤ a 2 on S2 × (0, δ], gi (·, u) ≥ m′i a 2 on Bp′i (ǫ′i ) ⊂ Di′
with points pi ∈ Di ⊂ S2 , p′i ∈ Di′ ⊂ S2 and constants mi , m′i ∈ (0, 1), ǫi , ǫ′i ∈ (0, π/2). With
these choices of initial data, conditions (13.6) and (13.7) in Theorem 13.1 are satisfied and the
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 85
corresponding conclusions for trapped surface and apparent horizon formation hold. And if we
have
′ ′ ′
0 + λ1 f1 + λ2 f2 + · · · + λk fk , ∂u φ(−1, 0, ω) + λ1 g1 + λ2 g2 + · · · + λk gk
′ ′ ′
= 0 + λ̃1 f1 + λ̃2 f2 + · · · + λ̃k fk , ∂u φ(−1, 0, ω) + λ̃1 g1 + λ̃2 g2 + · · · + λ̃k gk ,
since all the fi and gi with 1 ≤ i ≤ k are compacted supported in disjoint discs, it immediately
implies
λi = λ̃i and λ′i = λ̃i with 1 ≤ i ≤ k.
13.2. Instability Subject to Perturbations with Finite C 0 Norms. The above two theo-
rems in Section 13.1 generalize both the results in [2] by An-Luk and in [6] by An-Han with the
requirement:
1 1
fi (ω, u) ≥ mi a 2 on Bpi (ǫi ) ⊂ Di , gi (ω, u) ≥ m′i a 2 on Bp′i (ǫ′i ) ⊂ Di′ for u > 0,
and
fi (ω, u) = 0, gi (ω, u) = 0 for u ≤ 0.
In this subsection, we extend the instability theorems in Section 13.1 further by allowing choosing
continuous functions fi (ω, u) and gi (ω, u) for the u variable.
We start from reconsidering
2(1 − µ∗ )
Ω−1 trχ(u, u, ω) −
|u|
Z u Z u Z u
1
= −|χ̂|2 (u, u′ , ω)du′ + −|∇4 φ|2 (u, u′ , ω)du′ + − (trχ)2 (u, u′ , ω)du′ .
0 0 0 2
By the estimates derived above, we have
Z u
−1 2(1 − µ∗ ) |Ωχ̂|2 (−1, u′ , ω) + [Ωe4 φ(−1, u′ , ω) − Ωe4 φ(−1, 0, ω)]2 ′
|Ω trχ(u, u, ω) − − du |
|u| 0 |u|2 Ω2 (u, 0, ω)
1 1
ua ua 2 ua ua 2 µ∗ 1 u ua u2 a
≤ 2 2
· + 2 2 ·( + |u| 1−µ∗ a− 2 ) + 2 2 Ω4 (u, u, ω) + 2 2 · 2 (u, u, ω).
|u| Ω |u| |u| Ω |u| |u| Ω |u| Ω |u|
µ∗
Together with Ω(u, 0) = |u| 2(1−µ∗ ) , this implies
Z u
2(1 − µ∗ ) |Ωχ̂|2 (−1, u′ , ω) + [Ωe4 φ(−1, u′ , ω) − Ωe4 φ(−1, 0, ω)]2 ′
|Ω−1 trχ(u, u, ω) − − du |
|u| 0 |u|2 Ω2 (u, 0, ω)
1 1 5
ua 3 ua 2 Ω2 Ω2 ua 3
≤ 2 2
· Ω2 · [ 3 + 1 + ]≪ 2 2
· Ω2 .
|u| Ω |u|Ω 2 a2 a |u| Ω
For the f (ω, u) introduced in (13.6), with any constant C1 > 2 and Λ defined in Proposition
3.2, we further set
1 1 1
f (ω, u) = g(u)f˜(ω, u) with g(u) = (C1 Λ + C1 ) 2 · [ln(ln )]− 2 (13.8)
u
and f˜(ω, u) to satisfy
0 ≤ f˜ ≤ 1 on S2 × (0, δ] and f˜(·, u)≥m on Bp (ǫ) (13.9)
2
for a point p ∈ S and constants m ∈ (0, 1), ǫ ∈ (0, π/2). With this choice of g(u), we let
3 3 1 1 1
κ := ln[ g(u)−1 ] = ln {ln[(ln ) C1 Λ+C1 ]} 2 .
2 2 u
1 1
Now we have eκ = 3{ln[(ln 1/u) C1 Λ+C1 ]} 2 /2. And when 0 ≤ u ≪ 1, it holds
1 1
8Λκeκ ≤ e2κ = ln[(ln ) C1 Λ+C1 ] ≤ ln[l ln(u−1 )] = ln ln(u−l ) with 0 < l ≪ 1.
u
86 XINLIANG AN
with τ being a fixed number larger than 2.2. Back to (9.16), if requiring (13.10), we then verify
∗
4Λκ ¯ 1 1 −δ̃µ
K = max{(Ceκ )Ce , kφ̃kC 2 } ≤ u−l ≤ a 3 · Ω(u, 0)−δ̃ = a 3 · |u| 2(1−µ∗ ) . (13.13)
Now we revisit the arguments in Section 10, Section 11 and Section 12 of this article with
f (ω, u) replaced by g(u)f˜(ω, u). The arguments in Section 10 stay the same. For the arguments
in Lemma 11.3 of Section 11, we now consider the positive eigenfunction ψ1 associated with the
corresponding first eigenvalue µ1 > 0, i.e.,
1 − µ∗ − µ1∗
∂φ̌ S(φ̌)[ψ1 ] = ∆γ ψ1 − µ∗ ψ1 − f φ̌ ψ1 + ci ∇i ψ1 + cψ1
2
= −µ1 ψ1 + ci ∇i ψ1 + cψ1 .
With Ω being sufficiently small in the interior region, we have |ci | + |c| ≪ 1. To make sure that
1
∂φ̌ S(φ̌)[ψ1 ] < 0 on S2 , here we check the size of f φ̌− µ∗ . Using the facts f (ω, u) = g(u)f˜(ω, u)
1−µ∗
1 1
with 0 ≤ f˜ ≤ 1 and R = (ua)1−µ φ̌ µ∗ being the MOTS, we have uaφ̌ µ∗ = R 1−µ∗ ≈ uag(u),
∗
which is from the construction of the barriers. Here ≈ is independent of u. Hence it holds
1 1
f φ̌− µ∗ = g(u)f˜(u, ω)φ̌− µ∗ ≈ 1, being independent of u. Thus, the rest arguments stay the
same as in Section 11. We proceed similarly as in Lemma 11.4. With replacing f and f0 by
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 87
g(u)f˜ and g(u)f˜0 , we get the corresponding desired conclusion. All the other proofs for the
elliptic part can be obtained in the same way as in Section 11 and Section 12. In summary,
allowing the continuous perturbation for χ̂(−1, u, ω) and ∇4 φ(−1, u, ω) with respect to the u
variable, proceeding the same as in Section 13.1 we then prove
Theorem 13.3. Consider the characteristic initial value problem for the Einstein-scalar field
system (1.1). Assigning Christodoulou’s naked-singularity initial data in [18] along u = 0 with
−1 ≤ u ≤ 0 and prescribing perturbed initial data along u = −1 satisfying
X 1 1
uj a− 2 k(∂u )j ∇i χ̂kL∞ 2
u L (S−1,u )
+ uj a− 2 k(∂u )j ∇i ∇4 φkL∞ 2
u L (S−1,u )
≤ B, (13.14)
i≤5,j≤3
then for each B there exist a sufficient small δ = δ(B) and constant a = a(B)19, and the
Einstein-scalar field system admits a unique regular solution in the spacetime region (u, u, ω)
δ̃
satisfying 0 ≤ u ≤ δ and u ≤ |u|Ω2− 2 (u, 0) ≤ 1 with 0 < δ̃ ≪ 1.
Moreover, if we further assume
Z u
|χ̂|2 (−1, u′ , ω) + [∇4 φ(−1, u′ , ω) − ∇4 φ(−1, 0, ω)]2 du′ = g(u)f˜(ω, u)ua, (13.15)
0
with
1 1 1
g(u) = (C1 Λ + C1 ) 2 · [ln(ln )]− 2
u
where C1 > 2, Λ being a fixed positive constant, and f˜(ω, u) obeying
0 ≤ f˜ ≤ 1 on S2 × (0, δ] and f˜(·, u) ≥ m on Bp (ǫ) (13.16)
for a point p ∈ S2 and constants m ∈ (0, 1), ǫ ∈ (0, π/2), then within the solved hyperbolic
region, there exists a unique MOTS Mu on each H u with 0 < u ≤ δ and the collection of MOTS
{Mu } emerges and censors the central singularity. And they form an achronal apparent horizon.
We also have the corresponding nonlinear instability result with high co-dimensions. It is
summarized in Theorem 1.3. Here we prove it.
Proof. We proceed in the same way as in Theorem 13.2. Only we modify the requirement for
fi and gi as: For
1 1 1
1 ≤ i ≤ k, C1 > 2 and g(u) = (C1 Λ + C1 ) 2 · [ln(ln )]− 2 ,
u
the prescribed fi and gi obey
1 1 1 1
0 ≤ fi ≤ g(u) 2 a 2 on S2 × (0, δ], fi (·, u) ≥ mi g(u) 2 a 2 on Bpi (ǫi ) ⊂ Di ,
1 1 1 1
0 ≤ gi ≤ g(u) 2 a 2 on S2 × (0, δ], gi (·, u) ≥ m′i g(u) 2 a 2 on Bp′i (ǫ′i ) ⊂ Di′
with points pi ∈ Di ⊂ S2 , p′i ∈ Di′ ⊂ S2 and constants mi , m′i ∈ (0, 1), ǫi , ǫ′i ∈ (0, π/2). Note
that we can set a = 1. And as u → 0+ , we have
1
lim g(u) 2 = 0.
u→0+
Proof. Let fˆ(u) := [ln(ln u1 )]− 4 . Using the Sobolev-Slobodeckij space, we have
1
ZZ ZZ
|fˆ(x) − fˆ(y)|2 |fˆ(x) − fˆ(y)|2
kfˆk2 1 = 2
dxdy = 2 dxdy
Ḣ ([0,δ])
2 0≤x≤δ, |x − y| 0≤x≤y≤δ |x − y|2
0≤y≤δ
ZZ ZZ
|fˆ(x) − fˆ(y)|2 |fˆ(x) − fˆ(y)|2
=2 dxdy + 2 dxdy
0≤x≤y≤δ, |x − y|2 0≤x≤y≤δ, |x − y|2
x≤y≤2x y>2x
=Iˆ1 + Iˆ2 .
For Iˆ1 , applying the intermediate value theorem, with θ being a constant ∈ [0, 1], we get
ZZ 2
1 1 5 1 1
Iˆ =2 [ln(ln )]− 4 · 1 · dxdy
0≤x≤y≤δ, 4 x + θ(y − x) ln x+θ(y−x) x + θ(y − x)
x≤y≤2x
ZZ Z δ
1 1 1 1 1
.2 · dxdy . 1 2 · x dx . 1.
0≤x≤y≤δ, (ln 1 )2 y 2 0 (ln ) ln
x≤y≤2x x x δ
To bound Iˆ2 , we set x := eu and y := ev . The requirement 0 < 2x < y < δ implies
1 1
−∞ < u ≤ − ln − ln 2, −∞ < v ≤ − ln and v − u > ln 2 > 0.
δ δ
Utilizing (u, v) and (u, w) coordinates with w = v − u, we then obtain
ZZ 1
|[ln(ln y1 )]− 4 − [ln(ln x1 )]− 4 |2
1
Iˆ2 =2 dxdy
0≤x≤y≤δ, |y − x|2
y>2x
ZZ
1 1
.2 1
e−2v eu ev |[ln(−v)]− 4 − [ln(−u)]− 4 |2 dudv
−∞<u,v≤− ln δ,
v−u≥ln 2
Z +∞ Z − ln 1
δ −w 1 1
.2 e−w |[ln(−u − w)]− 4 − [ln(−u)]− 4 |2 dudw.
ln 2 −∞
1 1
of the above perturbation in (13.17) to be of size (ln 1δ )− 2 a 2 ≪ 1. And as δ → 0+ , we allow the
size of the scale-critical perturbation shrink to 0.
NAKED SINGULARITY CENSORING WITH ANISOTROPIC APPARENT HORIZON 89
∂x
×
(u, θ1 , θ2 )
∂u
From Chapter 2 in [21], we have
1 ∂ĝ 1 ∂
χ̂AB = φ2 = φ2 W 2 (θ1 , θ2 ) mAB (−1, u, θ1 , θ2 ) and
2 ∂u 2 ∂u
1 2 1 −1 CA ∂m ∂m
|χ̂| = (m ) ( )AB (m−1 )BD ( )DC .
2 g 8 ∂u ∂u
For our concrete construction, along u = −1 it hence holds
2 2
1 2 1 √ ∂x 1 1 0 1 ∂x
|χ̂| = · 2· (u, θ1 , θ2 ) · a 2 · tr = a· (u, θ1 , θ2 ) . (1.2)
2 g 8 ∂u 0 1 2 ∂u
With ∂x(u, θ1 , θ2 )/∂u as the freely prescribed scalar function on S−1,u , we therefore obtain the
desired localized-in-angle initial data for |χ̂|2g (−1, u, ω).
References
[1] X. An, Formation of trapped surfaces from past null infinity, preprint (2012), arXiv:1207.5271.
[2] X. An, J. Luk. Trapped surfaces in vacuum arising dynamically from mild incoming radiation, Adv. Theor.
Math. Phys. Volume 21, Number 1, 1-120 (2017).
[3] X. An, Emergence of apparent horizon in gravitational collapse, Ann. PDE 6, 10 (2020), 89 pages.
[4] X. An, A scale-critical trapped surface formation criterion: a new proof via signature for decay rates, Ann.
PDE 8, 3 (2022), 89 pages.
[5] X. An, N. Athanasiou, A scale-critical trapped surface formation criterion for the Einstein-Maxwell system,
arXiv: 2005.12699, preprint 105 pages, to appear in J. Math. Pures Appl.
[6] X. An, Q. Han, Anisotropic dynamical horizons arising in gravitational collapse, preprint (2020),
arXiv:2010:12524.
[7] X. An, T. He, Dynamics of apparent horizon and a null comparison principle, in preparation.
20In the below, we will prescribe Ψ(u, θ , θ ) in the north polar chart. And Ψ′ (u, θ ′ , θ ′ ) in the south polar
1 2 1 2
chart can be obtained accordingly.
90 XINLIANG AN
[8] X. An, Z. Lim, Trapped surface formation for spherically symmetric Einstein-Maxwell-charged scalar field
system with double null foliation, 34 pages, to appear in Ann. Henri Poincaré.
[9] X. An, H. Tan, in preparation.
[10] X. An, X. Zhang, Examples of naked singularity formation in higher-dimensional Einstein-vacuum space-
times, Ann. Henri Poincaré 19, 619-651 (2018).
[11] L. Andersson, M. Eichmair, J. Metzger, Jang’s equation and its applications to marginally trapped surfaces,
Complex analysis and dynamical systems IV. Part 2, 13-45, Contemp. Math., 554, Israel Math. Conf. Proc.,
Amer. Math. Soc., Providence, RI, 2011.
[12] L. Andersson, M. Mars, W. Simon, Local existence of dynamical and trapping horizons, Phys. Rev. Lett. 95
(2005), 111102.
[13] L. Andersson, M. Mars, W. Simon, Stability of marginally outer trapped surfaces and existence of marginally
outer trapped tubes, Adv. Theor. Math. Phys. 12 (2008), 853-888.
[14] L. Andersson, J. Metzger, The area of horizons and the trapped region, Comm. Math. Phys. 290 (2009),
no.3, 941-972.
[15] N. Athanasiou, P. Mondal, S. T. Yau, Formation of trapped surface in the Einstein-Yang-Mills system,
preprint (2023), arXiv: 2302.06915.
[16] D. Christodoulou, The formation of black holes and singularities in spherically symmetric gravitational
collapse, Comm. Pure Appl. Math. 44 (1991), no. 3, 339–373.
[17] D. Christodoulou, Bounded variation solutions of the spherically symmetric Einstein-scalar field equations,
Comm. Pure Appl. Math. 46 (1993), no. 8, 1131–1220.
[18] D. Christodoulou, Examples of naked singularity formation in the gravitational collapse of a scalar field,
Ann. of Math. (2) 140 (1994), no. 3, 607–653.
[19] D. Christodoulou, The instability of naked singularities in the gravitational collapse of a scalar field, Ann.
of Math. (2) 149 (1999), no. 1, 183–217.
[20] D. Christodoulou, On the global initial value problem and the issue of singularities, Class. Quantum Grav.
16 (1999), A23–A35.
[21] D. Christodoulou, The formation of black holes in general relativity, Monographs in Mathematics, European
Mathematical Soc. (2009).
[22] D. Christodoulou, S. Klainerman, The global nonlinear stability of the Minkowski space, Princeton mathe-
matical series 41 (1993).
[23] J. L. Costa, The formation of trapped surfaces in the gravitational collapse of spherically symmetric scalar
fields with a positive cosmological constant, Class. Quantum Grav. 37 (2020), 195022 (16pp).
[24] D. Gilbarg, N. S. Trudinger, Elliptic Partial Differential Equations of Second Order Grundlehren, Vol. 224,
Springer-Verlag, Berlin, (1983).
[25] Y. Guo, M. Hadzic, J. Jang, Naked singularities in the Einstein-Euler system, arXiv:2112.10826.
[26] Q. Han, F-H. Lin, Elliptic partial differential equations (2nd ed.), Courant lecture notes, Vol. 1, (2011).
[27] S. Klainerman, F. Nicolo, The evolution problem in General Relativity, Progress in Mathematical Physics,
Birkhaüser (2003).
[28] S. Klainerman, J. Luk, I. Rodnianski, A fully anisotropic mechanism for formation of trapped surfaces in
vacuum, Invent. math. 198, 1–26 (2014).
[29] S. Klainerman, I. Rodnianski, On the the formation of trapped surfaces, Acta Math. 208 (2012), no.2,
211–333.
[30] P. Le, The intersection of a hyperplane with a lightcone in the Minkowski spacetime, J. Differential Geom,
Vol. 109 (2018), no.3, 497-507.
[31] J. Li, J. Liu, A robust proof of the instability of naked singularities of a scalar field in spherical symmetry,
Commun. Math. Phys, Vol. 363, 561–578 (2018).
[32] J. Li, J. Liu, Instability of spherical naked singularities of a scalar field under gravitational perturbations,
to appear in J. Differential Geom.
[33] R. Penrose, Gravitational collapse and space-time singularities, Phys. Rev. Lett. 14 (1965), 57–59.
[34] R. Penrose, Gravitational collapse: The role of general relativity, Riv. Nuovo Cim., 1:252–276, 1969.
[35] I. Rodnianski, Y. Shlapentokh-Rothman, Naked Singularities for the Einstein Vacuum Equations: The
Exterior Solution, to appear in Ann. of Math.
[36] Y. Shlapentokh-Rothman, Naked Singularities for the Einstein Vacuum Equations: The Interior Solution,
arXiv:2204.09891.
[37] J. Singh, A construction of approximately self-similar naked singularities for the spherically symmetric
Einstein-scalar field system, arXiv:2210.11325.
[38] P. Yu, Dynamical Formation of black holes due to the condensation of matter field, preprint (2011),
arXiv:1105.5898.
[39] P. Zhao, D. Hilditch, J. Kroon, Trapped surface formation for the Einstein-scalar system, preprint (2023),
arXiv:2304.01695.