---
abstract: 'We use the inverse mean curvature flow to prove a sharp Alexandrov-Fenchel-type inequality for a class of hypersurfaces in certain locally hyperbolic manifolds. As an application we derive an optimal Penrose inequality for asymptotically locally hyperbolic graphs in any dimension $n\geq 3$. When the horizon has the topology of a compact surface of genus at least one, this provides an affirmative answer, for this class of initial data sets, to a question posed by Gibbons, Chruściel and Simon on the validity of a Penrose-type inequality for black hole solutions carrying a higher genus horizon.'
address: 'Federal University of Ceará, Department of Mathematics, Campus do Pici, Av. Humberto Monte, s/n, Bloco 914, 60455-760, Fortaleza/CE, Brazil.'
author:
- Levi Lopes de Lima
- Frederico Girão
title: A Penrose inequality for asymptotically locally hyperbolic graphs
---

[^1]

Introduction and statements of the results {#intro}
==========================================

Penrose-type inequalities express the widely accepted view in General Relativity that a black hole should contribute to the total mass of an isolated gravitational system through its size as measured by the area of its event horizon [@BC] [@Ma]. In the asymptotically flat case, which corresponds to the cosmological constant $\Lambda=0$, this inequality has been established for time-symmetric initial data sets by Huisken-Ilmanen [@HI] and Bray [@Br] in the physical dimension $n=3$ and by Bray-Lee [@BL] for $n\leq 7$. Recent contributions to the graph and conformally flat cases in all dimensions appear in [@L] [@dLG1] [@dLG2] [@FS] [@HW] [@MV] for the ADM mass and in [@GWW1] [@GWW2] [@LWX] for the Chern-Gauss-Bonnet mass.

Recently, there has been much interest in extending this type of result to asymptotically hyperbolic (AH) manifolds, corresponding to the case $\Lambda <0$. In particular, the conjectured Penrose inequality has been established when the manifold can be appropriately embedded as a graph in hyperbolic space [@dLG3] [@dLG4]; see also [@DGS] for a previous contribution. The purpose of this note is to extend the method in [@dLG4] to asymptotically locally hyperbolic (ALH) graphs, thus proving an optimal Penrose inequality for such manifolds in any dimension $n\geq 3$; see Theorem \[theomain\]. In the physical dimension $n=3$, when the horizon has the topology of a compact surface of genus $\gamma\geq 1$, our result provides, for this class of initial data sets, an affirmative answer (Corollary \[cormain\]) to a question posed by Chruściel-Simon [@CS] and Gibbons [@G] on the validity of a Penrose-type inequality for black hole solutions carrying a higher genus horizon; see also [@BC] [@Ma] for background on this problem. A key ingredient in the proof of Theorem \[theomain\] is a sharp Alexandrov-Fenchel-type inequality for a class of hypersurfaces in certain locally hyperbolic manifolds; see Theorem \[aftype\]. This result has an independent interest and its proof uses the inverse mean curvature flow and an integral inequality due to Brendle [@Br]. In this regard we mention that lately there has been considerable progress in establishing Alexandrov-Fenchel-type inequalities in hyperbolic space; see for instance [@WX] and the references therein.

We now proceed to a detailed account of our main results. A mass-type invariant for ALH manifolds has been introduced by Chruściel, Herzlich and Nagy [@CH] [@CN] [@H]. Their definition applies in particular to the class of ALH manifolds we consider here. We start by describing the corresponding [*locally hyperbolic*]{} (LH) reference metrics. Fix $\epsilon =-1, 0$ and let $(N^{n-1},h)$ be a closed space form of sectional curvature $\epsilon$. In the product manifold $P_\epsilon=I_\epsilon\times N$, where $I_{-1}=(1,+\infty)$ and $I_0=(0,+\infty)$, define the metric $$\label{half}
g_\epsilon=\frac{dr^2}{\rho_{\epsilon}(r)^2}+r^2h, \quad r\in I_\epsilon,$$ where $$\rho_\epsilon (r)=\sqrt{r^2+\epsilon}.$$ It is easy to check that $(P_\epsilon,g_\epsilon)$ is [locally hyperbolic]{} in the sense that its sectional curvature is constant and equal to $-1$. For further reference we also note that the manifold $Q_\epsilon=\mathbb R\times P_\epsilon$, endowed with the metric $$\label{aswe}
\overline g_\epsilon=\rho_\epsilon^2d\tau^2+g_\epsilon,$$ is locally hyperbolic as well. Moreover, if $\tau_0\in\mathbb R$ then the horizontal slice $P_{\epsilon}^{\tau_0}\subset Q_\epsilon$ given by $\tau=\tau_0$ is totally geodesic. This follows from the fact that the vertical vector field $\partial_t$ is Killing. In particular, each $P_\epsilon^{\tau_0}$ can be naturally identified to $P_\epsilon$.

Let $(M^n,g)$ be a complete $n$-dimensional Riemannian manifold, possibly carrying a compact inner boundary $\Sigma$. For simplicity, we assume that $M$ has a unique end, say $E$. Following [@CH] we say that $(M,g)$ is [*asymptotically locally hyperbolic*]{} (ALH) if there exists a chart $\Psi$ taking $E$ to the end of $P_\epsilon$ corresponding to $r=+\infty$ so that, as $r\to +\infty$, $$\label{alhcond}
\|\Psi_*g-g_{\epsilon}\|_{g_\epsilon}+ \|d\Psi_*g\|_{g_\epsilon}=O\left(r^{-\sigma}\right),$$ for some $\sigma>n/2$. We also assume that $R_{\Psi_*g}+n(n-1)$ is integrable, where $R$ denotes scalar curvature. For this kind of manifold, a mass-like invariant can be defined as $$\label{massl}
\mathfrak m_{(M,g)}=\lim_{r\to +\infty}c_n\int_{N_r}
  \left(\rho_\epsilon({\rm div}_{g_{\epsilon}}e-d{\rm tr}_{g_{\epsilon}}e)-i_{\overline \nabla\rho_\epsilon}e+
  ({\rm tr}_{g_{\epsilon}}e)d\rho_\epsilon\right)(\nu_r)dN_r,$$ where $e=\Psi_*g-g_{\epsilon}$, $N_r=\{r\}\times N$, $\overline\nabla$ is the connection of $g_\epsilon$, $\nu_r$ is the [*outward*]{} unit vector to $N_r$, which means that $\nu_r$ points toward the end $r=+\infty$, and $$c_n=\frac{1}{2(n-1)\vartheta_{n-1}},$$ where $\vartheta_{n-1}={\rm area}_{n-1}(N)$. We emphasize that a crucial point in justifying (\[massl\]) is precisely to check that the right-hand side does not depend on the chart at infinity $\Psi$ satisfying (\[alhcond\]). In any case, $\mathfrak m_{(M,g)}$ is an invariant of the asymptotic geometry of $(M,g)$ which measures the average rate of the convergence of $g$ to $g_\epsilon$ as $r\to +\infty$; details can be found in [@CH].

\[epsilon1\] [The above construction also applies to the case in which $\epsilon=1$ and $(N,h)$ is a round sphere, so that a mass invariant is also available. This is the asymptotically hyperbolic case, treated in [@dLG4].]{}

An ALH manifold $(M,g)$ as above can be thought of as the initial data set of a time-symmetric solution of the Einstein field equations with negative cosmological constant. The invariant $\mathfrak m_{(M,g)}$ is then interpreted as the total mass of the solution. Physical reasoning predicts that $\mathfrak m_{(M,g)}$ should have the appropriate sign under the relevant dominant energy condition, namely, $R_g\geq -n(n-1)$. When $M$ carries a compact minimal horizon $\Sigma$ one expects the invariant to satisfy a Penrose-type inequality in the sense that it should be bounded from below by a suitable expression involving the area $|\Sigma|$ of $\Sigma$. In order to figure out the correct form of this inequality, we consider the so-called [*Kottler black hole metrics*]{}, which are rather straightforward deformations of the LH metrics $g_\epsilon$ above.

Let us introduce a real parameter $m>0$ and consider the metric $$\label{bhmet}
g_{\epsilon,m}=\frac{dr^2}{\rho_{\epsilon,m}(r)^2}+r^2h,$$ where $$\rho_{\epsilon,m}(r)=\sqrt{r^2+\epsilon-\frac{2m}{r^{n-2}}}.$$ For each $m$ as above, it is easy to see that the function $$r\mapsto f_{\epsilon,m}(r)=r^n+\epsilon r^{n-2}-2m$$ is strictly positive for $r>r_{\epsilon,m}$, where $r_{\epsilon, m}$ is the unique positive zero of $f_{\epsilon,m}$. A further analysis shows that the function $m\mapsto r_{\epsilon,m}$ is smooth and strictly increasing with range $(0,+\infty)$ if $\epsilon=0$ and $(1,+\infty)$ if $\epsilon=-1$.

This discussion shows that the metric $g_{\epsilon,m}$ is well defined on the product $P_{\epsilon,m}=I_{\epsilon,m}\times N$, where $I_{\epsilon,m}=\{r;r>r_{\epsilon,m}\}$. Moreover, it extends smoothly to the slice $r=r_{\epsilon,m}$, the so-called [*horizon*]{}, denoted $\Hc_{\epsilon,m}$. This terminology can be justified as follows. The metric $g_{\epsilon,m}$ is [*static*]{} in the sense that $$\label{static}
(\overline\Delta \rho_{\epsilon,m})g_{\epsilon,m}-\overline\nabla^2
\rho_{\epsilon,m}+\rho_{\epsilon.m}{\rm Ric}_{g_{\epsilon.m}}=0,$$ where $\overline \nabla$ is the connection of $g_{\epsilon,m}$ and $\overline \Delta$ is the corresponding Laplacian. This means that the Lorentzian metric $$\tilde g_{\epsilon,m}=-\rho_{m,\epsilon}^2d\tau^2+g_{\epsilon,m},$$ defined on $Q_{\epsilon,m}=\mathbb R\times P_{\epsilon,m}$, is a static solution to the vacuum field equations with negative cosmological constant: $${\rm Ric}_{\tilde g_{\epsilon,m}}=-n\tilde g_{\epsilon,m}.$$ Moreover, the null hypersurface $r=r_{\epsilon,m}$ defines the event horizon surrounding the central singularity $r=0$. This justifies the horizon terminology and explains why $g_{\epsilon,m}$ is termed a black hole metric. We note in passing that $\overline g_\epsilon$ in (\[aswe\]) is obtained from $\tilde g_{\epsilon,0}$ by Wick rotation: $\tau\mapsto i\tau$.

A computation shows that if $(\theta_1,\cdots,\theta_{n-1})$ are coordinates in $N_r$ then the sectional curvatures of $g_{\epsilon,m}$ are $$K(\partial_r,\partial_{{\theta}_i})=-1-(n-2)\frac{m}{r^n}$$ and $$K(\partial_{{\theta}_i},\partial_{{\theta}_j})=-1+\frac{2m}{r^{n-2}}.$$ This not only shows that $g_{\epsilon,m}$ satisfies the appropriate dominant energy condition, since its scalar curvature is $$R_{g_{\epsilon,m}} =  -n(n-1),$$ but also suggests that $g_{\epsilon,m}$ is ALH in the sense described above. In fact, a straightforward computation gives $$\|g_{\epsilon,m}-g_{\epsilon}\|_{g_\epsilon}+ \|dg_{\epsilon,m}\|_{g_\epsilon}=O\left(mr^{-n}\right),$$ as expected. Using (\[massl\]) we finally conclude that $\mathfrak m_{(P_{\epsilon,m},g_{\epsilon,m})}=m$, which shows that $m$ should be interpreted as the total mass of $g_{\epsilon,m}$.

One immediately finds that the area $|\Hc_{\epsilon,m}|$ of the horizon $\Hc_{\epsilon,m}$ of $(P_{\epsilon,m},g_{\epsilon,m})$ relates to its mass $m$ by means of $$m=\frac{1}{2}\left(\left(\frac{|\Hc_{\epsilon,m}|}{\vartheta_{n-1}}\right)^{\frac{n}{n-1}}
+\epsilon \left(\frac{|\Hc_{\epsilon,m}|}{\vartheta_{n-1}}\right)^{\frac{n-2}{n-1}}\right).$$ Thus, in analogy with the standard Penrose inequality, it is natural to conjecture that if $(M,g)$ is an $n$-dimensional ALH manifold (with respect to the reference metric $g_\epsilon$) carrying an outermost minimal horizon $\Sigma$ and satisfying $R_g\geq -n(n-1)$ everywhere then there holds $$\label{propcon}
\mathfrak m_{(M,g)}\geq\frac{1}{2}\left(\left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n}{n-1}}
+\epsilon \left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n-2}{n-1}}\right),$$ with the equality occurring if and only if $(M,g)$ is isometric to the corresponding black hole metric.

