Stability of the Spacetime Positive Mass Theorem in Spherical Symmetry

TL;DR

This paper proves a stability version of the spacetime positive mass theorem under spherical symmetry, showing that small mass initial data approximate Minkowski space in a precise geometric sense, with implications for the Penrose inequality.

Contribution

It establishes the stability of the spacetime positive mass theorem in spherical symmetry, linking small ADM mass to convergence towards Minkowski space in a rigorous geometric framework.

Findings

Sequences with small ADM mass approximate Minkowski space

Convergence in the pointed volume preserving intrinsic flat sense

Second fundamental forms converge in $L^p$, $1 extless p extless 2$

Abstract

The rigidity statement of the positive mass theorem asserts that an asymptotically flat initial data set for the Einstein equations with zero ADM mass, and satisfying the dominant energy condition, must arise from an embedding into Minkowski space. In this paper we address the question of what happens when the mass is merely small. In particular, we formulate a conjecture for the stability statement associated with the spacetime version of the positive mass theorem, and give examples to show how it is basically sharp if true. This conjecture is then established under the assumption of spherical symmetry in all dimensions. More precisely, it is shown that a sequence of asymptotically flat initial data satisfying the dominant energy condition, without horizons except possibly at an inner boundary, and with ADM masses tending to zero must arise from isometric embeddings into a sequence of static spacetimes converging to Minkowski space in the pointed volume preserving intrinsic flat sense. The difference of second fundamental forms coming from the embeddings and initial data must converge to zero in $L^p$, $1\leq p<2$. In addition some minor tangential results are also given, including the spacetime version of the Penrose inequality with rigidity statement in all dimensions for spherically symmetric initial data, as well as symmetry inheritance properties for outermost apparent horizons.

Equations

$$16\pi \mu =R_g+({\text{Tr}}_{g}k{)}^2-|k{|}_g^2,8\pi J={\text{div}}_g\left(k-({\text{Tr}}_{g}k)g\right),$$

$$\mu \ge |J{|}_g.$$

$$\left|{\partial }^{{\beta }_1}(g_{ij}-{\delta }_{ij})\right|\le \frac{C}{|x{|}^{n-2+|{\beta }_1|}},\left|{\partial }^{{\beta }_2}{k}_{ij}\right|\le \frac{C}{|x{|}^{n-1+|{\beta }_2|}},$$

$$\left|R_g\right|\le \frac{C}{|x{|}^{n+1}},\left|{\text{Tr}}_{g}k\right|\le \frac{C}{|x{|}^n}.$$

$$E=\frac{1}{2(n-1){\omega }_{n-1}}\underset{r\to \infty }{\lim }{\int }_{S_r}(g_{ij,i}-g_{ii,j}){\nu }^j,$$

$$P_i=\frac{1}{2(n-1){\omega }_{n-1}}\underset{r\to \infty }{\lim }{\int }_{S_r}(k_{ij}-({\text{Tr}}_{g}k)g_{ij}){\nu }^j,$$

$$m=\sqrt{E^2-|P{|}^2}.$$

$$E\ge |P|.$$

$$E=|P|⟹E=|P|=0.$$

$$E=0⟹(M^n,g,k)\text{ embeds as initial data in Minkowski space}.$$

$$F:(M^n,g)\to ({\mathbb{R}}^{1,n},-d{t}^2+g_{\mathbb{E}}),F(x)=(f(x),x),$$

$$g=F^{\ast }(-d{t}^2+g_{\mathbb{E}})=-d{f}^2+g_{\mathbb{E}},$$

$$h=k.$$

$$g=g_{11}(r)d{r}^2+\rho (r{)}^2{g}_{S^{n-1}},k_{ij}=n_i{n}_j{k}_n(r)+(g_{ij}-n_i{n}_j)k_t(r),$$

$${\theta }_{\pm }:=H_{\Sigma }\pm T{r}_{\Sigma }k,$$

$${\Sigma }_A={\rho }^{-1}({\rho }_A)\text{ such that }{\mathrm{Vol}}_g({\Sigma }_A)=A={\omega }_{n-1}{\rho }_A^{n-1}.$$

$${\Omega }_A={\rho }^{-1}[0,{\rho }_A],{\Omega }_{A_1,A_2}={\rho }^{-1}[{\rho }_{A_1},{\rho }_{A_2}],$$

$$T_D({\Sigma }_A)=\{p\in M^n\mid \exists q\in {\Sigma }_{A}with{d}_g(p,q)<D\},$$

$$F_j:(M_j^n,g_j)\to (\mathbb{R}\times {\bar{M}}_j^n,-d{t}^2+{\bar{g}}_j),F_j(x)=(f_j(x),x),$$

$$g_j=F_j^{\ast }(-d{t}^2+{\bar{g}}_j)=-d{f}_j^2+{\bar{g}}_j,$$

$$(\mathbb{R}\times {\bar{M}}_j^n,-d{t}^2+{\bar{g}}_j))\text{ converge to Minkowski space }({\mathbb{R}}^{1,n},-d{t}^2+g_{\mathbb{E}})$$

$$\left(\,\Omega^{j}_{A}\cap T_{D}(\Sigma^{j}_{A})\,,\,\bar{g}_{j}\,\right)\,\,\stackrel{{\scriptstyle\mathcal{VF}}}{{\longrightarrow}}\,\,\left(\,B_{0}(\rho_{A})\,,^{{\color[rgb]{1,1,1}1}}{g_{\mathbb{E}}}\,\right).$$

$$\parallel h_j-k_j{\parallel }_{L^p\left({\Omega }_A^j\cap T_D({\Sigma }_A^j),{\bar{g}}_j\right)}\to 0.$$

$$d_{\mathcal{V}\mathcal{F}}({\Omega }_j,{\Omega }_j^{\prime })=d_{\mathcal{F}}({\Omega }_j,{\Omega }_j^{\prime })+|{\mathrm{Vol}}_j({\Omega }_j)-{\mathrm{Vol}}_{\infty }({\Omega }_j^{\prime })|.$$

$$d_{\mathcal{V}\mathcal{F}}({\Omega }_j,{\Omega }_{\infty })\to 0⟹\mathrm{Vol}(B_{p_j}(r))\to \mathrm{Vol}(B_{p_{\infty }}(r)),$$

$$\text{ the Jang deformations: }\left({\bar{M}}_j^n,{\bar{g}}_j=g_j+d{f}_j^2\right),$$

$$\text{ the conformal transformations: }\left({\tilde{M}}_j^n,{\tilde{g}}_j=u_j^{4/(n-2)}{\bar{g}}_j\right).$$

$$g=\left(1+z^{\prime }(r{)}^2\right)d{r}^2+r^2{g}_{S^{n-1}}$$

$$m_H(\Sigma )=\frac{1}{2}\sqrt{\frac{A}{{\omega }_2}}\left(1-\frac{1}{4\pi }{\int }_{\Sigma }{\left(\frac{H}{2}\right)}^2\right),$$

$$m(s)=\frac{{\rho }^{n-2}(s)}{2}(1-{\rho }^{\prime }(s{)}^2),$$

Full text

E. Bryden acknowledges the support of NSF Grant DMS-1612049. M. Khuri acknowledges the support of NSF Grant DMS-1708798. C. Sormani acknowledges the support of NSF Grant DMS-1612049.

1. Introduction

Let $(M^{n},g,k)$ be an initial data set for the Einstein equations. This means that $(M^{n},g)$ is a complete Riemannian manifold, possibly with boundary, and $k$ is a symmetric 2-tensor representing the second fundamental form of an embedding into spacetime. These satisfy the constraint equations

[TABLE]

where $\mu$ and $J$ are the energy and momentum density of the matter fields, and $R_{g}$ denotes scalar curvature. The dominant energy condition is satisfied if

[TABLE]

We will say that the initial data are asymptotically flat if there is an asymptotic end in the manifold $M^{n}$ that is diffeomorphic to the complement of a ball $\mathbb{R}^{n}\setminus B_{0}(\rho_{0})$, and there exists a constant $C$ such that in the coordinates $x$ provided by this asymptotic diffeomorphism

[TABLE]

for multi-indices $\beta_{1}\leq 2$, $\beta_{2}\leq 1$ and

[TABLE]

These fall-off conditions are modeled on those of the original Schoen-Yau proof of the positive mass theorem [34]. We believe our results should follow assuming the weaker asymptotic decay as in the work of Eichmair, Huang, Lee and Schoen [12, 13]; however, for the sake of simplicity of exposition this will not be done.

With the above setting, the ADM energy and linear momentum of the asymptotic end are finite, well-defined, and given by

[TABLE]

[TABLE]

where $S_{r}$ are coordinate spheres with unit outer normal $\nu$ and $\omega_{n-1}$ is the volume of the standard sphere ${\mathbb{S}}^{n-1}$. The ADM mass is then the Lorentz length of the energy-momentum 4-vector

[TABLE]

In this paper the main results will be concerned with spherically symmetric initial data. It turns out that in spherical symmetry, under the definition of asymptotic flatness in (1.3) and (1.4), the linear momentum vanishes $|P|=0$ and hence $m=E$ as is shown in Proposition 3.6.

The positive mass inequality asserts that an asymptotically flat complete initial data set satisfying the dominant energy condition has

[TABLE]

This was established by Eichmair, Huang, Lee, and Schoen in [13] for dimensions $3\leq n\leq 7$ by using stable marginally outer trapped surfaces (MOTS) in analogy with the minimal hypersurface technique deployed in the time-symmetric case, and in all dimensions $n\geq 3$ for spin manifolds by Bartnik [5] and Witten [40] (see also work of Parker and Taubes [30]). Earlier, the weaker inequality $E\geq 0$ was initially proven by Schoen and Yau [34] when $n=3$ with the help of Jang’s equation, and this reduction argument was later extended by Eichmair [12] to include dimensions $3\leq n\leq 7$.

The rigidity of the positive mass theorem may be broken into two statements. The first asserts:

[TABLE]

This was proven by Huang and Lee [17] for $3\leq n\leq 7$. Their approach only uses the positive mass inequality as input but not its proof, and thus can be extended to higher dimensions for spin manifolds. The spin case was previously treated by Beig and Chrusciel [6] for $n=3$ and Chrusciel and Maerten [10] for higher dimensions. The second statement is that

[TABLE]

As with the inequality, this was originally established by Schoen and Yau in [34] for three dimensions and extended by Eichmair in [12] to dimensions less than eight. Finally in the spin case this was treated for all dimensions in work of Beig, Chrusciel, and Maerten [6, 10]. Here we state the positive mass rigidity theorem in a particular way that allows us to propose a natural almost rigidity (or stability) conjecture.

