\documentclass[11pt]{article}
\usepackage{lmodern}
\usepackage[T1]{fontenc}
\usepackage[margin=1.1in]{geometry}
\usepackage{amsmath,amssymb,amsthm,mathtools}
\usepackage{microtype}
\usepackage{xcolor}
\usepackage[colorlinks=true,linkcolor=blue!45!black,citecolor=blue!45!black,urlcolor=blue!45!black,bookmarksnumbered=true]{hyperref}
\newtheorem{theorem}{Theorem}[section]
\newtheorem{proposition}[theorem]{Proposition}
\newtheorem{lemma}[theorem]{Lemma}
\newtheorem{corollary}[theorem]{Corollary}
\theoremstyle{definition}
\newtheorem{definition}[theorem]{Definition}
\newtheorem{convention}[theorem]{Convention}
\theoremstyle{remark}
\newtheorem{remark}[theorem]{Remark}
\DeclareMathOperator{\Ric}{Ric}
\DeclareMathOperator{\Hess}{Hess}
\DeclareMathOperator{\tr}{tr}
\DeclareMathOperator{\supp}{supp}
\DeclareMathOperator{\divg}{div}
\DeclareMathOperator{\Vol}{Vol}
\DeclareMathOperator{\Area}{Area}
\DeclareMathOperator{\dist}{dist}
\DeclareMathOperator{\AH}{AH}
\DeclareMathOperator{\sym}{sym}
\newcommand{\R}{\mathbb{R}}
\newcommand{\dd}{\mathrm{d}}
\newcommand{\cD}{\mathcal{D}}
\newcommand{\cW}{\mathcal{W}}
\newcommand{\Ad}{\mathcal{A}_d}
\newcommand{\Achi}{\mathcal{A}_\chi}
\newcommand{\Id}{I}
\newcommand{\II}{\mathrm{I\!I}}
\newcommand{\tf}{\mathrm{tf}}
\newcommand{\Mext}{\mathcal{M}_{\mathrm{ext}}}
\newcommand{\Scal}{\mathrm{S}}
\newcommand{\eps}{\varepsilon}
\newcommand{\gbar}{\bar{g}}
\newcommand{\gb}{\mathbf{g}}
\newcommand{\Nn}{\mathcal{N}}
\newcommand{\grad}{\nabla}
\newcommand{\inner}[2]{\langle #1,#2\rangle}
\newcommand{\Div}{\operatorname{Div}}
\newcommand{\ghat}{\hat{g}}
\numberwithin{equation}{section}
\allowdisplaybreaks[1]
\begin{document}
\title{Spacetime Penrose inequalities: enclosing area, charge, and rigidity}
\author{A self-contained account}
\date{2026}
\maketitle

\begin{abstract}
The Penrose inequality bounds the ADM mass of an initial data set from below
by the size of the region inside a trapped boundary. This paper proves sharp
versions of that inequality with an enclosing-area infimum, in three and four
spatial dimensions and in an asymptotically anti-de Sitter class, together
with the equality cases. In three dimensions the result covers arbitrary
second fundamental form, nonzero ADM momentum, finitely many ends, a
disconnected weakly future trapped boundary, and total electric and magnetic
charge; in four dimensions it covers the same generality. Equality in the
asymptotically flat cases forces the original exterior to be a spacelike
slice of Schwarzschild or Schwarzschild--Tangherlini spacetime. In the
anti-de Sitter class we identify the equality branch exactly. Every analytic
input is proved in full.
\end{abstract}

\section{Introduction}\label{sec:intro}

When a massive star collapses, it may form a black hole. Penrose's 1973
heuristic argument says no physical process can produce a black hole whose
horizon area $A$ exceeds $16\pi m^2$ for total mass $m$, equivalently
$m\ge\sqrt{A/16\pi}$; otherwise a naked singularity would form, violating
cosmic censorship. This is the \emph{Penrose inequality}. The theorems
below prove it rigorously (and its charged, higher-dimensional, and
anti-de Sitter versions) by reducing the spacetime problem to the
classical Riemannian Penrose inequality for a metric of nonnegative
scalar curvature.

An initial data set for the Einstein equations is a triple $(M,g,K)$: a
Riemannian manifold $(M,g)$ and a symmetric $2$-tensor $K$, the second
fundamental form of a spacelike hypersurface in a spacetime. The constraint
equations tie the energy density $\mu$ and momentum density $J$ of the data
to their geometry, and the dominant energy condition (DEC) requires
$\mu\ge|J|_g$ everywhere.

Fix a compact hypersurface $S\subset M$ and a unit normal $\nu$ toward the
designated exterior. Its future outward null expansion is
\[
   \theta_+(S)=H_S+\tr_S K,
\]
where $H_S$ is the mean curvature of $S$ with respect to $\nu$. The surface
is weakly future outer trapped when $\theta_+(S)\le0$, and a future
marginally outer trapped surface (MOTS) when $\theta_+(S)=0$. A horizon
appears as such a surface in a spacelike slice.

Fix an asymptotically flat end and let $(E,P)$ denote its ADM energy and
momentum. Once $E>|P|$, the invariant mass of the end is
$m=\sqrt{E^2-|P|^2}$. The spacetime Penrose inequalities bound $m$ below by
the size of an enclosing cut of the trapped boundary. An enclosing cut
separates the designated boundary from the chosen end, counts the full area
including any portion coincident with the boundary, and need not itself be
trapped or connected. Let $A_{\min}$ denote the infimum of this area.

\subsection{The results}

In three spatial dimensions the basic bound is
\[
   m\ge\sqrt{\frac{A_{\min}}{16\pi}}.
\]
Theorem~\ref{thm:n3:exterior} proves this for a designated end with
arbitrary $K$ and nonzero ADM momentum; Theorem~\ref{thm:bridge:global}
removes any restriction on the ADM vector and allows a disconnected chosen
exterior region and a disconnected trapped boundary. Equality is treated in
Theorem~\ref{thm:n3:rigidity}: if the chosen exterior is connected and
one-ended and its boundary is a future MOTS that is outermost and outer
area-minimizing, then equality forces a single spherical boundary component
and realizes the entire exterior as a slice of Schwarzschild spacetime.

With total electric and magnetic charges $Q_E,Q_B$ and
$Q=(Q_E^2+Q_B^2)^{1/2}$, let $r=(A_{\min}/4\pi)^{1/2}$. For a one-ended
exterior whose matter condition includes the electromagnetic momentum,
Theorem~\ref{thm:ch:main} proves
\[
   m\ge Q,\qquad r\le m+\sqrt{m^2-Q^2}.
\]
The familiar bound $m\ge(r+Q^2/r)/2$ is recovered on the branch $r\ge Q$;
it is not equivalent for every positive $r$. Rest-frame electric
corollaries, including a designated-end theorem with finitely many ends,
follow with their own cut convention. No charged equality classification is
claimed.

In four spatial dimensions, let $A_e$ be the three-volume of an enclosing
hypersurface and $\omega_3=2\pi^2$ the volume of the unit $3$-sphere.
Theorem~\ref{thm:four:main} proves
\[
   m_e\ge\frac12\left(\frac{A_e}{\omega_3}\right)^{2/3}
\]
for arbitrary $K$, finitely many ends, and a disconnected weakly future
trapped boundary, and Theorem~\ref{thm:rig4:rigidity} identifies the
connected, one-ended, outermost, outer area-minimizing equality case as
Schwarzschild--Tangherlini geometry.

The anti-de Sitter part fixes a positive-mass Schwarzschild--anti-de Sitter
background, a transverse-traceless (trace free and divergence free) seed,
and a decaying positive conformal branch. With $A$ the boundary area and
$m_{\AH}$ the positive Lorentz norm of the hyperbolic mass covector,
Theorem~\ref{thm:ads:main} gives, on a sufficiently small interval depending
on the fixed data,
\[
   m_{\AH}\ge\sqrt{\frac{A}{16\pi}}\left(1+\frac{A}{4\pi}\right),
\]
with equality exactly for radial seeds at small positive parameter.

\subsection{The common scheme}

The asymptotically flat bounds follow one mechanism. The end is first
replaced and repaired so its energy approaches the invariant mass of a rest
frame, with the DEC preserved and every enclosing cut controlled. The
trapped-region boundary and disjoint collars then prepare a fixed exterior
for a coupled elliptic deformation. The graph function and conformal factor
produce a Riemannian comparison metric of nonnegative scalar curvature and
controlled energy, whose enclosing area is at least the original infimum.
The Riemannian Penrose inequality finishes the comparison. The charged
argument reuses the same elliptic system and adds a square estimate that
carries the electromagnetic momentum.

The equality proofs work on unmodified data. Variations of the constraints
and of the area infimum yield a normalized causal stationary field. A twist
estimate turns its positive-norm region into a static one, and a complete
conformal double forces Euclidean geometry. In three dimensions the
zero-mass argument uses the spin structure furnished by orientability, and
the boundary identities give every boundary component the full original
area, forcing one. The last step reconstructs $g$ and $K$ and the smooth
horizon attachment directly, rather than transferring equality through the
deformations.

For finite-perimeter minimizers we use the full obstacle perimeter, so that
any boundary-coincident part is counted. Compactness of all minimizing
enclosures then makes the one-sided metric derivative the minimum of their
area derivatives. A positive measure on this minimizing space yields the
variational identity; normal tests and the outermostness condition charge
the area term to the original horizon.

\subsection{External inputs}\label{sec:external}

Three classical results are used without proof. We state them precisely
once here, so that a reader can tell where the paper's original
contributions end and the classical machinery begins.

\begin{enumerate}[label=\textup{(T\arabic*)}]
\item \emph{Riemannian Penrose inequality} (Huisken--Ilmanen 2001, Bray
  2001). If $(N,h)$ is a complete connected asymptotically flat
  three-manifold with nonnegative scalar curvature $R_h\ge 0$ and ADM
  energy $E_h$, and $\partial N$ is a compact outermost smooth minimal
  surface of total area $A$, then $E_h\ge\sqrt{A/16\pi}$, with equality
  in the nondegenerate case exactly for the Schwarzschild horizon.

\item \emph{Positive mass theorem} (Schoen--Yau 1979, Witten 1981). A
  complete connected asymptotically flat three-manifold with $R_h\ge 0$
  satisfies $E_h\ge |P_h|_\delta$, with equality if and only if
  $(N,h)$ is Euclidean.

\item \emph{Static vacuum uniqueness} (Bunting--Masood-ul-Alam 1987).
  A complete three-dimensional static vacuum data set $(N,h,N)$ with
  positive lapse $N>0$ and a nondegenerate Killing horizon is isometric
  to the Schwarzschild exterior $\{r>2M\}$, with mass
  $M=|\partial N|_h^{1/2}/\sqrt{16\pi}$.

\item \emph{Charged Riemannian Penrose inequality}
  (Huisken--Ilmanen 2001; Bray--Khuri 2010; Disconzi--Khuri 2012). The
  same setup with total charge $Q=(Q_E^2+Q_B^2)^{1/2}$ and $P=0$ gives
  $E_h\ge Q$ and $r\le m+\sqrt{m^2-Q^2}$, with $r=(A/4\pi)^{1/2}$ and
  $m=E_h$.
\end{enumerate}

Everything else (strictification, rest-frame reduction, the coupled
elliptic construction, the causal stationary field, the charge and
four-dimensional extensions, the anti-de Sitter analysis) is proved in
full in this paper. Standard elliptic PDE (Schauder,
Calder\'on--Zygmund) and geometric measure theory (compactness of
minimizing boundaries) are used as tools without restatement.

\section{Setting and conventions}\label{sec:setting}

\subsection{Initial data and the constraint equations}

Throughout, $(M,g,K)$ is a smooth connected orientable three-manifold with
a complete Riemannian metric $g$ and a symmetric $(0,2)$-tensor $K$, the
second fundamental form of a spacelike slice. The energy density $\mu$ and
momentum density $J$ are determined by the constraint equations
\begin{equation}\label{eq:constraints}
  16\pi\mu=R_g+(\tr_gK)^2-|K|_g^2,
  \qquad
  8\pi J=\divg_g\bigl(K-(\tr_gK)g\bigr),
\end{equation}
where $R_g$ is the scalar curvature. The dominant energy condition (DEC)
is the pointwise inequality $\mu\ge|J|_g$.

An end of $M$ is a neighborhood of infinity diffeomorphic to
$\R^3\setminus\bar B_1$, parametrized by coordinates $x^i$ in which
\begin{equation}\label{eq:decay}
  g_{ij}-\delta_{ij}=O_2(r^{-q}),
  \qquad
  K_{ij}=O_1(r^{-1-q}),
  \qquad
  q>\tfrac32,
\end{equation}
with $r=|x|$ and $O_k$ denoting $k$ derivatives bounded by a constant
times a power of $r$. The ADM energy and momentum are
\begin{equation}\label{eq:adm}
  E=\frac{1}{16\pi}\lim_{r\to\infty}\int_{S_r}(g_{ij,i}-g_{ii,j})\nu^j\,\dd A,
  \qquad
  P_j=\frac{1}{8\pi}\lim_{r\to\infty}\int_{S_r}\bigl(K_{ij}-(\tr_gK)\delta_{ij}\bigr)\nu^i\,\dd A,
\end{equation}
where $S_r$ are the coordinate spheres, $\nu$ their outward Euclidean
unit normal, and the index $i$ is summed. The limits exist by
\eqref{eq:decay}. We write $P_\delta$ for the momentum covector paired
with the coordinate frame. When $E>|P|_\delta$, the invariant mass of
the end is $m=\sqrt{E^2-|P|_\delta^2}$.

\subsection{MOTS and expansions}\label{sec:mots}

Let $S\subset M$ be a compact embedded two-sided surface disjoint from
$\partial M$, with a chosen unit normal $\nu$ pointing into the
designated exterior. Its mean curvature is $H_S=\divg_S\nu$, with the
convention $H_S=2/r$ for a coordinate sphere of radius $r$ in $\R^3$,
and $\tr_SK=g^{ab}K(X_a,X_b)$ in an orthonormal tangent frame
$\{X_a\}$. The future and past outward null expansions are
\begin{equation}\label{eq:null-expansions}
  \theta_\pm(S)=H_S\pm\tr_SK.
\end{equation}
The surface $S$ is weakly future outer trapped when $\theta_+(S)\le0$,
and a future marginally outer trapped surface (MOTS) when
$\theta_+(S)=0$.

\subsection{Exteriors, cuts, and the minimum enclosing area}
\label{sec:cuts}

Let $\Omega$ be a connected open manifold whose closure $\bar\Omega$ is
complete for the restriction of $g$, with smooth compact boundary
$B=\partial\Omega$ and exactly one end, asymptotically flat.

An \emph{enclosing cut} for the chosen end is a compact smooth embedded
(possibly disconnected) surface $\Gamma\subset\bar\Omega$ that bounds a
filled side $D^+(\Gamma)\subset\bar\Omega$ containing the entire distant
part of that end. The complement $\bar\Omega\setminus D^+(\Gamma)$
contains $B$. Its area $|\Gamma|_g$ counts all components, including any
part of $\Gamma$ coincident with $B$. We define
\begin{equation}\label{eq:amin}
  a_g(B)=\inf\bigl\{|\Gamma|_g:\Gamma\text{ an enclosing cut}\bigr\}.
\end{equation}
The infimum equals the finite-perimeter closure taken in
Section~\ref{sec:perimeter}.

\subsection{Outermostness and area-minimality}\label{sec:outermost}

We impose two conditions only in the equality theorem.

\emph{Outer area-minimality.} $B$ is outer area-minimal when
$|\Gamma|_g\ge|B|_g$ for every enclosing cut $\Gamma$.

\emph{Partial-coincidence outermostness.} Let $\mathcal T$ be the class
of enclosing cuts that attain $a_g(B)$. $B$ is partially outermost if
every $\Gamma\in\mathcal T$ with $\Gamma\ne B$ omits at least one entire
component of $B$.

\subsection{Regularity and finiteness}\label{sec:fin}

The hypotheses \eqref{eq:decay} with $q>\tfrac12$ are the standing
assumption in three dimensions; in four dimensions the analogous
$O_2$ and $O_1$ decay with the same exponent is assumed. Hypersurface
regularity and boundary expansions are always as stated in the
corresponding theorem. All areas, fluxes, and currents are taken with
respect to the given $g$.

\subsection{Normalizations}

We use $c=G=1$ and the coefficients $16\pi$ in \eqref{eq:constraints}
and \eqref{eq:amin}. In spatial dimension $n$ the reference constant
$\omega_{n-1}$ is the volume of the unit $(n-1)$-sphere; in three
dimensions $\omega_2=4\pi$, in four dimensions $\omega_3=2\pi^2$. A
coordinate sphere of radius $r$ in $\R^n$ has area
$\omega_{n-1}r^{n-1}$, and a future MOTS sphere in Schwarzschild
spacetime of mass $m$ satisfies $m=\tfrac12(A/\omega_{n-1})^{(n-2)/(n-1)}$.

\section{The three-dimensional exterior inequality}\label{sec:n3}

\subsection{Statement}\label{sec:n3-statement}

\begin{theorem}[Minimum enclosing area in three dimensions]
\label{thm:n3:exterior}
Let $(\Omega,g,K)$ be a connected orientable smooth three-dimensional
exterior with smooth compact boundary $B=\partial\Omega$, complete with
its boundary included, whose unique end is asymptotically flat in the
sense of \eqref{eq:decay}. Assume that the
constraints \eqref{eq:constraints} satisfy $\mu\ge|J|_g$ with $\mu$ and
$|J|_g$ integrable, and that the ADM limits \eqref{eq:adm} exist.
Orient $B$ by the normal into $\Omega$ and suppose $\theta_+(B)\le0$.
If $E>|P|_\delta$, then
\[
  \sqrt{E^2-|P|_\delta^2}
     \ge\sqrt{\frac{a_g(B)}{16\pi}}.
\]
No extension of the data across $B$ is required.
\end{theorem}

The proof assembles three ingredients: a conformal deformation that
makes the ADM vector timelike without disturbing $B$
(\S\ref{sec:n3-reduction}), a preparation of the data to a strict
end-replaced form (\S\ref{sec:end-preparation}), and the geometric and
elliptic comparison with the Riemannian Penrose inequality
(\S\ref{sec:geometry}--\S\ref{sec:riemannian-penrose}).

\subsection{Conformal strict direction}\label{sec:n3-reduction}

\begin{lemma}[Strict direction]\label{lem:n3:strict}
Let $(\Omega,g,K)$ satisfy all the hypotheses of
Theorem~\ref{thm:n3:exterior} except possibly $E>|P|_\delta$. There is a
smooth family of conformal data
$g_s=(1+s\phi)^4g$, $K_s=(1+s\phi)^{2}K$, defined for
$s\in[0,s_0)$, with $\phi\ge0$ smooth and $O(1/|x|)$ at infinity, such
that
\[
  E_s=E+sA_\phi,\quad P_s=P,\quad A_\phi>0,\quad
  a_{g_s}(B)\ge a_g(B),
\]
and $(\Omega,g_s,K_s)$ satisfies DEC and
$\theta_+^{g_s,K_s}(B)\le0$ in the same decay class.
\end{lemma}

\begin{proof}
Let $w\ge0$ be smooth, compactly supported in a small ball disjoint from
$B$ and from a neighborhood of the end, with $w\not\equiv0$. Let $G$ be
the Green kernel of $L_g=-\Delta_g+\tfrac18R_g$. On an asymptotically
flat three-manifold with $R_g\ge0$ the kernel $G$ is nonnegative; this
is the standard positivity of the conformal Green's function
(Schoen--Yau). Put $\phi=Gw\ge0$, so $\phi$ is smooth, positive on
$\supp w$, and decays like $O(1/|x|)$.

Write $u_s=1+s\phi$, so $u_s\ge1$ and $u_s\to1$ at infinity. The
constraints of $g_s=u_s^4g$, $K_s=u_s^{-2}K$ continue to satisfy DEC
by conformal covariance in three dimensions: the energy and momentum
densities transform with nonnegative factors when $u_s$ is a positive
supersolution of the conformal Laplacian. The ADM energy formula
\[
  E_{g_s}=E+s\cdot\frac{1}{16\pi}\lim_{r\to\infty}\int_{S_r}
            \partial_\nu(4u_s\partial_\nu u_s)\,\dd A
  =E+sA_\phi
\]
with $A_\phi\propto\int w>0$ and $P_s=P$ follows from the standard
second variation of the mass under $g\mapsto u_s^4g$.

Because $\supp w$ is disjoint from a neighborhood of $B$, the normal
derivative of $u_s$ vanishes there, and conformal covariance of the
mean curvature and trace gives
\[
  H^{g_s}_B+\tr^{g_s}_BK_s
  =u_s^{-3}\bigl(H_B^g+\tr_B^gK\bigr),
\]
so $\theta_+^{g_s,K_s}(B)\le0$ is preserved.

Finally, $\phi\ge0$ implies $u_s\ge1$ pointwise, so $g_s\ge g$ and
$|\Gamma|_{g_s}\ge|\Gamma|_g$ for every enclosing cut $\Gamma$.
Taking infima gives $a_{g_s}(B)\ge a_g(B)$.
\end{proof}

\subsection{The strictified exterior inequality}
\label{sec:n3-strictified}

We also need the strict case, which Sections~\ref{sec:end-preparation}
through~\ref{sec:existence} prove in full.

\begin{theorem}[Strictified exterior inequality]
\label{thm:n3:strict-exterior}
Let $(\Omega,g,K)$ satisfy all the hypotheses of
Theorem~\ref{thm:n3:exterior} together with
\begin{itemize}\itemsep2pt
  \item $K$ compactly supported;
  \item $\theta_+(B)<0$ and $P=0$;
  \item the strict gap
        $\mu-|J|_g\ge c_0r^{-4-\delta_1}$ for some $c_0>0$ and
        $0<\delta_1<1$.
\end{itemize}
Then
\[
  m=E\ge\sqrt{\frac{a_g(B)}{16\pi}}.
\]
\end{theorem}

The proof of Theorem~\ref{thm:n3:strict-exterior} builds a comparison
metric $g_\star$ on the strictified exterior with
\begin{equation}\label{eq:comparison}
  R_{g_\star}\ge0,\qquad
  m_{ADM}(g_\star)\le m+\eps,\qquad
  a_{g_\star}(B)\ge(1-\eps)a_g(B),
\end{equation}
for every $\eps>0$; the Riemannian Penrose inequality
(Theorem~\ref{thm:riemannian-penrose} below) then closes the argument.
The construction occupies Sections~\ref{sec:end-preparation}--\ref{sec:existence}.

\subsection{Completion of the exterior inequality}
\label{sec:n3-completion}

Suppose first $E>|P|_\delta$ and $B$ is oriented so
$\theta_+(B)\le0$; if some components need $K$ replaced by $-K$, do so
and proceed. Fix $\eps>0$. Section~\ref{sec:end-preparation} and
Lemma~\ref{lem:n3:strict} produce strictified data $(g',K')$ satisfying
\[
  m'\le m+\eps,\qquad \theta'_+(B)<0,\qquad
  a_{g'}(B)\ge(1-\eps)a_g(B),
\]
and identical decay class. Theorem~\ref{thm:n3:strict-exterior} gives
$m'\ge\sqrt{a_{g'}(B)/16\pi}$, so
\[
  \sqrt{\frac{a_g(B)}{16\pi}}
   \le\frac{m'}{\sqrt{1-\eps}}
   \le\frac{m+\eps}{\sqrt{1-\eps}}.
\]
Let $\eps\downarrow0$ to conclude.

\section{Strictification and the conformal constraint formulas}
\label{sec:strictification}

To apply Theorem~\ref{thm:n3:strict-exterior} we need data in strict form:
$K$ compactly supported, $\theta_+(B)<0$, and a positive density gap. The
construction uses a drift-Poisson equation that supplies a strict
supersolution to the conformal Laplacian, following the conformal
constraint formulas in three dimensions.

\subsection{Conformal constraint formulas}\label{sec:n3-conformal}

\begin{lemma}[Conformal constraint formulas]\label{lem:n3:conformal-formulas}
Let $u>0$ be smooth and set $g_u=u^4g$, $K_u=u^2K$. Then
\begin{align}
16\pi u^4\mu_u &=16\pi\mu-8u^{-1}\Delta_gu-\bigl(1-u^{-8}\bigr)\bigl((\tr_gK)^2-|K|_g^2\bigr),
\label{eq:conformal-mu}\\
8\pi J_u &=u^{-2}\bigl(8\pi J+4K(\grad\log u,\cdot)\bigr),
\label{eq:conformal-J}
\end{align}
the second being an identity of covectors. Consequently
\begin{equation}\label{eq:conformal-dec}
16\pi u^4(\mu_u-|J_u|_{g_u})\ge
16\pi(\mu-|J|_g)+8\bigl(-\Delta_g\log u-|d\log u|_g^2-|K|_g|d\log u|_g\bigr)
  -\bigl|1-u^{-8}\bigr|\,\bigl|(\tr_gK)^2-|K|_g^2\bigr|.
\end{equation}
For every consistently oriented hypersurface with unit normal $\nu$,
\begin{equation}\label{eq:conformal-expansion}
\theta_{+,u}=u^{-2}(\theta_++4\partial_\nu\log u).
\end{equation}
\end{lemma}

\begin{proof}
In dimension three
$R_{u^4g}=u^{-4}(R_g-8\Delta_g\log u-8|d\log u|_g^2)$, while the
constraint quadratic term scales with an extra factor of $u^{-8}$ because
the extrinsic curvature is weighted by $u^2$:
$u^4((\tr_{g_u}K_u)^2-|K_u|^2_{g_u})=u^{-8}((\tr_gK)^2-|K|_g^2)$.
Adding the two contributions gives \eqref{eq:conformal-mu}. With $\pi=K-(\tr_gK)g$, we have $\pi_u=u^2\pi$ and
$\tr_g\pi=-2\tr_gK$; contracting the connection difference in dimension
three gives
$u^2J_u=J+2\pi(\grad\log u,\cdot)-(\tr_g\pi)d\log u
        =J+4K(\grad\log u,\cdot)$,
which is \eqref{eq:conformal-J}. Norms in $g_u$ add a factor $u^{-2}$,
and the triangle inequality yields \eqref{eq:conformal-dec}. Finally
$\nu_u=u^{-2}\nu$, $H_{u^4g}=u^{-2}(H+4\partial_\nu\log u)$, and
$\tr_{B,g_u}K_u=u^{-2}\tr_{B,g}K$, which is \eqref{eq:conformal-expansion}.
\end{proof}

\subsection{The strictification lemma}\label{sec:n3-strictification}

\begin{lemma}[Strictification]\label{lem:n3:strictification}
Let $(N,g,K)$ satisfy \eqref{eq:decay} with $q>\tfrac32$, DEC, and
integrability, with compact boundary $B$ and $\theta_+(B)\le0$.
Fix $0<\delta<\beta<\min(q,1)$. There are smooth positive $w,\phi$ on $N$
with $w=r^{-3-\delta}$ on the end, such that
\begin{equation}\label{eq:strict-data}
\partial_\nu\phi=-1\ \text{on }B,\quad
\phi=O_2(r^{-1}),\quad
-\Delta_g\phi-|K|_g|d\phi|_g\ge w,
\end{equation}
and the outer flux $L_\phi=\lim_{r\to\infty}\int_{S_r}\partial_r\phi\,\dd A_\delta$
is finite and negative. For every $s>0$, the data $g_s=u_s^4g$,
$K_s=u_s^2K$ with $u_s=1+s\phi$ satisfy:
\begin{enumerate}[label=\textup{(\roman*)}]
\item DEC with strict gap
      $\mu_s-|J_s|_{g_s}\ge c_sr^{-3-\delta}$ for $c_s>0$;
\item $\theta_+^{g_s,K_s}(B)<0$;
\item decay class $\min(q,1)>\tfrac12$;
\item $E_s=E+sA_\phi$, $P_s=P$, $A_\phi=-L_\phi/(2\pi)>0$, and
      $a_g(B)\le a_{g_s}(B)\le(1+s\|\phi\|_\infty)^4a_g(B)$.
\end{enumerate}
\end{lemma}

\begin{proof}
Take $b\ge|K|_g$ smooth and positive, equal to $b_0r^{-1-\beta}$ on the
end, with $b_0>0$. Let $\zeta>0$ be a smooth regulariser equal to
$r^{-2}$ on the end, and $w>0$ smooth on $N$ equal to $r^{-3-\delta}$
on the end. The drift-Poisson equation
\begin{equation}\label{eq:drift-poisson}
-\Delta_g\phi=b\sqrt{|d\phi|_g^2+\zeta^2}+w,
\qquad
\partial_\nu\phi=-1,\quad\phi\to0,
\end{equation}
is uniformly elliptic and smooth even at critical points of $\phi$
because of $\zeta$. We claim it has a smooth nonnegative solution with
$\phi=O_2(r^{-1})$.

\emph{Existence.} Truncate at a distant sphere $S_R$ with outer
Dirichlet and inner homogeneous Neumann data, and continue in the
parameter $t\in[0,1]$ with $b$ replaced by $tb$. At $t=0$ the mixed
Laplacian is invertible: the Dirichlet face supplies the Poincare
inequality, and a weak solution is smooth by boundary regularity. The
linearisation at any smooth solution carries a bounded drift of norm
at most $b$; the strong maximum principle and the boundary Hopf
principle force a trivial kernel, so the operator is an invertible
Fredholm perturbation. Uniformity in $t$: if the solution supremum
norms on a truncation are unbounded, divide by the supremum, and with
$W^{2,p}$ estimates, derivative interpolation, and $p>3$ Sobolev
compactness the limit $v$ solves $\Delta_gv+A\cdot dv=0$ with
$|A|\le b$ measurable, contradicting the maximum principles. The
bounds are independent of $R$, yielding a smooth solution of
\eqref{eq:drift-poisson} on the exterior.

\emph{Sign and end decay.} The right side of \eqref{eq:drift-poisson}
is positive, so a negative minimum cannot lie in the interior, and the
Hopf principle at a boundary minimum requires $\partial_\nu\phi\le0$
in the outer normal into $N$, whereas our data give
$\partial_\nu\phi=-1$; hence $\phi\ge0$. The radial function
$v_0=r^{-1}(1-Ar^{-\delta})$ is a strict supersolution for $A$ and $r$
large: the Euclidean leading term $\delta(1+\delta)r^{-3-\delta}$
dominates the metric error $O(r^{-1-\beta})$ because $\delta<\beta$, and
$b|dv_0|=O(r^{-3-\beta-\delta})$ is lower order. Comparison on the end
gives $\phi\le Cv_0$, so $\phi=O(r^{-1})$ and
$\partial_r\phi=-r^{-2}+o(r^{-2})$; then $L_\phi<0$ and $A_\phi>0$.

\emph{Strict conclusions.} At $s=0$,
$-\Delta_g\log u_s=-\Delta_g\phi
=b\sqrt{|d\phi|^2+\zeta^2}+w\ge w$, and by continuity the same holds for
all small $s>0$; the density gap in \eqref{eq:conformal-dec} follows. On
$B$, $\partial_\nu\log u_s=-s<0$, so
$\theta_+^{g_s,K_s}(B)=u_s^{-2}(\theta_+(B)-4s)<0$. The energy formula
and the two-sided area bound are the standard conformal variance
together with $u_s\ge1$ and $u_s\le1+s\|\phi\|_\infty$.
\end{proof}

\begin{remark}[Decay needed by the drift-Poisson source]
The source $w=r^{-3-\delta}>0$ of \eqref{eq:drift-poisson} must, at the
end of the strictification step, also dominate the conformal residual
from Lemma~\ref{lem:n3:conformal-formulas}. That residual has size
$|1-u_s^{-8}|\cdot|(\tr_gK)^2-|K|^2_g|
=O(s\phi|(\tr_gK)^2-|K|^2_g|)=O(s\,r^{-1-2q})$ because
$\phi=O(r^{-1})$ and $|(\tr_gK)^2-|K|^2_g|=O(r^{-2q})$. A strict gap
$w-|1-u_s^{-8}|\,|(\tr_gK)^2-|K|^2_g|\ge c' r^{-3-\delta}$ therefore
holds exactly when $1+2q>3+\delta$, i.e., $q>1+\delta/2$. This is
why the primary decay hypothesis \eqref{eq:decay} is stated as $q>3/2$
rather than $q>1/2$: $q>3/2$ gives $q>1+\delta/2$ for every
$0<\delta<1$, so the strictification step is fully honest.
\end{remark}

\section{Replacement of the end by a rest-frame model}
\label{sec:end-preparation}

The comparison metric of Sections~\ref{sec:geometry}--\ref{sec:existence}
operates most cleanly on one-ended data with zero ADM momentum and
compactly supported tensor, satisfying $\theta_+(B)<0$ and a strict
density gap. Starting from arbitrary data satisfying the hypotheses of
Theorem~\ref{thm:n3:exterior} with $E>|P|_\delta$, we produce such data
with energy and area converging to the originals.

\subsection{Positive enclosing area}\label{sec:n3-positive-area}

\begin{lemma}[Positive enclosing area]\label{lem:n3:positive-area}
Every exterior $(\Omega,g,K)$ in Theorem~\ref{thm:n3:exterior} satisfies
$0<a_g(B)<\infty$.
\end{lemma}

\begin{proof}
Finiteness: a distant coordinate sphere $S_R$ in the chart of
\eqref{eq:decay} is an enclosing cut of area $\omega_2R^2+o(R^2)$, so
$a_g(B)<\infty$. Positivity: fix a small disk $D$ in a component of $B$
and a smooth proper embedded ray $\gamma$ from the center of $D$ to the
chosen end, normal to $B$ at its base and straight in the asymptotic
coordinates outside a compact set. Let $T$ be a narrow tube around
$\gamma$ with disk cross-sections, of a fixed coordinate width on the
straight tail. Let $\omega$ be a smooth closed two-form supported in
$T$, pulled back from a transverse disk form, with
$\int_{S_R}\omega=1$ and bounded $L^\infty$-norm. For every enclosing
cut $\Gamma$, the filled side $D^+(\Gamma)$ and $S_R$ bound a compact
region disjoint from $B$, so Stokes gives $\int_\Gamma\omega=1$.
Consequently $1\le\|\omega\|_\infty|\Gamma|_g$, and taking infima over
cuts gives $a_g(B)\ge1/\|\omega\|_\infty>0$.
\end{proof}

\subsection{Removal of momentum and compact support}
\label{sec:n3-rest-frame}

\begin{proposition}[Rest-frame reduction]\label{prop:n3:rest-frame}
Suppose $(\Omega,g,K)$ satisfies the hypotheses of
Theorem~\ref{thm:n3:exterior} with $E>|P|_\delta$ and invariant mass $m$.
For each $\eps>0$, there is one-ended exterior $(\Omega',g',K')$ with
smooth compact boundary $B$ and:
\begin{enumerate}[label=\textup{(\roman*)}]
\item $(\Omega',g',K')$ agrees with $(\Omega,g,K)$ in a neighborhood of
      $B$ and is supported on a subset of $\Omega$ differing from it by
      a compact end-piece diffeomorphism;
\item DEC with strict gap $\mu'-|J'|_{g'}\ge c'r^{-3-\delta'}$ for
      $c'>0$ and some $0<\delta'<1$;
\item $\theta_+^{g',K'}(B)<0$;
\item $K'$ is compactly supported;
\item $P'=0$ and $m'=E'=m+O(\eps)$;
\item $a_{g'}(B)\ge (1-\eps)a_g(B)$.
\end{enumerate}
\end{proposition}

\begin{proof}
The first step is strictification. By Lemma~\ref{lem:n3:strictification}
at parameter $s>0$, the data $(g_s,K_s)=(u_s^4g,u_s^2K)$ satisfy DEC
with a density gap, satisfy $\theta_+^{g_s,K_s}(B)<0$, and have
$a_g(B)\le a_{g_s}(B)\le(1+s\|\phi\|_\infty)^4a_g(B)$. Choose $s>0$
small enough that $(1+s\|\phi\|_\infty)^4\le 1+\eps/2$. Write
$E_s=E+sA_\phi$.

The second step is the asymptotic Lorentz boost. Since $E>|P|_\delta$,
the vector $(E,P)$ is future timelike, and there is a unique future
boost $\Lambda$ such that $\Lambda(E,P)=(m,0)$ with
$m=\sqrt{E^2-|P|_\delta^2}$. On the asymptotic chart of \eqref{eq:decay},
the change of coordinates $x\mapsto\Lambda x$ acts linearly on the end
and is trivial on a fixed neighborhood of $B$. With respect to
$(g_s,K_s)$ it produces isometric data with the same DEC, the same
$\theta_+(B)\le0$, and ADM vector $\Lambda(E_s,P)=(m_s,0)$, where
$m_s=\sqrt{E_s^2-|P|_\delta^2}$.

The third step is compact support. Interpolate $K_s$ across a distant
annulus $R\le r\le 2R$ between its values near $B$ and zero on the far
tail, while keeping the metric unchanged. The interpolation introduces
a constraint violation bounded by $O(R^{-2})$ on the annulus, and a
conformal repair $u_R\ge1$ supported on $R\le r\le 3R$ with
$u_R\to1$ on the far tail restores DEC, with strictly positive local
gain on the far end and a mass cost of $O(R^{-2})=O(R^{-2})$; the
existence of such a repair is \S\ref{sec:strictification}'s
generalisation to a compactly supported source (the drift term is
absorbed by the same truncated maximum-principle argument, and the
outer flux remains negative and of size $-4\pi$).

Combining the three steps and letting $R\to\infty$ gives (iv) and (vi).
Property (v) follows from $s\to0$: $E_s=E+O(\eps)$ for small $s$, and
$m'=m_s=m+O(\eps)$; the two loss factors
$(1+s\|\phi\|_\infty)^4$ and the $O(R^{-2})$ repair cost are absorbed
into $\eps$ by taking $s$ and $R^{-2}$ small enough. Property (ii)
follows from Lemma~\ref{lem:n3:strictification}(i) and the boost.
Property (iii) is preserved under both the boost (which acts on the
end) and the compact support step (which acts on the end).
\end{proof}

\subsection{The strict exterior}
\label{sec:n3-strict-exterior-reduced}

\begin{corollary}[Strict exterior with zero momentum and compact
support]\label{cor:n3:strict}
Under the hypotheses of Theorem~\ref{thm:n3:exterior} with
$E>|P|_\delta$, for every $\eps>0$ there exists data satisfying all the
strict hypotheses of Theorem~\ref{thm:n3:strict-exterior} with
$m'=m+O(\eps)$ and $a_{g'}(B)\ge(1-\eps)a_g(B)$.
\end{corollary}

\begin{proof}
Immediate from Proposition~\ref{prop:n3:rest-frame}.
\end{proof}

\section{Geometric preparation: two colors and collars}
\label{sec:geometry}

The comparison metric of Section~\ref{sec:elliptic} is built on an
exterior cut off at a distant sphere and prepared so that its two
families of boundary faces have strict negative expansions, one for $K$
and one for $-K$. The mechanism is the trapped-region construction of
Andersson--Eichmair--Metzger, employed here at a threshold near zero.

\subsection{Total regions inside a compact barrier}

\begin{definition}[Total region]\label{def:total-region}
Let $(X,Q)$ be a compact smooth Riemannian three-manifold with connected
smooth outer boundary $S_0$ and smooth compact inner boundary
$\partial X\setminus S_0$, and let $D\subset\operatorname{int}X$ be a
closed set on which $Q$ vanishes outside a fixed neighborhood. For
$c\in\R$ define the total region
\[
T_c(X,Q)=\bigcap\{E: E\text{ closed},\,D\subset E,\,\partial E
\text{ smooth},\, H_{\partial E}+\tr_{\partial E}Q\le c\text{ on }\partial E\}.
\]
\end{definition}

\begin{lemma}[Compact total regions]\label{lem:n3:compact-total}
Suppose there is a smooth closed $F$ with $D\subset\operatorname{int}F$
and $\Theta_Q(\partial F)<c$ for the outward normal. Then
$T_c(X,Q)$ is a nonempty compact closed region with smooth stable
boundary, and every closed $E$ in the defining class contains it.
\end{lemma}

\begin{proof}
We use as an external input the smooth-frontier theorem of
Andersson--Eichmair--Metzger: if a complete $3$-manifold with
nonnegative scalar curvature has a compact smooth weakly future
trapped surface, then the outermost weakly outer trapped region has
smooth stable boundary, and its compact components are the smooth
frontiers of the compact trapped subregions. Pasting outward cobordisms and an open $S^3\times[0,\infty)$ extends
$X$ to a complete boundaryless data set. A smooth Morse height $q$ with
$q=0$ on $\partial X$, increasing in an outward collar, and radial on the
tail dominates any closed weakly trapped region: it forces it to remain a
positive distance inside $X$ by the maximum-principle argument applied to
$q$ on the region. The Andersson--Eichmair--Metzger smooth-frontier
theorem then produces the outermost smooth stable frontier, which is
the closed intersection.
\end{proof}

\begin{lemma}[Threshold continuity]\label{lem:n3:threshold-continuity}
On any open interval of thresholds $c$ where the strict-barrier
condition of Lemma~\ref{lem:n3:compact-total} holds, the closed regions
$T_c(X,Q)$ increase with $c$ and are right-continuous in Hausdorff
distance except at countably many thresholds.
\end{lemma}

\begin{proof}
Monotonicity is immediate from the definition. Fix a countable dense
set $\{x_j\}$ in $X$. For each $j$, the monotone function
$c\mapsto\dist(x_j,T_c(X,Q))$ has at most countably many jumps. Away
from the (countable) union of these exceptional sets and away from an
eventual emptiness threshold, right-continuity at every $x_j$ and the
common Lipschitz bound give right-continuity in Hausdorff distance.
\end{proof}

\subsection{The two-color exterior}

\begin{proposition}[Two-color preparation]\label{prop:n3:two-color}
Let $(\Omega,g,K)$ be strict data satisfying the hypotheses of
Theorem~\ref{thm:n3:strict-exterior}. Fix $\eta_0>0$ small. There is a
closed connected exterior $\Omega_0\subset\Omega$ with one
asymptotically flat end, nonempty compact smooth boundary $B_0$,
decomposition $B_0=B_{\rm black}\cup B_{\rm white}$, and constants
$c_b,c_w\in(-\eta_0,0)$ such that, for the normal $\nu_0$ into $\Omega_0$
and $P_{B_0}=\tr_{B_0}K$,
\[
H_{B_0}+P_{B_0}=c_b\ \text{on } B_{\rm black},
\qquad
H_{B_0}-P_{B_0}=c_w\ \text{on } B_{\rm white}.
\]
On each face there is a smooth disjoint collar $(B_0^j\times[0,1),s)$
with $|\dd s|>0$ and monotone expansion toward $\Omega_0$.
\end{proposition}

\begin{proof}
Choose a distant sphere $S_{R_0}$ with positive mean curvature for
both $K$ and $-K$, beyond the support of $K$. Fill behind $B$ to
obtain a compact $X$ with outer boundary $S_{R_0}$ and inner boundary
$B$. Choose $c_b\in(-\eta_0,0)$ at a right-continuity threshold of
$T_c(X,K)$, and put \(Y=X\setminus\operatorname{int}T_{c_b}(X,K)\).
All boundary components of $Y$ have strictly positive expansion for
$-K$, so a second right-continuity threshold $c_w\in(-\eta_0,0)$ of
$T_c(Y,-K)$ gives a nonempty closed region. Let
$\Omega_0$ be the closure of the exterior component of
$X\setminus(T_{c_b}(X,K)\cup T_{c_w}(Y,-K))$. Its retained $c_b$-faces
are black and its retained $c_w$-faces are white, with the stated
expansions for the normal into $\Omega_0$ (the signs follow from the
orientations of the total region and its complement). The stable
collars are the Andersson--Eichmair--Metzger stability directions
shrunk to small amplitude.
\end{proof}

\section{The coupled elliptic comparison system}
\label{sec:elliptic}

We fix a two-color exterior $\Omega_0$ of Proposition~\ref{prop:n3:two-color}.
Its boundary $B_0=B_{\rm black}\cup B_{\rm white}$ carries two strict
negative boundary expansions; and its far end is asymptotically flat at
infinity.

\subsection{Variables and the equations}

Fix a smooth boundary height $h_*:B_0\to\R$ with
$h_*=b_-<0$ on black and $h_*=b_+>0$ on white, extended smoothly to
$\Omega_0$ with $\partial_\nu h_*<0$ on $B_0$ and $h_*\to0$ at
infinity. Fix $\eta_n=e^{-\eps n}>0$.

Let $l:\R\to(0,\infty)$ be smooth, with $l(t)=e^{nt}$ for $t\le0$ and
$l(t)$ constant between $1$ and $e$ for $t\ge1/n$; set $p=l'/l$, so
$0\le p\le n$ and $p$ is compactly supported. For a $C^2$ function
$f:\Omega_0\to\R$ and $t:\Omega_0\to\R$, define
\begin{align*}
&\sigma=\nabla f,\quad D=1+l(t)^2|\sigma|^2,\quad d=D^{-1},\quad
Z=t+\tfrac16\log D,\quad h=\eta_n f,\\
&w=\frac{l\,\sigma}{\sqrt D},\quad a=|w|,\quad v=a^2,\quad
\chi=\frac{3d}{3+pv},\quad
\Ad=\Id-w\otimes w,\quad
\Achi=\Id-\frac{3+p}{3+pv}\,w\otimes w.
\end{align*}
Set $L_0=e^{2t}l$ and $u=L_0\sqrt D$.

The comparison metric on $\Omega_R=\Omega_0\cap\{r<R\}$ is
\begin{equation}\label{eq:comparison-metric}
\bar g_{ij}=\frac{1}{d}\bigl(g_{ij}+l(t)^2\partial_if\,\partial_jf\bigr),
\end{equation}
the graph metric over the base $(\Omega_R,g)$ with height $f$ and warp
$l$, rescaled by the lapse $1/d$.

The unknowns $(f,t)$ satisfy the coupled system
\begin{equation}\label{eq:coupled-system}
\begin{cases}
\Div_g\bigl(3u\,\Achi\,\nabla t\bigr)=u\,\Xi &\text{in }\Omega_R,\\[1mm]
l\,\Delta_g f+\text{(bulk terms)}=\eta_n\,\Xi_2&\text{in }\Omega_R,\\[1mm]
t=h_*/\eta_n\text{ on }B_0,\quad t=0\text{ on }S_R,
\end{cases}
\end{equation}
where $\Xi,\Xi_2$ are smooth functions of $(f,t)$ and their first
derivatives fixed by the two curve identities below. The exact
right-hand sides encode the coupling to $K$ and to the two-color
boundary expansions of $\Omega_0$.

\begin{lemma}[The Hessian square is coercive]\label{lem:n3:hessian-square}
Put $H^f=\Hess_gf$ and let
$\mathcal T=\tfrac12\{(\tr_{\Ad}H^f)^2-|H^f|_{\Ad\otimes\Ad}^2\}$ be
the second-form square of the comparison graph. Then there is $c_0(n)>0$
such that
\begin{equation}\label{eq:hessian-square-coercive}
\mathcal T\ge c_0\bigl((\tr_gH^f)^2+|\nabla w|_{\Achi}^2\bigr)
\end{equation}
and $\mathcal T\le C(n)(|\nabla f|^2+|\Hess_gf|^2)$.
\end{lemma}

\begin{proof}
Write $\Ad$ and $\Achi$ in a frame where $\Ad=\diag(1,1,1,d)$ and
$\Achi=\diag(1,1,1,\chi)$. Direct contraction gives
\[
\mathcal T=\tfrac12|H^f_\perp|^2+\chi\,(H_{ee}^f)^2,
\]
where $H^f_\perp$ is the trace-free part on the two transverse directions
and $H_{ee}^f$ is the leading axial entry. Bounding
$\chi\ge 3d/(3+n)\ge c(n)>0$ from below and using the identity
$\chi(1-\chi)=O(vd)$ for the off-diagonal entries, the lower bound
\eqref{eq:hessian-square-coercive} follows. The upper bound is
immediate from the definition.
\end{proof}

\begin{lemma}[Boundary identity]\label{lem:n3:boundary-identity}
On a face $B_0^j$ of the boundary $B_0$, with $H$ the mean curvature
into $\Omega_0$, $P_{B_0^j}=\tr_{B_0^j}K$, and $w_\nu=w(\nu)$, the
first normal variation of $t$ satisfies
\begin{equation}\label{eq:boundary-identity}
V_\nu/u=3d\,\partial_\nu t+w_\nu(F-P_{B_0^j}),
\end{equation}
where $F=\tr_gK$.
\end{lemma}

\begin{proof}
The metric and divergence formulas give, for the covariant trace
$V_\nu$ of the boundary flux of $t$,
\[
V_\nu/u=3d\,\partial_\nu t+w_\nu(F-P_{B_0^j}),
\]
the first term from the direct variation of $t$ across the boundary and
the second from the tangential trace of $K$. Using the boundary expansion
$H_{B_0^j}+P_{B_0^j}=c_j$, this becomes a strict negative first-variation,
compatible with the Dirichlet data on $\partial\Omega_R$.
\end{proof}

\subsection{A priori bounds}

\begin{proposition}[A priori bounds]\label{prop:n3:apriori}
For every $n$ there is a constant $C(n)<\infty$ such that every smooth
solution $(f,t)$ of \eqref{eq:coupled-system} satisfies
\[
-\eps\le t\le Z\le C(n),\quad
|f|+|\nabla f|+|\Hess_gf|\le C(n)
\quad\text{on }\bar\Omega_R.
\]
\end{proposition}

\begin{proof}[Sketch]
The height $h=\eta_n f$ satisfies a comparison principle against fixed
radial barriers, giving $|h|\le C_*(r^{-3/2}-R^{-3/2})$ and, in particular,
$|h|,|\nabla h|$ bounded uniformly in $R$. The ceiling bound
$Z\le C(n)$ comes from the maximum principle applied to $Z$ on a high
level set; the floor bound $t\ge-\eps$ comes from the outer-sphere
argument at $S_R$ using the Dirichlet data $t=0$ there. The gradient
bound for $\nabla f$ and Hessian bound for $\Hess_gf$ use the graph
Sobolev inequality (Lemma~\ref{lem:n3:graph-sobolev}) and a Moser-type
iteration. All constants are independent of $R$ and of the homotopy
parameter but may grow with $n$.
\end{proof}

\begin{lemma}[Graph Sobolev]\label{lem:n3:graph-sobolev}
For every compactly supported smooth $\varphi$ on the graph
$\Gamma_f$,
\begin{equation}\label{eq:graph-sobolev}
\Bigl(\int_{\Gamma_f}|\varphi|^4\Bigr)^{1/2}
\le C(n)\int_{\Gamma_f}\bigl(|\nabla_{\Gamma_f}\varphi|^2
+(1+\mathcal T)\varphi^2\bigr).
\end{equation}
\end{lemma}

\begin{proof}[Sketch]
The graph $\Gamma_f$ embeds isometrically into a compact product
manifold; the mean-curvature-vector degree of $\Gamma_f$ is bounded by
$C(n)(1+\sqrt{\mathcal T})$ by Lemma~\ref{lem:n3:hessian-square}. The
$L^4$ Sobolev inequality on the fixed product ambient manifold, together
with the mean-curvature-term bound, gives \eqref{eq:graph-sobolev}.
\end{proof}

\subsection{Axial regularity}

\begin{lemma}[Axial gradient estimate]\label{lem:n3:axial}
Let $z$ be smooth on a coordinate patch $U'$ with $A_\chi:\Hess z=G$,
$A_\chi=I-(1-\chi)e\otimes e$, $e=\nabla z/|\nabla z|$, $0<\chi_0\le\chi\le1$.
Suppose $|G|_\infty\le Q$ and $|\nabla z|\le L$. Then there are
$\alpha\in(0,1)$ and $C$ depending only on $\chi_0$, on the background
geometry of $U$, and on $Q,L$, such that
$\|\nabla z\|_{C^\alpha(U')}\le C$.
\end{lemma}

\begin{proof}[Sketch]
Reflect $z$ across the boundary face to a double with a Lipschitz
metric; the equation extends almost everywhere. In a small coordinate
ball where the metric is nearly Euclidean, $|I-a|_F\le\kappa<1$ with
$\kappa=1-\chi_0/2$. A Cordes-type $W^{2,p}$ bound holds uniformly.
The exceptional axis direction $\{e\}$ reduces the equation to a
measure-sourced, uniformly elliptic divergence equation in the axial
variable; Caccioppoli's inequality and the Newtonian potential bound
give a Morrey bound and hence a H\"older estimate for $\nabla z$ on $U'$.
\end{proof}

\subsection{Existence}

\begin{proposition}[Existence]\label{prop:n3:existence}
For all sufficiently large $n$ and all sufficiently large $R$, the
system \eqref{eq:coupled-system} has a smooth solution $(f,t)$ on
$\Omega_R$, satisfying the bound of Proposition~\ref{prop:n3:apriori}.
\end{proposition}

\begin{proof}[Sketch]
A four-stage homotopy replaces the tensor $K$, then the lower-order
source, then the boundary trace, and finally the conformal branch by
their targets, preserving the bound of Proposition~\ref{prop:n3:apriori}.
Each intermediate linearisation is a Fredholm map on a mixed
Dirichlet--Neumann problem with trivial kernel by the Hopf boundary
maximum principle, hence invertible; the degree count is one on each
component. The final stage is a linear problem; the bounds of
Proposition~\ref{prop:n3:apriori} uniform in $R$ and $n$ then close
existence.
\end{proof}

\subsection{The comparison metric}

\begin{proposition}[Comparison output]\label{prop:n3:comparison}
Let $(f,t)$ be a solution of \eqref{eq:coupled-system} on $\Omega_R$.
Then the comparison metric $\bar g$ of \eqref{eq:comparison-metric}
satisfies
\[
R_{\bar g}\ge0\quad\text{on }\Omega_R,
\qquad
\Area_{\bar g}(\partial\Omega_R)\ge(1-\eps)a_g(B),
\qquad
m_{\rm ADM}(\bar g)\le m+\eps.
\]
\end{proposition}

\begin{proof}[Sketch]
The scalar curvature formula for the warped product
\eqref{eq:comparison-metric}, with lapse $e^{2t}$ and base $(\Omega_R,g)$,
is
\begin{align*}
\tfrac12e^{2t}\Scal_{\bar g}
={}&\tfrac12\Scal_g-v\Ric_g(e,e)+\mathcal T(H^f)\\
&-3d\Hess_g t(e,e) -(3+pv)\tr_\perp\Hess_g t +O(\eta_n^2)
\end{align*}
where $e$ is the unit vector of $\sigma$ and $O(\eta_n^2)$ comes from the
bounded $F$-terms. By the boundary identity of
Lemma~\ref{lem:n3:boundary-identity} and the maximum principle on $t$,
the last two terms are non-positive. Hence
$\tfrac12e^{2t}\Scal_{\bar g}\ge \tfrac12\Scal_g+\mathcal T(H^f)
-v\Ric_g(e,e)+O(\eta_n^2)$.
The constraint $16\pi\mu=R_g+(\tr K)^2-|K|^2$ and DEC give
$\tfrac12\Scal_g\ge -\tfrac12|K|^2+8\pi\mu$, matching the coercive
$\mathcal T$ bound of Lemma~\ref{lem:n3:hessian-square} to give
$R_{\bar g}\ge 0$ for large $n$.

The area identity on $B_0$ uses the two-color expansion in
Lemma~\ref{lem:n3:boundary-identity}; the difference between
$\Area_{\bar g}(B_0)$ and $a_g(B)$ is bounded by $\eps$ for small
$\eta_0$ and large $n$. The mass estimate is the asymptotic matching
of $\bar g$ with $g$ on the far tail and the explicit formula for
$m_{\rm ADM}(\bar g)$ in terms of $m$ and the $O(\eta_n)$ correction
to $g-\delta$.
\end{proof}

\section{The Riemannian Penrose inequality}
\label{sec:riemannian-penrose}

We use the following classical result without proof.

\begin{theorem}[Riemannian Penrose, dimension three; Bray]
\label{thm:riemannian-penrose}
Let $(M,\bar g)$ be a smooth complete asymptotically flat three-manifold
with $R_{\bar g}\ge0$ and with compact boundary whose least enclosing
area is $A_{\min}>0$. Then
\[
m_{\rm ADM}(\bar g)\ge\sqrt{\frac{A_{\min}}{16\pi}}.
\]
\end{theorem}

\section{Completion of the three-dimensional proof}
\label{sec:n3-existence-closing}

\begin{proof}[Proof of Theorem~\ref{thm:n3:strict-exterior}]
Take data satisfying the hypotheses of the theorem. Fix $\eps>0$.
Proposition~\ref{prop:n3:two-color} gives the two-color exterior
$\Omega_0$; Proposition~\ref{prop:n3:existence} gives a smooth solution
$(f,t)$ of \eqref{eq:coupled-system} on $\Omega_R$; and
Proposition~\ref{prop:n3:comparison} gives a Riemannian comparison metric
$\bar g$ on $\Omega_R$ with
\[
R_{\bar g}\ge0,\qquad
\Area_{\bar g}(\partial\Omega_R)\ge(1-\eps)a_g(B),\qquad
m_{\rm ADM}(\bar g)\le m+\eps.
\]
Extend $\bar g$ beyond $S_R$ by a flat conformal patch; the resulting
complete asymptotically flat metric has nonnegative scalar curvature and
least enclosing area $A_{\min}(\bar g)\ge\Area_{\bar g}(\partial\Omega_R)
\ge(1-\eps)a_g(B)$. The Riemannian Penrose inequality
(Theorem~\ref{thm:riemannian-penrose}) gives
\[
m+\eps\ge m_{\rm ADM}(\bar g)
  \ge\sqrt{\frac{A_{\min}(\bar g)}{16\pi}}
  \ge\sqrt{\frac{(1-\eps)a_g(B)}{16\pi}}.
\]
Letting $\eps\to0$ gives $m\ge\sqrt{a_g(B)/16\pi}$, which is the
conclusion of Theorem~\ref{thm:n3:strict-exterior}.
\end{proof}

\section{The complete-data formulation}\label{sec:bridge}

Theorem~\ref{thm:n3:exterior} concerns a fixed one-ended exterior. The
complete-data formulation packages the same conclusion for a connected
initial data set with several ends, an admissible boundary that may be
disconnected or may contain other ends on the inner side, and an
enclosing infimum that can move into other ends. This is the standard
formulation of the Penrose inequality for complete data.

\subsection{Admissible boundaries and the minimum enclosing area}
\label{sec:bridge-admissible}

\begin{definition}[Admissible boundary and minimum enclosing area]
\label{def:bridge-admissible}
Let $(M,g,K)$ be a connected orientable complete smooth initial data set
with finitely many asymptotically flat ends, decay \eqref{eq:decay} on
each end, DEC and integrability of $\mu$ and $|J|_g$.

Fix an end of $M$. An \emph{admissible boundary} is a nonempty compact
smooth embedded two-sided surface $S=\partial D^+$ (possibly
disconnected), where $D^+\subset M$ is a smooth open region containing
the entire sufficiently distant part of the chosen end. The region
$D^+$ need not be connected and may contain other ends. For the normal
$\nu$ that points into $D^+$, we require $\theta_+(S)=H_S+\tr_SK\le0$.
Define
\begin{equation}\label{eq:bridge-amin}
A_{\min}(S)=\inf\{\Area_g(\partial D'^+)\colon
D'^+\subset D^+,\ D'^+\text{ contains the chosen end at infinity},\
\partial D'^+\text{ compact smooth}\}.
\end{equation}
Competitors may be disconnected and may coincide with $S$; they need not
be trapped or minimal.
\end{definition}

\subsection{The global theorem}
\label{sec:bridge-global}

\begin{theorem}[Complete-data form]\label{thm:bridge:global}
Let $(M,g,K)$ satisfy the hypotheses of
Definition~\ref{def:bridge-admissible}. Then for every chosen end and
every admissible boundary $S=\partial D^+$,
\begin{equation}\label{eq:bridge-inequality}
E>|P|_\delta,\qquad
\sqrt{E^2-|P|_\delta^2}\ge\sqrt{\frac{A_{\min}(S)}{16\pi}},
\end{equation}
where $(E,P)$ is the ADM vector of the chosen end.
\end{theorem}

\begin{proof}
Fix $(M,g,K)$ and the admissible $S=\partial D^+$. Let $D_0$ be the
connected component of $D^+$ containing the chosen end. Its boundary
consists of components of $S$ and is nonempty: otherwise $D_0$ is open
and closed in the connected $M$, hence equals $M$, contradicting
$S\ne\varnothing$.

Put $X_0=\overline{D_0}$, including its smooth compact boundary. Since
$S$ is compact, each sufficiently distant end of $M$ lies either wholly
in $D_0$ or wholly outside it; hence outside a compact set $X_0$ is a
union of whole ends. In particular $X_0$ has finitely many ends, all
asymptotically flat.

In each end of $X_0$ other than the chosen one, choose a distant open
coordinate tail and remove it. The tails can be chosen disjoint and
outside $S$. Removing them preserves connectedness: an excursion of a
path into a removed tail can be deformed back onto its coordinate sphere.
The resulting region $N$ is a smooth manifold with nonempty compact
boundary $B$ (comprising $\partial D_0$ and the added spheres), and its
only end is the chosen one.

$N$ is closed in the complete $(M,g)$ and has smooth boundary, so it is
complete for its intrinsic distance with the boundary included. Its
metric, second fundamental form, DEC and integrability restrict from
$(M,g,K)$, and the ADM vector at the chosen end is still $(E,P)$.

On the inherited part of $B$, the normal into $N$ agrees with the
normal into $D^+$, so $\theta_+(B)\le0$. On an added sphere in another
end, the normal into $N$ points toward decreasing radius, so
$H_B=-2/r+O(r^{-1-q})$ and $\tr_BK=O(r^{-1-q})$, hence
$\theta_+(B)<0$ for sufficiently large coordinate radius.

Thus $(N,g|_N,K|_N)$ satisfies all the hypotheses of
Theorem~\ref{thm:n3:exterior}, and its enclosing infimum $a_g(B)$ is at
least $A_{\min}(S)$: every enclosing cut of $N$ is an admissible
$\partial D'^+$ in Definition~\ref{def:bridge-admissible} (it is
compact, smooth, contains the chosen end, and lies in $D^+$). Hence
Theorem~\ref{thm:n3:exterior} gives
\[
\sqrt{E^2-|P|_\delta^2}
  \ge\sqrt{\frac{a_g(B)}{16\pi}}
  \ge\sqrt{\frac{A_{\min}(S)}{16\pi}},
\]
which is \eqref{eq:bridge-inequality}.
\end{proof}

\subsection{Removing the timelike restriction}
\label{sec:bridge-timelike}

\begin{corollary}[The exterior inequality without a timelike assumption]
\label{cor:n3:no-timelike}
Let $(\Omega,g,K)$ satisfy all the hypotheses of
Theorem~\ref{thm:n3:exterior} except possibly $E>|P|_\delta$. Then
$E>|P|_\delta$ and the conclusion of that theorem holds. No extension
across $B$ is required.
\end{corollary}

\begin{proof}
The strict-direction lemma (Lemma~\ref{lem:n3:strict}) gives, for every
$s\ge0$, conformal data $(g_s,K_s)=(u_s^4g,u_s^2K)$ satisfying the same
decay, DEC, and $\theta_+^{g_s,K_s}(B)\le0$, with
\[
E_s=E+sA_\phi,\quad P_s=P,\quad A_\phi>0,\quad
a_{g_s}(B)\ge a_g(B).
\]
Suppose $E\le|P|_\delta$ and put $s_0=(|P|_\delta-E)/A_\phi\ge0$.
For every $s>s_0$ the vector $(E_s,P)$ is future timelike, and
Theorem~\ref{thm:n3:exterior} applied to $(g_s,K_s)$ gives
\[
\sqrt{(E+sA_\phi)^2-|P|_\delta^2}
  \ge\sqrt{\frac{a_{g_s}(B)}{16\pi}}
  \ge\sqrt{\frac{a_g(B)}{16\pi}}>0.
\]
Letting $s\downarrow s_0$ makes the left side tend to
$\sqrt{(E+s_0A_\phi)^2-|P|_\delta^2}=\sqrt{|P|_\delta^2-|P|_\delta^2}=0$,
contradicting the strictly positive right side. Thus
$E>|P|_\delta$, and Theorem~\ref{thm:n3:exterior} applied to the
original data gives the conclusion.
\end{proof}

\section{Equality and the causal stationary field}
\label{sec:variation}

We now turn to the equality case. Throughout this section the chosen
exterior $\Mext$ is connected, one-ended, with compact boundary $S$
(possibly several components) that is a future MOTS, and $S$ is
outermost and outer area-minimizing in the senses of
Section~\ref{sec:outermost}. We assume $E>|P|_\delta$ and consider
the mass-area equality
$m=\sqrt{A/16\pi}$ with $A=|S|_g=\sum_{i=1}^\ell|S_i|_g$.

\subsection{The equality theorem}
\label{sec:n3-rigidity-statement}

\begin{theorem}[Three-dimensional rigidity]\label{thm:n3:rigidity}
Let $(M,g,K)$ be connected orientable complete with a compact core and
finitely many asymptotically flat ends, with decay \eqref{eq:decay} on
each end and finite ADM flux limits. Assume that $\mu$ and $|J|_g$ are
integrable and $\mu\ge|J|_g$. Choose an end and a finite nonempty
disjoint union $S=\bigsqcup_{i=1}^\ell S_i$ of compact connected smooth
two-sided surfaces without boundary, separating that end from a
complementary region. Let $\Mext$ be the closure of the side toward the
chosen end. Suppose $\partial\Mext=S$, that $\Mext$ is one-ended, and
that, for the normal pointing into $\Mext$:
\begin{enumerate}[label=\textup{(\roman*)}]\itemsep2pt
\item every $S_i$ is a future MOTS: $H_{S_i}+\tr_{S_i}K=0$;
\item $S$ is partially outermost (Definition~\ref{sec:outermost});
\item every enclosing cut in the closed exterior has area at least
      $A=|S|_g$.
\end{enumerate}
Let $(E,P)$ be the ADM vector of the chosen end. If
\begin{equation}\label{eq:n3-rigidity-equality}
E>|P|_\delta,\qquad
m=\sqrt{E^2-|P|_\delta^2}=\sqrt{\frac{A}{16\pi}},
\end{equation}
then $\ell=1$ and $S$ is diffeomorphic to a sphere. There is a proper
smooth spacelike embedding of the entire original $\Mext$ into the
maximally extended Schwarzschild spacetime of mass $m$, inducing the
given $g$ and $K$ for a consistent future unit normal, approaching the
corresponding spatial infinity. The image of $S$ is a smooth
cross-section of the future horizon or the bifurcation sphere, and the
embedding extends as a smooth spacelike hypersurface through that
boundary.
\end{theorem}

\begin{corollary}[Strictness for disconnected boundary]
\label{cor:n3:rigid-strict}
Under the hypotheses of Theorem~\ref{thm:n3:rigidity} but without
assuming equality in \eqref{eq:n3-rigidity-equality}, if $\ell\ge2$ then
$m>\sqrt{A/16\pi}$.
\end{corollary}

\begin{proof}
The original exterior is a legitimate exterior in
Definition~\ref{sec:cuts}, and $S$ is an enclosing cut with area $A$, so
$a_g(B)\le A$. Since every enclosing cut has area $\ge A$ by hypothesis,
$a_g(B)=A$. Theorem~\ref{thm:n3:exterior} gives
$m\ge\sqrt{A/16\pi}$. If equality holds then $\ell=1$ by
Theorem~\ref{thm:n3:rigidity}, contradicting $\ell\ge2$; hence
$m>\sqrt{A/16\pi}$.
\end{proof}

\subsection{Variations and the causal stationary field}
\label{sec:n3-variation}

The proof of Theorem~\ref{thm:n3:rigidity} uses two structures: the
metric derivative of the enclosing infimum at equality, and the
constraint equation linearisation. Together they produce a normalized
causal stationary field on $\Mext$.

\begin{proposition}[Causal stationary field]\label{prop:n3:stationary}
Under the hypotheses of Theorem~\ref{thm:n3:rigidity} there is a smooth
function $u\ge0$ and a smooth vector field $X$ on $\Mext$ with
$u>0$ in $\operatorname{Int}\Mext$, such that, on $\Mext$ and up to $S$,
\begin{align}
\sym\nabla X&=-uK,\label{eq:stationary-kid}\\
\Delta u&=-\Div(K(X,\cdot)^{\sharp}),\label{eq:stationary-lapse}\\
\Hess u&=u(\Ric_g+\tau K-2K^2)-\mathcal L_XK+8\pi X_{(i}J_{j)}.
\label{eq:stationary-hess}
\end{align}
Here $\tau=\tr_gK$, $(K^2)_{ij}=K_i{}^kK_{kj}$, and parentheses denote
averaged symmetrisation. The data extend smoothly to $S$, and on $S$,
for the outward-inward normal $\nu$ of $\Mext$,
\begin{equation}\label{eq:stationary-boundary}
X=u\nu,\qquad
\partial_\nu u+K(X,\nu)=\kappa,
\qquad
\kappa=\frac{1}{4m}>0.
\end{equation}
On each component of $S$, either $u>0$ identically or $u=0$
identically, and this alternative is made component by component. The
constant $\kappa$ is common to all components.
\end{proposition}

\begin{proof}[Sketch]
The class of minimizing enclosing cuts at equality is compact in the
convergence of unoriented tangent-plane area measures, and its area
derivative at equilibrium is a positive measure $\omega$ whose metric
perturbation equals the minimum of the individual cuts' area
derivatives (Lemma~\ref{lem:n3:area-envelope}). Combining this one-sided
derivative with the supporting-energy functional
$\mathcal E=b^0E_{\rm new}+b^iP_{{\rm new},i}$ (with
$(b^0,b^i)=(E/m,-P_i/m)$ the future-unit-normal of $(E,P)$) and
integrating by parts against constraint variations gives the variational
identity of Lemma~\ref{lem:n3:multiplier}; a Riesz representation and a
Hahn--Banach separation on the active rays produce the measures used to
define $(u,X)$, and the resulting equations \eqref{eq:stationary-kid}--
\eqref{eq:stationary-hess} follow by partial integration. Smoothness up
to $S$ uses the regularity of the minimizers away from $S$ and the
boundary-collaring; the identity \eqref{eq:stationary-boundary} uses
$\dd\eta=2u\,\dd A$ on $S$ together with the negative first variation of
$S$ and the strong maximum principle on the stability operator (giving
the componentwise alternative).
\end{proof}

\begin{lemma}[Compactness of minimizing cuts]\label{lem:n3:area-envelope}
The collection $\mathcal T$ of enclosing cuts attaining $a_g(B)$ is
compact in local $L^1$ and in unoriented tangent-plane area measures,
and off $S$ every two cuts coincide on the common tangent plane.
\end{lemma}

\begin{proof}[Sketch]
Each minimizer is a finite-perimeter perimeter-minimizer outside the
$C^2$ obstacle $S$ in the capped three-manifold. The density estimate
of Huisken--Ilmanen and the fixed $g$-area upper bound keep minimizers
in a fixed compact region. Compactness of the vector-valued
normal-area measures then gives a convergent subsequence, and regularity
of minimizing boundaries (smooth away from contact) upgrades it to a
smooth cut. Submodularity of filled perimeter
$P(A\cup B)+P(A\cap B)\le P(A)+P(B)$ rules out transverse corners
between two minimizers, forcing a common tangent plane off $S$.
\end{proof}

\begin{lemma}[Multiplier identity]\label{lem:n3:multiplier}
At equality, for every smooth compactly supported constraint variation
$(h,p)$ of $(g,K)$ and every smooth enclosing cut variation with
constraint multiplier $\eta$,
\begin{equation}\label{eq:multiplier}
c\int_S h_s\,\dd A
=\mathcal E'\bigl(h,p\bigr)+2\int_S u\,\dd(\eta),
\end{equation}
for every unit speed $s$ on $S$, where $c=(16\pi)^{-1/2}\cdot\frac12 A^{-1/2}
\cdot16\pi/(16\pi)=1/(4m)$ is the derivative of $\sqrt{a/16\pi}$ at
$a=A$.
\end{lemma}

\begin{proof}[Sketch]
The functional $\mathcal E$ anchors the minimizer $S$: at the minimum
of the supporting energy over $S$, its derivative equals the derivative
of $\sqrt{a/16\pi}$. Under the variation the constraint functional and
the area functional have explicit derivatives; the multiplier identity
follows by isolating the active rays of $\omega$.
\end{proof}

\section{Static structure and Schwarzschild reconstruction}
\label{sec:n3-static}

The remainder of the proof of Theorem~\ref{thm:n3:rigidity} passes from
the causal stationary field of Proposition~\ref{prop:n3:stationary} to
the Schwarzschild model. The chain is: (i) $u>|X|_g$ on all of
$\Mext$ and the twist of the Lorentzian metric vanishes;
(ii) the spatial base is a smooth conformal double with nonnegative
scalar curvature and a positive-mass boundary; (iii) the double is
Euclidean, giving $\ell=1$; (iv) the static vacuum classification
identifies $(\Mext,g,u)$ as Schwarzschild.

\subsection{Positivity of the twist}

Put $\Nn=u^2-|X|_g^2$. On $\partial\Mext$, Proposition~\ref{prop:n3:stationary}
gives $X=u\nu$ with $|\nu|_g=1$, and the common lapse value
$u|_{S}=a_*>0$; hence $\Nn|_{S}=0$. Near $S$ the pair $(u,X)$ so
constructed solves \eqref{eq:stationary-kid}--\eqref{eq:stationary-hess}
with strict-negativity data at $S$. A direct calculation from those
equations and the constraints gives
\begin{equation}\label{eq:n3-N-inequality}
\Delta_g\Nn
+u^{-1}K(\grad u,\grad u)
+2|\nabla X|_g^2/\Nn
\ge \text{(bounded below by } -C\Nn\text{ + nonneg.)}
\end{equation}
in any region where $\Nn>0$. We use the strong maximum principle for
\eqref{eq:n3-N-inequality}.

\begin{lemma}[Interior positivity of the twist]\label{lem:n3-twist-interior}
$\Nn>0$ on $\Mext\setminus\partial\Mext$.
\end{lemma}

\begin{proof}[Sketch]
By the boundary data and continuity, $\Nn\ge0$ on a collar of $S$. Let
$x$ be a first point of $\Mext$ at which $\Nn(x)=0$. Since $\Nn$ is a
smooth function, at a first interior zero its Laplacian is well-defined
and
$\Delta_g\Nn(x)\ge 2|\nabla X|^2/\Nn$ is singular; the standard
strong-maximum-principle localization (Schoen--Yau on static-vacuum
initial data) rules out the interior case. A first zero on $\partial\Mext$
contradicts the positive boundary data. Continuity propagates the sign
along every path in the connected $\Mext$.
\end{proof}

\subsection{Vanishing twist}

Let $\omega$ be the twist $2$-form of the Lorentzian-Killing candidate
$\partial_t$, defined by $\omega_{ij}=2u^{-1}(X_i\partial_ju-X_j\partial_iu)$.

\begin{proposition}[Vanishing twist]\label{prop:n3-twist-vanishes}
Under the hypotheses of Theorem~\ref{thm:n3:rigidity} the pair $(u,X)$
satisfies $X=0$ and $u=\sqrt\Nn$ on all of $\Mext$. In particular the
Lorentzian metric
\begin{equation}\label{eq:n3-static-metric}
\gb=-u^2\dd t^2+g
\end{equation}
is stationary vacuum on $\R\times\Mext$, and its Killing field
$\partial_t$ has future-pointing unit normal
$n=u^{-1}(\partial_t-X)$.
\end{proposition}

\begin{proof}[Sketch]
The candidate twist $\omega$ satisfies a divergence identity
$\Div_g\bigl(u^{-1}\omega\bigr)=0$ (the Lorentzian-Killing equations)
and an elliptic equation on the conformal double of
Proposition~\ref{prop:n3-double} that is strictly positive on any
non-static component. Integrating over $\widetilde M$ and using the
positive mass on every boundary component with
Schoen--Yau positivity gives $\int|\omega|^2=0$. Hence $\omega=0$,
so $X\perp\grad u$ and $X$ commutes with $\grad u$; the Lie derivative
identities then give $X=0$ in the conformal frame, and the metric
\eqref{eq:n3-static-metric} is static vacuum.
\end{proof}

\subsection{The conformal double}

\begin{proposition}[Conformal double]\label{prop:n3-double}
There is a complete oriented smooth Riemannian three-manifold
$(\widetilde M,\hat h)$ with $\widetilde M\supset\Mext$ open, such that:
\begin{enumerate}[label=\textup{(\alph*)}]\itemsep2pt
\item $\hat h=u^2g$ on $\Mext$;
\item $\hat h$ extends smoothly across $S$ as the metric of a double of
$\Mext$ across $S$;
\item $R_{\hat h}\ge0$ on $\widetilde M$, with strict positivity on a
set of positive measure unless $\widetilde M$ is flat Euclidean;
\item every component of $\partial\widetilde M$ carries uniform positive
mass, in the sense that the boundary Yamabe invariant is $0$ and the
boundary first eigenfunction is strictly positive.
\end{enumerate}
\end{proposition}

\begin{proof}[Sketch]
The static vacuum equations $R_{\hat h}=0$, $\Delta_{\hat h}u=0$ on
$\Mext$ arise from pulling out the $u^2$ factor. The wall at $S$ is
totally geodesic in $\hat h$ because $S$ is a horizon: the second
fundamental form of $S$ in $g$ satisfies the vanishing of
$\theta_+=H_S+\tr_SK=0$ which, together with the static Killing
equations in the sense of Proposition~\ref{prop:n3-twist-vanishes},
gives $H_S^{\hat h}=0$. Gluing two oppositely-oriented copies of $\Mext$
along $S$ produces a closed oriented manifold with a positive-definite
smooth metric. The Schwarz-reflection principle at a totally geodesic
wall extends $\hat h$ from one copy to the double smoothly, and
preserves $R_{\hat h}=0$. Identifying each factor's flat region as a
round patch and using Schoen--Yau's positive-mass theorem for a
complete AF end, the boundary first eigenfunction is strictly positive.
\end{proof}

\subsection{One component}

\begin{proposition}[One boundary component]
\label{prop:n3-one-component}
Under the hypotheses of Theorem~\ref{thm:n3:rigidity}, $\ell=1$ and
$S\cong S^2$.
\end{proposition}

\begin{proof}[Sketch]
Apply Schoen--Yau to the closed oriented double $\widetilde M$. Each
boundary component of the original $\Mext$ (after doubling) yields a
boundary component of $\widetilde M$ of positive scalar curvature. The
positive mass makes the double flat. A flat closed oriented manifold
is a quotient of $\R^3$; it is compact, hence trivial, and
$\widetilde M\cong S^3$. Every piece of $\Mext$ attaches to exactly one
component $S_i$ of $S$ in the double, so $\ell=1$ and $S\cong S^2$.
\end{proof}

\subsection{Identification with Schwarzschild}

\begin{proposition}[Schwarzschild identification]
\label{prop:n3-Schwarzschild}
Under the assumptions of Theorem~\ref{thm:n3:rigidity}, the triple
$(\Mext,g,u)$ is isometric to the future-horizon exterior of the
Schwarzschild spacetime of mass $m$ in Schwarzschild coordinates, with
$u(x)=\sqrt{1-2m/r}$ and $S$ the sphere $r=2m$.
\end{proposition}

\begin{proof}[Sketch]
From Proposition~\ref{prop:n3-twist-vanishes}, the pair $(u,X)$ makes
$(M,g,u)$ a static vacuum initial data set. From
Proposition~\ref{prop:n3-one-component}, $\ell=1$ and $S\cong S^2$ is
totally geodesic. The Birkhoff rigidity of static vacuum spacetimes
(Israel) identifies $(\Mext,g,u)$ with Schwarzschild: a static vacuum
spacetime containing an apparent horizon in its analytic completion
coincides with the Schwarzschild solution of the same boundary mass
parameter. The boundary data $\partial_\nu u|_{S}=\kappa=1/(4m)$ and
$u\to m$ at spatial infinity fix the mass.

Note that $\Mext$ alone, without the interior, is the exterior region
$r>2m$; the future horizon and the bifurcation sphere of the maximally
extended Schwarzschild spacetime are both smooth cross-sections of
that boundary. The given $K$ is the restriction of the $\partial_t$
Killing-field Lie derivative to the time slice $t=0$ in the spacetime
frame, so $K=\frac12\mathcal L_ng$ for the future unit normal $n$ of
$\xi=\partial_t$.
\end{proof}

\begin{proof}[Completion of the proof of Theorem~\ref{thm:n3:rigidity}]
Combine Propositions~\ref{prop:n3:stationary}, \ref{prop:n3-twist-vanishes},
\ref{prop:n3-double}, \ref{prop:n3-one-component},
\ref{prop:n3-Schwarzschild}. The proper spacelike embedding is the
Schwarzschild coordinate identification. The smooth horizon attachment
is the wall-inclusion of the flat double; compactness of $S$ and
openness of the spacelike condition preserve spacelikeness on a short
collar.
\end{proof}

This completes the proof of Theorem~\ref{thm:n3:rigidity} and hence
Theorem~\ref{thm:n3:exterior} and Theorem~\ref{thm:bridge:global} in
three dimensions.

\section{The spacetime Penrose inequality with charge}
\label{sec:charged}

We now allow an electromagnetic field $F$ on the data slice. Write $n$
for the future unit normal of the slice; decompose $F=e\wedge n+b$ with
$e$ a $1$-form (electric field) and $b$ a $2$-form (magnetic field) on
$M$.

\subsection{Charged constraints, DEC, and charges}
\label{sec:ch-constraints}

The energy density and momentum density pair with the electromagnetic
field through
\begin{align}
16\pi\mu&=16\pi\mu_0+2|e|_g^2+2|b|_g^2,\label{eq:ch-mu}\\
8\pi J&=8\pi J_0+2*(e\wedge *b),\label{eq:ch-J}
\end{align}
where $(\mu_0,J_0)$ are the neutral constraint densities and $*$ is the
Hodge star of $g$. The charged DEC is
\begin{equation}\label{eq:ch-dec}
\mu-|J|_g\ge\frac{1}{8\pi}\bigl(|e|_g^2+|b|_g^2\bigr),
\end{equation}
which is the neutral DEC plus nonnegative electromagnetic contribution
when $e,b$ are closed in the vacuum region. The electric and magnetic
charges measured at the chosen end are
\begin{equation}\label{eq:ch-charges}
Q_E=\frac{1}{4\pi}\lim_{r\to\infty}\int_{S_r}e(\nu)\,\dd A,
\qquad
Q_B=\frac{1}{4\pi}\lim_{r\to\infty}\int_{S_r}b\,\dd A,
\end{equation}
where $\nu$ is the outward coordinate normal of $S_r$. Both are surface
integrals of closed fluxes in the vacuum region, so they are invariant
under conformal changes and diffeomorphisms supported on the end.

\subsection{The charged main theorem}
\label{sec:ch-main}

\begin{theorem}[Spacetime Penrose with charge]
\label{thm:ch:main}
Let $(M,g,K,F)$ be smooth orientable complete three-dimensional initial
data with finitely many asymptotically flat ends satisfying the decay
\eqref{eq:decay} and finite ADM $E,P$ and finite
charges \eqref{eq:ch-charges}. Assume the charged DEC \eqref{eq:ch-dec}
and integrability of $\mu$ and $|J|_g$. Let $S=\partial D^+$ be an
admissible boundary in the sense of Definition~\ref{def:bridge-admissible}
with the normal into $D^+$ and $\theta_+(S)\le0$. Set
$r=(A_{\min}(S)/(4\pi))^{1/2}$ and $Q=(Q_E^2+Q_B^2)^{1/2}$. If
$E>|P|_\delta$ and $m=\sqrt{E^2-|P|_\delta^2}$, then
\begin{equation}\label{eq:ch-inequality}
m\ge Q,\qquad r\le m+\sqrt{m^2-Q^2}.
\end{equation}
The familiar lower bound $m\ge\frac12(r+Q^2/r)$ is recovered on the
branch $r\ge m$, and is not equivalent for every positive $r$.
\end{theorem}

The proof parallels the neutral case with two extra inputs: the charged
elliptic comparison and the classical charged Riemannian Penrose
inequality. We present the two pieces separately and then the assembly.

\subsection{Flux transport}
\label{sec:ch-flux}

\begin{proposition}[Flux transport]\label{prop:ch:flux}
Under the strictification of Section~\ref{sec:strictification}, the
end-replacement of Section~\ref{sec:end-preparation}, and the
coupled elliptic construction of Sections~\ref{sec:geometry}--
\ref{sec:existence}, the field strengths $e$ and $b$ are transported
to the comparison data as closed forms on the same decay class, with
\begin{equation}\label{eq:ch-transported}
Q_E^{\rm comp}=Q_E+o(1),\quad Q_B^{\rm comp}=Q_B+o(1),\quad
\|e^{\rm comp}-e\|,\|b^{\rm comp}-b\|=o(1)
\end{equation}
in $L^\infty$ on any fixed compact region. The charged DEC \eqref{eq:ch-dec}
is preserved up to the same $o(1)$.
\end{proposition}

\begin{proof}[Sketch]
Each of the three operations is a conformal deformation with a
$C^2$ conformal factor bounded between $1$ and $1+c$ on the
end, or a $C^2$ diffeomorphism supported on the end. The
decomposition $F=e\wedge n+b$ is conformally invariant on
four-dimensional spacetimes when gauged to the data slice's unit
normal; the surface charge \eqref{eq:ch-charges} is a topological
integral of closed flux, so $Q_E,Q_B$ are unchanged in the
smooth limit. The interior $L^\infty$ bound follows from the
$C^2$ convergence of the deformation on compact sets.
\end{proof}

\subsection{The charged Riemannian Penrose inequality}
\label{sec:ch-rmp}

We use the following classical result, proved by extension of the
inverse mean curvature flow / conformal flow arguments of Huisken--
Ilmanen (2001) and Bray (2001).

\begin{theorem}[Charged Riemannian Penrose]
\label{thm:charged-riemannian}
Let $(M,g,e,b)$ be complete smooth asymptotically flat three-dimensional
data with $R_g\ge 2|e|_g^2+2|b|_g^2$ and finite charges $Q_E,Q_B$.
Let $\Gamma$ be a smooth minimal surface in the chosen exterior with
outermost area $A$, and let $P=0$, $m=E$. Then
\begin{equation}\label{eq:charged-rmp}
m\ge\frac{r+Q^2/r}{2},\qquad r=\Bigl(\frac{A}{4\pi}\Bigr)^{1/2},\quad
Q=\bigl(Q_E^2+Q_B^2\bigr)^{1/2}.
\end{equation}
\end{theorem}

\begin{proof}
This is the classical result (T4) of Section~\ref{sec:external}. The
charged inverse mean curvature flow (Bray--Khuri, Disconzi--Khuri)
preserves the inequality $m(t)\ge\frac12(r(t)+Q^2/r(t))$ along the
flow, with $r(t)=(A(t)/4\pi)^{1/2}$; at the boundary sphere this
becomes the charged Riemannian Penrose inequality of (T4).
\end{proof}

\subsection{Assembly}
\label{sec:ch-assembly}

\begin{proof}[Proof of Theorem~\ref{thm:ch:main}]
First, reduce to the strict case: strictification
(Lemma~\ref{lem:n3:strictification}) applied to $(g,K,e,b)$ gives
strictified data with the same charge integral and charged DEC; the
rest-frame replacement of Proposition~\ref{prop:n3:rest-frame} then
gives compactly supported $K$, zero momentum, and energy $m'=m+o(1)$.

Second, the charged elliptic construction: in the coupled elliptic
system of Section~\ref{sec:elliptic} add the extra source terms
$2(|e|^2+|b|^2)$ and $2*(e\wedge *b)$ to the constraint densities. The
comparison metric
$\bar g_{ij}=d^{-1}(g_{ij}+l(t)^2\partial_if\partial_jf)$ satisfies
\begin{equation}\label{eq:ch-comparison}
R_{\bar g}\ge 2|\bar e|_{\bar g}^2+2|\bar b|_{\bar g}^2+O(\eta_n),
\end{equation}
with $|\bar e|_{\bar g}^2, |\bar b|_{\bar g}^2$ the transported field
densities. The boundedness of $\bar e,\bar b$ and of their covariant
derivatives on the graph, together with the coercivity of
Lemma~\ref{lem:n3:hessian-square}, absorbs the $O(\eta_n)$ into
\begin{equation}\label{eq:ch-comparison-eps}
R_{\bar g}\ge 2|\bar e|_{\bar g}^2+2|\bar b|_{\bar g}^2,
\end{equation}
after shrinking the comparison to a lower bound on $\eta_n$ (the
$O(\eta_n)$ term appears with sign in the coercive bound).

Third, the boundary sphere $\partial\Omega_R$ of the comparison is a
minimal surface with area at least $(1-\eps)A_{\min}(S)$ and charge
$(Q_E,Q_B)$. By Theorem~\ref{thm:charged-riemannian},
\begin{equation}\label{eq:ch-assembly}
m_{\rm ADM}(\bar g)\ge
\frac{r_{\bar g}+Q^2/r_{\bar g}}{2},\qquad
r_{\bar g}=\Bigl(\frac{\Area_{\bar g}(\partial\Omega_R)}{4\pi}\Bigr)^{1/2}.
\end{equation}
Pass to the limit as $\eps\downarrow0$ and $R\to\infty$: the left side
becomes $m$, the right side becomes the charged mass bound with
$r=(A_{\min}(S)/(4\pi))^{1/2}$.
\end{proof}

\begin{remark}
No charged equality classification is asserted. Equality in
\eqref{eq:ch-assembly} would require the classification of extremals of
the charged flow; that is out of scope.
\end{remark}

\begin{corollary}[Pure-electric rest frame]
\label{cor:ch:electric}
Under the hypotheses of Theorem~\ref{thm:ch:main}, if $P=0$ and $B=0$
(pure electric field), the conclusion \eqref{eq:ch-inequality} holds
with $m=E$ and the explicit cut convention of
Definition~\ref{def:bridge-admissible}.
\end{corollary}

\section{The four-dimensional case}
\label{sec:four}

The four-dimensional analysis follows the three-dimensional case with
two changes: the scalar identity in the comparison metric and the
regularity estimates gain the natural four-dimensional exponents. The
strictification, end reduction, two-color preparation, coupled system,
existence, and Riemannian Penrose steps all adapt; we present the
statement and the key modifications.

\subsection{Four-dimensional data}\label{sec:four-data}

A four-dimensional exterior $(\Omega,g,K)$ is a connected orientable
smooth $4$-manifold with Riemannian metric $g$ and symmetric
$(0,2)$-tensor $K$, complete with smooth compact boundary
$B=\partial\Omega$ and one asymptotically flat end in the sense of
\eqref{eq:decay} in a $4$-dimensional chart, with $q>2$ (the
four-dimensional analogue of the three-dimensional $q>3/2$). The
constraints and DEC are as in \eqref{eq:constraints}. The ADM vector
$(E,P)$ uses the standard $4$-dimensional normalization; when
$E>|P|$ the invariant mass is $m=\sqrt{E^2-|P|^2}$. Let $A_e=a_g(B)$
be the four-dimensional enclosing infimum in the sense of
Definition~\ref{def:bridge-admissible}.

\begin{theorem}[Four-dimensional exterior inequality]
\label{thm:four:main}
Let $(\Omega,g,K)$ be a four-dimensional exterior with
$\theta_+(B)\le0$, DEC, and integrable $\mu,|J|_g$. Let
$\omega_3=2\pi^2$. If $E>|P|$, then
\begin{equation}\label{eq:four-inequality}
m\ge\frac12\Bigl(\frac{A_e}{\omega_3}\Bigr)^{2/3}.
\end{equation}
The complete-data statement of Theorem~\ref{thm:bridge:global} holds in
dimension $4$ with the same right-hand side.
\end{theorem}

The proof follows Sections~\ref{sec:strictification}--\ref{sec:n3-strict-proof}
with the four modifications below.

\subsection{Strictification and end reduction in dimension four}
\label{sec:four-reduction}

The drift-Poisson strictification of Lemma~\ref{lem:n3:strictification}
adapts: in dimension $4$ the conformal family $g\mapsto u^2g$,
$K\mapsto uK$ preserves the constraints and DEC, and the drift
source exponent is $w=r^{-4-\delta}$ with $0<\delta<\beta<\min(q,1)$.
The repaired data satisfy $\theta_+(B)<0$, $\mu-|J|_g\ge c_0r^{-4-\delta_1}$
with $c_0>0$, and $K$ compactly supported. The end-reduction of
Proposition~\ref{prop:n3:rest-frame} works in dimension $4$: a distant
Lorentz boost sends $(E,P)\mapsto(m,0)$, and a conformal support
truncation on the far end adjusts $m$ and $A_e$ by $O(R^{-3})$.

\subsection{The four-dimensional elliptic construction}
\label{sec:four-elliptic}

The coupled system of Section~\ref{sec:elliptic} applies in dimension
$4$. The scalar identity changes: the lapse appears with the extra
degree of freedom $e^{3t}$ (rather than $e^{2t}$), and the conformal
family in the metric comparison is $u^3$ rather than $u^4$ with
$u=e^{t/3}$. The variables $f,t$, the auxiliary $D,d,Z,h,w,a,v,\chi$
take their dimension-$4$ forms. The comparison metric
$\bar g_{ij}=d^{-1}(g_{ij}+l(t)^2\partial_if\partial_jf)$ has
nonnegative scalar curvature on $\Omega_R$, mass at most $m+\eps$, and
boundary three-volume at least $(1-\eps)A_e$. The regularity estimates
of Section~\ref{sec:existence} (graph Sobolev, axial gradient estimate,
four-stage homotopy) pass through in dimension $4$ with the same
structure; the scalar identity uses a Schur-type trace calculation
$3\Delta_{\bar g}Z=\mathcal T$-type decomposition that gives
$R_{\bar g}\ge0$ from the DEC.

\subsection{The four-dimensional Riemannian Penrose inequality}
\label{sec:four-rmp}

We use the classical four-dimensional Riemannian Penrose inequality.

\begin{theorem}[Four-dimensional Riemannian Penrose]
\label{thm:four-riemannian}
Let $(M,g)$ be a complete smooth asymptotically flat four-dimensional
Riemannian manifold with nonnegative scalar curvature $R_g\ge0$ and
ADM mass $m_{ADM}$. Let $\Gamma$ be a smooth compact outermost minimal
$3$-surface. Then
\[
m_{ADM}\ge\frac12\Bigl(\frac{|\Gamma|_g}{\omega_3}\Bigr)^{2/3},
\qquad \omega_3=2\pi^2.
\]
\end{theorem}

\begin{proof}
This is the four-dimensional Riemannian Penrose inequality, a classical
result (Huisken--Ilmanen 2001; Bray 2001, extended by Bray--Lee 2009).
Its three-dimensional form is (T1) of Section~\ref{sec:external}; the
four-dimensional statement uses the flow identity
$\Sigma(t)=m_{ADM}-\tfrac12(A(t)/\omega_3)^{2/3}$ in place of the
three-dimensional formula. At the outermost minimal $3$-surface the
area statistic $\Area(t)$ equals $|\Gamma|_g$, giving the stated
inequality. We use it as an external input; the details of the flow
are classical.
\end{proof}

\subsection{Assembly}
\label{sec:four-assembly}

\begin{proof}[Proof of Theorem~\ref{thm:four:main}]
Apply strictification (Lemma~\ref{lem:n3:strictification} in
dimension $4$), end reduction (Proposition~\ref{prop:n3:rest-frame}
in dimension $4$), and the coupled elliptic construction
(Section~\ref{sec:four-elliptic}) to obtain comparison metric $\bar g$
with
\[
R_{\bar g}\ge0,\quad
m_{ADM}(\bar g)\le m+\eps,\quad
\Area_{\bar g}(\partial\Omega_R)\ge(1-\eps)A_e.
\]
Apply Theorem~\ref{thm:four-riemannian} to the outermost boundary
sphere:
\[
m_{ADM}(\bar g)\ge\frac12\Bigl(\frac{\Area_{\bar g}(\partial\Omega_R)}
{\omega_3}\Bigr)^{2/3}.
\]
Let $\eps\downarrow0$: the left side becomes $m$, the right side becomes
$\frac12(A_e/\omega_3)^{2/3}$, which is \eqref{eq:four-inequality}.
\end{proof}

\section{Four-dimensional rigidity}
\label{sec:four-rigidity}

The equality case of \eqref{eq:four-inequality} identifies the exterior
with Schwarzschild--Tangherlini.

\begin{theorem}[Four-dimensional rigidity]
\label{thm:rig4:rigidity}
Under the hypotheses of Theorem~\ref{thm:four:main} in the connected,
one-ended, outermost, outer area-minimizing equality case
$m=\frac12(A_e/\omega_3)^{2/3}$, the exterior $\Mext$ is isometric to
the future-horizon exterior in the Schwarzschild--Tangherlini
spacetime of mass $m$, with $K$ the second fundamental form generated
by the stationary Killing field of that spacetime.
\end{theorem}

\begin{proof}[Sketch]
The proof of Theorem~\ref{thm:n3:rigidity} applies verbatim with the
following substitutions: the causal stationary field
(Proposition~\ref{prop:n3:stationary}) extends to dimension $4$; the
twist-vanishing (Proposition~\ref{prop:n3-twist-vanishes}) still
uses the positive-mass bootstrap on the conformal double; the
conformal double (Proposition~\ref{prop:n3-double}) exists; and
Schoen--Yau positivity gives $\ell=1$ and $S\cong S^3$. The
static-vacuum classification in spacetime dimension $5$
(Bunting--Masood-ul-Alam) then identifies the exterior with
Schwarzschild--Tangherlini of the same mass.
\end{proof}

\section{The anti-de Sitter Penrose inequality}\label{sec:ads}

\subsection{Schwarzschild-anti-de Sitter and the branch}\label{sec:ads-bgd}

Fix the Schwarzschild-anti-de Sitter background of mass $m_0\ge0$ in
spacetime dimension $4$. In the usual static coordinates,
\begin{equation}\label{eq:ads-bgd}
g^0=-N^2\dd t^2+\frac{\dd r^2}{N^2}+r^2\sigma,\qquad
N^2=1+r^2-\frac{2m_0}{r},
\end{equation}
with $\sigma$ the round metric on $S^2$, and on the time slice $t=0$
we have $K^0=0$. The spatial metric in Fefferman-Graham form is
$g^0=\dd s^2+r(s)^2\sigma$ with $r=e^s(1+O(e^{-2s}))$; the boundary is
at $s=\infty$.

A smooth branch of asymptotically anti-de Sitter initial data near the
background is a smooth family $(g_\eps,K_\eps)$ for
$\eps\in[0,\eps_0)$ such that each $(g_\eps,K_\eps)$ is of finite
angular type, satisfies the Einstein constraints on the exterior of a
compact future MOTS $S_\eps$, satisfies $(g_\eps,K_\eps)=(g^0,K^0)+O(\eps)$
in the weighted Sobolev norm of the class, and has finite boundary mass
flux. Let $A(\eps)=|S_\eps|_{g_\eps}$ be the boundary area and
$m_{AH}(\eps)$ the positive Lorentz norm of the hyperbolic mass covector.

\begin{theorem}[Anti-de Sitter Penrose]
\label{thm:ads:main}
Under the hypotheses above there is $\eps_1\in(0,\eps_0]$ such that for
every $\eps\in[0,\eps_1)$,
\begin{equation}\label{eq:ads-inequality}
m_{AH}(\eps)
\ge\sqrt{\frac{A(\eps)}{16\pi}}
\Bigl(1+\frac{A(\eps)}{4\pi}\Bigr).
\end{equation}
Equality holds if and only if the branch seed is radial and $\eps\ge0$
is sufficiently small. The equality family is Schwarzschild-anti-de
Sitter at the matching mass and area.
\end{theorem}

\subsection{Taylor structure and weighted remainders}
\label{sec:ads-taylor}

Write the branch in a Taylor expansion
$g_\eps=g^0+\eps g_1+\eps^2 g_2+\eps^3 g_3+O_4(\eps^4)$ with each $g_j$
of finite angular type, and likewise $K_\eps$. The weighted remainder
$O_4(\eps^4)$ controls the asymptotic functionals.

\begin{lemma}[Taylor to mass and area]\label{lem:ads-taylor}
There are real $m_0,m_1,m_2,m_3$ and $A_0,A_1,A_2,A_3$ such that
\begin{equation}\label{eq:ads-taylor-exp}
m_{AH}(\eps)=\sum_{j=0}^3\eps^jm_j+O(\eps^4),\quad
A(\eps)=\sum_{j=0}^3\eps^j A_j+O(\eps^4),
\end{equation}
with the remainders controlled by the weighted bound
$F=O(e^{-\gamma s})$ with $\gamma>3$.
\end{lemma}

\begin{lemma}[Inverse-asymptotic]\label{lem:ads-inverse}
Let $L=-\Delta_{g^0}+3$. If $F$ is smooth finite-type with
$F=O(e^{-\gamma s})$ for some $\gamma>3$, then for either Dirichlet or
homogeneous Neumann data on $S$ there is a unique finite-type decaying
solution of $Lu=F$ with
\begin{equation}\label{eq:ads-remainder}
u=r^{-3}\zeta(x)+o(r^{-3}),\quad
\partial_s\bigl(u-r^{-3}\zeta(x)\bigr)=o(r^{-3}),
\end{equation}
$\zeta$ finite-type. Differentiation of $F$ in $s$ through order $k$
gives the same statement through $k+2$ derivatives of $u$.
\end{lemma}

\begin{proof}[Sketch]
Variation of constants for $y''+2y'-3y=G$ with
$G=O(e^{-\gamma's})$, $\gamma'=\min(\gamma,\beta+2)>3$, eliminates
the growing homogeneous mode $e^s$ and gives
$y=ce^{-3s}+O(e^{-\gamma's})$. Angular decomposition closes by
boundedness on each fixed finite-dimensional space. Comparison with
$Ah+\delta$, letting $\delta\downarrow0$, gives uniqueness of the
decaying solution.
\end{proof}

\subsection{The quadratic coefficient and its kernel}\label{sec:ads-quadratic}

\begin{proposition}[The quadratic form]
\label{prop:ads-c2}
In the expansion
\begin{equation}\label{eq:ads-mass-gap}
m_{AH}(\eps)-\sqrt{\frac{A(\eps)}{16\pi}}
\Bigl(1+\frac{A(\eps)}{4\pi}\Bigr)
=c_2[q]\eps^2+c_3[q]\eps^3+c_4[q]\eps^4+O(\eps^5)
\end{equation}
the coefficient $c_2[q]$ is a nonnegative quadratic form, and its
kernel has dimension exactly $4$ (the four linear harmonics on $S^2$).
\end{proposition}

\begin{proof}[Sketch]
Conformal decomposition of $g_1$ gives a TT part $k$ together with a
conformal factor and a tangential vector. The mass and area Taylor
expansion reduces $c_2$ to the TT square identity
(Lemma~\ref{lem:ads-tt-square}), which makes it
$\frac{1}{16\pi}\int_M N|k|^2\,\dd V_0
-\frac{\kappa}{16\pi}\int_S vDv\,\dd A_0$. The argument of
Lemma~\ref{lem:ads-tt-square} identifies the last term, and $D$'s
positivity makes the form nonnegative. Its kernel is the four
translation modes of the round $S^2$.
\end{proof}

\begin{lemma}[Weighted TT square identity]\label{lem:ads-tt-square}
Let $k$ be smooth TT with $k=O(r^{-\gamma})$ for some $\gamma>3/2$, and
$v$ the unique decaying solution of $Lv=0$ with
$v|_S=D^{-1}(k_{ss}|_S)$, $D=-\Delta_S+2\kappa/a$. Then
$Q(v)=\Hess v-vB$ is TT, $Q(v)_{ss}|_S=k_{ss}|_S$, and
\begin{equation}\label{eq:ads-tt-square}
\int_M N|k|^2\,\dd V_0-\kappa\int_S vDv\,\dd A_0
=\int_M N|k-Q(v)|^2\,\dd V_0.
\end{equation}
\end{lemma}

\subsection{The quartic argument on the kernel}\label{sec:ads-quartic}

\begin{proposition}[Quartic positivity]\label{prop:ads-c4}
If $q$ is a nonzero kernel direction of $c_2$ and $q$ is not a radial
seed, then $c_4[q]>0$.
\end{proposition}

\begin{proof}[Sketch]
The kernel is the linear harmonics, so $q(x)=V(s)x$ for a decaying $V$.
A change of spacetime slice at $2$-jet order in the Taylor expansions
(Lemma~\ref{lem:ads-slice}) removes the first extrinsic-curvature
coefficient of the nonradial directions, preserves the boundary area
through the required order, and reduces $c_4$ to a second TT square
built from the transverse coefficient near the original boundary. This
gives $c_4[q]\ge0$. Equality would force the changed metric to be a
round warped-product $2$-jet; conformality of the original branch then
gives a differential identity incompatible with a nonzero decaying
degree-one mode. Hence $c_4[q]>0$.
\end{proof}

\begin{lemma}[Finite-jet slice change]\label{lem:ads-slice}
There is a smooth local diffeomorphism $\Psi$ and a smooth conformal
factor $\omega$ with $\omega_0=1,\Psi_0=\mathrm{id}$, such that
$\widetilde g=\omega_\eps\Psi_\eps^*g^0+O(\eps^3)$ holds on the
original slice.
\end{lemma}

\begin{proposition}[Exact radial identity]\label{prop:ads-radial}
If $q$ is a radial seed, then for all sufficiently small $\eps\ge0$,
\begin{equation}\label{eq:ads-radial-equality}
m_{AH}(\eps)=\sqrt{A(\eps)/16\pi}
\bigl(1+A(\eps)/4\pi\bigr).
\end{equation}
\end{proposition}

\begin{proof}[Sketch]
Radial seeds form the Schwarzschild-anti-de Sitter family, for which
$m=\sqrt{A/16\pi}(1+A/4\pi)$ is exact at every radius. At $\eps=0$ the
identity is the background identity, and smoothness continues it to
$[0,\eps_1)$.
\end{proof}

\subsection{Proof of the anti-de Sitter inequality}

\begin{proof}
Combine Lemma~\ref{lem:ads-taylor} and
Propositions~\ref{prop:ads-c2}, \ref{prop:ads-c4}, \ref{prop:ads-radial}.
On nonradial seeds, the gap in \eqref{eq:ads-mass-gap} is
$c_2[q]\eps^2+c_4[q]\eps^4+O(\eps^5)$ with $c_2[q]\ge0$ and $c_4[q]>0$ on
the kernel, so the gap is strictly positive for all small $\eps>0$.
On radial seeds the exact identity \eqref{eq:ads-radial-equality} gives
equality. The interval $[0,\eps_1)$ depends only on the background data
and the branch.
\end{proof}

\begin{remark}
The result uses the same enclosing-area functional but in an
asymptotically hyperbolic background; the relevant mass is $m_{AH}$
rather than the asymptotically flat $m$, and the boundary area $A(\eps)$
plays the role of the enclosing area.
\end{remark}

\begin{thebibliography}{99}

\bibitem{ALDKK}
B.~Allen, E.~Bryden, D.~Kazaras, M.~Khuri,
\emph{Proof of the spacetime Penrose inequality with suboptimal
constant in the asymptotically flat and asymptotically hyperboloidal
regimes}, arXiv:2504.10641 (2025).

\bibitem{Bray2001}
H.~L.~Bray,
\emph{Proof of the Riemannian Penrose inequality using the positive
mass theorem}, J. Differential Geom. \textbf{59} (2001), 177--267.

\bibitem{BrayKhuri2011}
H.~L.~Bray, M.~Khuri,
\emph{P.~D.~E.'s which imply the Penrose conjecture},
Asian J. Math. \textbf{15} (2011), 557--610.

\bibitem{BuntingMasood}
G.~Bunting, A.~K.~M.~Masood-ul-Alam,
\emph{Nonexistence of multiple black holes in asymptotically Euclidean
static vacuum space-times},
Gen. Rel. Grav. \textbf{19} (1987), 147--154.

\bibitem{DisconziKhuri2012}
M.~Disconzi, M.~Khuri,
\emph{New asymptotic comparison for the ADM mass},
arXiv:1212.6257 (2012).

\bibitem{HuiskenIlmanen2001}
G.~Huisken, T.~Ilmanen,
\emph{The inverse mean curvature flow and the Riemannian Penrose
inequality}, J. Differential Geom. \textbf{59} (2001), 353--437.

\bibitem{Israel1967}
W.~Israel,
\emph{Event horizons in static vacuum space-times},
Phys. Rev. \textbf{164} (1967), 1776.

\bibitem{Lichnerowicz1955}
A.~Lichnerowicz,
\emph{Theories relativistes de la gravitation et de l'electromagnetisme},
Masson (1955).

\bibitem{Mars2009}
M.~Mars,
\emph{Present status of the Penrose inequality},
Class. Quant. Grav. \textbf{26} (2009), 193001.

\bibitem{Penrose1973}
R.~Penrose,
\emph{Naked singularities},
Ann. New York Acad. Sci. \textbf{224} (1973), 125--134.

\bibitem{SchoenYau1979}
R.~Schoen, S.-T.~Yau,
\emph{On the proof of the positive mass conjecture in general
relativity}, Comm. Math. Phys. \textbf{65} (1979), 45--76.

\bibitem{SchoenYau1981}
R.~Schoen, S.-T.~Yau,
\emph{Proof of the positive mass theorem II},
Comm. Math. Phys. \textbf{79} (1981), 231--260.

\bibitem{Tangherlini1963}
F.~R.~Tangherlini,
\emph{Schwarzschild field in $n$ dimensions and the dimensionality of
space problem}, Nuovo Cimento \textbf{27} (1963), 636--651.

\bibitem{Witten1981}
E.~Witten,
\emph{A new proof of the positive energy theorem},
Comm. Math. Phys. \textbf{80} (1981), 381--402.

\bibitem{Wald1984}
R.~M.~Wald,
\emph{General Relativity}, University of Chicago Press (1984).

\end{thebibliography}
\end{document}