In the physical dimension $n=3$, (\[propcon\]) first appears in [@CS] as a conjectured inequality whose veracity would follow in case the use of the so-called Geroch’s monotonicity of the Hawking mass under the inverse mean curvature flow, as envisaged by Gibbons [@G], could be justified. Contrary to this initial expectation, Neves [@N] has shown that, at least in the AH case, the convergence properties of the flow at infinity are insufficient to implement Geroch’s scheme. Similar remarks should also apply in the general ALH context, even though Lee and Neves [@LN1] [@LN2] have recently established that Geroch’s strategy works in the so-called non-positive mass range; see Remark \[negmass\]. Despite these negative results, we show here that the method developed in [@dLG3] [@dLG4] can be adapted to handle the special case of graphs in [*any*]{} dimension $n\geq 3$. To explain this we observe that each $(P_{\epsilon,m},g_{\epsilon,m})$ can be isometrically immersed as a radially symmetric vertical graph inside $(Q_{\epsilon},\overline g_\epsilon)$: the defining function $u_{\epsilon,m}:I_{\epsilon,m}\to \mathbb R$ satisfies $u_{\epsilon,m}(r_{\epsilon,m})=0$ and $$\label{embed}
\rho_{\epsilon}(r)^2\left(\frac{du_{\epsilon,m}}{dr^2}\right)^2
=\frac{1}{\rho_{\epsilon,m}(r)^2}-\frac{1}{\rho_{\epsilon}(r)^2}.$$ It is clear from this that the horizon $\Hc_{\epsilon,m}$ lies on the totally geodesic horizontal slice $P_\epsilon^{0}$, with the intersection $M\cap P_\epsilon^{0}$ being orthogonal along $\Hc_{\epsilon,m}$. This motivates us to consider a more general class of hypersurfaces in $(Q_{\epsilon},\overline g_\epsilon)$.