Theorem 1.1** (Positive Mass Rigidity Theorem [12, 17, 34]).**

Let $(M^{n},g,k)$ be a complete asymptotically flat initial data set, with $3\leq n\leq 7$, and satisfying the dominant energy condition. If the ADM mass vanishes $m=0$, then $M^{n}$ is diffeomorphic to $\mathbb{R}^{n}$ and $(M^{n},g)$ can be isometrically embedded as a graph in Minkowski space. That is

[TABLE]

where $g_{\mathbb{E}}$ is the Euclidean metric and

[TABLE]

and the second fundamental form, $h$, of the embedding agrees with that of the initial data

[TABLE]

The purpose of this paper is to establish an almost rigidity or stability version of this theorem in the spherically symmetric setting. We will say that the initial data are spherically symmetric if $M^{n}$ is diffeomorphic to $\mathbb{R}^{n}\setminus B_{0}(r_{0})$ or $\mathbb{R}^{n}$ and the metric and second fundamental form may be expressed by

[TABLE]

for some radial functions $g_{11}$, $\rho$, $k_{n}$, and $k_{t}$, where $n=\sqrt{g^{11}}\partial_{r}$ is the unit normal to coordinate spheres. This decomposition for $k$ exhibits its normal and tangential components with respect to the coordinate spheres, and is motivated by the implicit assumption that the initial data come from a spherically symmetric spacetime in which $k$ is the ‘time derivative’ of $g$ which already has this structure.

The boundary, if nonempty, of the initial data will consist of apparent horizons. Recall that the strength of the gravitational field around a hypersurface $\Sigma\subset M^{n}$ may be measured by the null expansions (null mean curvatures) given by

[TABLE]

where $H_{\Sigma}$ is the mean curvature with respect to the unit normal pointing towards spatial infinity. These quantities can be interpreted as the rate at which the area of a shell of light changes as it moves away from the surface in the outward future/past direction ($+$/$-$). Future or past trapped surfaces are defined by the inequalities $\theta_{+}<0$ or $\theta_{-}<0$, respectively, and may be thought of as lying in a region of strong gravity. If $\theta_{+}=0$ or $\theta_{-}=0$, then $\Sigma$ is called a future or past apparent horizon; these naturally arise as boundaries of future or past trapped regions. Furthermore, such surfaces will be referred to as an outermost apparent horizon if it is not enclosed by any other apparent horizon. In Lemma 3.1 it is shown that the outermost apparent horizon inherits the symmetry of its ambient space. In this text the abbreviated term horizon will often be used for these objects.

We will consider asymptotically flat $(M^{n},g,k)$ that have either no horizons or only a horizon on an inner boundary, in which case the boundary is an outermost apparent horizon. Under these conditions for spherically symmetric initial data, it is shown in Lemma 3.4 that the areas (or $n-1$ dimensional volumes) of the level sets of $\rho$ are increasing. Thus we may define the level set

[TABLE]

We will study regions within and between these level sets

[TABLE]

as well as their tubular neighborhoods

[TABLE]

where $d_{g}(p,q)$ denotes the distance between $p$ and $q$.

In order to state the main theorem we need the notion of uniform asymptotic flatness. A sequence of initial data $(M^{n}_{j},g_{j},k_{j})$ will be referred to as uniformly asymptotically flat if each member of the sequence is asymptotically flat according to (1.3) and (1.4), and the constants $\rho_{0}$ and $C$ in the definition are independent of $j$.

Theorem 1.2**.**

Fix $A>0$ and $D>\rho_{A}$, and consider a sequence of uniformly asymptoticaly flat spherically symmetric initial data sets $(M_{j}^{n},g_{j},k_{j})$ satisfying the dominant energy condition and with no closed horizons except possibly the inner boundary. If their ADM masses converge to zero $m_{j}\rightarrow 0$, then there exist Riemannian manifolds $(\bar{M}^{n}_{j},\bar{g}_{j})$ diffeomorphic to $(M^{n}_{j},g_{j})$ with graphical isometric embeddings

[TABLE]

[TABLE]

such that the static spacetimes

[TABLE]

in that the base manifolds converge in the pointed volume preserving intrinsic flat sense to Euclidean space. More precisely, regions within $\Sigma_{A}$ in $(\bar{M}^{n}_{j},\bar{g}_{j})$ converge to balls in Euclidean space

[TABLE]

Furthermore if there is a uniform constant $C$ such that $\parallel k_{j}\parallel_{L^{2}(M^{n}_{j})}\leq C$, then for any $1\leq p<2$ the second fundamental forms $h_{j}$ of the graphs satisfy

[TABLE]

The intrinsic flat distance $d_{\mathcal{F}}(\Omega_{j},\Omega^{\prime}_{j})$ between pairs of compact oriented Riemannian manifolds with boundary was first introduced by the third author with Wenger in [38]. Intuitively it measures the filling volume between the given manifolds. It is [math] if and only if there is an orientation preserving isometry between the manifolds $\Omega_{j}$ and $\Omega^{\prime}_{j}$ [38]. The volume preserving intrinsic flat distance was introduced in [36] and includes an extra term involving the global difference of volumes

[TABLE]

This has been studied by Portegies in [32] and by Jauregui-Lee in [20]. In particular they have shown that

[TABLE]

for sequences of points $p_{j}\in\Omega_{j}$ converging to $p_{\infty}\in\Omega_{\infty}$.

Theorem 1.2 has been proven for time-symmetric initial data sets by Lee and the third author in [23]. In that setting $k_{j}=0$ and $f_{j}$ can be taken to be constant so that $h_{j}=0$ and (1.23) follows trivially. The intrinsic flat convergence is proven in [23] by constructing an explicit filling manifold. In fact, LeFloch and the third author have proven in [24] that the metric tensors converge in the $H^{1}_{loc}$ sense. Note that examples in [23] demonstrate that even in the spherically symmetric time-symmetric setting one can have sequences with ADM mass converging to [math] which do not converge to regions in Euclidean space in the smooth or Gromov-Hausdorff sense. Applying techniques from [23] in our Example 2.4, it is shown why one needs a tubular neighborhood in (1.22). Furthermore, in Example 2.1 and Example 2.2 we demonstrate the need to assume that there are no interior horizons.

Without time symmetry, when $k_{j}\neq 0$, Theorem 1.2 makes no claim as to the convergence of the original sequence of Riemannian manifolds $(M_{j}^{n},g_{j})$. Example 2.7 illustrates why the initial sequence $(M_{j}^{n},g_{j})$ need not converge in any reasonable sense. There we construct sequences of initial data sets of zero mass lying in Minkowski space which become increasingly null on large regions, so that volumes disappear instead of converging.

Conjecture 1.3**.**

Theorem 1.2 holds without requiring spherical symmetry when suitable definitions are made for the regions $\Omega_{A}$. To achieve the conclusion exactly as stated we expect that $E\to 0$ should replace $m\to 0$ in the hypotheses for the general case. It is possible that a similar statement holds for $m\to 0$, but the approach would have to be different from the one used here in light of examples with boost. In the outline below, we clarify which steps strongly use spherical symmetry and which hold more generally.

The corresponding almost rigidity or stability conjecture in the time-symmetric case was stated and proven in the spherically symmetric setting by Lee and the third author in [23]. It has been confirmed in the graph setting by Huang, Lee, and the third author in [18] and for geometrostatic manifolds by the third author with Stavrov in [37]. Initial controls on the metric tensor towards proving the time-symmetric conjecture have been found by Allen for regions covered by smooth inverse mean curvature flow in [1], and by the first author for axisymmetric manifolds [8].

With the definition of asymptotic flatness used here the ADM mass $m$ and ADM energy $E$ agree in spherical symmetry since the linear momentum $P$ vanishes (Proposition 3.6). In general when mass and energy differ, Conjecture 1.3 could be quite subtle in the case of large linear momentum, as the construction of the base manifolds $(\bar{M}^{n},\bar{g})$ presented here does not behave well in this setting. The methods used here and based on the Jang equation are tailored to the situation when $E$ is small, which will not be the case if $|P|$ stays uniformly away from zero.

We now give an outline of the proof of Theorem 1.2, which is modeled on the Schoen-Yau approach to the positive mass theorem [34]. It should be noted that some of the arguments do not require spherical symmetry. The first step is to solve the 2nd order quasi-linear elliptic Jang equation for each $(M_{j}^{n},g_{j},k_{j})$ to obtain solutions $f_{j}$ with asymptotically cylindrical blow-up at the outermost apparent horizon. See Section 4 for details. The original study of this equation in [34] observed that the solution only blows-up at apparent horizons. Prescribed blow-up at the outermost horizon was obtained in work of Eichmair, Han, the second author, and Metzger for low dimensions [12, 14, 27]. In the spherically symmetric case, the equation can be reduced to a 1st order ODE and the desired solutions can be produced in any dimension [Theorem 4.1]. From the solutions a sequence of Riemannian manifolds,

[TABLE]

can be constructed which serve as the base for the ambient static spacetimes of Theorem 1.2. Schoen-Yau showed that the scalar curvature of the Jang metric is nonnegative modulo a divergence term as stated in (4.8). The Jang manifolds are also uniformly asymptotically flat [Lemma 4.3], and have the same ADM masses as the original initial data [Corollary 4.2]. A primary difference is that they have a cylindrical end where previously there was a boundary. In Example 2.5 we explicitly solve the Jang equation for a constant time slice of the Schwarzschild spacetime so that one can see precisely how this step may be implemented constructively.

The nonnegativity property of the scalar curvature of $(\bar{M}_{j}^{n},\bar{g}_{j})$ allows one to further conformally transform the Jang deformations to Riemannian manifolds of zero scalar curvature,

[TABLE]

Along the cylindrical ends $u_{j}$ decays exponentially fast to zero and hence conformally closes this end (see [34] and Proposition 5.1). In general the masses of the conformal deformations converge to zero. Thus if no horizons are present or one restricts attention to domains outside the outermost minimal surface, it is expected (by the time-symmetric almost rigidity conjecture) that regions in $(\tilde{M}^{n}_{j},\tilde{g}_{j})$ converge to balls in Euclidean space in the intrinsic flat sense. The hope is then to prove that the conformal factors $u_{j}$ are sufficiently close to $1$ in order to establish that regions in $(\bar{M}^{n}_{j},\bar{g}_{j})$ converge to balls in Euclidean space as well. See Remark 5.2.