\[alshyp\] We say that a complete hypersurface $(M,g)\subset (Q_{\epsilon},\overline g_\epsilon)$, possibly carrying a compact inner boundary $\Sigma$, is [*asymptotically locally hyperbolic*]{} if there exists a compact set $K\subset M$ so that $M-K$ can be written as a graph over the end $E_0$ of the horizontal slice $P_{\epsilon}^{0}\subset Q_{\epsilon}$, with the graph being associated to a smooth function $u$ such that the asymptotic condition (\[alhcond\]) holds for the nonparametric chart $\Psi_u(x,u(x))=x$, $x\in E_0$. Moreover, we assume that $R_{{\Psi_u}_*g}+n(n-1)$ is integrable.

Under these conditions, the mass $\mathfrak m_{(M,g)}$ can be computed by taking $\Psi=\Psi_u$ in (\[massl\]). More precisely, if we assume that the inner boundary $\Sigma$ lies on some totally geodesic, horizontal slice $P_{\epsilon}^{\tau_0}$, which we of course identify with $P_\epsilon$, and moreover that the intersection $M\cap P^{\tau_0}_\epsilon$ is orthogonal along $\Sigma$, so that $\Sigma\subset M$ is minimal and hence a horizon indeed, then the computations in [@dLG2] actually give the following integral formula for the mass: $$\label{massform}
\mathfrak m_{(M,g)}=c_n\int_M\Theta\left(R_g+n(n-1)\right) dM+
c_n\int_\Sigma \rho_{\epsilon} Hd\Sigma,$$ where $\Theta=\langle \partial/\partial t,N\rangle$, $N$ is the unit normal to $M$, which we choose so as to point upward at infinity, and $H$ is the mean curvature of $\Sigma\subset P_{\epsilon}$ with respect to its [*inward*]{} pointing unit normal, which means that the normal points in the direction [*opposite*]{} to the end of $P_{\epsilon}$ given by $r=+\infty$. In particular, if $R_g\geq -n(n-1)$ and $M$ is a graph ($\Theta>0$) then $$\label{massineq}
\mathfrak m_{(M,g)}\geq c_n\int_\Sigma \rho_{\epsilon} Hd\Sigma.$$

We are now in a position to use the following Alexandrov-Fenchel-type inequality.

\[aftype\] Let $\Sigma\subset P_\epsilon$ be a compact embedded hypersurface which is [*star-shaped*]{} in the sense that it can be written as a radial graph over a slice $N_r$ and [*strictly mean convex*]{} in the sense that its mean curvature satisfies $H>0$. Then there holds $$\label{aftype2}
c_n\int_\Sigma \rho_{\epsilon} Hd\Sigma\geq \frac{1}{2}\left(
\left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n}{n-1}}
+\epsilon \left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n-2}{n-1}}\right),$$ with the equality occurring if and only if $\Sigma$ is a slice.

By combining (\[massineq\]) and (\[aftype2\]) we immediately obtain the first part of our main result.

\[theomain\] If $M\subset Q_{0,\epsilon}$ is an ALH graph as above, so that its horizon $\Sigma\subset P^{\tau_0}_\epsilon$ is star-shaped and mean convex in the sense that $H\geq 0$, then $$\label{penrose}
\mathfrak m_{(M,g)}\geq \frac{1}{2}\left(
\left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n}{n-1}}
+\epsilon\left(\frac{|\Sigma|}{\vartheta_{n-1}}\right)^{\frac{n-2}{n-1}}\right),$$ with the equality holding if and only if $(M,g)$ is (congruent to) the graph realization of the corresponding black hole solution.

\[effec\] [Under the conditions of Theorem \[theomain\], the mass $\mathfrak m_{(M,g)}$ is always positive due to (\[massineq\]). In particular, if $\epsilon=-1$ the lower bounds (\[aftype2\]) and (\[penrose\]) only are effective if we further assume that $|\Sigma|>\vartheta_{n-1}$.]{}

As already observed, the following corollary provides a positive answer to a question posed by Gibbons [@G] and Chruściel-Simon [@CS] for the class of initial data sets we consider.

\[cormain\] If the horizon $\Sigma$ is a surface of genus $\gamma\geq 1$ then $$\label{haw}
\mathfrak m_{(M,g)}\geq \left(\frac{4\pi}{\vartheta_2}\right)^{{3}/{2}}
\sqrt{\frac{|\Sigma|}{16\pi}}\left(1-\gamma+\frac{|\Sigma|}{4\pi}\right),$$ where for $\gamma=1$ we assume the normalization $\vartheta_2=4\pi$. Moreover, the equality holds if and only if $(M,g)$ is (congruent to) the graph realization of the corresponding black hole solution.

If $\gamma\geq 2$ this follows by taking $n=3$ and $\epsilon=-1$ in the theorem and observing that Gauss-Bonnet gives $\vartheta_2=4\pi (\gamma-1)$. If $\gamma=1$ we take $\epsilon=0$ and use the normalization.

\[negmass\] [*As discussed in [@CS], when $\epsilon=-1$ the Kottler metrics (\[bhmet\]) also describe static black hole solutions when the parameter $m$ becomes negative up to a certain critical value, namely, $$m_{\rm crit}=-\frac{(n-2)^{\frac{n-2}{2}}}{n^{\frac{n}{2}}}.$$ In this regard we mention that Lee and Neves [@LN1] [@LN2] used the Huisken-Ilmanen’s formulation of the inverse mean curvature flow to establish a Penrose-type inequality for conformally compact ALH $3$-manifolds in this mass range. More precisely, they prove (\[haw\]) with the mass replaced by the supremum of the mass aspect function, which is assumed to be non-positive along the boundary at infinity. In particular, their manifolds always have [*non-positive*]{} mass while our graphs necessarily satisfy $\mathfrak m_{(M,g)}>0$; see Remark \[effec\]. Thus, their result and Corollary \[cormain\] are in a sense complementary to each other.*]{}

The proofs of Theorems \[aftype\] and \[theomain\] are presented in Section \[proof\]. This uses the asymptotics of solutions of the mean curvature flow established in Section \[imcflh\].