In the spherically symmetric setting the conformal deformation is Euclidean space $\tilde{g}_{j}=g_{\mathbb{E}}$, so the mass is [math] and there are no horizons (cf. Lemma 5.3). Therefore $\bar{g}_{j}$ is related to the Euclidean metric via the conformal factor $u_{j}$, and establishing (1.22) is reduced to controlling $u_{j}$. A global $L^{2}$ gradient bound in terms of the mass is obtained from the stability property associated with the Jang surface in Lemma 5.6, and this is parlayed into $C^{0,\frac{1}{2}}$ control away from the center of the manifold in Proposition 5.10 by using the uniform asymptotically flat assumption and Lemma 5.9. We then have that $u_{j}\rightarrow 1$ uniformly on appropriate subdomains avoiding the center. This work is completed in Section 5.

The volume preserving intrinsic flat convergence (1.22) is proven in Section 6 in two main steps. The overall approach is to apply a result of Lakzian and the third author [22], as stated in Proposition 6.1, to achieve intrinsic flat convergence. First, control on $u_{j}$ is used to show that regions avoiding the center are smoothly close to annuli in Euclidean space (see Lemma 6.2). Secondly, we prove that the volumes of the regions closer to the center are small using area monotonicity and the coarea formula in Lemma 6.5. The remaining required terms are estimated in various lemmas throughout the section. Moreover, volume convergence follows from the above and is stated in Lemma 6.8. Without spherical symmetry one might imagine doing something similar, cutting out many wells rather than just the center as in joint work of the third author with Stavrov in [37], or using a completely different approach as in joint work of the third author with Huang and Lee in [18].

In Section 7 convergence of second fundamental forms (1.23) is established, where the proof relies on nonnegativity of the spacetime Hawking mass. Control on $|h_{j}-k_{j}|_{\bar{g}_{j}}$ away from the center is given in Proposition 7.1 using estimates for the conformal factors $u_{j}$, as well as the stability property associated with the Jang surface. While $|h_{j}-k_{j}|_{\bar{g}_{j}}$ might be large near the center, with the additional hypothesis on $k_{j}$ and the small volume inside in Proposition 7.2, convergence in the desired tubular neighborhood is achieved in Theorem 7.3. Finally, in Section 8 we prove Theorem 1.2 using all the above.

Acknowledgements: The authors would like to thank Dan Lee for helpful conversations, and Walter Simon for comments. Christina Sormani gratefully acknowledges office space in the Simons Center for Geometry and Physics, Stony Brook University at which most of the research for this paper was performed, and also support from a CUNY Fellowship Leave. Edward Bryden would like to thank the CUNY Graduate Center for allowing him to visit in Spring 2019 during which this paper was completed.

2. Examples

In this section we provide some examples which illustrate the importance of various hypotheses in Theorem 1.2, and some intuition as to what is happening in the proof. In the time-symmetric setting, where the objects of study are manifolds with nonnegative scalar curvature, examples are given with closed interior horizons (Examples 2.1 and  2.2) that fail to have volume preserving intrinsic flat convergence to Euclidean space.

The additional assumption of no closed interior horizons and no boundary is also considered. In this setting, a time-symmetric initial data set is a graph over itself and the solution to Jang’s equation is constant. An example within this context is provided that contains a deep well, demonstrating why tubular neighborhoods are introduced in order to obtain volume preserving intrinsic flat convergence (Example 2.4).

The proof of Theorem 1.2 will then be applied to slices of the Schwarzschild spacetime, so one can see what happens when there is a boundary horizon. Difficulties arise even in this time-symmetric example because solutions to Jang’s equation blow-up near the horizon so that the base spaces $(\bar{M}^{n},\bar{g})$ possess an asymptotically cylindrical end. In Example 2.5 we see exactly how Schwarzschild slices arise as graphs over base spaces which are close to Euclidean space in the volume preserving intrinsic flat sense. This clarifies why the proof of Theorem 1.2 is delicate when the manifolds have boundary.

Finally we consider examples which are not time-symmetric. In Example 2.7 it is demonstrated how even the restriction to sequences of spacelike graphs $(M_{j}^{n},g_{j},k_{j})$ in Minkowski space does not allow for proper control over the original sequence of Riemannian manifolds $(M_{j}^{n},g_{j})$. This justifies why Theorem 1.2 only deals with convergence of the base manifolds, and does not address convergence of the given sequence of initial data.

2.1. Horizons in Time-Symmetric Examples

The assumption of no closed interior horizons is necessary to avoid the formation of bubbles and other phenomena which may occur behind a horizon. Since the inner boundary is allowed to be a horizon, these hypotheses mean that the main theorem applies within the domain of outer communication. This is consistent with the basic intuition that the ADM mass cannot effectively ‘see’ within a black hole. The following well-known examples explain why one cannot hope to obtain volume preserving intrinsic flat convergence without these assumptions.

Example 2.1**.**

Riemannian Schwarzschild space is a constant time slice of the Schwarzschild spacetime, with a metric that can be written as

[TABLE]

where $z^{\prime}(r)=\sqrt{2m/(r-2m)}$. It has a horizon (minimal surface) at $r=2m$, and can be extended smoothly past the horizon by writing $r$ as a function of $z$. In fact the graph is that of a parabola. If we take a sequence of Riemannian Schwarzschild spaces of smaller and smaller mass $m$, this parabola becomes more vertical and the vertex decreases to the origin. See Figure  1. This sequence of Riemannian Schwarzschild spaces converges smoothly to Euclidean space on compact sets that avoid the increasingly thin necks, and by any weak notion of convergence is seen to converge to a double sheeted Euclidean space. The volumes of balls centered around points on the horizon converge to twice the volume of a Euclidean ball. It is only by removing the part behind the horizon that one may consider the limit to be a single Euclidean space, and obtain volume preserving convergence to Euclidean space.

Example 2.2**.**

Start with the Riemannian Schwarzschild initial data viewed as a parabola using the function $r(z)$ described in the previous example. Keep the region outside the horizon exactly isometric to Riemannian Schwarzschild, but behind the black hole attach a round sphere in a $C^{1}$ way. This is achieved by ensuring that the induced metrics and mean curvatures of the interface surfaces agree from both sides. This guarantees that the scalar curvature is distributionally nonnegative across the interface surface. Furthermore since Riemannian Schwarzschild is scalar flat and the sphere has positive scalar curvature, applying Ricci flow for a very short time, one obtains a smooth metric which is $C^{1}$ close to the original and has positive scalar curvature everywhere. The almost spherical region of the resulting manifold is called a bubble. See Figure 2. Now perform this construction with a sequence of Schwarzschild spaces having masses converging to zero, while keeping the bubble the same size throughout the sequence. By any notion of weak convergence this sequence converges to a Euclidean space with a sphere attached to it. In analogy with the previous example, balls of a fixed radius centered at points on the horizon have volumes converging to the sum of the volume of a ball of the same radius in Euclidean space plus a ball of the same radius in the sphere. Again we do not obtain volume preserving intrinsic flat convergence to Euclidean space, unless the part inside the horizon is cut out.

2.2. Deep Wells in Time-Symmetric Examples

Let us now consider Theorem 1.2 in the time-symmetric case where there are no horizons and no boundary. In such a setting the solution to Jang’s equation is constant, so the theorem states that volume preserving intrinsic flat convergence occurs within the initial data themselves, as opposed to convergence of ambient spacetimes in which the data embed. This was established by D. Lee and the third author in [23]. In this subsection we first recall in Lemma 2.3 an example construction technique from [23]. We then review an example, Example 2.4, with deep wells demonstrating the need to use tubular neighborhoods to obtain volume preserving intrinsic flat convergence.

Recall the definition of Hawking mass for a surface $\Sigma$ in a Riemannian 3-manifold

[TABLE]

where $A$ and $H$ are the area and mean curvature of $\Sigma$. As described in Section 3.3 this may be generalized in spherical symmetry to higher dimensions $n$ by

[TABLE]

where the metric is expressed in a radial arclength coordinate $s$ and with area radius function $\rho(s)$. The first variation of Hawking mass becomes

[TABLE]

which is nonnegative for nondecreasing area radius functions and nonnegative scalar curvature.

Lemma 2.3**.**

Let $\mathcal{M}$ denote the collection of asymptotically flat spherically symmetric manifolds

[TABLE]

with nonnegative scalar curvature that have no closed interior minimal surfaces and either no boundary, or minimal surface boundary $r^{-1}(r_{min})$. Let $\mathcal{H}$ be the collection of admissible Hawking mass functions, that is increasing functions $m:[r_{min},\infty)\to\mathbb{R}$ such that

[TABLE]

and

[TABLE]

for $r>r_{\textrm{min}}\geq 0$. There is a constructive bijection between $\mathcal{M}$ and $\mathcal{H}$ such that

[TABLE]

In [23] this result was used to construct an example with an arbitrarily deep well. Here we also describe the volumes in this example, justifying the necessity of cutting off the region using tubular neighborhoods of fixed size $D$ to obtain volume preserving intrinsic flat convergence in Theorem 1.2.

Example 2.4**.**

Given $A>0$, $L>0$, and $\delta>0$ there exists $(M^{n},g)\in\mathcal{M}$ with ADM mass $m<\delta$ such that the distance $d(\Sigma_{min},\Sigma_{A})>L$, where $\Sigma_{A}$ is a symmetry sphere with area $A$ and $\Sigma_{min}$ is either the boundary $\partial M^{n}$ or the pole. See Figure 3. In fact the example is constructed by first choosing a radius $r_{\epsilon}$ depending on $\delta$, and then constructing the admissible Hawking mass function so that the distance between the levels $r_{\epsilon}$ and $r_{\epsilon}/2$ is an arbitrary value $L$. Thus the volume between these level sets, computed with the coarea formula, provides a lower bound

[TABLE]

where $A_{\delta}$ is the area of the level set $r^{-1}(r_{\epsilon})$ which depends on $\delta$ but not $L$. Taking a sequence with

[TABLE]

we obtain a sequence of spherically symmetric examples which are increasingly deep and have

[TABLE]

By (1.22) of Theorem 1.2 it holds that for any fixed $D>0$,

[TABLE]

This fixed distance $D>0$ of the tubular neighborhood is needed to cut off the arbitrarily large volumes in the arbitrarily deep wells.

2.3. Riemannian Schwarzschild Space

Let us now consider Theorem 1.2 in the time-symmetric case where there is a boundary horizon; for simplicity we restrict the discussion in this subsection to dimension $n=3$. Even though this setting is time-symmetric the solutions of Jang’s equation will not be constant, rather they will blow-up at the horizon boundary. These solutions are then used to embed the initial data $(M^{3},g,0)$ into $({\mathbb{R}}\times\bar{M}^{3},-dt^{2}+\bar{g})$ as a graph over the base space $(\bar{M}^{3},\bar{g})$. This base will have different properties than the original data, in particular it will have an asymptotically cylindrical end. Next an appropriate conformal factor $u$ is found so that the new metric $\tilde{g}=u^{4}\bar{g}$ is scalar flat. These are some of the main steps in the proof of Theorem 1.2, and will in this subsection be computed explicitly for Riemannian Schwarzschild initial data.