[**Acknowledgements:**]{} The authors would like to thank Prof. A. Neves for making available the preprint [@LN1].

The inverse mean curvature flow in LH manifolds {#imcflh}
===============================================

In this section we study the asymptotics of the inverse mean curvature for a class of hypersurfaces in $P_\epsilon$. Our main result, Theorem \[longtime\], follows from standard arguments in the theory of such flows; see for instance [@Ge1], [@Ge2] , [@BHW] and [@D]. For convenience we include all the relevant estimates here, even though our treatment is necessarily brief. The interested reader will find a more detailed account in the references above.

Hypersurfaces in LH manifolds {#hyperlh}
-----------------------------

It is convenient to consider a parameter $s$ so that $$ds=\frac{dr}{\rho_{\epsilon}(r)},$$ which gives $$s=\left\{
\begin{array}{cc}
\log r, & \epsilon =0\\
\log (2\sqrt{r^2-1}+2r), & \epsilon=-1
\end{array}
\right.$$ Thus, $s\in\mathbb R$ if $\epsilon=0$ and $s>\log 2$ if $\epsilon=-1$. Using this coordinate, it follows from (\[half\]) that $$\label{change}
g_{\epsilon}=ds^2+\lambda_\epsilon(s)^2h,$$ where $$\lambda_\epsilon(s)=\left\{
\begin{array}{cc}
e^s, & \epsilon =0\\
\frac{e^{s}}{4}+e^{-s}, & \epsilon=-1
\end{array}
\right.$$ We note the useful identities $$\label{usef}
\dot\lambda^2_\epsilon=\lambda^2_\epsilon+\epsilon, \quad \ddot\lambda_\epsilon=\lambda_\epsilon.$$ Also, in terms of $s$ we have $$\rho_\epsilon(s)=\left\{
\begin{array}{cc}
e^s, & \epsilon =0\\
\sqrt{\left(\frac{e^s}{4}+e^{-s}\right)^2-1}, & \epsilon=-1
\end{array}
\right.$$ This shows that $$\label{fund}
\dot\lambda_\epsilon=\rho_\epsilon.$$

We consider compact hypersurfaces $\Sigma\subset P_{\epsilon}$ which are [*star-shaped*]{} in the sense that they can be written as radial graphs over some slice $N_s=\{s\}\times N$. We choose the [*inward*]{} unit normal $\xi$ to $\Sigma$, which means that $\xi$ points toward the direction [*opposite*]{} to the end of $P_\epsilon$ defined by $s=+\infty$. Thus, if $\Sigma$ is graphically represented as $$\theta\in N\mapsto (u(\theta),\theta)\in P_\epsilon,$$ then $$\xi=\frac{1}{W}\left(\frac{u^i}{\lambda_\epsilon(u)^2}\partial_{\theta_i}
-\partial_s\right),\quad W=\sqrt{1+|D v|_h^2},$$ where $u_i=\partial_{\theta_i}u$, $u^i=h^{ij}u_j$, $D$ is the covariant derivative of $h$ and $$\label{change2}
v=\varphi(u), \quad \dot \varphi=1/\lambda_\epsilon.$$ We note that the tangent space to the graph is spanned by $$Z_i=\lambda_\epsilon v_i\partial_s+\partial_{\theta_i},\quad i=1,\cdots,n-1,$$ so that the induced metrics are $$\label{met1}
\eta_{ij}=\lambda_\epsilon^2(h_{ij}+v_iv_j)$$ and $$\eta^{ij}=\frac{1}{\lambda_{\epsilon}^{2}}\left(h^{ij}-\frac{v^iv^j}{W^2}\right).$$ The second fundamental form of the graph is $$c_{ij}=\frac{\lambda_\epsilon}{W}\left(\dot\lambda_\epsilon(h_{ij}+v_iv_j)-
v_{ij}\right),$$ so the shape operator is $$a^i_j=\eta^{ik}c_{kj}=\frac{\dot\lambda_\epsilon}{W\lambda_\epsilon}\delta^i_j-
\frac{1}{W\lambda_\epsilon}\tilde h^{ik}v_{kj},$$ where $$\tilde h^{ij}=h^{ij}-\frac{v^iv^j}{W^2}.$$ Thus, the mean curvature is $$\label{meancur}
H=a^i_i=(n-1)\frac{\dot\lambda_\epsilon}{W\lambda_\epsilon}-\frac{\tilde h^{ik}}{W\lambda_\epsilon}v_{ki}.$$ We recall that $$\label{cauchy}
|a|^2\geq \frac{H^2}{n-1},$$ with the equality holding if and only if $\Sigma$ is totally umbilical.

We also consider the [*support function*]{} $$\label{supp}
p=-\langle \overline \nabla\rho_\epsilon,\xi\rangle=\frac{\lambda_\epsilon}{W},$$ which relates to $\rho_\epsilon$ and $H$ by means of the following Minkowski-type identity: $$\label{minko}
\Delta \rho_\epsilon=(n-1)\rho_\epsilon -Hp,$$ where $\Delta={\rm div}\circ \nabla$ is the Laplacian of $\Sigma$. This is an easy consequence of the fact that, for $m=0$, (\[static\]) reduces to $$\label{confo}
\overline\nabla^2\rho_\epsilon=\rho_\epsilon g_\epsilon.$$

The inverse mean curvature flow {#subimcf}
-------------------------------

We consider a one-parameter family $$t\in [0,\delta)\mapsto \Sigma_t\subset P_\epsilon$$ of star-shaped, strictly mean convex hypersurfaces evolving under the inverse mean curvature flow. This means that the corresponding parametrization $X:[0,\delta)\times N\to P_\epsilon$ satisfies $$\label{imcfpar}
\frac{\partial X}{\partial t}=-\frac{\xi}{H}.$$ The following evolution equations are well-known [@D] [@dLG4].

\[evolequations\] Under the flow (\[imcfpar\]) we have:

1.  The area element $d\Sigma$ evolves as $$\label{evolvelemvol}
    \frac{\partial}{\partial t}d\Sigma=d\Sigma.$$ In particular, if $A$ is the area of $\Sigma$ then $$\label{evolvarea}
    \frac{dA}{dt}=A;$$

2.  The shape operator evolves as $$\label{evolvshape}
    \frac{\partial a}{\partial t}=\frac{\nabla^2H}{H^2}-2\frac{\nabla H\cdot\nabla H}{H^3}-\frac{a^2}{H}+\frac{I}{H},$$ where $\cdot$ is the symmetric product and $I$ is the identity map. As a consequence, the mean curvature evolves as $$\label{evolvmean}
    \frac{\partial H}{\partial t}=\frac{\Delta H}{H^2}-2
    \frac{|\nabla H|^2}{H^3} -\frac{|a|^2}{H}+\frac{n-1}{H},$$ and the tensor $b=Ha$ evolves as $$\label{simmons}
    \frac{\partial b}{\partial t} \approx \frac{\Delta b}{H^2}
     -2\frac{b^2}{H^2}+2(n-1)\frac{b}{H^2},$$ where $\approx$ means that first order terms in $b$ have been omitted in the right-hand side.

3.  The static potential evolves as $$\label{evolvstat}
     \frac{\partial \rho_\epsilon}{\partial t}=\frac{p}{H}.$$

4.  The support function evolves as $$\label{evolvsupp}
     \frac{\partial p}{\partial t}=\frac{\rho_\epsilon}{H}+\frac{1}{H^2}
     \langle \nabla \rho_\epsilon,\nabla H\rangle.$$ As a consequence of (\[minko\]), $$\label{evolvsuppint}
     \frac{d}{dt}\int_\Sigma p=n\int_\Sigma\frac{\rho_\epsilon}{H}d\Sigma.$$

Up to tangential diffeomorphisms, (\[imcfpar\]) means that $$\label{par}
\frac{\partial u}{\partial t}=\frac{W}{H},$$ where $u=u(t,\cdot)$ is the evolving graphing function. Equivalently, $$\label{imcfnonpar}
\frac{\partial v}{\partial t}=\frac{W}{\lambda_\epsilon H},$$ where $u$ and $v$ are related by (\[change2\]). It follows from (\[meancur\]) that (\[imcfnonpar\]) is a parabolic equation for $v$, so that a smooth solution always exists for any strictly mean convex initial hypersurface, at least for small $\delta$. We now develop the required apriori estimates to make sure that any such solution is defined for all time $t>0$. Our final result, Theorem \[longtime\], is far from providing the optimal asymptotic behavior of solutions but it suffices for our purposes. Our presentation follows [@BHW] and [@Ge1] closely; see also [@D]. As it is usually the case, a key ingredient here is the Parabolic Maximum Principle.

In the following, whenever $U$ is a time-dependent function we use the notation $U\sim e^{\frac{t}{n-1}}$ to mean that $U$ is bounded from above and below by a positive constant times $e^{\frac{t}{n-1}}$. Also, we set $$\underline U(t)=\inf_N U(t,\cdot), \quad \overline U(t)=\sup_N U(t,\cdot).$$ The next $C^0$ estimate is an immediate consequence of the fact that, in view of (\[meancur\]), the mean curvature is comparable to $
(n-1){\dot\lambda_\epsilon}/{W\lambda_\epsilon}
$ at an extremal point of $v$.

\[c0est\] As long as the solution $u$ of (\[par\]) exists, one has $$\label{c0estpre}
\lambda_\epsilon(u)\sim e^{\frac{t}{n-1}}.$$ More precisely, $$\label{c0est2}
e^{\frac{t}{n-1}}\lambda_\epsilon(\underline u(0))\leq \lambda_{\epsilon}(\underline u(t))\leq \lambda_{\epsilon}(\overline u(t)) \leq e^{\frac{t}{n-1}}\lambda_\epsilon(\overline u(0)).$$

We now turn to $C^1$ estimates. For this we first need to control the quantity determining the speed of the flow, namely, the mean curvature.

\[conth\] As long as the mean curvature remains positive, it is bounded from above as follows: $$\label{conth2}
H\leq n-1 + O(e^{-\frac{2t}{n-1}}).$$

It follows from (\[evolvmean\]) and (\[cauchy\]) that $$\frac{\partial H^2}{\partial t}\leq 2\frac{\Delta H}{H}-4
\frac{|\nabla H|^2}{H^2} -\frac{2}{n-1}{H^2}+2(n-1),$$ which implies that $\varrho=H^2-(n-1)^2$ satisfies $$\frac{d \overline \varrho}{dt}\leq - \frac{2}{n-1}\overline \varrho.$$ The result follows.

If we differentiate (\[imcfnonpar\]) with respect to $v^kD_k$ then we obtain that the quantity $$\label{quantomeg}
\omega=\frac{|D v|_h^2}{2}$$ evolves as $$\begin{aligned}
\frac{\partial \omega}{\partial t} & = & \frac{\tilde h^{ij}}{W^2F^2}\omega_{ij}-
  \frac{1}{F^2}\frac{\partial F}{\partial v_i}\omega_i-2(n-2)\epsilon\frac{\omega}{W^2F^2}\nonumber\\
  & & \quad -\frac{1}{W^2F^2}\tilde h_{ij}h^{kl}v_{ij}v_{jl}-2(n-1) H^{-2}\omega,
\label{evolvomega}\end{aligned}$$ where $$F=\frac{\lambda_\epsilon H}{W}.$$

\[apc1\] We have $$\label{apc12}
|D v|_h=O(e^{-\frac{t}{n-1}}),$$ or equivalently, $$\label{apc123}
|D u|_h=O(1).$$ In particular, $$\label{supernew}
W=1+O(e^{-\frac{2t}{n-1}}).$$

From (\[evolvomega\]), (\[conth2\]) and (\[c0est2\]), $$\frac{d\overline \omega}{dt}
   \leq  -\left(\frac{2}{n-1}+O(e^{-\frac{2t}{n-1}})\right)\overline \omega,$$ and the assertion follows.

The next proposition shows not only that strict mean convexity is preserved but also that the mean curvature remains bounded away from zero.

The mean curvature is uniformly bounded from below by a positive constant.

First note that (\[c0est2\]) and (\[apc12\]) imply $p \sim e^{\frac{t}{n-1}}$. Also, we have that $\kappa=\log (1/pH)$ satisfies $$\frac{\partial \kappa}{\partial t}=\frac{\Delta \kappa}{H^2}-\frac{|\nabla \log p^{-1}|^2}{H^2}+
\frac{|\nabla \log H|^2}{H^2} -\frac{n-1}{H^2},$$ so (\[conth2\]) gives $$\frac{d\overline \kappa}{dt}\leq -\frac{1}{n-1}+O(e^{-\frac{2t}{n-1}}).$$ Hence, $e^\kappa\leq O(e^{-\frac{t}{n-1}})$, which gives $H^{-1}\leq O(1)$.

Now we are ready to establish $C^2$ apriori estimates for the flow.

\[c2est\] The shape operator remains uniformly bounded as long as the flow exists.

Let $\mu$ be the maximal eigenvalue of $b=Ha$. Clearly, $|b|\leq C(\mu+1)$ for some constant $C>0$. On the other hand, (\[simmons\]) gives $$\frac{d\overline \mu}{dt}\leq -2\frac{\overline \mu^2}{H^2}+2(n-1)\frac{\overline \mu}{H^2},$$ so $\mu$ is uniformly bounded from above. Thus, $b$ is bounded and since $H$ is uniformly bounded, $a$ is bounded as well.

The following result establishes the longtime existence for solutions of the inverse mean curvature flow (\[imcfpar\]) and describes the asymptotics of the graphing function $u$ in a way that suffices for our purposes.

\[longtime\] For any initial hypersurface which is star-shaped and strictly mean convex, the solution of (\[imcfpar\]) exists for all time $t>0$. Also, the evolving hypersurface expands toward infinity, remaining star-shaped and strictly mean convex. Moreover, there exists $\alpha\in \mathbb R$ so that the rescaled function $$\tilde u=u-\frac{t}{n-1}$$ converges to $\alpha$ in the sense that $$\label{sharpest}
|D \tilde u|_h+|D^2 \tilde u|_h=o(1).$$

The longtime existence follows by applying Parabolic Regularity Theory [@Ge2] [@Li] in the usual way to the above apriori estimates. Also, these same estimates guarantee that star-shapedness and strict mean convexity are preserved along the flow. We now complete the asymptotic analysis by showing that the first and second order derivatives of $\tilde u$ actually vanish at infinity as indicated in (\[sharpest\]). A key ingredient here is a sharper control on the mean curvature.

\[vanish\] We have $$\label{vanishwh}
\frac{W}{H}= \frac{1}{n-1} + O(te^{-\frac{2t}{n-1}}).$$ In particular, $$\label{vanish2}
H=n-1 +O(te^{-\frac{2t}{n-1}}).$$

As explained in [@BHW], it follows from the evolution equation for the uniformly bounded quantity $\tau=W/H=\lambda_\epsilon/pH$ that $$\frac{d}{dt}\log \overline\tau\leq \frac{2}{n-1}-2\overline \tau +O(e^{-\frac{2t}{n-1}}),$$ so that $$\begin{aligned}
\frac{d\overline \tau}{dt} & \leq & \frac{2\overline \tau}{n-1}-2\overline \tau^2+O(e^{-\frac{2t}{n-1}})\\
&\leq & O(e^{-\frac{2t}{n-1}}),\end{aligned}$$ whenever $\tau\geq 1/(n-1)$. This gives $$\frac{W}{H}\leq \frac{1}{n-1}+ O(te^{-\frac{2t}{n-1}}),$$ and the reverse inequality follows from (\[conth2\]) and (\[supernew\]).

\[huwrem\] [Using the arguments in [@BHW] it can be shown that (\[vanish2\]) leads to a rather precise control on the shape operator, namely, $$|a^i_j-\delta^i_j|\leq O(t^2e^{-\frac{2t}{n-1}}).$$ In words, the evolving hypersurface becomes totally umbilical at an exponential rate.]{}

So far we have essentially followed [@BHW] and [@D]. Inspired by [@Ge1] we now introduce the uniformly bounded function $$\label{defw}
w=(e^u+1)e^{-\frac{t}{n-1}}=e^{\tilde u}+e^{-\frac{t}{n-1}}.$$

\[limit\] The functions $w$ and $\tilde u$ converge pointwisely as $t\to +\infty$.

By (\[defw\]) it suffices to prove the assertion for $w$. By (\[par\]) we get $$\label{evolvew}
\frac{\partial w}{\partial t}  =  \frac{W}{H}{e^u}e^{-\frac{t}{n-1}}-\frac{w}{n-1},$$ so that (\[vanishwh\]) and (\[c0est2\]) give $$\frac{\partial w}{\partial t}  =
   \frac{e^u}{n-1}e^{-\frac{t}{n-1}}-\frac{w}{n-1}+ O(te^{-\frac{2t}{n-1}}),$$ and hence, $$\label{evolvv2}
\frac{\partial w}{\partial t} =
 -\frac{1}{n-1}e^{-\frac{t}{n-1}}+O(te^{-\frac{2t}{n-1}}).$$ Thus, for all $t$ large enough there holds $$\frac{\partial w}{\partial t}\geq -\frac{2}{n-1}e^{-\frac{t}{n-1}},$$ or $$\frac{\partial}{\partial t}\left(w-e^{-\frac{2t}{n-1}}\right)\geq 0.$$ The assertion follows since $w=O(1)$.

We now prove that the limiting functions above are actually constant. This is the content of the next result, which also proves (\[sharpest\]) and concludes the proof of Theorem \[longtime\].

\[final\] There holds $$\label{estfinal}
|D w|_h+|D^2 w|_h=o(1).$$

First note that (\[change2\]), (\[meancur\]) and our previous estimates imply that (\[evolvew\]) is a [*uniformly*]{} parabolic equation for $w$. Also, (\[apc123\]) and (\[evolvv2\]) imply that both $\partial w/\partial t$ and $|Dw|_h$ decay exponentially. Since $w$ remains uniformly bounded in time, parabolic regularity [@Ge2] [@Li] implies that all higher order derivatives of $w$ are uniformly bounded as well. Standard interpolation inequalities then guarantee that $|D^2w|_h$ decays exponentially, which proves (\[estfinal\]).

The proofs of Theorems \[aftype\] and \[theomain\] {#proof}
==================================================

As remarked in the Introduction, the Penrose inequality (\[penrose\]) in Theorem \[theomain\] is an immediate consequence of the Alexandrov-Fenchel-type inequality (\[aftype2\]) in Theorem \[aftype\]. We now explain how the above results on the asymptotics of the inverse mean curvature flow may be used to prove (\[aftype2\]). A first ingredient is a Heintze-Karcher-type inequality due to Brendle [@Br].

\[bren\] If $\Sigma\subset P_\epsilon$ is star-shaped and strictly mean convex then $$\label{bren2}
(n-1)\int_\Sigma\frac{\dot\lambda_\epsilon}{H}d\Sigma\geq \int_\Sigma pd\Sigma.$$ Moreover, if the equality holds then $\Sigma$ is totally umbilical.

The static warped product $(P_\epsilon,g_\epsilon)$, with $g_\epsilon$ given by (\[change\]), has $\rho_\epsilon=\dot\lambda_\epsilon$ as its static potential and satisfies the general structure conditions in [@Br].

We now determine the asymptotic behavior of certain integral quantities associated to solutions of the flow (\[imcfpar\]). If $\Sigma\subset P_\epsilon$ is closed we define $$\Ic(\Sigma)=\int_\Sigma \rho_\epsilon Hd\Sigma=\int_\Sigma \dot\lambda_\epsilon Hd\Sigma,$$ $$\Jc(\Sigma)=\int_\Sigma pd\Sigma),$$ $$\Kc(\Sigma)=\vartheta_{n-1}\Ac(\Sigma)^{\frac{n}{n-1}},$$ where $\Ac(\Sigma)={|\Sigma|}/{\vartheta_{n-1}}$ is the normalized area of $\Sigma$, and $$\Lc(\Sigma)=\frac{\Ic(\Sigma)-(n-1)\Kc(\Sigma)}{\Ac(\Sigma)^{\frac{n-2}{n-1}}}.$$

\[asymbe\] If $\Sigma_t$ is a solution of (\[imcfpar\]) with $\Sigma_0$ star-shaped and strictly mean convex then $$\label{asymbe1}
\lim_{t\to +\infty}\frac{\Jc}{\Ac^{\frac{n}{n-1}}}(\Sigma_t)=\vartheta_{n-1},$$ and $$\label{asymbe2}
\liminf_{t\to +\infty}\Lc(\Sigma_t)\geq (n-1)\vartheta_{n-1}\epsilon.$$

First note that (\[sharpest\]) gives $$|D u|_h+|D^2u|_h=o(1),$$ and hence $$|D v|_h+|D^2v|_h=o(e^{-\frac{t}{n-1}}).$$ Thus, $$W^{-1}=1-\frac{1}{2}|D v|_h^2 +o(e^{-\frac{4t}{n-1}}),$$ so that (\[supp\]) gives $$p=\lambda_\epsilon-\frac{\lambda_\epsilon}{2}|D v|_h^2+o(e^{-\frac{3t}{n-1}}).$$ On the other hand, (\[met1\]) gives $$\label{expvol}
\sqrt{\det \eta}=\lambda_\epsilon^{n-1}\sqrt{\det h}\left(1+\frac{1}{2}|D v|_h^2+o(e^{-\frac{4t}{n-1}})\right).$$ We thus obtain the expansions $$\Jc(\Sigma_t)=\vartheta_{n-1}\fint\lambda_\epsilon^n+o(e^{\frac{(n-4)t}{n-1}}),$$ $$\label{exp2}
\Ac(\Sigma_t)=\fint \lambda_\epsilon^{n-1}+o(e^{\frac{(n-3)t}{n-1}})$$ and $$\Ac(\Sigma_t)^{\frac{n}{n-1}}=\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}+
o(e^{\frac{(n-2)t}{n-1}}),$$ where $
\fint
=\frac{1}{\vartheta_{n-1}}\int
$ and $\int$ means integration over $N$ with respect to $h$. In particular, $$\label{inp}
\frac{\Jc}{\Ac^{\frac{n}{n-1}}}=\frac{\vartheta_{n-1}\fint\lambda_\epsilon^n}
  {\left(\fint\lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}+o(e^{\frac{(n-2)t}{n-1}})}
  +o(e^{-\frac{4t}{n-1}}).$$ Using standard mean value theorems it is easy to check that $$\left|\frac{\fint\lambda_\epsilon^n}{\left(
 \fint\lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}}-1\right|\leq C\sup_N|D v|_h=o(1),$$ so that $$\label{same}
\lim_{t\to +\infty}\frac{\fint\lambda_\epsilon^n}{\left(
 \fint\lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}}=1,$$ and this together with (\[inp\]) proves (\[asymbe1\]).

We now prove (\[asymbe2\]). From (\[meancur\]) we get $$\dot\lambda_\epsilon H=(n-1)\frac{\dot\lambda_\epsilon^2}{\lambda_\epsilon}+o(e^{-\frac{t}{n-1}}),$$ and recalling (\[expvol\]) we obtain $$\int_\Sigma \rho_\epsilon H d\Sigma=(n-1)\int\dot\lambda_\epsilon^2\lambda_\epsilon^{n-2} +o(e^{\frac{(n-2)t}{n-1}}).$$ On the other hand, from (\[exp2\]) we have $$\Ac(\Sigma_t)^{\frac{n-2}{n-1}}=\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n-2}{n-1}}+o(e^{\frac{(n-4)t}{n-1}}).$$ Hence, from (\[usef\]) we get $$\begin{aligned}
\liminf_{t\to +\infty} \Lc(\Sigma_t) & = & (n-1)\vartheta_{n-1}\liminf_{t\to +\infty}
   \frac{\fint\dot\lambda_\epsilon^2\lambda_\epsilon^{n-2}-
     \left(\fint\lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}+o( e^{\frac{(n-2)t}{n-1}})}{\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n-2}{n-1}}+o(e^{\frac{(n-4)t}{n-1}})}\\
& \geq & (n-1)\vartheta_{n-1}\epsilon\liminf_{t\to +\infty}
    \frac{\fint\lambda_\epsilon^{n-2}}{\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n-2}{n-1}}+o(e^{\frac{(n-4)t}{n-1}})}+\\
    & & \quad + (n-1)\vartheta_{n-1}\liminf_{t\to +\infty}\frac{\fint \lambda_\epsilon^n-
     \left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}}
     {\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n-2}{n-1}}+o(e^{\frac{(n-4)t}{n-1}})}+\\
     & & \quad\quad +(n-1)\vartheta_{n-1}\liminf_{t\to +\infty}
     \frac{o(e^{\frac{n-2}{n-1}})}{\left(\fint \lambda_\epsilon^{n-1}\right)^{\frac{n-2}{n-1}}+o(e^{\frac{(n-4)t}{n-1}})}.\end{aligned}$$ Since $\lambda_\epsilon(u)\sim e^{\frac{t}{n-1}}$, the last limit on the right-hand side vanishes. Also, the argument leading to (\[same\]) shows that the first one is $1$. Finally, by Jensen’s inequality, $$\left(\fint\lambda_\epsilon^{n-1}\right)^{\frac{n}{n-1}}\leq\fint (\lambda_\epsilon^{n-1})^{\frac{n}{n-1}}=\fint\lambda_\epsilon^n,$$ which implies that the second limit is nonnegative. This proves (\[asymbe2\]).