Example 2.5**.**

Recall that the induced metric on a time slice $M^{3}=r^{-1}[2m,\infty)$ of the Schwarzschild spacetime of mass $m$ can be written in the form

[TABLE]

The Jang equation may be solved explicitly in this case for a blow-up solution. To see this observe that from [7] and the discussion in Section 4.2, Jang’s equation may be reduced to a first order ODE by setting

[TABLE]

where $g^{11}=1-2m/r$. Namely, the Jang equation in this case becomes simply

[TABLE]

and the blow-up solution is $v=\left(\frac{2m}{r}\right)^{2}$. It follows that

[TABLE]

so that

[TABLE]

Therefore the Jang metric is

[TABLE]

with

[TABLE]

This clearly has an asymptotically cylindrical end as $r\rightarrow 2m$. Figure 4 illustrates how the Riemannian Schwarzschild geometry embeds as a graph over this Jang deformation. Note that it becomes increasingly null upon approach to the horizon.

The conformal deformation to zero scalar curvature can also be given explicitly for this Schwarzschild example. To do this let $r=r(\tilde{r})$ be such that

[TABLE]

then

[TABLE]

We may solve for $\tilde{r}$ in terms of $r$ by

[TABLE]

Now set $u^{-4}=(r/\tilde{r})^{2}$ so that ${g_{\mathbb{E}}}=u^{4}\bar{g}$. Thus, $u=\sqrt{\tilde{r}/r}$ serves as the desired conformal factor yielding a scalar flat deformation. The fact that this conformal change resulted in a Euclidean metric is not special to the Schwarzschild example, as will be seen in Section 5.

Figure 5 depicts a sequence of Riemannian Schwarzschild manifolds $(M^{3}_{j},g_{j})$ with masses $m_{j}\rightarrow 0$, embedded as graphs over a sequence of base Jang deformations $(\bar{M}^{3}_{j},\bar{g}_{j})$ with asymptotically cylindrical ends that converge in the pointed volume preserving intrinsic flat sense to Euclidean space.

2.4. Graphs in Minkowski Space

If $(M_{j}^{n},g_{j},k_{j})$ arise as spacelike graphs $t=f_{j}(x)$ in Minkowski space, then the Jang metric $\bar{g}_{j}$ obtained by solving Jang’s equation is exactly the Euclidean metric and $k_{j}$ is the second fundamental form of the graph. In the notation of Theorem 1.2 we have

[TABLE]

Theorem 1.2 is trivially true and there is nothing to prove. On the other hand, such examples can exhibit pathological behavior from the point of view of establishing volume preserving intrinsic flat convergence of $(M_{j}^{n},g_{j})$. In this subsection an example is presented to demonstrate why we say nothing about the limiting behavior of $(M_{j}^{n},g_{j})$ in Theorem 1.2. In particular, even for sequences of spacelike graphs in Minkowski space one cannot hope for more than subsequential convergence, and the limiting space need not be well-behaved.

Remark 2.6**.**

If $(M_{j}^{n},g_{j},k_{j})$ arise as spacelike graphs $t=f_{j}(x)$ in Minkowski space, then $g_{j}=-df_{j}^{2}+{g_{\mathbb{E}}}$ is positive definite, so the Lipschitz norm satisfies $Lip_{g_{\mathbb{E}}}(f_{j})<1$. Thus by Arzela-Ascoli a subsequence of the $f_{j}$ converge to a Lipschitz function $f_{\infty}$ with $Lip_{g_{\mathbb{E}}}(f_{j})\leq 1$. However, the graph of $f_{\infty}$ need not be spacelike!

Example 2.7**.**

Consider a sequence of graphs $t=f_{j}(x)$ in Minkowski space which are spherically symmetric, with $Lip_{g_{\mathbb{E}}}(f_{j})<1$, and converging to $f_{\infty}$ that satisfies

[TABLE]

as in Figure 6, where $r$ is the radial distance function for $g_{\mathbb{E}}$. It can be arranged that $g_{j}$ converge in the $C^{0}$ sense. The limit is a semidefinite metric $g_{\infty}$ with

[TABLE]

By the coarea formula

[TABLE]

since

[TABLE]

In contrast

[TABLE]

Thus we find nice behavior of the base spaces $(\bar{M}_{j}^{n},\bar{g}_{j})$ as described in Theorem 1.2, with pathological limiting behavior for the original sequence $(M_{j}^{n},g_{j})$.

3. Manifolds with Spherical Symmetry

In this section we prove that outermost apparent horizons inherit the symmetries of the asymptotically flat initial data sets in which they lie (Lemma 3.1), that the areas of symmetry spheres are monotonic in spherically symmetric initial data sets without horizons or with outermost apparent horizon boundary (Lemma 3.4), and establish the spacetime Penrose Inequality in all dimensions under the assumption of spherically symmetry (Theorem 3.5). These results are of use to us when proving Theorem 1.2. Prior work in these directions is reviewed within.

3.1. Horizons in Initial Data With Symmetry

Lemma 3.1**.**

Let $(M^{n},g,k)$ be an asymptotically flat initial data set which admits a continuous symmetry with generator $\eta$. This means that $\eta$ is a Killing field which leaves $k$ invariant, and thus the following Lie derivatives vanish $\mathfrak{L}_{\eta}g=\mathfrak{L}_{\eta}k=0$. If the outermost apparent horizon is smooth, then $\eta$ must be tangential to it.

In particular, if $(M^{n},g,k)$ is spherically symmetric with smooth outermost apparent horizon then this surface is also spherically symmetric.

Remark 3.2**.**

The existence of outermost apparent horizons due to appropriate trapping is proven in [4, 11]. Like minimal surfaces, they are shown to have a singular set which is no larger than codimension 7. Thus for $2\leq n\leq 7$ the outermost apparent horizon is smooth.

Remark 3.3**.**

The symmetry inheritance property for stable MOTS was already observed in Theorem 8.1 of [3] (see also [29]) when $n=3$. There a spacetime perspective was taken, as opposed to the initial data point of view used here.

Proof.

The following argument is a generalization of that in [9] for outermost minimal surfaces in axisymmetry. Suppose that the outermost apparent horizon $\Sigma$ does not admit the stated symmetry. Then the Killing field $\eta$ is not tangential to $\Sigma$ at all points. Thus, if $\varphi_{t}$ denotes the flow of this Killing field so that $\partial_{t}\varphi_{t}=\eta\circ\varphi_{t}$, then there is a nonzero $t_{0}$ near zero such that a domain within $\varphi_{t_{0}}(\Sigma)$ lies outside of $\Sigma$. Furthermore, observe that since $(M^{n},g,k)$ is invariant under the action of $\varphi_{t}$ the surface $\varphi_{t_{0}}(\Sigma)$ is an apparent horizon of the same type.

Consider now the compact set $\mathcal{U}$ which is the union of all smooth compact embedded apparent horizons within $M^{n}$, and define the trapped region $\mathcal{T}$ to be the union of $\mathcal{U}$ with all the bounded components of $M^{n}\setminus\mathcal{U}$. As described in [2, Theorem 3.3] the outermost apparent horizon arises as the boundary $\partial\mathcal{T}$, moreover it is embedded and smooth away from a singular set of Hausdorff codimension at most 7; in fact it will be smooth by the assumptions of this lemma. Because $\partial\mathcal{T}$ must enclose both $\Sigma$ and $\varphi_{t_{0}}(\Sigma)$, it cannot agree with $\Sigma$ at all points. This, however, contradicts the outermost assumption for $\Sigma$. ∎

3.2. Monotonicity of Area

In this and the following subsection, the discussion is relevant to the spacetime Penrose inequality in spherical symmetry. The next result may be derived from the arguments in [26, Section 4].

Lemma 3.4**.**

Let $(M^{n},g,k)$ be a spherically symmetric asymptotically flat initial data set as in (1.14). Outside of the outermost apparent horizon, the area of symmetry spheres in $(M^{n},g)$ is a strictly increasing function of the radial coordinate. In particular, the warping function $\rho$ defining $g$ is also an increasing function.

Proof.

Let $S_{r}$ denote the level sets of $r$. Since

[TABLE]

the null expansions (null mean curvatures) are given by

[TABLE]

Since the null expansions are both positive near infinity, as the mean curvature $H$ dominates $\text{Tr}_{S_{r}}k$ according to decay rates, when moving inwards from infinity they must remain positive outside of the outermost apparent horizon where $\theta_{+}=0$ or $\theta_{-}=0$. Note that here we are using Lemma 3.1 which asserts that the outermost apparent horizon (if present) is one of the spheres $S_{r_{0}}$. Therefore

[TABLE]

To finish the proof simply add the two null expansions, and observe that since both are positive we obtain

[TABLE]

Hence, the warping function $\rho$ is an increasing function. ∎

3.3. The Spacetime Penrose Inequality in All Dimensions

The spherically symmetric Penrose inequality without the maximal assumption was established in dimension $n=3$ in [15], although the case of equality was not treated. A similar result is stated in [16] for all dimensions, but the hypotheses are too strong for our purposes and they also do not address the case of equality. The full result including the case of equality was given in [7] for $n=3$, and here we easily extend it to all dimensions.

Theorem 3.5**.**

Let $(M^{n},g,k)$, $n\geq 3$ be an asymptotically flat spherically symmetric initial data set satisfying the dominant energy condition, and let $\mathcal{A}_{0}$ denote the area of the outermost apparent horizon. Then

[TABLE]

and equality holds if and only if the initial data outside the outermost apparent horizon arise from an embedding into the Schwarzschild spacetime. In particular, for a sequence of initial data with $m_{j}\rightarrow 0$ we have $\mathcal{A}_{j}\rightarrow 0$.

Proof.

We will follow and generalize the arguments of [7] to higher dimensions. The proof is based on the generalized Jang equation introduced in [7]. The corresponding Jang deformation is similar to that of the original with the addition of an extra function $\phi$ that plays the role of warping factor for embeddings into a static spacetime, namely $\bar{g}=g+\phi^{2}df^{2}$ where $f$ satisfies equation (3) of [7] for a canonical choice of $\phi$. The scalar curvature of the spherically symmetric generalized Jang metric

[TABLE]

is given by

[TABLE]

Here $\bar{s}$ denotes radial distance from the boundary so that $\bar{s}=0$ corresponds to the outermost apparent horizon. The fact that the outermost apparent horizon is a level set of $\bar{s}$ is a consequence of Lemma 3.1.

By comparing arbitrary spherically symmetric metrics to that of Schwarzschild we may derive and generalize the Hawking mass (Misner-Sharp mass in spherical symmetry [28]) to higher dimensions

[TABLE]

where

[TABLE]

A direct computation yields

[TABLE]

Therefore integrating produces

[TABLE]

We may now choose $\phi=\rho_{\bar{s}}$ and follow the arguments in [7, page 750]. This allows one to integrate away the divergence term appearing in $\bar{R}$, leaving only nonnegative terms on the right-hand side of (3.11). It follows that $\bar{m}(\infty)\geq\bar{m}(0)$, and this gives the desired inequality

[TABLE]

since the ADM mass of the Jang metric agrees with that of the given initial data in addition to the fact that the area of the apparent horizon agrees in both metrics as well. The case of equality follows directly from the arguments of [7]. ∎

3.4. Decay of the Second Fundamental Form

We now prove that under the definition of asymptotic flatness for spherical symmetry given in Section 1, the ADM linear momentum vanishes $|P|=0$, and hence the ADM mass coincides with the ADM energy $m=E$.

Proposition 3.6**.**

Under the asymptotic decay conditions (1.3) and (1.4), a spherically symmetric initial data set $(M^{n},g,k)$ satisfies the stronger decay

[TABLE]

In particular, the ADM linear momentum vanishes $|P|=0$ and the ADM mass agrees with the ADM energy $m=E$.

Proof.

Recall that the spherically symmetric initial data $(M^{n},g,k)$ may be expressed by

[TABLE]

where $n=\partial_{s}$. Consider the divergence constraint

[TABLE]

Observe that

[TABLE]

and therefore

[TABLE]

It follows that

[TABLE]

Furthermore

[TABLE]

and hence

[TABLE]

The only nonzero component is in the $\partial_{s}$ direction, and from this we find that

[TABLE]

Since

[TABLE]

this may be rewritten as

[TABLE]

A priori the assumed decay for $k$ leads to $k_{t}=O(s^{1-n})$. However the ODE (3.23) shows that the fall-off is stronger. Indeed, from the dominant energy condition and the assumed decay

[TABLE]

and thus

[TABLE]

It then follows from the trace decay in (3.24) that $k_{n}$ also satisfies the fall-off in (3.25), and hence

[TABLE]

which implies that the ADM linear momentum vanishes $|P|=0$. ∎

4. Solving Jang’s Equation to Obtain the Base Manifolds

4.1. Review of Jang’s Equation Without Symmetry

Given an initial data set $(M^{n},g,k)$, the following quantities may be used to measure how far away it is from being realized as a graph $t=f(x)$ in Minkowski space

[TABLE]

In particular, such an embedding exists if and only if $\bar{g}=g_{\mathbb{E}}$ and $\bar{k}=0$. A necessary condition for this to occur is the Jang equation [2]

[TABLE]

where $f^{a}=g^{ab}\partial_{b}f$. Thus, one may think of Jang’s equation as an attempt to find a candidate graph for an embedding into Minkowski space.

Even though such an embedding may not exist for the given initial data, we may use these ideas to construct an isometric embedding into a relevant static spacetime. Namely consider the map

[TABLE]

defined by

[TABLE]

Then

[TABLE]

and the second fundamental form is

[TABLE]

Another motivation for the Jang equation which is pertinent to the positive mass theorem, is to consider it as a method for deforming initial data to obtain weakly nonnegative scalar curvature. In this setting one may view the Jang graph as a submanifold of the $(n+1)$-dimensional dual Riemannian manifold $(\mathbb{R}\times M^{n},dt^{2}+g)$. The Jang metric $\bar{g}$ is then the induced metric on the graph and $h$ is again the second fundamental form. If $k$ is extended trivially off the $t=0$ slice to the whole $(n+1)$-dimensional ambient space, then the Jang equation simply states that the Jang surface $\bar{M}^{n}$ satisfies the apparent horizon equation

[TABLE]

where $H_{\bar{M}^{n}}$ is the mean curvature of the Jang surface. A computation [34] then shows that the scalar curvature of the Jang metric is nonnegative modulo a divergence term whenever the dominant energy condition is satisfied, that is

[TABLE]

where

[TABLE]

This positivity property for $R_{\bar{g}}$ allows one to conformally transform $\bar{g}$ to zero scalar curvature. Hence through the Jang deformation combined with a conformal transformation, the initial data is taken into the time-symmetric setting.

4.2. Existence of Solutions to Jang’s Equation

The Jang deformation preserves uniform asymptotic flatness as will be shown below, and preserves the mass so that $\bar{m}=m$. An interesting feature of the Jang equation’s existence theory is its ability to detect apparent horizons. That is, it can only blow-up at apparent horizons in which case it approximates a cylinder over these surfaces [34]. In dimension $n=3$ it has been shown that this cylindrical blow-up behavior can in fact be prescribed at the outermost apparent horizon [12, 14, 27]. Suppose that the boundary is decomposed into a disjoint union of future ($+$) and past ($-$) apparent horizon components $\partial M^{3}=\partial_{+}M^{3}\cup\partial_{-}M^{3}$, and that there are no other apparent horizons present. If in a neighborhood of $\partial_{\pm}M^{3}$ there are constants $l\geq 1$ and $c>0$ such that

[TABLE]

where $\tau(x)=\operatorname{dist}(x,\partial M^{3})$ and $S_{\tau}$ are surfaces of constant distance to the boundary, then there exists a smooth solution $f$ of Jang’s equation with the property that $f(x)\rightarrow\pm\infty$ as $x\rightarrow\partial_{\pm}M^{3}$. Furthermore, the asymptotics for this blow-up are given by

[TABLE]

for some positive constants $c_{1}$, $c_{2}$. In the spherically symmetric case, this type of existence result may be established in all dimensions.

Recall the form of the spherically symmetric initial data

[TABLE]

defined on the compliment of a ball $M^{n}=\mathbb{R}^{n}\setminus B_{0}(r_{0})$. It is assumed that $\partial M^{n}=S_{r_{0}}$ is the only apparent horizon, which means that the null expansions satisfy

[TABLE]

and that either $\theta_{+}(r_{0})=0$, $\theta_{-}(r_{0})=0$, or $\theta_{+}(r_{0})=\theta_{-}(r_{0})=0$ depending on whether $S_{r_{0}}$ is a future horizon, past horizon, or both respectively. As observed in [25] the Jang equation in spherical symmetry may be reduced to a first order ODE by setting

[TABLE]

The equation (4.2) then becomes

[TABLE]

Observe that $|v|\leq 1$ and blow-up occurs precisely when $v=\pm 1$. A maximum principle type argument shows that the outermost horizon condition (4.13) ensures that $|v|<1$ away from $S_{r_{0}}$. Building upon this estimate, existence and uniqueness for the spherically symmetric Jang equation may be established following the arguments of [7, Theorem 2]; this prior results was stated for dimension three but the proof carries over to higher dimensions. The result may be stated as follows, under the hypothesis that the initial data satisfy the following fall-off conditions in the asymptotic end

[TABLE]

Theorem 4.1**.**

Assume that the initial data set is spherically symmetric, smooth, either complete or with outermost apparent horizon boundary, and satisfies the asymptotics (4.16). Then there exists a unique solution $v\in C^{\infty}((r_{0},\infty))\cap C^{1}([r_{0},\infty))$ of (4.15) (the spherically symmetric Jang equation) such that $-1<v(r)<1$, $r>r_{0}$, with $v(r_{0})$ taking the value [math] or $\pm 1$ depending on whether $r_{0}=0$ and the manifold is complete, or $S_{r_{0}}$ is a past (future) horizon, respectively. Furthermore, in the asymptotic end the decay is of the form

[TABLE]

In all cases this gives rise to a spherically symmetric solution $f\in C^{\infty}(M^{n})$ of the Jang equation (4.2) satisfying

[TABLE]

Corollary 4.2**.**

Under the hypotheses of Theorem 4.1 the Jang deformed initial data $(\bar{M}^{n},\bar{g})$ is smooth, spherically symmetric, and asymptotically flat with the same mass $\bar{m}=m$. If $(M^{n},g,k)$ has a boundary then the Jang deformation, in addition, has an asymptotically cylindrical end satisfying the asymptotics (4.11).

4.3. Uniform Asymptotics for the Solution of Jang’s Equation

The asymptotic fall-off for solutions of Jang’s equation given in the previous theorem depend on spherical symmetry and are not necessarily uniform. By allowing for a slightly weaker fall-off we may obtain uniform fall-off in the general case independent of any symmetry.

Lemma 4.3**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of uniformly asymptotically flat initial data, and let $f_{j}\in C^{\infty}(M^{n})$ be corresponding solutions of Jang’s equation with $f_{j}(x)\rightarrow 0$ as $|x|\rightarrow\infty$ where $x\in\mathbb{R}^{n}\setminus B_{\rho_{A}}$ are coordinates given by the asymptotic diffeomorphisms. Then for any small $\varsigma>0$ there exist unform constants $\mathcal{C}$ and $\bar{r}$, depending only on $\varsigma$, such that

[TABLE]

where $\bar{r}>\rho_{A}$. In particular, the sequence of Jang deformations $(\bar{M}^{n}_{j},\bar{g}_{j})$ is uniformly asymptotically flat.

Proof.

We shall adapt to our purposes an argument of Schoen-Yau which can be found in [34, pages 248-9] (see also [12, Proposition 4]). Let $r(x)=|x|$ and for $0<p<n-2$, $\lambda>0$, $r>\lambda^{\frac{1}{p+1}}$ define the radial function

[TABLE]

Observe that there is a constant $c_{1}=c_{1}(p)$ such that

[TABLE]

and

[TABLE]

A computation shows that the Jang operator evaluated at this radial function yields

[TABLE]

where $c_{2}$ is a uniform constant arising from the uniformly asymptotically flat condition. We may then choose a uniform $\lambda$ large enough to ensure that the right-hand side of (4.23) is nonpositive for $r>\lambda^{\frac{1}{p+1}}$, and thus $\bar{f}$ is a super-solution. Similarly, $-\bar{f}$ is a sub-solution on this domain. Since $f$ and $\bar{f}$ both vanish at spatial infinity, and the derivative (4.22) is infinity, a maximum principle argument guarantees that $-\bar{f}\leq f\leq\bar{f}$. Therefore

[TABLE]