We now identify monotone quantities along solutions of the inverse mean curvature flow.

\[mono\] Along solutions of (\[imcfpar\]) we have $$\label{mono1}
\frac{d}{dt}\left(\frac{\Jc-\Kc}{\Ac^{\frac{n}{n-1}}}\right)\geq 0,$$ and $$\label{mono2}
\frac{d\Lc}{dt}\leq 0.$$ Moreover, if the equality holds in any of these inequalities then $\Sigma$ is totally umbilical.

From (\[evolvarea\]), (\[evolvsuppint\]) and (\[bren2\]) we get $$\frac{d}{dt}\left(\Jc-\Kc\right)\geq \frac{n}{n-1}\left(\Jc-\Kc\right),$$ which immediately implies (\[mono1\]). In particular, if we let a star-shaped and strictly mean convex hypersurface $\Sigma\subset P_\epsilon$ flow under (\[imcfpar\]) and use (\[asymbe1\]) we deduce that $$\label{indep}
\Jc(\Sigma)\leq \Kc(\Sigma).$$ Moreover, if the equality holds then $\Sigma$ is totally umbilical.

On the other hand, a straightforward computation using Proposition \[evolequations\] gives $$\frac{d\Ic}{dt}=2\int_\Sigma\frac{\rho K}{H}d\Sigma +2\Jc,$$ where $K$ is the [*extrinsic*]{} scalar curvature of $\Sigma$, i.e. the sum of products of its principal curvatures. It then follows from Newton-McLaurin’s inequality that $$\frac{d\Ic}{dt}\leq \frac{n-2}{n-1}\Ic +2\Jc,$$ so that $$\begin{aligned}
\frac{d}{dt}\left(\Ic-(n-1)\Kc\right) & \leq & \frac{n-2}{n-1}(\Ic-(n-1)\Kc) +2(\Jc-\Kc)\\
  &\leq & \frac{n-2}{n-1}(\Ic-(n-1)\Kc),\end{aligned}$$ where we used (\[indep\]) in the last step. This easily yields (\[mono2\]) and completes the proof.