From this, higher order fall-off follows by rescaling combined with the Schauder estimates as in Proposition 3 of [34]. Lastly, we may set $\bar{r}=\lambda^{\frac{1}{p+1}}$ to obtain the statement of this lemma. ∎

5. The Conformal Transformations

In this section we construct the conformal transformations and control the conformal factor as described in the introduction.

5.1. Review of Conformal Change without Symmetry

In the previous section we have obtained, from the given initial data $(M^{n},g,k)$, a Jang deformation $(\bar{M}^{n},\bar{g})$ which is complete, asymptotically flat, and with an additional asymptotically cylindrical end if the original data possessed a boundary. The positivity property (4.8) for the scalar curvature of the Jang metric leads to a stability-type inequality via integration by parts combined with Cauchy-Schwarz

[TABLE]

for all $\phi\in C^{\infty}_{c}(\bar{M}^{n})$, where $c_{n}=\frac{n-2}{4(n-1)}$. The left-hand side arises from the basic quadratic form associated with the conformal Laplacian $L_{\bar{g}}=\Delta_{\bar{g}}-c_{n}R_{\bar{g}}$, and asserts that on compact subsets this operator has nonnegative spectrum (for the Dirichlet problem). In fact the spectrum is strictly positive, since if the principal eigenvalue is zero each term on the right-hand side of (5.1) would vanish, implying that $R_{\bar{g}}=0$ and the principal eigenfunction is harmonic, which is impossible. Thus, a standard exhaustion argument together with asymptotic analysis [12, 34] shows that there is a positive solution $u>0$ of the zero scalar curvature equation

[TABLE]

for some constant $\alpha$. This allows a conformal transformation $\tilde{g}=u^{\frac{4}{n-2}}\bar{g}$ to zero scalar curvature in which the relation between the masses is given by $\tilde{m}=\bar{m}+2\alpha$ [34, page 259].

Moreover in the case that $(\bar{M}^{n},\bar{g})$ possesses an additional cylindrical end, the solution $u$ tends to zero in the limit along that end. In fact the decay along the cylindrical end is exponentially fast $u\sim e^{-\gamma s}$, where $s$ is an arclength parameter along the cylindrical end and $\gamma$ is the principal eigenvalue of $\Delta_{\partial M^{n}}-c_{n}R_{\partial M^{n}}$. In spherical symmetry additional assumptions are not required to obtain $\gamma>0$ since the scalar curvature of the outermost apparent horizon $R_{\partial M^{n}}$ is positive, although in the general case a sufficient condition is for the dominant energy condition to be strict near the horizon as was used in [34]. The next proposition records these observations.

Proposition 5.1**.**

Given a smooth Jang deformation $(\bar{M}^{n},\bar{g})$ there exists a positive solution $u\in C^{\infty}(\bar{M}^{n})$ of (5.2), so that the conformal metric $\tilde{g}=u^{\frac{4}{n-2}}\bar{g}$ has zero scalar curvature $R_{\tilde{g}}=0$ and $(\tilde{M}^{n},\tilde{g})$ is asymptotically flat with mass

[TABLE]

If an asymptotically cylindrical end is present in the Jang deformation, then the conformal factor is asymptotic to $e^{-\gamma s}$ along this end with an arclength parameter $s$. Furthermore, if the initial data are spherically symmetric then the function $u$ and hence metric $\tilde{g}$ are also spherically symmetric.

Remark 5.2**.**

Without any symmetry in dimension $n=3$, it follows from a slightly generalized positive mass inequality [34] that $\tilde{m}\geq 0$, and in addition that $\alpha\leq 0$ (see (5.8) below). Thus if $m=\bar{m}\rightarrow 0$ then $\tilde{m}\rightarrow 0$. The almost rigidity conjecture in the time-symmetric setting then suggests, modulo horizon issues, that $(\tilde{M}^{3},\tilde{g})$ is close in the intrinsic flat sense to Euclidean space.

5.2. Spherically Symmetry Gives $\tilde{g}=g_{E}$

We point out that in the spherically symmetric case the existence of a conformal transformation to zero scalar curvature may be obtained from the alternate observation that all spherically symmetric metrics are conformally flat. This is related to a rigidity phenomena associated to zero scalar curvature in spherical symmetry. The following result, which is similar to Birkhoff’s Theorem [39] in general relativity, is well-known although the authors do not know of a proper reference in the literature and thus include it here. It should be noted that a proof may be obtained from the arguments of [23], although here we give another approach.

Lemma 5.3**.**

Let $(M^{n},g)$, $n\geq 3$ be spherically symmetric, complete, scalar flat, and asymptotically flat. Then either it is isometric to flat Euclidean space $(\mathbb{R}^{n},\delta)$ or the constant time slice of a Schwarzschild spacetime, that is $M^{n}\cong\mathbb{R}^{n}\setminus\{0\}$ and

[TABLE]

Proof.

This conclusion may be derived from the observation that the Hawking mass of radial spheres is constant under the assumption of zero scalar curvature. The inverse mean curvature flow proof of the Penrose inequality [19] then guarantees the desired result where the parameter $m$ is the value of the constant Hawking mass.

An alternative proof is to directly compute the scalar curvature of the spherically symmetric metric in polar form as in (3.7), and analyze the ODE as is done in [31, page 70]. It is found that $\rho^{\prime}=1+c\rho^{2-n}$ for some constant $c<0$, the area radius function $\rho(r)>0$ for all $r\in(-\infty,\infty)$, and it has a unique minimum (corresponding to a minimal surface) at $\rho=|c|^{\frac{1}{n-2}}$. It follows that $\rho$ may be treated as a radial coordinate and so

[TABLE]

Since $\rho(r)$ has a reflection symmetry across the minimal surface, this may be doubled to obtain the stated conclusion. ∎

This rigidity result suggests that the conformal transformation obtained in Proposition 5.1, in the case of spherical symmetry, gives rise to Euclidean space as we will now see.

Corollary 5.4**.**

Let $(M^{n},g,k)$ be a spherically symmetric, asymptotically flat initial data set satisfying the dominant energy condition which is either complete or has an outermost apparent horizon boundary. Then the conformally transformed Jang deformation of Proposition 5.1 is isometric to Euclidean space $(\tilde{M}^{n},\tilde{g})\cong(\mathbb{R}^{n},g_{\mathbb{E}})$.

Proof.

We only treat the case with boundary, as the case without boundary is similar. Let $\tau$ be the radial distance function from the boundary for $(M^{n},g)$. With the help of (4.11) the Jang metric takes the form

[TABLE]

Set $s=\tfrac{2\sqrt{c}}{l}\tau^{-l/2}$ and use the expansion $\rho=\rho(0)+O(\tau)$ to obtain

[TABLE]

which illustrates the cylindrical asymptotics. Since $u\sim e^{-\gamma s}$, a straightforward computation shows that the Hawking mass (3.8) of the radial spheres $m(s)\rightarrow 0$ as $s\rightarrow\infty$. Moreover, since $R_{\tilde{g}}=0$, as in the proof of Lemma 5.3 the Hawking mass of these spheres must be constant, and hence zero. The desired conclusion now follows. ∎

Remark 5.5**.**

The results of this subsection rely heavily on the spherically symmetric assumption. Without this hypothesis, the conformally changed manifold is only known to be scalar flat rather than isometric to Euclidean space. Any progress in showing that this manifold is close to Euclidean space will most likely be obtained in conjunction with advances toward the Riemannian version of the stability conjecture.

5.3. Controlling the Conformal Factor $u$

For the remainder of this section we will examine how the mass may be used to control the conformal factor $u$ and the Jang deformation. Since $\tilde{g}$ is flat its mass vanishes $\tilde{m}=0$, and therefore the formula of Proposition 5.1 relating the masses of each deformation yields $-2\alpha=m$; where we have also used that the Jang transformation preserves mass. Furthermore, multiplying equation (5.2) through by $u$ and integrating by parts, and using the divergence structure present in $R_{\bar{g}}$ yields

[TABLE]

It follows that $L^{2}$ gradient bounds for the conformal factor are given in terms of the mass.

Lemma 5.6**.**

Assume that the given spherically symmetric initial data $(M^{n},g,k)$ satisfies the dominant energy condition. If $u$ is the solution of (5.2) given by Proposition 5.1, then

[TABLE]

Remark 5.7**.**

We remark that a version of Lemma 5.6 is likely to hold without the assumption of spherical symmetry, when a smooth Jang deformation exists. Namely, if the positive mass theorem is valid for $(\tilde{M}^{n},\tilde{g})$ then $-2\alpha\leq m$, and (5.9) again follows from (5.8). The issue is that the $n$-dimensional Riemannian positive mass theorem [33, 35] is not immediately applicable, as $(\tilde{M}^{n},\tilde{g})$ may not be a smooth manifold, geometrically or topologically. The later pathology arises from the fact that in higher dimensions, although the outermost apparent horizon is of positive Yamabe invariant, it may not be of spherical topology. Nonetheless, Eichmair [12] has found an effective way to deal with these concerns.

Next observe that the $L^{2}$ gradient bound for $u$ on the Jang surface may be translated into a similar bound for $\log u$ on Euclidean space. This follows from Corollary 5.4 since $g_{\mathbb{E}}=u^{-\frac{4}{n-2}}\bar{g}$, and in particular

[TABLE]

Corollary 5.8**.**

Under the hypotheses of Lemma 5.6

[TABLE]

The global Sobolev bounds for the conformal factor obtained from the stability inequality can be parlayed into $C^{0}$ and even Hölder estimates away from the central fixed point of the spherical symmetry. To accomplish we will first need to obtain uniform control for $u$ in the asymptotically flat end.

Lemma 5.9**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of spherically symmetric, uniformly asymptotically flat initial data satisfying the dominant energy condition and with either outermost apparent horizon boundary or no boundary. Let $(\bar{M}_{j}^{n},\bar{g}_{j})$ be the corresponding Jang deformations conformally related to $(\tilde{M}^{n}_{j},\tilde{g}_{j})$ via conformal factors $u_{j}$ solving (5.2). Then there exist uniform constants $c$ and $\bar{r}$ such that

[TABLE]

where $r=|x|$ is the radial coordinate from (1.3) and $\bar{r}$ is as in Proposition 4.3.

Proof.

For convenience, within the proof the subscript $j$ will be suppressed. Coordinate spheres in the asymptotic region will be denoted by $S_{r}$; note that they are distinct from the surfaces $\Sigma_{A}$ used in other parts of the manuscript. Let $\bar{M}^{n}_{r}$ denote the component of $\bar{M}^{n}$ lying inside $S_{r}$. Multiply equation (5.2) through by $u$ and integrate by parts up to a coordinate sphere $S_{r}$, $r\geq\bar{r}$, and use the divergence structure present in $R_{\bar{g}}$ as in (5.8) to find

[TABLE]

where $\bar{\nu}$ is the unit outer normal with respect to $\bar{g}$. Since all quantities are spherically symmetric and $r$ is arbitrary, we obtain the differential inequality

[TABLE]

According to the uniform fall-off (4.16) and (4.19), $q(\partial_{r})$ may be estimated to yield

[TABLE]

where $\bar{c}>0$ is a uniform constant. Integrating on the interval $[r,r_{1}]$ produces

[TABLE]

Now let $r_{1}\rightarrow\infty$ and use that $u(r_{1})\rightarrow 1$ to obtain

[TABLE]

with $c=\bar{c}(2n-2-2\varsigma)^{-1}$. ∎

Let $\tilde{r}$ denote the radial distance function for $g_{\mathbb{E}}=\tilde{g}$, so that

[TABLE]

Consider annular domains $\Omega_{\tilde{A}_{1},\tilde{A}_{2}}\subset\mathbb{E}^{n}$ whose boundary consists of two coordinate spheres having areas $\tilde{A}_{1}<\tilde{A}_{2}$. The next result shows that the conformal factors defining $\tilde{g}$ are uniformly close to $1$ in Hölder space away from the center of the spherical symmetry. Note that in light of Example 2.4 we see that it is not possible to have such $C^{0}$ control near the center.

Proposition 5.10**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of spherically symmetric, uniformly asymptotically flat initial data satisfying the dominant energy condition and with either outermost apparent horizon boundary or no boundary. Let $(\bar{M}_{j}^{n},\bar{g}_{j})$ be the corresponding Jang deformations conformally related to $(\tilde{M}^{n}_{j},\tilde{g}_{j})$ via conformal factors $u_{j}$ solving (5.2). Let $\tilde{A}_{0}>0$ be fixed, assume $m_{j}<\tilde{A}_{0}^{-2}$ is uniformly small, and set $\tilde{A}_{j}=\frac{1}{\sqrt{m_{j}}}$. Then there exists a uniform constant $C$ such that

[TABLE]

Proof.

For convenience the subscript $j$ will be suppressed in the proof. Estimate (5.11) shows that for any $\tilde{r}_{1}<\tilde{r}_{2}$ we have

[TABLE]

where $\tilde{c}_{n}=(n-2)^{2}c_{n}^{-1}/8$. With the help of Hölder’s inequality it follows that

[TABLE]

If

[TABLE]

denotes the area radius for the inner boundary of $\Omega_{\tilde{A}_{0},\tilde{A}_{j}}$ then

[TABLE]

for $\tilde{r}_{1}\geq\tilde{r}_{0}$, which yields one half of the desired Hölder estimate.

The next goal is to obtain $C^{0}$ bounds. Note that (5.21) implies

[TABLE]

In order to control $u(\tilde{r}_{2})$ uniformly we will utilize Lemma 5.9. The estimate there, however, is given in terms of the radial coordinate $r$ associated with the uniform asymptotic coordinates of $\bar{g}$. The two coordinates $r$ and $\tilde{r}$ may be compared in the asymptotic end by relating the volumes of coordinate spheres. Let $S_{r}$ denote a coordinate sphere of radius $r$, then uniform asymptotic flatness shows that its volume with respect to $\bar{g}$ satisfies

[TABLE]

for some uniform constant $c_{1}$ (uniform constants will be denoted by $c_{i}$, $i=1,2,\ldots$). With the help of (5.12), the volume of this same sphere computed with respect to $g_{\mathbb{E}}=u^{\frac{4}{n-2}}\bar{g}$ may be estimated by

[TABLE]

On the other hand

[TABLE]

and hence

[TABLE]

We may now combine (5.12), (5.24), and (5.28) to find for large $\tilde{r}_{2}$ that

[TABLE]

By choosing $\tilde{r}_{2}=m^{-\frac{1}{2}}$ the desired result follows. ∎

6. Intrinsic Flat Convergence of the Base Manifolds

In this section we prove the following proposition which implies (1.22) of Theorem 1.2. Recall that

[TABLE]

is the region within the level set $\Sigma^{j}_{A}$ of area $A$ with respect to $g_{j}$ and $\bar{g}_{j}$. A priori we do not know the area of the level set $\Sigma^{j}_{A}$ with respect to the metric $\tilde{g}_{j}=u_{j}^{\frac{4}{n-2}}\bar{g}_{j}=g_{\mathbb{E}}$. This section strongly uses spherical symmetry and the fact that the metrics $\bar{g}$ and $g$ are monotone as was proven in Lemma 3.4.

Proposition 6.1**.**

Given any $A>0$, $D>0$, and $\epsilon>0$ there exists $\delta=\delta(A,D,\epsilon)>0$ such that if mass $m_{j}<\delta$ then

[TABLE]

where

[TABLE]

and

[TABLE]

Furthermore for fixed $A>0$ and $D>0$, if $m_{j}\to 0$ then

[TABLE]

To prove the volume preserving intrinsic flat convergence statement (6.2), we will show they have diffeomorphic subregions $W_{j}$ and $W^{\prime}_{j}$ where the metric tensors are $C^{0}$ close and the volumes not covered by these subregions are small. This will be stated more precisely below. Based on the examples of Section 2, we cannot expect the region near the asymptotically cylindrical end of $(\Omega_{j},\bar{g_{j}})$ to be $C^{0}$ close to $(\Omega^{\prime}_{j},g_{\mathbb{E}})$. We will therefore choose an

[TABLE]

and then cut out the annular domain $\Omega^{j}_{A_{\epsilon},A}$ from inside $\Omega_{j}$ and $\Omega^{\prime}_{j}$ to create regions where we will have the appropriate control. Set

[TABLE]

and

[TABLE]

The precise choice of $A_{\epsilon}$ will be made later in Subsection 6.5.

6.1. Estimates on $W_{j}$

To begin the proof of Proposition 6.1 we apply the results of the previous section to establish uniform closeness between $\bar{g}_{j}$ and $g_{\mathbb{E}}$ on diffeomorphic domains. This is a primary step towards proving the intrinsic flat distance estimate (6.2). Secondly we show the smooth convergence to a ball in $\mathbb{E}^{n}$ (6.5).

Lemma 6.2**.**

Given any $\varepsilon>0$ and fixed $A>A_{\epsilon}>0$ there exists a $\delta=\delta(\varepsilon,A_{\epsilon})>0$ such that

[TABLE]

whenever the mass $m_{j}<\delta$.

Proof.

Choose $\tilde{A}_{0}<A_{\epsilon}$ and $\tilde{A}>A$. Then, according to Proposition 5.10, there is a uniform constant $C$ such that

[TABLE]

where $\Omega_{\tilde{A}_{0},\tilde{A}}\subset\mathbb{E}^{n}$ is the annular domain whose boundary consists of two coordinate spheres having areas $\tilde{A}_{0}<\tilde{A}$. Thus, for all $m_{j}$ sufficiently small we find that $\Omega_{A_{\epsilon},A}^{j}\subset\Omega_{\tilde{A}_{0},\tilde{A}}$ and

[TABLE]

from which the desired result follows. ∎

Lemma 6.3**.**

Fix $A>0$ and $D>0$, and let the symmetry spheres $\Sigma^{j}_{A}\subset\bar{M}_{j}^{n}$ be defined by

[TABLE]

If $m_{j}\to 0$ then

[TABLE]

and so the spheres

[TABLE]

and the balls

[TABLE]

Proof.

Observe that (6.13) follows from the $C^{0}$ convergence of the metric tensors in Lemma 6.2. Moreover (6.14) and (6.15) follow immediately from (6.13). ∎

6.2. Applying the Lakzian-Sormani Theorem

In work of Lakzian and the third author [22], the following theorem was proven providing a concrete means to estimate the intrinsic flat distance. Intuitively this theorem observes that two manifolds $\Omega_{j}$ and $\Omega^{\prime}_{j}$, as in (6.3) and (6.4), are close in the intrinsic flat sense if they have diffeomorphic subregions $W_{j}$ and $W^{\prime}_{j}$, as in (6.7) and (6.8) where the metric tensors are $C^{0}$ close and the volumes not covered by these subregions are small. See Figure 8. When using this result for the current problem note that a dictionary between the notation is $\Omega=\Omega_{j}$, $\Omega^{\prime}=\Omega^{\prime}_{j}$, $W=W_{j}$, $W^{\prime}=W^{\prime}_{j}$, $g=\bar{g}_{j}$, and $g^{\prime}=g_{\mathbb{E}}$.

Theorem 6.4**.**

[22]* Suppose that $(\Omega,g)$ and $(\Omega^{\prime},g^{\prime})$ are oriented precompact Riemannian manifolds with diffeomorphic subregions $W\subset\Omega$ and $W^{\prime}\subset\Omega^{\prime}$. Identifying $W=W^{\prime}$ assume that*

[TABLE]

Taking the extrinsic diameters

[TABLE]

define a hemispherical width

[TABLE]

Taking the difference in distances with respect to the outside manifolds, set

[TABLE]

and define the height

[TABLE]

Then

[TABLE]

6.3. Estimating the Volumes

In this subsection all volumes and areas appearing in the first line of (6.21) will be shown to be uniformly bounded, and those in the second line will be shown to be arbitrarily small. In addition, it will be established that the volumes $\operatorname{Vol}_{\bar{g}_{j}}(\Omega_{j})$ converges to that of a ball in Euclidean space.

Lemma 6.5**.**

Volumes outside the diffeomorphic subregions may be estimated by

[TABLE]

and

[TABLE]

In particular, both are less than $2^{n-1}DA_{\epsilon}$.

Proof.

First observe that

[TABLE]

By Lemma 3.4, the largest area of the radial levels in this set is $A_{\epsilon}$. Applying the coarea formula to these levels and the fact that the depth of the set is $D$, we have

[TABLE]

Applying the coarea formula again and using (6.16) yields

[TABLE]

∎

Lemma 6.6**.**

Volumes of the diffeomorphic subregions may be estimated by

[TABLE]

and

[TABLE]

In particular, both are less than $2^{n-1}DA$.

Proof.

This follows from Lemma 3.4 and the coarea formula

[TABLE]

Furthermore

[TABLE]

∎

Lemma 6.7**.**

Boundary areas of the diffeomorphic subregions may be estimated by

[TABLE]

and

[TABLE]

In particular, both are less than $2^{n}A$.

Proof.