The proof of Theorem \[aftype\] now follows straightforwardly if we let $\Sigma$ flow under (\[imcfpar\]). Combining (\[mono2\]) and (\[asymbe2\]) we obtain $\Lc(\Sigma)\geq (n-1)\vartheta_{n-1}\epsilon$, which is just a rewriting of the Alexandrov-Fenchel-type inequality (\[aftype2\]). Moreover, if the equality holds then $\Sigma$ is totally umbilical. In particular, $\Sigma$ has constant mean curvature and results in [@Mo] imply that $\Sigma$ is a slice. This proves Theorem \[aftype\]. Finally, since any mean convex hypersurface can be arbitrarily approximated by strictly mean convex ones, the Penrose-type inequality (\[penrose\]) holds. In case of equality, it follows from (\[massform\]) that $\Sigma$ satisfies $$\label{ellip}
R_g=-n(n-1)$$ everywhere. As explained in [@dLG4], the argument in [@HW] can be adapted to show that $\Sigma$ is an elliptic solution of (\[ellip\]). Thus, Alexandrov reflection can be used in the vertical direction to make sure it is congruent to the graph realization of $(P_{\epsilon,m},g_{\epsilon,m})$ described in (\[embed\]), which is an elliptic solution as well. This completes the proof of Theorem \[theomain\].