We know that $\partial W_{j}$ has at most two components, both of which are radial levels of area less than $A$ by monotonicity. The same holds for $\partial W^{\prime}_{j}$, except that as above the upper bound on the outer area is $(1+\varepsilon)^{n-1}A$. ∎

The next lemma will be used to obtain volume preserving intrinsic flat convergence. It follows from the last few lemmas.

Lemma 6.8**.**

The difference of total volumes may be estimated by

[TABLE]

Proof.

From Lemma 6.5 we have

[TABLE]

and

[TABLE]

Moreover by (6.16)

[TABLE]

∎

6.4. Estimating Distances and Diameters

Lemma 6.9**.**

Let $D\geq\rho_{A}$, then

[TABLE]

where $D_{0}\leq 4\pi D$.

Proof.

This follows because the depth of the tubular neighborhood is $D$, and the largest symmetry sphere satisfies

[TABLE]

So by the triangle inequality the diameter of $\Omega_{j}$ is no larger than $2D+\pi D\leq 4\pi D$. In addition

[TABLE]

So by the triangle inequality the diameter of $\Omega^{\prime}_{j}$ is no larger than $2D+(1+\varepsilon)\pi D\leq 4\pi D$. ∎

Lemma 6.10**.**

The difference of distances satisfies

[TABLE]

where $W_{j}$ is identified with $W^{\prime}_{j}$.

Proof.

Let $\sigma=\sigma_{x,y}$ be a line segment from $\sigma(0)=x$ to $\sigma(1)=y$ in $\Omega^{\prime}_{j}$ so that

[TABLE]

Then there is a first time that $\sigma(t_{1})\in\Sigma_{A_{\epsilon}}^{j}$ and a last time that $\sigma(t_{2})\in\Sigma_{A_{\epsilon}}^{j}$. See Figure 9.

By triangle inequality we have

[TABLE]

Here we have used Lemma 6.9 and the fact that $\Sigma_{A_{\epsilon}}^{j}$ is connected.

Let $\gamma=\gamma_{x,y}$ be a $\bar{g}_{j}$-length minimizing geodesic from $\gamma(0)=x$ to $\gamma(1)=y$ in $\Omega_{j}$ so that

[TABLE]

Then there is a first time that $\gamma(t_{1})\in\Sigma_{A_{\epsilon}}^{j}$ and a last time that $\gamma(t_{2})\in\Sigma_{A_{\epsilon}}^{j}$. By the triangle inequality we have

[TABLE]

Combining these with $D_{0}\leq 4\pi D$ from Lemma 6.9 we obtain the desired result. ∎

6.5. Proof of Proposition 6.1

First note that by Lemmas 6.5-6.9 applied to Theorem 6.4 we have

[TABLE]

where

[TABLE]

Furthermore by Lemma 6.8

[TABLE]

Given $A>0$, $D>0$, and $\epsilon>0$ we must choose $\delta=\delta(A,D,\epsilon)>0$ sufficiently small so that

[TABLE]

whenever $\delta<m_{j}$. Therefore choose $\varepsilon$ and $A_{\epsilon}$ appropriately in order to satisfy

[TABLE]

[TABLE]

and

[TABLE]

where

[TABLE]

Observe that these together imply (6.48). It remains only to show that (6.51) may be obtained from choosing $\varepsilon$ and $A_{\epsilon}$ small.

According to the definition of $\Lambda_{j}$ and $\omega$, (6.51) will be valid if

[TABLE]

[TABLE]

and

[TABLE]

Moreover, from Lemma 6.10 it follows that (6.54) holds if

[TABLE]

and

[TABLE]

Thus we choose the area of the inner level set $A_{\epsilon}=A_{\epsilon}(A,D)>0$ small enough so that (6.49), and (6.57) hold. Recall that $\rho_{A_{\epsilon}}$ is the radius of the $n-1$-sphere whose area is $A_{\epsilon}$. In addition, we choose $\varepsilon=\varepsilon(A,D,\epsilon)$ small enough so that (6.50), (6.53), (6.55), and (6.56) hold. It follows that there exists $\delta=\delta(\varepsilon,A_{\epsilon})>0$ such that $m_{j}<\delta$ implies $d_{\mathcal{VF}}(\Omega_{j},\Omega^{\prime}_{j})<\epsilon$. Finally (6.5) was established in Lemma 6.3.

7. Convergence of the Second Fundamental Form Differences

In Theorem 1.2 a statement (1.23) is given concerning the convergence of the difference of the second fundamental forms of the initial data and Jang graphs. In this section we give a proof of that statement.

7.1. Outside Control

The pointwise control on the conformal factor gives rise to $L^{2}$ control of the difference of second fundamental forms. Recall that $\Omega_{A}$ is the domain inside the symmetry sphere of volume $A$ in $(\bar{M}^{n},\bar{g})$ and $(M^{n},g)$. By combining Proposition 5.10 and (5.8) we obtain the following.

Proposition 7.1**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of spherically symmetric, uniformly asymptotically flat initial data satisfying the dominant energy condition and with either outermost apparent horizon boundary or no boundary. For every $A_{0}>0$ there exists a uniform constant $C$ such that

[TABLE]

Proof.

According to Proposition 5.10 and the asymptotics of $u_{j}$, we find by choosing $0<\tilde{A}_{0}<A_{0}$ that

[TABLE]

for some uniform constant $c_{1}$. Thus for large $j$, the inequality (5.8) implies

[TABLE]

∎

7.2. Inside Control

The conformal factors $u_{j}$ decay exponentially fast along the cylindrical ends of $(\bar{M}^{n}_{j},\bar{g}_{j})$, and thus the techniques used above to obtain estimates outside a fixed level surface do not extend down the cylindrical end. Nevertheless, with an extra hypothesis of uniform $L^{2}$ control for the initial data second fundamental forms $k_{j}$, we are able to obtain uniform $L^{2}$ bounds for the difference of second fundamental forms on tubular neighborhoods of the anchor surface. It may be possible to weaken this hypothesis, however we expect that some condition on $k_{j}$ is necessary in order to control the difference of second fundamental forms.

Proposition 7.2**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of spherically symmetric, uniformly asymptotically flat initial data satisfying the dominant energy condition and with either outermost apparent horizon boundary or no boundary. Further, assume the uniform bound $\parallel k_{j}\parallel_{L^{2}(M_{j}^{n})}\leq B$. Given $A>0$ and $D>0$ there exists $C$ depending only on $A$, $B$, and $D$ such that

[TABLE]

Proof.

For convenience in this proof we will drop the subscript $j$. Recall from (5.1) that

[TABLE]

for all $\phi\in C^{\infty}_{c}(\bar{M}^{n})$. Choose a particular test function such that $0\leq\phi\leq 1$, $\phi\equiv 1$ on $T_{D+1}(\Sigma_{A})$, and $|\nabla\phi|_{\bar{g}}\leq 2$ then

[TABLE]

Write

[TABLE]

such that the interval $0<\bar{s}<D+2$ covers $T_{D+1}(\Sigma_{A})$, and notice that

[TABLE]

Integrating by parts produces

[TABLE]

Since the mean curvature of coordinate spheres is $\bar{H}=(n-1)\rho^{-1}\dot{\rho}$, this may be rewritten as

[TABLE]

We now estimate each term of (7.10) separately. By uniform asymptotic flatness of $\bar{g}$ in Lemma 4.3 and asymptotic control (7.2) of $u$, we have that

[TABLE]

is uniformly bounded and thus

[TABLE]

Next according to Lemma 6.8 it holds that $\text{Vol}_{\bar{g}}(T_{D+1}(\Sigma))$ is uniformly bounded, and hence by Hölder’s inequality

[TABLE]

The remaining term of (7.10) may also be estimated by the $L^{2}$ norm of $\bar{H}$.

To complete the proof we will use the extra hypothesis concerning $k$ to bound $\bar{H}$ in $L^{2}$. First observe that

[TABLE]

so that

[TABLE]

where $H$ is the mean curvature with respect to $g$. Furthermore, let $\psi:T_{D+1}(\Sigma_{A})\rightarrow M^{n}$ be the diffeomorphism onto its image associated with the coordinate change $\bar{s}\rightarrow s$. Then

[TABLE]

where we have used that the spacetime Hawking mass

[TABLE]

is nonnegative [21], and the areas of coordinate spheres agree with respect to both metrics $g$ and $\bar{g}$, that is $A(s)=A(\bar{s})$. Again using that $\rho(D+2)$ is uniformly bounded produces

[TABLE]

from which the desired result follows. ∎

7.3. Global Control

Here we show that the difference of second fundamental forms tends to zero in $L^{p}$, $1\leq p<2$ on the tubular neighborhoods $\Omega_{A}\cap T_{D}(\Sigma_{A})\subset\bar{M}^{n}$, when the masses go to zero. This result, and in fact even a uniform bound, cannot be extended to the entire domain $\Omega_{A}$. The reason is due to the cylindrical asymptotics contained within this domain, and the observation that in the limit along the cylinder $h-k$ converges to $II-i^{*}k$, where $II$ is the second fundamental form of $\partial M^{n}$ in $M^{n}$ and $i:\partial M^{n}\hookrightarrow M^{n}$ is inclusion. Since $II-i^{*}k$ does not necessarily vanish, it follows that $h-k$ is not necessarily in $L^{p}(\Omega_{A})$.

Theorem 7.3**.**

Let $(M^{n}_{j},g_{j},k_{j})$ be a sequence of spherically symmetric, uniformly asymptotically flat initial data satisfying the dominant energy condition and with either outermost apparent horizon boundary or no boundary. Further, assume the uniform bound $\parallel k_{j}\parallel_{L^{2}(M_{j}^{n})}\leq B$ and that $m_{j}\rightarrow 0$. Then for any $1\leq p<2$, $A>0$, and $D>0$

[TABLE]

Proof.

Let $A_{\varepsilon}>0$ be fixed. Then by Hölder’s inequality

[TABLE]

Using the notation and results of Lemma 6.5, observe that

[TABLE]

Furthermore by Lemma 6.8

[TABLE]

by Proposition 7.1

[TABLE]

and by Proposition 7.2

[TABLE]

It follows that

[TABLE]

This can be made arbitrarily small by choosing $A_{\varepsilon}$ small and $j$ appropriately large. ∎

8. Proof of Theorem 1.2

The proof of the main theorem follows quickly from Proposition 6.1 and Theorem 7.3. Recall that

[TABLE]

and

[TABLE]

By the triangle inequality

[TABLE]

As $m_{j}\rightarrow 0$ we have that (6.2) produces

[TABLE]

and (6.5) yields

[TABLE]

since smooth convergence implies convergence in the volume preserving intrinsic flat sense. Finally, Theorem 7.3 gives (1.23).