[9999999]{}

H. L. Bray, Proof of the Riemannian Penrose inequality using the positive mass theorem. [*J. Differential Geom.*]{} 59 (2001), no. 2, 177-267.

H. L. Bray and D. L. Lee, On the Riemannian Penrose inequality in dimensions less than eight. [*Duke Math. J.*]{} 148 (2009), no. 1, 81-106.

S. Brendle, Constant mean curvature surfaces in warped product manifolds, [*arXiv:1105.4273*]{}, to appear in [Publ. Math. IHES]{}.

S. Brendle, P.-K. Hung, M.-T. Wang, A Minkowski-type inequality for hypersurfaces in the Anti-deSitter-Schwarzschild manifold, [*arXiv:1209.0669*]{}.

P. T. Chruściel, M. Herzlich, The mass of asymptotically hyperbolic Riemannian manifolds, [*Pacific J. Math.*]{} 212 (2003), 2, 231-26

H. L. Bray, P. T. Chruściel, The Penrose inequality, [*The Einstein equations and the large scale behavior of gravitational fields*]{}, 39-70, Birkhäuser, Basel, 2004.

Chruściel, P. T., Nagy, G., The mass of spacelike hypersurfaces in asymptotically anti-de Sitter space-times, [*Adv. Theor. Math. Phys.*]{} 5 (2001), 4, 697-754.

P. T. Chruściel, W. Simon, Towards the classification of static vacuum spacetimes with negative cosmological constant, [*J. Math. Phys.*]{} 42 (2001), no. 4, 1779-1817.

M. Dahl, R. Gicquaud, A. Sakovich, Penrose type inequalities for asymptotically hyperbolic graphs, [*arXiv:1201.3321*]{}.

L. L. de Lima, F. Girão, The ADM mass of asymptotically flat hypersurfaces, [*arXiv:1108.5474*]{}, to appear in Transactions of the AMS.

L. L. de Lima, F. Girão, A rigidity result for the graph case of the Penrose inequality, [*arXiv:1205.1132*]{}.

L. L. de Lima, F. Girão, Positive mass and Penrose type inequalities for asymptotically hyperbolic hypersurfaces, [*arXiv:1201.4991*]{}.

L. L. de Lima, F. Girão, An Alexandrov-Fenchel-type inequality in hyperbolic space with an application to a Penrose inequality, [*arXiv:1209.0438*]{}.

Q. Ding, The inverse mean curvature flow in rotationally symmetric spaces, [*Chin. Ann. Math. Ser. B*]{}, 32 (2011), no. 1, 27-44.

A. Freire, F. Scwartz, Mass-capacity inequalities for conformally flat manifolds with boundary, [*arXiv:1107.1407*]{}.

C. Gerhardt, Inverse curvature flows in hyperbolic space, [*J. Differential Geom.*]{} 89 (2011), no. 3, 487-527

C. Gerhardt, [*Curvature problems*]{}, Series in Geometry and Topology, 39. International Press, Somerville, MA, 2006.

Y. Ge, G. Wang, J. Wu, A new mass for asymptotically flat manifolds, [*arXiv:1211.3645*]{}.

Y. Ge, G. Wang, J. Wu, The Gauss-Bonnet-Chern mass of conformally flat manifolds, [*arXiv:1212.3213*]{}.

G. W. Gibbons, Some comments on gravitational entropy and the inverse mean curvature flow, [*Classical Quantum Gravity*]{} 16 (1999), no. 6, 1677-1687.

M. Herzlich, Mass formulae for asymptotically hyperbolic manifolds, [*AdS/CFT correspondence: Einstein metrics and their conformal boundaries*]{}, 103-121, Eur. Math. Soc., Zürich, 2005.

L.-H. Huang, D. Wu, The equality case of the Penrose inequality for asymptotically flat graphs, [*arXiv:1205.2061*]{}.

G. Huisken and T. Ilmanen, The inverse mean curvature flow and the Riemannian Penrose inequality, [*J. Differential Geom.*]{} 59 (2001), no. 3, 353-437.

M.-K. G. Lam, The Graphs Cases of the Riemannian Positive Mass and Penrose Inequalities in All Dimensions, [*arXiv:1010.4256*]{}.

D. Lee, A. Neves, Penrose inequality for higher genus Kottler metrics of nonpositive mass, [*preprint*]{}.

D. Lee, A. Neves, A static uniqueness theorem for higher genus Kottler metrics, slides of a talk delivered at Tsinghua Sanya International Mathematical Forum.

G. M. Lieberman, [*Second order parabolic differential equations*]{}, World Scientific Publishing Co., Inc., River Edge, NJ, 1996.

H. Li, Y. Wei, C. Xiong, The Gauss-Bonnet-Chern mass for graphic manifolds, [*arXiv:1304.6561*]{}.

M. Mars, Present status of the Penrose inequality, [*Classical Quantum Gravity*]{} 26 (2009), no. 19, 193001.

H. Mirandola, F. Vitório, The positive mass Theorem and Penrose inequality for graphical manifolds, [*arXiv:1304.3504*]{}.

S. Montiel, Unicity of constant mean curvature hypersurfaces in some Riemannian manifolds, [*Indiana Univ. Math. J.*]{} 48 (1999), no. 2, 711-748.

A. Neves, Insuficcient convergence of inverse mean curvature flow on asymptotically hyperbolic manifolds, [*J. Differential Geom.*]{} 84 (2010), no. 1, 191-229.

G. Wang, C. Xia, Isoperimetric type problems and Alexandrov-Fenchel type inequalities in the hyperbolic space, [*arXiv:1304.1674*]{}.

[^1]: The first author was partially supported by Fapesp/SP Grant 2012/09925-4
