On center of mass and foliations by constant spacetime mean curvature surfaces for isolated systems in General Relativity

TL;DR

This paper introduces a new foliation by constant spacetime mean curvature surfaces in asymptotically flat spacetimes, enabling a covariant definition of total center of mass that behaves like a point particle under Einstein evolution.

Contribution

It constructs a unique STCMC-foliation near infinity, providing a new, covariant way to define the total center of mass in General Relativity, improving upon previous CMC-based methods.

Findings

Existence of a unique STCMC-foliation near infinity

Definition of a new covariant total center of mass

Center of mass transforms under Poincaré group and evolves like a point particle

Abstract

We propose a new foliation of asymptotically Euclidean initial data sets by 2-spheres of constant spacetime mean curvature (STCMC). The leaves of the foliation have the STCMC-property regardless of the initial data set in which the foliation is constructed which asserts that there is a plethora of STCMC 2-spheres in a neighborhood of spatial infinity of any asymptotically flat spacetime. The STCMC-foliation can be understood as a covariant relativistic generalization of the CMC-foliation suggested by Huisken and Yau. We show that a unique STCMC-foliation exists near infinity of any asymptotically Euclidean initial data set with non-vanishing energy which allows for the definition of a new notion of total center of mass for isolated systems. This STCMC-center of mass transforms equivariantly under the asymptotic Poincar\'e group of the ambient spacetime and in particular evolves under the Einstein evolution equations like a point particle in Special Relativity. The new definition also remedies subtle deficiencies in the CMC-approach to defining the total center of mass suggested by Huisken and Yau which were described by Cederbaum and Nerz.

Equations

$$\mathrm{Scal}-|K{|}^2+(\mathrm{tr}K{)}^2$$

$$\mathrm{div}(K-(\mathrm{tr}K)g)$$

$$|g_{ij}-{\delta }_{ij}|+|\vec{x}||{\partial }_k{g}_{ij}|+|\vec{x}{|}^2|{\partial }_k{\partial }_l{g}_{ij}|$$

$$|K_{ij}|+|\vec{x}||{\partial }_k{K}_{ij}|$$

$$|\mu |+|J_i|$$

$$E$$

$$P^j$$

$$m$$

$$C_{\text{BOˊM}}^l$$

$$\vec{c}(\Sigma )$$

$$\vec{C}\text{raise}0.40903pt:=\underset{\sigma \to \infty }{\lim }\vec{c}({\Sigma }^{\sigma }),$$

$$\left|g_{ij}^{\mathrm{odd}}\right|+|\vec{x}|\left|{\partial }_k(g_{ij}^{\mathrm{odd}})\right|+|\vec{x}{|}^2\left|{\partial }_k{\partial }_l(g_{ij}^{\mathrm{odd}})\right|$$

$$|{\pi }_{ij}^{\mathrm{even}}|+|\vec{x}||{\partial }_k({\pi }_{ij}^{\mathrm{even}})|$$

$$|{\mu }^{\mathrm{odd}}|+|(J^j{)}^{\mathrm{odd}}|$$

$$f^{\mathrm{odd}}(\vec{x})\text{raise}0.40903pt:=\frac{1}{2}(f(\vec{x})-f(-x)),f^{\mathrm{even}}(\vec{x})\text{raise}0.40903pt:=\frac{1}{2}(f(\vec{x})+f(-\vec{x})).$$

$$\Vert f{\Vert }_{L^p(\Sigma )}\text{raise}0.40903pt:={\left({\int }_{\Sigma }|f{|}^{p}d\mu \right)}^{\frac{1}{p}}$$

$$\Vert f{\Vert }_{W^{0,p}(\Sigma )}\text{raise}0.40903pt:=\Vert f{\Vert }_{L^p(\Sigma )},\Vert f{\Vert }_{W^{k+1,p}(\Sigma )}\text{raise}0.40903pt:=\Vert f{\Vert }_{L^p(\Sigma )}+r{\Vert {\nabla }^{\Sigma }f\Vert }_{W^{k,p}(\Sigma )},k=0,1,\dots$$

$$\mathcal{H}=\sqrt{H^2-P^2}.$$

$$|K|\le C_{\mathcal{I}}|\vec{x}{|}^{-2}$$

$$Z^i\text{raise}0.40903pt:=\frac{1}{32\pi E}\underset{r\to \infty }{\lim }{\int }_{{\mathbb{S}}_r^2}\frac{x^i{\left({\sum }_{k,l}{\pi }_{kl}{x}^k{x}^l\right)}^2}{r^3}d{\mu }^{\delta }$$

$${\vec{C}}_{\text{STCMC}}={\vec{C}}_{\text{BOˊM}}+\vec{Z},$$

$$\begin{array}{r}{C}_{\text{STCMC}}^i=\frac{1}{16\pi E}\underset{r\to \infty }{\lim }[{\int }_{{\mathbb{S}}_r^2}\left(x^i\sum _{k,l}({\partial }_k{g}_{kl}-{\partial }_l{g}_{kk})\frac{x^l}{r}-\sum _k\left(g_{ki}\frac{x^k}{r}-g_{kk}\frac{x^i}{r}\right)\right)d{\mu }^{\delta }\\ +{\int }_{{\mathbb{S}}_r^2}\frac{x^i{\left({\sum }_{k,l}{\pi }_{kl}{x}^k{x}^l\right)}^2}{2r^3}d{\mu }^{\delta }],i=1,2,3.\end{array}$$

$$\mathcal{I}(t)=(\{t\}\times M^3,g(t),K(t),\mu (t),J(t))\hookrightarrow (\mathbb{R}\times M^3,\mathfrak{g})$$

$${\frac{d}{dt}|}_{t=0}{\vec{C}}_{STCMC}=\frac{\vec{P}}{E}.$$

$$2{\mu }_{\tau }$$

$$r\text{raise}0.40903pt:=\sqrt{\frac{|\Sigma {|}_g}{4\pi }},\text{ and }{z}^i\text{raise}0.40903pt:=\frac{1}{|\Sigma {|}_{\delta }}{\int }_{\Sigma }{x}^{i}d{\mu }^{\delta },i=1,2,3,$$

$$\Sigma \in \mathcal{A}(a,b,\eta ),$$

$$|\vec{z}|\le ar+b{r}^{1-\eta },r^{2+\eta }\le \underset{\Sigma }{\min }|\vec{x}{|}^{\frac{5}{2}+\epsilon },{\int }_{\Sigma }{H}^{2}d\mu -16\pi (1-\gamma )\le \frac{b}{r^{\eta }},$$

$$r^{-1}\Vert \mathring{A}{\Vert }_{W^{1,2}(\Sigma )}+\Vert \mathring{A}{\Vert }_{L^{\infty }(\Sigma )}\le C{r}^{-\frac{3}{2}-\epsilon }.$$

$$\Vert f{\Vert }_{W^{2,\infty }({\mathbb{S}}_r^2(\vec{z}))}\le C{r}^{\frac{1}{2}-\epsilon },$$

$$\Vert \psi -\mathrm{Id}{\Vert }_{W^{2,2}({\mathbb{S}}_r^2(\vec{z}))}\le C{r}^{\frac{3}{2}-\epsilon },$$

$$\Vert u^2-1{\Vert }_{W^{2,2}({\mathbb{S}}_r^2(\vec{z}))}\le C{r}^{\frac{1}{2}-\epsilon }.$$

$$(1-\overline{a})r-C{r}^{\max \{\frac{1}{2}-\epsilon ,1-\overline{\eta }\}}\le |\vec{x}|$$

Full text

We show that a unique STCMC-foliation exists near infinity of any asymptotically Euclidean initial data set with non-vanishing energy which allows for the definition of a new notion of total center of mass for isolated systems. This STCMC-center of mass transforms equivariantly under the asymptotic Poincaré group of the ambient spacetime and in particular evolves under the Einstein evolution equations like a point particle in Special Relativity. The new definition also remedies subtle deficiencies in the CMC-approach to defining the total center of mass suggested by Huisken and Yau [26] which were described by Cederbaum and Nerz [10].

1. Introduction and goals

In General Relativity, isolated (gravitating) systems are individual or clusters of stars, black holes, or galaxies that do not interact with any matter or gravitational radiation outside the system under consideration. Intuitively, they should have a total center of mass which should in a suitable sense behave as a point particle in Special Relativity. In this paper, we suggest a definition of total center of mass for suitably isolated systems and argue that this center of mass notion indeed behaves as a point particle in Special Relativity in a suitable sense (meaning it transforms equivariantly under the asymptotic Poincaré group of the ambient spacetime). In particular, we will show that the center of mass notion we suggest evolves in time under the Einstein evolution equations like a point particle in Special Relativity.

The main idea of our approach is to modify the definition of center of mass given by Huisken and Yau [26] for asymptotically Euclidean Riemannian manifolds — using an asymptotic foliation by $2$-spheres of constant mean curvature (CMC), see Section 2 — by staging it in a Lorentzian (spacetime) setting or in other words by staging it in asymptotically Euclidean initial data sets. More specifically, we will prove existence and uniqueness of an asymptotic foliation by $2$-spheres of constant spacetime mean curvature under optimal asymptotic decay assumptions. Here, “spacetime constant mean curvature (STCMC)” means that the co-dimension $2$ mean curvature vector $\vec{\mathcal{H}}$ of each $2$-sphere has constant Lorentzian length $\mathcal{H}$.

It is well known that this STCMC-condition can be reformulated in terms of initial data sets, namely as the product of the inner and outer “expansions” (or “null mean curvatures”, see Remark 5.2) with respect to any given null frame along a $2$-surface. On the other hand, the STCMC-condition is naturally independent of the initial data set in which the foliation is constructed. Our result thus asserts that there is a plethora of STCMC-surfaces in a neighborhood of spatial infinity of any asymptotically flat spacetime.

Furthermore, the new construction of a center of mass will be shown to remedy the subtle deficiencies of the Huisken and Yau approach [26] described by Cederbaum and Nerz [10]. Last but not least, we will provide an asymptotic flux integral formula for the center of mass extending that of Beig and Ó Murchadha [4]. The analytic techniques in our proofs rely on those developed by Metzger [31] and Nerz [36, 37].

Concluding this introduction, we would like to point out that the notion of spacetime mean curvature of $2$-surfaces in initial data sets has independently been considered in other contexts, both before and after the results of this paper had been announced. For example, the inverse spacetime mean curvature flow has been studied by Bray, Hayward, Mars, and Simon in [6] and by Hangjun Xu in his thesis [42].

The STCMC-condition is (trivially) satisfied by marginally outer/inner trapped surfaces (MOTS/MITS), extremal surfaces (see e.g. [21]), and generalized apparent horizons (see e.g. [30], [7]), with spacetime mean curvature $\mathcal{H}=0$ in all those cases. More generally, 2-surfaces with constant spacetime mean curvature are critical points for the area functional inside the future-directed null-cone, with mean curvature vector pointing in the direction in which the expansion of the surface is extremal. The aforementioned generalized apparent horizons have outer area minimizing property which is appealing in the view of spacetime Penrose Inequality. We would like to point the reader to the interesting work by Carrasco and Mars [9] giving insights into the (over-)generality of $\mathcal{H}=0$ as a condition for a horizon. In a recent paper of Cha and Khuri [11], the area $A$ of the outermost STCMC-surface with $\mathcal{H}=2$ appears in the conjectured Penrose Inequality $m\geq\sqrt{\nicefrac{{A}}{{16\pi}}}$ expected to hold for an asymptotically anti-de Sitter initial data set of mass $m$ satisfying the dominant energy condition.

Because of the spacetime geometry nature of the STCMC-condition, we expect that STCMC-surfaces and STCMC-foliations will have a number of applications beyond the definition of a center of mass of an isolated system as well as beyond the setting of asymptotically Euclidean initial data sets. For example, a special subfamily of STCMC-surfaces foliating a null hypersurface implicitly appears in recent work by Klainerman and Szeftel [28], where they arise as surfaces with both constant outer and constant inner expansion.

Structure of the paper.

In Section 2, we will summarize the necessary definitions and notations as well as more details on the background and existing work on the total center of mass of isolated systems. In Section 3, we will state our main results and very briefly explain the strategy of our proofs. The following sections will be dedicated to the more technical components of the proof with Section 4 focusing on a priori estimates for STCMC-surfaces, Section 5 discussing the linearization of spacetime mean curvature, Section 6 asserting existence of the STCMC-foliation, Section 7 introducing the coordinate expression of the center of mass associated with the STCMC-foliation, and Section 8 proving the claimed law of time evolution under the Einstein evolution equations. Appendix A collects results such as Sobolev Inequalities on $2$-surfaces, while Appendix B studies STCMC-surfaces in normal geodesic coordinates. Finally, in Section 9, we will discuss an exemplary initial data set highlighting the differences between the newly suggested notion of center of mass and the existing one suggested by Huisken and Yau.

Acknowledgements

We would like to thank Julien Cortier who was involved in this project at early stages and contributed with many insightful discussions.

The authors would like to extend thanks to the Erwin Schrödinger Institute, the Institut Henri Poincaré, the Mathematisches Forschungsinstitut Oberwolfach, and the Max Planck Institute for Mathematics for allowing us to collaborate in stimulating environments.

CC is indebted to the Baden-Württemberg Stiftung for the financial support of this research project by the Eliteprogramme for Postdocs. The work of CC is supported by the Institutional Strategy of the University of Tübingen (Deutsche Forschungsgemeinschaft, ZUK 63). CC also thanks the Fondation des Sciences Mathématiques de Paris for generous support. The work of AS was supported by Knut and Alice Wallenberg Foundation and Swedish Research Council (Vetenskapsrådet). We also thank the Crafoord Foundation for generous support.

2. Preliminaries

Recall that an initial data set for the Einstein equations is a tuple $(M^{3},g,K,\mu,J)$ where $(M^{3},g)$ is a smooth Riemannian manifold and $K$ is a smooth symmetric $(0,2)$-tensor field on $M^{3}$ playing the role of the second fundamental form of $M^{3}$ in an ambient Lorentzian spacetime. The (scalar) local energy density $\mu$ and the ($1$-form) local momentum density $J$ defined on $M^{3}$ can be read off from the constraint equations

[TABLE]

Here, $\operatorname{tr}$, $\operatorname{div}$, and $|\cdot|$ denote the trace, the divergence, and the tensor norm with respect to $g$, respectively, and $\mathrm{Scal}$ denotes its scalar curvature. Sometimes we will find it convenient to use the conjugate momentum tensor $\pi\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}(\operatorname{tr}K)g-K$.

The constraint equations (1) arise as a consequence of the Gauss–Codazzi–Mainardi equations from the Einstein equations $\mathfrak{{Ric}}-\tfrac{1}{2}\,\mathfrak{Scal}\,\mathfrak{g}=\mathfrak{T}$ satisfied by a given spacetime $(\mathfrak{M}^{1,3},\mathfrak{g})$ with energy-momentum tensor $\mathfrak{T}$, where $g$ is the Riemannian metric induced by the Lorentzian metric $\mathfrak{g}$ on the spacelike hypersurface $M^{3}$ and $K$ is the induced second fundamental form. Letting $\eta$ denote the timelike future unit normal to the initial data set $(M^{3},g,K,\mu,J)$, the energy and momentum density are derived from $\mathfrak{T}$ via $\mu=\mathfrak{T}(\eta,\eta)$, and $J=\mathfrak{T}(\eta,\cdot)$, and the stress tensor $S$ on $M^{3}$ is defined by $S=\mathfrak{T}(\cdot,\cdot)$. The constraint equations (1) thus necessarily hold on any spacelike hypersurface (or ”initial data set”) $(M^{3},g,K,\mu,J)$ in the spacetime $(\mathfrak{M}^{1,3},\mathfrak{g})$.

In order to model an ”isolated system”, we will assume that the ambient spacetime $(\mathfrak{M}^{1,3},\mathfrak{g})$ with its energy-momentum tensor $\mathfrak{T}$ and the choice of initial data set $(M^{3},g,K,\mu,J)$ are such that the initial data set is ”asymptotically Euclidean”, a notion made precise in the following standard definition.

Definition 2.1** (Asymptotically Euclidean initial data sets).**

Let $\varepsilon\in(0,\tfrac{1}{2}]$ and let $(M^{3},g,K,\mu,J)$ be a smooth initial data set. Assume there is a smooth coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ defined in the region exterior to a compact set $\mathcal{B}\subset M^{3}$. We say that $\mathcal{I}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}(M^{3},g,K,\mu,J)$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set (with respect to $\vec{x}$) if there is a constant $C=C(\mathcal{I},{\vec{x}})$ such that, in the coordinates $\vec{x}=(x^{1},x^{2},x^{3})\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$, we have the pointwise estimates

[TABLE]

for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and for all $i,j,k,l\in\{1,2,3\}$. Here, slightly abusing notation, we have silently pushed forward all tensor fields on $M^{3}$ (including scalars) and written $g_{ij}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}(\vec{x}_{*}g)_{ij}$ as well as $K_{ij}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}(\vec{x}_{*}K)_{ij}$, etc. The Kronecker delta $\delta_{ij}$ denotes the components of the Euclidean metric with respect to the coordinates $\vec{x}$. By another slight abuse of notation, we will refer to the above constant $C$ as $C_{\mathcal{I}}$, suppressing the dependence on the chart $\vec{x}$.

Asymptotically Euclidean initial data sets are well-known to possess well-defined total energy, linear momentum, and mass. More precisely, if $\mathcal{I}=(M,g,K,\mu,J)$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set for any $\varepsilon>0$ (naturally extending the definition to $\varepsilon>\tfrac{1}{2}$), its (ADM-)energy $E$ and its (ADM-)linear momentum $\vec{P}=(P^{1},P^{2},P^{3})$ are given by

[TABLE]

respectively, where $d\mu^{\delta}$ denotes the area measure induced on the coordinate sphere $\{{|\vec{x}\,|=r\}}$ by the Euclidean metric $\delta$ and ADM stands for Arnowitt–Deser–Misner [1]. The quantities $E$ and $\vec{P}$ are well-defined under the asymptotic conditions imposed here for arbitrary $\varepsilon>0$ [3, 16] — meaning the expressions converge and $E$ is asymptotically independent of the chart $\vec{x}$ while $\vec{P}$ is asymptotically covariant under chart deformations in a suitable way. From them, one defines the (ADM-)mass by

[TABLE]

whenever this expression makes sense, that is whenever the energy-momentum $4$-vector $(E,\vec{P})$ is causal with respect to the Minkowski metric of Special Relativity.

Remark 2.2* (Bounds on $\varepsilon$).*

For $\varepsilon\leq 0$ in the above definition, one can find an asymptotic chart $\vec{x}$ (meaning a coordinate transformation outside a compact set) on the canonical Euclidean initial data set $\mathcal{I}_{\text{Eucl.}}=(\mathbb{R}^{3},\delta,K\equiv 0,\mu\equiv 0,J\equiv 0)$ with respect to which $\mathcal{I}_{\text{Eucl.}}$ is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean but the expression $E$ does not vanish as it should for Euclidean space, see Denisov and Soloviev [19]. This explains the suggestive notation of the decay order as $\nicefrac{{1}}{{2}}+\varepsilon$.

On the other hand, if $\varepsilon>\tfrac{1}{2}$ for an initial data set $\mathcal{I}$, a simple computation shows that it has $E=\vec{P}=0$ which is non-desirable in the context of discussing the center of mass and asymptotic foliations by constant mean curvature. This explains why we exclude this case in Definition 2.1.

Remark 2.3* (Asymptotically Euclidean Riemannian manifolds).*

If a Riemannian manifold $(M^{3},g)$ (with asymptotic chart $\vec{x}$) satisfies (2a) and its scalar curvature satisfies $|\mathrm{Scal}|\leq C|\vec{x}\,|^{-3-\varepsilon}$ for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$, we say that $(M^{3},g)$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean manifold. This is called the “Riemannian case”, the reason being that one can reinterpret this as saying that the “trivially extended initial data set” $(M^{3},g,K\equiv 0,\mu=\tfrac{1}{2}\,\mathrm{Scal},J\equiv 0)$ satisfies (2a)-(2c). In the Riemannian case, the notions “mass” $m$ and “energy” $E$ can be and are used interchangeably.

2.1. Center of mass

We now proceed to discussing the total center of mass of an asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ with energy $E\neq 0$. The assumption $E\neq 0$ is both technical (as many definitions of center of mass explicitly divide by $E$) and physically reasonable when considering the center of mass.

First, let us remark that our field knows many definitions of center of mass for isolated systems. The first definitions were given in terms of asymptotic flux integral expressions in coordinates, similar to those of energy and linear momentum above, see (6) below and the text surrounding it. In 1996, Huisken and Yau [26] proved existence and uniqueness of a foliation by constant mean curvature $2$-spheres near infinity of an asymptotically Euclidean Riemannian manifold of positive energy $E>0$ and related it to a definition of center of mass in a way described below and in more detail in Section 7. More recently, Chen, Wang, and Yau [13] suggested a new definition of center of mass for isolated systems which is constructed from optimal isometric embeddings into the flat Minkowski spacetime of Special Relativity. For a brief, non-complete summary of other definitions of center of mass, please see [10].

Flux integral definition. The most prominent flux integral notion of center of mass $\vec{C}_{\text{B\'{O}M}}=(C^{1}_{\text{B\'{O}M}},C^{2}_{\text{B\'{O}M}},C^{3}_{\text{B\'{O}M}})$ for asymptotically Euclidean initial data sets was introduced by Beig and Ó Murchadha [4] as the asymptotic flux integral

[TABLE]

a definition going back in parts to Regge and Teitelboim [39]. See Szabados [41] for valuable critical comments on this definition, and see Section 7 for a covariant generalization of this formula following from our work.

The center of mass integral $\vec{C}_{\text{B\'{O}M}}$ will in general not converge for initial data sets $\mathcal{I}=(M^{3},g,K,\mu,J)$ which are merely $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean with respect to some chart $\vec{x}$ and have $E\neq 0$. It will however converge once one assumes that the initial data set satisfies certain asymptotic symmetry conditions in the given chart $\vec{x}$, as for example the Regge–Teitelboim conditions introduced in [39], see [4, 17, 24] and Definition 2.5 below. We also point out that the expression for $\vec{C}_{\text{B\'{O}M}}$ does not explicitly depend on the second fundamental form $K$ of the initial data set.

Definitions via foliations. Several authors define the center of mass of an initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ via a foliation by $2$-spheres near infinity. Following Cederbaum and Nerz [10], we will call such definitions “abstract” in contrast to the more explicit “coordinate definitions” of center of mass, see below.

The first abstract definition of center of mass was given in 1996 by Huisken and Yau [26], who proved existence and uniqueness of a foliation near the asymptotic end of an asymptotically Euclidean Riemannian manifold by closed, stable $2$-spheres of constant mean curvature, the CMC-foliation. This goes back to an idea of Christodoulou and Yau [15]. In 2006, Metzger [31] considered a foliation by $2$-spheres of constant null mean curvature (also called constant expansion) and concluded that this foliation is not fully suitable for defining a center of mass. For a more detailed review of foliations suggested to study in this context and of recent progress in terms of necessary and sufficient asymptotic decay conditions, please see [10].

Huisken and Yau [26] also assign a coordinate center to their foliation. It is constructed from the abstract CMC-center as a “Euclidean center” of the CMC-foliation as follows: First, any closed, oriented $2$-surface $\Sigma\hookrightarrow\mathbb{R}^{3}$ has a Euclidean coordinate center $\vec{c}\,(\Sigma)$ defined by

[TABLE]

Picking a fixed asymptotically flat coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$, this definition can naturally be extended to closed, oriented $2$-surfaces $\Sigma\hookrightarrow M^{3}\setminus\mathcal{B}$ by pushing $\Sigma$ forward to $\mathbb{R}^{3}$ and identifying $\vec{c}\,(\Sigma)\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\vec{c}\,(\vec{x}(\Sigma))$, slightly abusing notation. We will also call this center Euclidean center of $\Sigma$ (with respect to $\vec{x}$). This naturally extends to asymptotic foliations:

Definition 2.4** (Coordinate center of a foliation).**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+{\varepsilon}}$-asymptotically Euclidean initial data set for a chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$. Let $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ be a foliation of the asymptotic end $M^{3}\setminus\mathcal{B}$ of $M^{3}$ with area radius $r(\Sigma^{\sigma})=\sqrt{\nicefrac{{|\Sigma^{\sigma}|}}{{4\pi}}}$ of $\Sigma^{\sigma}$ diverging to $\infty$ as $\sigma\to\infty$. Denote by $\vec{c}\,(\Sigma^{\sigma})$ the Euclidean coordinate center of the leaf $\Sigma^{\sigma}$ with respect to $\vec{x}$. Then the (Euclidean) coordinate center $\vec{C}=(C^{1},C^{2},C^{3})$ of the foliation $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ (with respect to the asymptotic chart $\vec{x}$) is given by

[TABLE]

in case the limit exists. Otherwise, we say that the coordinate center of the foliation $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ diverges (with respect to the asymptotic chart $\vec{x}$).

The vector $\vec{C}$ can be pictured to describe a point in the target $\mathbb{R}^{3}$ of the asymptotically flat coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$, but it need not lie in the image of $\vec{x}$, and indeed will often lie inside $\overline{B_{R}(0)}$. This means it cannot necessarily be pulled back into $M^{3}$. Moreover, $\vec{C}$ depends on the choice of asymptotic chart $\vec{x}$ — at least a priori.

Coming back to the CMC-foliation constructed by Huisken and Yau [26], it is well-known that the coordinate center $\vec{C}_{\text{HY}}$ of the CMC-foliation of a suitably asymptotically flat Riemannian manifold $(M^{3},g)$ or initial data set $(M^{3},g,K,\mu,J)$ of non-vanishing energy $E$ with respect to a given asymptotic chart $\vec{x}$ coincides with the Beig–Ó Murchadha center of mass vector $\vec{C}_{\text{B\'{O}M}}$ defined by (6), provided that some additional symmetry assumptions are satisfied, see Huang [25], Eichmair and Metzger [20], and Nerz [34]. The most optimal result to date [36, Theorem 6.3] states that for $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean Riemannian manifolds with $E\neq 0$ satisfying the $C^{2}_{1+\varepsilon}$-Regge–Teitelboim condition (see Definition 2.5 below), we have $\vec{C}_{\text{HY}}=\vec{C}_{\text{B\'{O}M}}$ whenever both definitions converge, and that divergence of one implies divergence of the other. Again, let us point out that the construction of $\vec{C}_{HY}$ does not explicitly depend on the second fundamental form $K$ of the initial data set under consideration. We furthermore note that the product $E\,\vec{C}_{\text{B\'{O}M}}$ is sometimes referred to as the “center of mass charge” in the literature, even when $E=0$. We will not follow this usage here.

In this paper, we construct a novel geometric foliation $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ of the asymptotically flat end $M^{3}\setminus\mathcal{B}$ of a given $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $(M^{3},g,K,\mu,J)$ with non-vanishing energy $E\neq 0$, namely a foliation with “constant spacetime mean curvature (STCMC)”-leaves, see Section 3. The general approach to define the coordinate center of a foliation $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ described above will then be applied to this new foliation to obtain a new definition of the coordinate center of mass of an initial data set as well as a coordinate expression analogous to and extending (6), see Section 7.

2.2. Miscellannea

Here we collect some other definitions for future reference.

Regge–Teitelboim condition for initial data sets. With the exception of the later part of Section 7, we will not assume that the initial data sets under consideration satisfy any asymptotic symmetry assumptions, in particular the Regge–Teitelboim conditions. However, it will be useful in our discussion to refer to those conditions which is why we define them here.

Definition 2.5** (Regge–Teitelboim conditions for initial data sets).**

We say that a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ satisfies the $C^{2}_{\gamma+\varepsilon}$-Regge–Teitelboim conditions for $\gamma>\tfrac{1}{2}$ (with respect to the given chart $\vec{x}$ with respect to which it is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean) if there is a constant $C=C(\mathcal{I},{\vec{x}},\gamma)$ such that

[TABLE]

holds for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and for all $i,j,k,l\in\{1,2,3\}$. Here, as usual, we have denoted the even and odd parts of any continuous function $f\colon\mathbb{R}^{3}\setminus\overline{B_{R}(0)}\to\mathbb{R}$ by

[TABLE]

Remark 2.6* (Regge–Teitelboim conditions for Riemannian manifolds).*

We say that a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean Riemannian manifold $(M^{3},g)$ satisfies the $C^{2}_{\gamma+\varepsilon}$-(Riemannian) Regge–Teitelboim conditions on $\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ for $\gamma>\tfrac{1}{2}$ if the above inequalities are satisfied for $\pi\equiv K\equiv 0$, i.e. if (9a) holds and if $|\mathrm{Scal}^{\mathrm{odd}}|\leq C|\vec{x}\,|^{-\frac{5}{2}-\gamma-\varepsilon}$ for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$.

Weighted Sobolev spaces. In this paper, we use the following definition of Sobolev spaces, which is well-suited for keeping track of fall-off rates of different quantities associated with our foliation. Suppose that $(\Sigma,g^{\Sigma})$ is a closed (compact without boundary), oriented $2$-surface in an asymptotically Euclidean $3$-manifold $(M^{3},g)$ of suitable regularity. For $p\in[1,\infty)$, the Lebesgue space $L^{p}(\Sigma)$ is defined as the set of all measurable functions $f\colon\Sigma\to\mathbb{R}$ such that their $L^{p}$-norm

[TABLE]

is finite. Recall also that the $L^{\infty}$-norm of a measurable $f\colon\Sigma\to\mathbb{R}$ is defined by $\|f\|_{L^{\infty}(\Sigma)}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathrm{ess}\sup_{\Sigma}|f|$. Then for $p\in[1,\infty]$, the Sobolev norms are defined as follows:

[TABLE]

where $r\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\sqrt{\nicefrac{{|\Sigma|}}{{4\pi}}}$ is the area radius of $\Sigma$. The Sobolev space $W^{k,p}(\Sigma)$ is the set of all functions with finite $W^{k,p}$-norm. This definition naturally extends to the case of tensor fields on $\Sigma$. Appendix A in particular collects some Sobolev Inequalities for functions on $2$-surfaces $(\Sigma,g^{\Sigma})$ embedded in Euclidean space.

3. Main results, motivation, and the strategy of the proof

Given a $2$-dimensional surface $\Sigma$ in an initial data set $(M^{3},g,K,\mu,J)$, we denote its mean curvature inside the Riemannian manifold $(M^{3},g)$ with respect to the outward pointing unit normal111Please note that we use the convention for the sign of the second fundamental form ensuring that $H=2$ with respect to the outward pointing unit normal for the unit round sphere in $\mathbb{R}^{3}$. by $H$ and set $P\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\operatorname{tr}_{\Sigma}K$, as usual. The spacetime mean curvature (STMC) of $\Sigma$ is defined by the length of the spacetime mean curvature vector $\vec{\mathcal{H}}$

[TABLE]

We will suggestively write $\mathcal{H}^{\sigma}_{\tau}$ to denote the spacetime mean curvature of a surface called $\Sigma^{\sigma}_{\tau}$ etc, $\widetilde{\mathcal{H}}$ to denote the spacetime mean curvature of a surface called $\widetilde{\Sigma}$ etc., whenever the initial data set inducing the intrinsic and extrinsic geometry on the surface is clear from context.

In this paper we prove the following theorems.

Theorem 6.2** (Existence of STCMC-foliation).**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}}>0$ depending only on $\varepsilon$, $C_{\mathcal{I}}$, and $E$, a compact set $\mathcal{K}\subset M^{3}$, and a bijective $C^{1}$-map $\Psi\colon(\sigma_{\mathcal{I}},\infty)\times\mathbb{S}^{2}\to M^{3}\setminus\mathcal{K}$ such that each of the surfaces $\Sigma^{\sigma}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Psi(\sigma,\mathbb{S}^{2})$ has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma})\equiv\nicefrac{{2}}{{\sigma}}$ provided that $\sigma>\sigma_{\mathcal{I}}$.

Theorem 6.12** (Uniqueness of STCMC-foliation).**

Let $a\in[0,1)$, $b\geq 0$, and $\eta\in(0,1]$ be constants and let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}}$ depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$, and $E$, such that for all $\sigma>\sigma_{\mathcal{I}}$, there is a unique surface $\Sigma^{\sigma}\in\mathcal{A}(a,b,\eta)$ with constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to $\mathcal{I}$.

Here, $\mathcal{A}(a,b,\eta)$ is an a priori class of “asymptotically centered” spheres introduced in Section 4. It has been shown in particular by Brendle and Eichmair [8] that such an a priori condition is necessary to obtain uniqueness of CMC-surfaces in general, see the discussion in Section 6.3. As STCMC-surfaces generalize CMC-surfaces, their observation applies here, too.

We also obtain a coordinate expression $\vec{C}_{\text{STCMC}}$ for the STCMC-center of mass, see below. It differs from the Beig–Ó Murchadha formula $\vec{C}_{\text{B\'{O}M}}$ given in (6) by a term $\vec{Z}$, as stated in the following theorem.

Theorem 7.5** (STCMC-coordinate expression).**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with respect to an asymptotic coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and decay constant $C_{\mathcal{I}}$, with non-vanishing energy $E\neq 0$. Assume in addition that

[TABLE]

for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and that $g$ satisfies the Riemannian $C^{2}_{\nicefrac{{3}}{{2}}+\varepsilon}$-Regge–Teitelboim condition. Then the coordinate center $\vec{C}_{\text{STCMC}}$ of the unique foliation by surfaces of constant spacetime mean curvature is well-defined if and only if the correction term

[TABLE]

limits exist for $i=1,2,3$. In this case, we have

[TABLE]

where $\vec{C}_{\text{B\'{O}M}}$ is the Beig–Ó Murchadha center of mass and $\vec{Z}=(Z^{1},Z^{2},Z^{3})$, or equivalently

[TABLE]

An example of an initial data set with $\vec{C}_{\text{STCMC}}\neq\vec{C}_{\text{B\'{O}M}}$ or in other words with $\vec{Z}\neq 0$ will be discussed in Section 9. The above formula for $\vec{C}_{\text{STCMC}}$ allows to compute the STCMC-center of mass of an initial data set explicitly, once an asymptotic chart $\vec{x}$ has been picked. However, as the assumptions of Theorem 7.5 suggest, this formula cannot be expected to always converge. See Conjecture 7.9 and the text above of it for a discussion of when the coordinate expression for $\vec{C}_{\text{STCMC}}$ should converge, without reference to $\vec{C}_{\text{B\'{O}M}}$ and $\vec{Z}$ and without any Regge–Teitelboim conditions nor additional decay assumptions on $K$.

We get the following theorem on the time evolution of the STCMC-foliation and center of mass. The full covariance of the STCMC-foliation under the Poincaré group is discussed in Section 8.2.

Theorem 8.1** (Time evolution of STCMC-foliation).**

Let $(\mathbb{R}\times M^{3},\mathfrak{g})$ be a smooth, globally hyperbolic Lorentzian spacetime satisfying the Einstein equations with energy momentum tensor $\mathfrak{T}$. Suppose that, outside a set of the form $\mathbb{R}\times\mathcal{K}$, $\mathcal{K}\subset M^{3}$ compact, there is a diffeomorphism $\operatorname{Id}_{\mathbb{R}}\times\,\vec{x}\colon\mathbb{R}\times(M^{3}\setminus\mathcal{K})\to\mathbb{R}\times(\mathbb{R}^{3}\setminus\overline{B_{R}(0)})$ which gives rise to asymptotic coordinates $(t,\vec{x}\,)$ on $\mathbb{R}\times(M^{3}\setminus\mathcal{K})$.

Assume that $\mathcal{I}_{0}=(\{0\}\times M^{3},g,K,\mu,J)\hookrightarrow(\mathbb{R}\times M^{3},\mathfrak{g})$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with respect to the coordinate chart $\vec{x}$ and with $E\neq 0$, and suppose additionally that $K=O_{1}(|\vec{x}\,|^{-2})$ with constant $C_{\mathcal{I}}$ as $|\vec{x}\,|\to\infty$. Now consider the $C^{1}$-parametrized family of $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data sets

[TABLE]

with respect to $\vec{x}$ which starts from $\mathcal{I}(0)=\mathcal{I}_{0}$, and which exists for all $t\in(-t_{*},t_{*})$ for some $t_{*}>0$. Assume furthermore that the constants $C_{\mathcal{I}(t)}$ are uniformly bounded on $(-t_{*},t_{*})$, without loss of generality such that $C_{\mathcal{I}(t)}\leq C_{\mathcal{I}_{0}}$.

Assume the foliation $\mathcal{I}(t)$ has initial lapse $N=1+O_{2}(|\vec{x}\,|^{-\frac{1}{2}-\varepsilon})$ as $|\vec{x}\,|\to\infty$ with decay measuring constant denoted by $C_{N}$ and initial shift $X=0$, and suppose furthermore that the initial stress tensor $S$ of $\mathcal{I}_{0}$ satisfies $S=O(|\vec{x}\,|^{-\frac{5}{2}-\varepsilon})$ as $|\vec{x}\,|\to\infty$. There is a constant $\overline{t}>0$, depending only on $\varepsilon$, $C_{\mathcal{I}_{0}}$, $C_{N}$, and $E(0)$ such that the following holds: If the initial data set $\mathcal{I}_{0}$ has well-defined STCMC-center of mass $\vec{C}_{\text{STCMC}}\,(0)$ then the STCMC-center of mass $\vec{C}_{\text{STCMC}}\,(t)$ of $\mathcal{I}(t)$ is also well-defined for $|t|<\overline{t}$. Furthermore, the initial velocity at $t=0$ is given by

[TABLE]

Moreover, we have that $\left.\frac{d}{dt}\right|_{t=0}E=0$ and $\left.\frac{d}{dt}\right|_{t=0}\vec{P}=\vec{0}$.

3.1. Strategy of the proofs of Theorems 6.2 and 6.12

The underlying structure of the proofs of Theorems 6.2 and 6.12 presented in Section 6 and several of the lemmas proved in the same section is a method of continuity inspired by Metzger [31, 32]. Given an initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$, we will consider the one-parameter family of initial data sets $\mathcal{I}_{\tau}=(M^{3},g,\tau K,\mu_{\tau},\tau J)$, $\tau\in[0,1]$, with $\mu_{\tau}$ given through the constraint equations (1) as

[TABLE]

For $\tau=0$, we thus consider the Riemannian manifold $(M^{3},g)$ with $2\mu=2\mu_{0}=\mathrm{Scal}$ while for $\tau=1$, we study the original initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ with $\mu=\mu_{1}$. It is straightforward to see that if the original initial data set $\mathcal{I}$ is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean with respect to an asymptotical chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ then all initial data sets $\mathcal{I}_{\tau}$ are also $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean with respect to the same chart and comparable constants. In particular, the Riemannian manifold $(M^{3},g)$ is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean in this chart. This is what will allows us to drop the explicit mention of the chart in the proofs. Moreover, we note that the energy $E_{\tau}$ computed for the initial data set $\mathcal{I}_{\tau}=(M^{3},g,\tau K,\mu_{\tau},\tau J)$ does in fact not depend on $\tau$ and can and will thus be called $E$. We globally assume in this paper that $E\neq 0$ and we will fix the background Riemannian manifold $(M^{3},g)$ once and for all.

For second fundamental form $K=0$, the desired STCMC-foliation coincides with the classical CMC-foliation. From Nerz’ work [36], we thus know that the theorems and lemmas we will prove for initial data sets hold in the Riemannian setting under the Riemannian version of our assumptions, see also Remark 2.3. In other words, we know that our claims hold for $\tau=0$ in the method of continuity approach described above. In Section 6, we will recall Nerz’ corresponding theorems in our notation.

As usual, we will appeal to the Implicit Function Theorem in order to show openness of the interval in the method of continuity. Closedness follows from a standard convergence argument.

4. A priori estimates on STCMC-surfaces

When deforming the foliation by $2$-surfaces of constant mean curvature to the foliation by $2$-surfaces of constant spacetime mean curvature, we need to keep track of how the geometry of the leaves changes. For this, following [32] and [36], we will now introduce an a priori class of closed, oriented $2$-surfaces having the properties that their “area radius”, “coordinate radius”, and “mean curvature radius” as defined below are comparable in a certain sense.

In this section, we will not make explicit reference to the asymptotic coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ in most estimates, however the asymptotic coordinates $\vec{x}$ will be used in order to compute the coordinate radius and the center of mass of a given $2$-surface $\Sigma\hookrightarrow M^{3}$ (or “$\,\Sigma\hookrightarrow\mathcal{I}\,$”). We will always and mostly tacitly assume that $\Sigma\hookrightarrow M^{3}\setminus\mathcal{B}$ so that it lies in the domain of the asymptotic coordinate chart.

Definition 4.1**.**

Let $(M^{3},g)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically flat manifold with asymptotic coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$. Given any closed, oriented $2$-surface $\Sigma\hookrightarrow M^{3}\setminus\mathcal{B}$, we define its area radius $r=r(\Sigma)$ and (Euclidean) coordinate center $\vec{z}=\vec{z}\,(\Sigma)$, $\vec{z}=(z^{1},z^{2},z^{3})$, by

[TABLE]

respectively, where $d\mu^{\delta}$ denotes the area element on $\Sigma$ induced by the Euclidean metric $\delta$. Given constants $a\in[0,1)$, $b\geq 0$, and $\eta\in(0,1]$, we say that $\Sigma$ belongs to the a priori class of $(M^{3},g)$-asymptotically centered surfaces,

[TABLE]

if its area radius $r$, coordinate center $\vec{z}$, coordinate radius $|\vec{x}\,|$, and mean curvature $H$ satisfy the following estimates

[TABLE]

where $\gamma$ denotes the genus of $\Sigma$.

Remark 4.2*.*

We will use the same a priori classes in the context of asymptotically Euclidean initial data sets $\mathcal{I}=(M^{3},g,K,\mu,J)$, where the definition of $\mathcal{A}(a,b,\eta)$ only depends on the Riemannian manifold part $(M^{3},g)$. This will later be important when we consider families of initial data sets of the form $\mathcal{I}_{\tau}=(M^{3},g,\tau K,\mu_{\tau},\tau J)$, see Section 6 and (12).

Remark 4.3*.*

Note that for $r>1$, $0\leq a\leq\overline{a}<1$, $0\leq b\leq\overline{b}$, and $0<\overline{\eta}\leq\eta\leq 1$, we have $\mathcal{A}(a,b,\eta)\subseteq\mathcal{A}(\overline{a},\overline{b},\overline{\eta})$.

Example 4.4*.*

Let $(M^{3},g)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean manifold with non-vanishing energy $E\neq 0$. Then the unique leaves of the constant mean curvature foliation $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{0}}$ constructed in [36] are asymptotically centered in this sense. More specifically, there are constants $b>0$ and $\sigma_{0}>0$ depending only on $C_{\mathcal{I}}$ such that $\Sigma^{\sigma}\in\mathcal{A}(a=0,b,\eta=\varepsilon)$ for $\sigma>\sigma_{0}$. See [36, Section 5] for details.

Proposition 4.5**.**

Suppose that $\overline{a}\in[0,1)$, $\overline{b}\geq 0$, $\overline{\eta}\in(0,1]$, and assume that $0\leq a\leq\overline{a}$, $0\leq b\leq\overline{b}$, and $\overline{\eta}\leq\eta\leq 1$. Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set. Then there exist constants $\overline{\sigma}$ and $C$ depending only on $\varepsilon$, $\overline{a}$, $\overline{b}$, $\overline{\eta}$, and $C_{\mathcal{I}}$ such that the following a priori conclusions hold for any closed, oriented $2$-surface $\Sigma\hookrightarrow\mathcal{I}$ with $\Sigma\in\mathcal{A}(a,b,\eta)$: Suppose that $\Sigma$ has constant spacetime mean curvature $\mathcal{H}\equiv\nicefrac{{2}}{{\sigma}}$ in $\mathcal{I}$ for some $\sigma>\overline{\sigma}$. Then $\Sigma$ is a topological sphere and the tracefree part $\mathring{A}$ of its second fundamental form satisfies

[TABLE]

Furthermore, there exists a function $f\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\mathbb{R}$ such that $\Sigma$ is the graph of $f$ and

[TABLE]

as well as a conformal parametrization $\psi\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\Sigma$ which satisfies

[TABLE]

where $\operatorname{Id}$ denotes the trivial embedding $(\mathbb{S}^{2}_{r}(\vec{z}\,),g^{\mathbb{S}^{2}_{r}(\vec{z}\,)})\hookrightarrow(\mathbb{R}^{3},\delta)$. The conformal factor $u\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\mathbb{R}^{+}$ such that $\psi^{*}g^{\Sigma}=u^{2}g^{\mathbb{S}^{2}_{r}(\vec{z}\,)}$ satisfies

[TABLE]

Finally, the Euclidean distance to the coordinate origin $|\vec{x}\,|$ (on $\Sigma$), the area radius $r$, and the spacetime mean curvature radius $\sigma$ are comparable in the following sense:

[TABLE]

Remark 4.6*.*

The conclusions of this theorem are mostly the same as those in [36, Proposition 4.4], only for STCMC- rather than CMC-surfaces. However, we cannot directly refer to this result because, roughly speaking, it assumes that the mean curvature $H$ of $\Sigma$ falls off like $H-\tfrac{2}{\sigma}=O(r^{-\frac{3}{2}-\varepsilon})$, whereas the relation

[TABLE]

recalling $P=\operatorname{tr}_{\Sigma}K$, and the definition of the a priori class $\mathcal{A}(a,b,\eta)$ — which coincides with that in [36] —, only ensure via the second inequality in (15) that

[TABLE]

which does not a priori give us $H-\tfrac{2}{\sigma}=O(r^{-\frac{3}{2}-\varepsilon})$. We will thus need to extend the result and its proof to our setting.

Proof.

Within this proof, $C$ will always be a generic constant depending only on $\overline{\sigma}$, $\overline{a}$, $\overline{b}$, $\overline{\eta}$, and $C_{\mathcal{I}}$. With Remark 4.6 in mind, we need to improve the estimate in (22). For this purpose, we first note that by the definition of $(M^{3},g)$ being $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean and by the second inequality in (15), we have

[TABLE]

which implies $\|\mathrm{Scal}\|_{L^{1}(\Sigma)}\leq Cr^{-\eta}$. Similarly, with $\nu$ denoting the unit normal of $\Sigma$ in $(M^{3},g)$, we get $\|\mathrm{Ric}(\nu,\nu)\|_{L^{1}(\Sigma)}\leq Cr^{-\eta}$. Combining this with the last inequality of (15), we conclude by the Gauss equation and the Gauss–Bonnet Theorem that

[TABLE]

hence $\|\mathring{A}\|_{L^{2}(\Sigma)}\leq Cr^{-\frac{\eta}{2}}$. Then by Lemma A.1 we also have $\|\mathring{A}^{\delta}\|_{L^{2}(\Sigma,\delta^{\Sigma})}\leq Cr^{-\frac{\eta}{2}}$, where $\delta^{\Sigma}$ is the induced metric of the embedding $(\Sigma,\delta^{\Sigma})\hookrightarrow(\mathbb{R}^{3},\delta)$. We are now in a position to apply the result of De Lellis and Müller [18, Theorem 1.1] (see also [32, Section 2.3] where this result is reformulated in a scale invariant form) to conclude that $\Sigma$ is a topological sphere, with a conformal parametrization $\psi\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\Sigma$ and the conformal factor $u\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\mathbb{R}^{+}$ such that $\psi^{*}\delta^{\Sigma}=u^{2}g^{\mathbb{S}^{2}_{r}(\vec{z}\,)}$ satisfying

[TABLE]

In order to prove that $\sigma$ and $r$ are comparable, we estimate

[TABLE]

where we have used (2) and the second inequality in (15) in the last line. Here, we have

[TABLE]

by [18, Theorem 1.1] (see also (2.4) in [32]), and

[TABLE]

by Lemma A.1 combined with (15), and

[TABLE]

by (22).

Summing up, we conclude that

[TABLE]

To prove that $|\vec{x}\,|$ and $r$ are comparable, note that by the first inequality in (15) and because $0<1-\overline{a}\leq 1-a$, the coordinate origin $\vec{0}$ lies inside $\mathbb{S}^{2}_{r}(\vec{z}\,)$ for $r>\left(\frac{\overline{b}}{1-\overline{a}}\right)^{\frac{1}{\overline{\eta}}}$. For such large radii, we thus elementarily find

[TABLE]

with the help of the first inequality in (15). By (23a) and the Sobolev Inequality in the form of Lemma A.2, it follows that $|\psi-\operatorname{Id}|\leq C|\vec{x}\,|^{1-\frac{\eta}{2}}$. Combining this with the above inequalities, we conclude that, on $\Sigma$, we have

[TABLE]

provided that the area radius $r$ of $\Sigma$ is sufficiently large. Via (24), we can alternatively state that (25) holds if the spacetime mean curvature radius $\sigma$ satsifies $\sigma>\overline{\sigma}$ for a suitably large $\overline{\sigma}$ only depending on $\varepsilon$, $\overline{a}$, $\overline{b}$, $\overline{\eta}$, and $C_{\mathcal{I}}$.

Bootstrapping.

With these new bounds (24) and (25) at hand, we can apply [36, Proposition 4.1] with $\kappa$ chosen as $\tfrac{3}{2}+\varepsilon>1$, $\eta$ chosen as our $\tfrac{\eta}{2}>0$, and $c_{1},c_{2}$ chosen as our generic constant $C$. As all the estimates going into verifying the assumptions from [36, Proposition 4.1] hold pointwise in our case, the assumptions are indeed satisfied for any $p>2$. Note that the existence of the uniform Sobolev Inequality assumed in [36, Proposition 4.1] is well-established in our setting, and goes back to [26, Proposition 5.4] which holds for surfaces in asymptotically Euclidean manifolds with general asymptotics as described in Section 2. Again via (24), this gives us (16) for $\sigma>\overline{\sigma}$, with suitably enlarged $\overline{\sigma}$ only depending on $\varepsilon$, $\overline{a}$, $\overline{b}$, $\overline{\eta}$, and $C_{\mathcal{I}}$.

As a consequence of (16), the estimates (23a)-(23b) improve, and we get (18) and (19). Similarly, repeating the above argument that we used to derive (24) and (25), we obtain the improved radius comparison (20) and (21).

Finally, now that we have a pointwise bound on the tracefree part of the second fundamental form $\mathring{A}$ accompanying the pointwise estimate (22) for the mean curvature $H$, it follows that $\Sigma$ is the graph of a function $f\in W^{2,\infty}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ such that (17) holds for $\sigma>\overline{\sigma}$, for again suitably enlarged $\overline{\sigma}$ only depending on $\varepsilon$, $\overline{a}$, $\overline{b}$, $\overline{\eta}$, and $C_{\mathcal{I}}$, see e.g. [36, Corollary E.1], which adapts [18, Theorem 1.1] to our setting. To be more precise, [36, Corollary E.1] is only not stated invariantly under scaling but with $|\Sigma|=4\pi$, but it is straightforward to adapt it to include the area radius for our purposes. This finishes the proof of Proposition 4.5. ∎

5. The linearization of spacetime mean curvature

In this section, we will introduce the spacetime mean curvature map $\mathcal{H}$ in a given initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$. We will analyze its properties in a neighborhood of a given $2$-surface $\Sigma$ having constant spacetime mean curvature. We will show that the linearization of the map $\mathcal{H}$ is invertible when the linearization is computed with respect to normal variations within the given initial data set $\mathcal{I}$. This will later be used to ensure that the CMC-foliation of $(M^{3},g)$ constructed in [36] can be pushed via a method of continuity to an STCMC-foliation of $\mathcal{I}$.

Throughout this section, we will assume that $\mathcal{I}=(M^{3},g,K,\mu,J)$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$ and with fixed asymptotic coordinates $\vec{x}$. Furthermore, it will be assumed that $\Sigma$ is a fixed $2$-surface of constant spacetime mean curvature $\mathcal{H}(\Sigma)\equiv\nicefrac{{2}}{{\sigma}}$ which has sufficiently large mean curvature radius $\sigma$ and which for some fixed $a\in[0,1)$, $b\geq 0$, and $\eta\in(0,1]$ belongs to the a priori class $\mathcal{A}(a,b,\eta)$, see Definition 4.1. In this setting, we know from Proposition 4.5 that $\Sigma$ is a topological sphere, and that its coordinate radius, area radius, and mean curvature radius are comparable as stated in (20), (21). This in particular implies that

[TABLE]

for $\sigma>\overline{\sigma}$, where $\overline{\sigma}$ and the constants hidden in the $O$-notation only depend on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$.

5.1. Stability operators associated with prescribed (spacetime) mean curvature surfaces

In a neighborhood of $\Sigma$, we introduce normal geodesic coordinates $y\colon\Sigma\times(-\xi,\xi)\to M^{3}$ for some $\xi>0$, such that $y(\cdot,0)=\operatorname{Id}_{\Sigma}$, and $\tfrac{\partial y}{\partial t}=\nu^{\Sigma_{t}}$, with $\nu^{\Sigma_{t}}$ being the outward unit normal to $\Sigma_{t}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}y(\Sigma,t)$. For a function $f\in C^{\infty}(\Sigma)$ with $|f|<\xi$, we define the graph of $f$ over $\Sigma$ as

[TABLE]

Then, slightly abusing notation, let $\mathcal{H}\colon C^{\infty}(\Sigma)\to C^{\infty}(\Sigma)$ be the operator which assigns to a function $f$ the spacetime mean curvature $\mathcal{H}(f)$ of $\operatorname{graph}f$ (with respect to the fixed initial data set $\mathcal{I}$). The linearization of this map $\mathcal{H}$ is computed in the following lemma.

Lemma 5.1**.**

Let $\Sigma\hookrightarrow M^{3}$ be a closed, oriented $2$-surface. Let $\mathcal{V}\colon\Sigma\times(-\xi,\xi)\to M^{3}$ be the normal variation with $\mathcal{V}(\cdot,0)=\operatorname{Id}_{\Sigma}$ and $\left.\frac{\partial\mathcal{V}}{\partial t}\right|_{t=0}=f\nu$ for $f\in C^{\infty}(\Sigma)$. Then the linearization $L^{\mathcal{H}}$ of the spacetime mean curvature map at $\Sigma$ is given by

[TABLE]

where $\triangle^{\Sigma}$, $\nabla^{\Sigma}$ denote the Laplacian and covariant gradient on $(\Sigma,g^{\Sigma})$, respectively.

Proof.

This follows from the definition of spacetime mean curvature $\mathcal{H}=\sqrt{H^{2}-P^{2}}$ and the well-known formulas for $\left.\frac{\partial H(\mathcal{V}(\cdot,t))}{\partial t}\right|_{t=0}$ and $\left.\frac{\partial P(\mathcal{V}(\cdot,t))}{\partial t}\right|_{t=0}$, see Metzger [32, Lemma 5.1]. ∎

The map $L^{\mathcal{H}}$ naturally extends to a bounded mapping $L^{\mathcal{H}}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$. In Section 5.3, we will prove that this mapping has a bounded inverse, for which it is convenient to rewrite the above expression for $L^{\mathcal{H}}$ in the form

[TABLE]

where

[TABLE]

Since the denominator is clearly bounded and bounded away from zero by our assumptions on $\Sigma$, the (bounded) invertibility of $L^{\mathcal{H}}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$ will follow once we show that $\mathcal{L}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$ is invertible with bounded inverse.

Remark 5.2*.*

Recall that the $H\pm P$-stability operator $L^{H\pm P}$ of the map $H\pm P$ (surfaces of constant expansion or null mean curvature) is given by

[TABLE]

see [32]. As it turns out, the analytic properties of $L^{H\pm P}$ imply that constant expansion foliations do not provide an adequate notion of center of mass, in contrast to the STCMC-foliation studied here. The main difference is that the contribution of the second fundamental form $K$ in the $H\pm P$-stability operator is large, while it is rescaled by a factor $\nicefrac{{P}}{{H}}$ in the STCMC-stability operator. The largeness of the contribution of $K$ in the $H\pm P$-stability operator will cause the geometric centers of the surfaces of the foliation to drift away in the direction of the linear momentum $\vec{P}$ in general, see Metzger [32, Section 7]. This can only be avoided by imposing very fast fall-off conditions on $K$ to ensure that $\vec{P}=0$. Furthermore, a certain smallness assumption on $K$ is also directly required to ensure the invertibility of $L^{H\pm P}$, and hence the existence of the constant expansion foliation, see [37, Theorem 3.1].

As a consequence of the factor $\nicefrac{{P}}{{H}}$ in the STCMC-stability operator, no smallness assumption on $K$ will be needed to ensure the existence of the foliation by surfaces of constant spacetime mean curvature. Furthermore, we will see that the leaves of this foliation do not translate as their spacetime mean curvature approaches zero, provided that the standard asymptotic symmetry conditions are imposed.

As we will see, the operator $\mathcal{L}$, and consequently the operator $L^{\mathcal{H}}$, is in many respects similar to the standard (CMC-)stability operator of $\Sigma$, namely to

[TABLE]

This operator has been intensively studied, see e.g. [2].

5.2. Eigenvalues and eigenfunctions of $-\Delta^{\Sigma}$

In preparation for proving the invertibility of the operator $\mathcal{L}$, we summarize the spectral properties of the operator $-\Delta^{\Sigma}$, the Laplacian calculated with respect to the metric $g^{\Sigma}$ induced by $(\Sigma,g^{\Sigma})\hookrightarrow(M^{3},g)$. For this, let us first consider the operator $-\Delta^{\mathbb{S}^{2}_{r}}$, the Laplacian calculated with respect to the standard round metric $\delta^{\mathbb{S}^{2}_{r}}$ on $\mathbb{S}^{2}_{r}=\mathbb{S}^{2}_{r}(\vec{0}\,)\hookrightarrow(\mathbb{R}^{3},\delta)$. The eigenvalues of $-\Delta^{\mathbb{S}^{2}_{r}}$ are $\nicefrac{{l(l+1)}}{{r^{2}}}$ for $l\geq 0$, and the eigenspace corresponding to $\nicefrac{{l(l+1)}}{{r^{2}}}$ is given by the space of homogeneous harmonic polynomials of degree $l$ restricted to $\mathbb{S}^{2}_{r}$, see e.g. [12, Chapter II.4]. In particular, the first non-zero eigenvalue of $-\Delta^{\mathbb{S}^{2}_{r}}$ is $\nicefrac{{2}}{{r^{2}}}$, the corresponding eigenspace is spanned by the restrictions to $\mathbb{S}^{2}_{r}$ of the coordinate functions $x^{1}$, $x^{2}$, $x^{3}$ on $\mathbb{R}^{3}$. In the following, we enumerate the eigenvalues of $-\Delta^{\mathbb{S}^{2}_{r}}$ counting their multiplicity by

[TABLE]

and we denote the associated complete $L^{2}(\mathbb{S}^{2}_{r})$-orthonormal system of eigenfunctions by $\{f_{i}^{\delta}\}_{i=0}^{\infty}$. Without loss of generality, we may assume that the chosen enumeration is such that

[TABLE]

Note that the tracefree part of the Hessian of each of these functions vanishes, and that we have

[TABLE]

where $\nabla^{\mathbb{S}^{2}_{r}}$ denotes the gradient with respect to $\delta^{\mathbb{S}^{2}_{r}}$.

In order to describe the eigenvalues and eigenfunctions of the operator $-\Delta^{\Sigma}$, note that by Proposition 4.5 there is a vector $\vec{z}\in\mathbb{R}^{3}$ and a conformal parametrization $\psi\colon\mathbb{S}^{2}_{r}(\vec{z}\,)\to\Sigma$ such that

[TABLE]

where $r$ is the area radius of $(\Sigma,g^{\Sigma})$. As all spheres of radius $r$ in Euclidean space are isometric, we can easily “translate” $\psi$ to a conformal parametrization $\overline{\psi}\colon\mathbb{S}^{2}_{r}\to\Sigma$ such that

[TABLE]

where $r$ still denotes the area radius of $(\Sigma,g^{\Sigma})$.

We will now describe a complete orthonormal system in $L^{2}(\Sigma)$ consisting of the eigenfunctions $\{f_{i}\}_{i=0}^{\infty}$ such that $-\Delta^{\Sigma}f_{i}=\lambda_{i}f_{i}$, with $0=\lambda_{0}<\lambda_{i}\leq\lambda_{i+1}$, $i=1,2,\ldots$, again counted with multiplicity. The eigenfunctions $f_{i}$ will be chosen so that $\overline{\psi}^{\,*}f_{i}$ is asymptotic to $f_{i}^{\delta}$ for each $i=1,2,\dots$. For simplicity of notation, in what follows we will identify $f_{i}\colon\Sigma\to\mathbb{R}$ with its pullback $\overline{\psi}^{\,*}f_{i}$ without further ado. This enumeration and identification will also allow us to prove useful estimates for the eigenvalues and eigenfunctions of $-\Delta^{\Sigma}$.

Lemma 5.3**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set. Let $a\in[0,1)$, $b\geq 0$, $\eta\in(0,1]$, and consider a $2$-surface $\Sigma\hookrightarrow\mathcal{I}$ such that $\Sigma\in\mathcal{A}(a,b,\eta)$ with respect to $\mathcal{I}$. Then there exist constants $C>0$ and $\overline{\sigma}>0$ depending only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$ such that if $\Sigma$ has constant spacetime mean curvature $\mathcal{H}\equiv\nicefrac{{2}}{{\sigma}}$ for $\sigma>\overline{\sigma}$, then there is a complete orthonormal system in $L^{2}(\Sigma)$ consisting of the eigenfunctions $\{f_{i}\}_{i=0}^{\infty}$ such that

[TABLE]

with $0=\lambda_{0}<\lambda_{i}\leq\lambda_{i+1}$, $i=1,2,\ldots$, counted with multiplicity,

and such that for $i=1,2,3$ the following estimates hold

[TABLE]

Furthermore, $\lambda_{0}=0$ and

[TABLE]

Proof.

By (31) and Lemma A.2 we have

[TABLE]

Applying the Rayleigh Theorem (see e.g. [12, Chapter II.5]), we see from the above estimate and (21) that

[TABLE]

where the infimum is taken over all $f\in W^{1,2}(\Sigma)$ with $\int_{\Sigma}f\,d\mu=0$ and $\|f\|_{L^{2}(\Sigma)}=1$. Of course, the $O$-term constant and the lower bound on $\sigma$ coming from this calculation only depend on $\varepsilon$, $a$, $b$, $\mu$, and $C_{\mathcal{I}}$. We will now construct the respective eigenfunctions $f_{i}$, for which we will use the fact that these functions are solutions to the equation

[TABLE]

where $\lambda^{\delta}_{i}=\nicefrac{{2}}{{r^{2}}}$ for $i=1,2,3$. Noting that the right hand side of the equation equals $-\Delta^{\mathbb{S}^{2}_{r}}f_{i}-\lambda_{i}^{\delta}f_{i}$, and using integration by parts it is straightforward to check that it is orthogonal in $L^{2}(\mathbb{S}^{2}_{r})$ to any element in the kernel of the self-adjoint differential operator in the left hand side. Thus, by the Fredholm Alternative [5, Appendix I, Theorem 31], for every $i=1,2,3$ there is a unique solution $f_{i}-f_{i}^{\delta}\in W^{2,2}(\mathbb{S}^{2}_{r})$ orthogonal in $L^{2}(\mathbb{S}^{2}_{r})$ to the linear space spanned by $f_{i}^{\delta}$, $i=1,2,3$. Note that we may without loss of generality assume that $\|f_{i}\|_{L^{2}(\Sigma)}=1$ so that $\|\Delta^{\Sigma}f_{i}\|_{L^{2}(\Sigma)}=\lambda_{i}=O(r^{-2})$. Since $\mathrm{Scal}^{\Sigma}=\tfrac{2}{\sigma^{2}}+O(\sigma^{-\tfrac{5}{2}-\varepsilon})$ as a consequence of the Gauss equation (see e.g. (43) below), in view of (21) and Lemma A.3 we have $\|f_{i}\|_{W^{2,2}(\Sigma)}=O(1)$ as $\overline{f}_{i}=0$. With the above estimates at hand, it is now straightforward to check that

[TABLE]

so by standard elliptic regularity (see e.g. [5, Appendix H, Theorem 27]222We cite this result for the unit sphere and apply rescaling to extend it to spheres of radius $r>0$.) applied to the operator on the left hand side of (37), we have

[TABLE]

whenever $\sigma>\overline{\sigma}$, for suitably large $C>0$ and $\overline{\sigma}>0$ depending only on $\varepsilon$, $a$, $b$, $\nu$, and $C_{\mathcal{I}}$. This defines the eigenfunctions $f_{i}$, $i=1,2,3$, up to applying the Gram-Schmidt process in the case of multiple eigenvalues. Note that (38) with (21), (31), and the properties of the functions $f_{i}^{\delta}$ implies (34) and (35).

We have $\lambda_{0}=0$ by definition. Using again the Rayleigh Theorem and (21), it is also straightforward to check that $\lambda_{i}>\nicefrac{{5}}{{\sigma^{2}}}$ for $i=4,5,\ldots$ , since the respective eigenvalues of $-\Delta^{\mathbb{S}^{2}_{r}}$ satisfy $\lambda_{i}^{\delta}\geq\nicefrac{{6}}{{r^{2}}}$, whenever $\sigma>\overline{\sigma}$ for suitably large $\overline{\sigma}$ only depending on $\varepsilon$, $a$, $b$, $\mu$, and $C_{\mathcal{I}}$. This concludes the proof. ∎

We can now give a more detailed characterization of the lowest eigenvalues $\lambda_{i}$, $i=1,2,3$. More specifically, in the following lemma we show that these eigenvalues are computed in terms of the Hawking mass

[TABLE]

of $\Sigma$ in the initial data set $\mathcal{I}$. We will drop the explicit reference to $\Sigma$ later and will write $m_{\mathcal{H}}$ instead of $m_{\mathcal{H}}(\Sigma)$. This lemma and its proof are very similar to [36, Lemma 4.5], but rephrased in the spacetime setting.

Lemma 5.4**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with the energy $E$. Suppose that $a\in[0,1)$, $b\geq 0$, $\eta\in(0,1]$, and that $\Sigma\in\mathcal{A}(a,b,\eta)$ with respect to $\mathcal{I}$ is a surface with Hawking mass $m_{\mathcal{H}}(\Sigma)$. Then there exist constants $C>0$ and $\overline{\sigma}>0$, depending only on $\varepsilon$, $a$, $b$, $\eta$, $|E|$, and $C_{\mathcal{I}}$ such that if $\Sigma$ has constant spacetime mean curvature $\mathcal{H}\equiv\nicefrac{{2}}{{\sigma}}$ for $\sigma>\overline{\sigma}$ then the following estimates hold

[TABLE]

and

[TABLE]

Proof.

A polarized version of the standard Bochner formula (see e.g. Proposition 33(3) in [38, Chapter 3]) in dimension 2 applied to the eigenfunctions $f_{i}$ and $f_{j}$ for $i,j=1,2,3$ reads

[TABLE]

Integrating this identity, using the Divergence Theorem on the closed surface $\Sigma$, integrating by parts, and recalling (34), we obtain

[TABLE]

for some constant $C>0$ and all $\sigma>\overline{\sigma}>0$, with $C$ and $\overline{\sigma}$ only depending on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$. Next, the Gauss equation combined with (16) and (26) gives us

[TABLE]

possibly enlarging $C>0$ and $\overline{\sigma}>0$ without introducing new dependencies.

Substituting this into (42) and using (35), (32) together with the fact that our initial data set is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean we conclude that, by partial integration, we get

[TABLE]

When $i\neq j$, $i,j=1,2,3$, this gives us (41) once we recall that $\mathrm{Scal}=O(\sigma^{-3-\varepsilon})$ as a consequence of Definition 2.1 with possibly enlarged $C>0$ and $\overline{\sigma}>0$. In the case $i=j$, $i,j=1,2,3$, one arrives at (40) by combining (44), (21), (32), and the fact that our initial data set is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean, as well as using

[TABLE]

This last inequality follows from (43), the Gauss-Bonnet Theorem, and the definition of $r$, $\sigma$, and $m_{\mathcal{H}}(\Sigma)$ with possibly enlarged $C>0$ and $\overline{\sigma}>0$. This proves the claims of the lemma. ∎

Remark 5.5*.*

Since $P^{2}=O(\sigma^{-3-2\varepsilon})$ in (5.2), this lemma remains valid if we replace the Hawking mass $m_{\mathcal{H}}(\Sigma)$ by the Geroch mass $m_{H}(\Sigma)$ (also sometimes referred to as “(Riemannian) Hawking mass”) given by

[TABLE]

The same remark will hold true for the subsequent results. However, we choose to use $m_{\mathcal{H}}(\Sigma)$, and not $m_{H}(\Sigma)$, throughout to emphasize the spacetime nature of our result.

5.3. Invertibility of the operator $\mathcal{L}$

Section 5.2 above provides the following description of the eigenvalues of the Laplacian $-\Delta^{\Sigma}$:

• •

$\lambda_{0}=0$,

• •

for $i=1,2,3$ the eigenvalues $\lambda_{i}$ are characterized by formula (40),

• •

for $i=4,5,\ldots$ we have $\lambda_{i}>\nicefrac{{5}}{{\sigma^{2}}}$.

It turns out to be useful to decompose functions $h\in L^{2}(\Sigma)$ with respect to the $L^{2}(\Sigma)$-complete orthonormal system $\{f_{0},f_{1},f_{2},f_{3},\dots\}$ of eigenfunctions corresponding to $-\Delta^{\Sigma}$ when analyzing $\mathcal{L}h$ and the $L^{2}(\Sigma)$-adjoint $\mathcal{L}^{*}h$ (the latter being of interest as we are aiming for a Fredholm Alternative argument). More specifically, it is useful to split any given function $h\in L^{2}(\Sigma)$ into its mean value

[TABLE]

its translational part

[TABLE]

and the difference part333We do not call $h^{d}$ the “deformational part” as some other authors do, because $h^{d}$ also contains the mean value information. In other words, we primarily use this splitting to distinguish between the eigenfunctions corresponding to eigenvalues of different magnitude.

[TABLE]

Proposition 5.6**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with energy $E$. Suppose that $a\in[0,1)$, $b\geq 0$, $\eta\in(0,1]$, and that $\Sigma\in\mathcal{A}(a,b,\eta)$ with respect to $\mathcal{I}$ is a surface with non-vanishing Hawking mass $m_{\mathcal{H}}(\Sigma)\neq 0$. Then there exist constants $C>0$ and $\overline{\sigma}>0$, with $C$ depending only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$ and $\overline{\sigma}$ in addition depending on $|E|$ in (54), such that if $\Sigma$ has constant spacetime mean curvature $\mathcal{H}\equiv\nicefrac{{2}}{{\sigma}}$ for $\sigma>\overline{\sigma}$, the following estimates

[TABLE]

hold for any $h,h_{1},h_{2}\in W^{2,2}(\Sigma)$. The same estimates apply to the $L^{2}(\Sigma)$-adjoint $\mathcal{L}^{*}$. Moreover, the Hawking mass $m_{\mathcal{H}}(\Sigma)$ and the energy $E$ are related by

[TABLE]

In particular, the operator $\mathcal{L}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$ is invertible as long as the energy $E$ of the initial data set does not vanish and $\overline{\sigma}$ is sufficiently large, depending only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$.

Proof.

In this proof, $C>0$ and $\overline{\sigma}>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$, and, in the case of (54), also on $|E|$.

Proving (50).

By definition of $\mathcal{L}$ in (28), we have

[TABLE]

It follows from Proposition 4.5 that $|\mathring{A}|^{2}=O(\sigma^{-3-{\varepsilon}})$, and we know from (26) that $H^{2}=4\sigma^{-2}+O(\sigma^{-3-2{\varepsilon}})$, and that $\frac{P}{H}=O(\sigma^{-\frac{1}{2}-\varepsilon})$. Furthermore, the definition of $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data sets implies that $\mathrm{Ric}(\nu,\nu)=O(\sigma^{-\frac{5}{2}-{\varepsilon}})$ and that $\nabla_{\nu}\operatorname{tr}K-\nabla_{\nu}K(\nu,\nu)=O(\sigma^{-\frac{5}{2}-{\varepsilon}})$. Hence,

[TABLE]

By (32), we have

[TABLE]

whereas (35) implies by a Cauchy–Schwarz Inequality that

[TABLE]

recalling that $h^{t}\in\operatorname{Span}(f_{1},f_{2},f_{3})$. This proves (50).

Proving (51).

Arguing as above, by Proposition 4.5, Lemma 5.3, and our decay assumptions on the initial data set, we have that

[TABLE]

for any $i,j=1,2,3$. It then follows by Lemma 5.4 that

[TABLE]

and

[TABLE]

In particular, we see that (51) holds for $h_{1}^{t},h_{2}^{t}\in\{f_{1},f_{2},f_{3}\}$. The general case follows by bilinearity and by the Cauchy–Schwarz Inequality on $\mathbb{R}^{3}$.

Proving (52).

By definition of $\mathcal{L}$, we have that

[TABLE]

Next, by (26), we have

[TABLE]

while the $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotic decay assumptions as well as (26), (35), and the Cauchy–Schwarz Inequality lead to

[TABLE]

Proposition 4.5 gives us that $|\mathring{A}|^{2}=O(\sigma^{-3-2\varepsilon})$. Summarizing, we found

[TABLE]

recalling that $u^{d}$ is $L^{2}(\Sigma)$-orthogonal to $f_{i}$ for $i=1,2,3$.

Proving (53).

We will use the following manifest relation for the linear operator $\mathcal{L}$

[TABLE]

Arguing similarly to how we argued above, we now integrate by parts and use Proposition 4.5, (26), and (36) giving $\lambda_{i}>\nicefrac{{5}}{{\sigma^{2}}}$ for $i=4,5,\dots$, and the asymptotic decay conditions on $\mathcal{I}$ to estimate from below the expression

[TABLE]

as $h^{d}-h^{0}\in\operatorname{Span}(f_{4},f_{5},\dots)$. Here, the factor $\tfrac{7}{4}<3$ is chosen for later convenience. Hence by a Cauchy–Schwarz Inequality on $\int_{\Sigma}(h^{d}-h^{0})\mathcal{L}(h^{d}-h^{0})\,d\mu$, we obtain

[TABLE]

Note also that $h_{0}$ is a constant, so that

[TABLE]

and thus

[TABLE]

where again, the factor $\tfrac{7}{4}<2$ is chosen for later convenience. Using (58), integration by parts and finally Young’s Inequality, one can also check with the same decay arguments as above that

[TABLE]

[TABLE]

Combing this estimate with (57) and (59), (53) follows from (56) once we recall that $h^{d}-h^{0}$ is $L^{2}(\Sigma)$-orthogonal to $h^{0}$.

Proving (55).

To see that $E$ and $m_{\mathcal{H}}(\Sigma)$ are as close as claimed, we recall the well-known fact that the Geroch mass $m_{H}(\Sigma)$ of sufficiently round large surfaces in a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically flat initial data set $\mathcal{I}$ is close to the energy $E$ of $\mathcal{I}$. More specifically, Lemma A.1 in [36] (relying on [23] and [33]) and (5.2) imply that

[TABLE]

Thus $m_{\mathcal{H}}(\Sigma)\neq 0$ if $E\neq 0$ as long as $\sigma>\overline{\sigma}$ with $C>0$ and $\overline{\sigma}>0$ sufficiently large, depending only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$.

Proving (54).

We pick a number $\kappa>0$ such that $1-4\varepsilon<2\kappa<1$ and consider two cases. In this part, we will abbreviate $m_{\mathcal{H}}(\Sigma)\mathrel{=\mkern-5.2mu\raise 0.40903pt\hbox{\rm:}}m_{\mathcal{H}}$.

Case 1.

Suppose that $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\geq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$. As a consequence, using (50), (53), and Young’s Inequality, we obtain for any $0<\alpha<1$, e.g. $\alpha=\tfrac{1}{2}$, that

[TABLE]

provided that $\sigma>\overline{\sigma}$, where now $\overline{\sigma}$ may actually depend on $E$ as we used (55) in the last step. Thus (54) holds in case $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\geq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$.

Case 2.

Now assume that $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\leq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$. In this case,

[TABLE]

where by (51), (55), and because $h^{t}$ and $h^{d}$ are $L^{2}(\Sigma)$-orthogonal, we see that for $\sigma>\overline{\sigma}$, $\overline{\sigma}$ suitably large depending now in addition on $E$, we have

[TABLE]

where we used that $\varepsilon\leq\tfrac{1}{2}$ by definition. Moreover, using the Cauchy–Schwarz Inequality and the assumption $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\leq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$, we estimate

[TABLE]

Further, arguing once more as above with the explicit form of $\mathcal{L}$ in (28), using the asymptotic decay conditions of $\mathcal{I}$, (16), (26), and integration by parts, one confirms that

[TABLE]

using again the assumption $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\leq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$ as well as $\|h^{t}\|_{L^{2}(\Sigma)}\leq\|h\|_{L^{2}(\Sigma)}$ in the last step. It then follows from (60), the Cauchy–Schwarz Inequality, and the bounds on $\kappa$ that (54) also holds when $\|h^{d}\|^{2}_{L^{2}(\Sigma)}\leq\sigma^{-\frac{1}{2}-\kappa}\|h\|^{2}_{L^{2}(\Sigma)}$.

Combining Case 1 and Case 2, we have thus shown (54). To conclude the proof, it only remains to show that $\mathcal{L}^{*}$ also satisfies the estimates (50)–(54) and that $\mathcal{L}$ is invertible provided that the initial data set $\mathcal{I}$ has non-vanishing energy $E\neq 0$ and the spacetime mean curvature radius $\sigma$ of $\Sigma$ is sufficiently large, $\sigma>\overline{\sigma}$, $\overline{\sigma}$ depending on $E$.

Invertibility of $\mathcal{L}$ and estimates on $\mathcal{L}^{*}$.

The operator $\mathcal{L}$ is not self-adjoint, but its $L^{2}(\Sigma)$-adjoint $\mathcal{L}^{*}$ has very similar structure, differing only in the last term. In $\mathcal{L}h$, this term reads $\frac{P}{H}K(\nabla^{\Sigma}h,\nu)$, while in $\mathcal{L}^{*}h$, this term becomes $-\operatorname{div}^{\Sigma}\left(\frac{P}{H}K(\cdot,\nu)\right)h$. Recall that

[TABLE]

Going back to the proofs of (50)–(54), we see that all of them work if we replace $\mathcal{L}$ by $\mathcal{L}^{*}$ modulo exchanging the performance of partial integration with the decay estimate (61) and vice versa. This, in particular (54), implies that the $L^{2}(\Sigma)$-kernel of $\mathcal{L}^{*}$ is trivial, and hence $\mathcal{L}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$ is invertible by the Fredholm Alternative [5, Appendix I, Theorem 31], as long as $m_{\mathcal{H}}(\Sigma)\neq 0$ which is guaranteed from $E\neq 0$ and (55). The Fredholm Alternative applies as $\mathcal{L}$ is clearly a linear elliptic operator as its symbol is that of the Laplacian $-\Delta^{\Sigma}$ and because $\Sigma$ is compact. ∎

Corollary 5.7**.**

For every $h\in W^{2,2}(\Sigma)$, we have

[TABLE]

for $\sigma>\overline{\sigma}$ and with $C>0$, $\overline{\sigma}>0$ only depending on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$, and $E$.

Proof.

Note by (43) we have $\mathrm{Scal}^{\Sigma}=\tfrac{2}{\sigma^{2}}+O(\sigma^{-\tfrac{5}{2}-\varepsilon})$, hence, in the view of (21), Lemma A.3 applies to $\Sigma$. Combined with the Cauchy-Schwartz Inequality and (21), this result gives us

[TABLE]

Recalling the definition of the operator $\mathcal{L}$ (see (28)) and the fall-off properties of the initial data set, we further find that

[TABLE]

hence

[TABLE]

Combining (62) with this estimate and (54) we thereby obtain

[TABLE]

This proves the estimate for $h$. The estimate for $h^{d}$ is proven similarly, using (53) instead of (54). ∎

6. Existence and uniqueness of the STCMC-foliation

In this section, we will prove that any $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ is foliated (i.e., roughly speaking, covered without gaps or overlaps), outside a compact set, by $2$-surfaces of constant spacetime mean curvature (STCMC). We also prove a uniqueness result for STCMC-surfaces.

6.1. Existence of the STCMC-foliation

In [36], Nerz proved the following result, rephrased here in our notation. Note that because of time symmetry $K\equiv 0$, the CMC-foliation constructed by Nerz can be viewed as a special case of the STCMC-foliation under consideration here.

Theorem 6.1** (Nerz 2015).**

Let $(M^{3},g)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean Riemannian manifold viewed as a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}_{0}=(M^{3},g,K\equiv 0,\mu=\tfrac{1}{2}\mathrm{Scal},J\equiv 0)$, with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}_{0}}>0$ depending only on $\varepsilon$, $C_{\mathcal{I}_{0}}$, and $E$, a compact set $\mathcal{K}_{0}\subset M^{3}$, and a bijective $C^{1}$-map $\Psi_{0}\colon(\sigma_{\mathcal{I}_{0}},\infty)\times\mathbb{S}^{2}\to M^{3}\setminus\mathcal{K}_{0}$ such that each of the surfaces $\Sigma_{0}^{\sigma}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Psi_{0}(\sigma,\mathbb{S}^{2})$ has constant mean curvature $H(\Sigma_{0}^{\sigma})\equiv\nicefrac{{2}}{{\sigma}}$ provided that $\sigma>\sigma_{\mathcal{I}_{0}}$.

This result is a starting point for proving the following theorem, which is essentially the main result of this paper. For the sake of clarity of exposition, we provide the proof of the following theorem right away, saving the verification of some preliminary lemmas for later. We state Theorem 6.2 here in a notation convenient for its proof.

Theorem 6.2** (Existence of STCMC-foliation).**

Let $\mathcal{I}_{1}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}_{1}}>0$ depending only on $\varepsilon$, $C_{\mathcal{I}_{1}}$, and $E$, a compact set $\mathcal{K}_{1}\subset M^{3}$, and a bijective $C^{1}$-map $\Psi_{1}\colon(\sigma_{\mathcal{I}_{1}},\infty)\times\mathbb{S}^{2}\to M^{3}\setminus\mathcal{K}_{1}$ such that each of the surfaces $\Sigma_{1}^{\sigma}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Psi_{1}(\sigma,\mathbb{S}^{2})$ has constant spacetime mean curvature $\mathcal{H}(\Sigma_{1}^{\sigma})\equiv\nicefrac{{2}}{{\sigma_{1}}}$ provided that $\sigma>\sigma_{\mathcal{I}_{1}}$.

Remark 6.3*.*

As the proof of Theorem 6.2 will show, the surfaces $\Sigma_{1}^{\sigma}$ are in fact asymptotically centered in the sense of Definition 4.1, more specifically, $\Sigma_{1}^{\sigma}\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ for all $\sigma>\sigma_{\mathcal{I}_{1}}$, with constants $b_{\mathcal{I}_{1}}>0$, $\eta_{\mathcal{I}_{1}}\in(0,1]$, and $\sigma_{\mathcal{I}_{1}}>0$ defined in the proof of Theorem 6.2, and depending only on $\varepsilon$, $C_{\mathcal{I}_{1}}$, and $E$.

Proof.

The family of closed, oriented $2$-surfaces $\{\Sigma_{1}^{\sigma}\}_{\sigma>\sigma_{\mathcal{I}_{1}}}$ will be constructed via a method of continuity, see also Section 3. Roughly speaking, we will deform the constant (automatically spacetime) mean curvature foliation $\{\Sigma_{0}^{\sigma}\}_{\sigma>\sigma_{\mathcal{I}_{0}}}$ of the initial data set $\mathcal{I}_{0}$ from Theorem 6.1 along the curve of initial data sets $\{\mathcal{I}_{\tau}\}_{\tau\in[0,1]}$, where $\mathcal{I}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}(M^{3},g,\tau K,\mu_{\tau},\tau J)$ is as described in Section 3.1, arriving at the foliation of the initial data set $\mathcal{I}_{1}$ by constant spacetime mean curvature surfaces $\{\Sigma^{\sigma}_{1}\}_{\sigma>\sigma_{\mathcal{I}_{1}}}$. In order to make this idea more precise, we introduce the following construction.

By Theorem 6.1, we know that for every $\sigma>\sigma_{\mathcal{I}_{0}}$ there is a closed, oriented $2$-surface $\Sigma^{\sigma}_{0}\hookrightarrow M^{3}$ with constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{0})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{0}$. Furthermore, the proof of this result in [36] shows that there are constants $b_{\mathcal{I}_{0}}\geq 0$ and $1\geq\eta_{\mathcal{I}_{0}}>0$ such that $\Sigma^{\sigma}_{0}\in\mathcal{A}(0,b_{\mathcal{I}_{0}},\eta_{\mathcal{I}_{0}})$ for all $\sigma>\sigma_{\mathcal{I}_{0}}$. We recall from [36] that $b_{\mathcal{I}_{0}}$ and $\eta_{\mathcal{I}_{0}}$ only depend on $\varepsilon$, $C_{\mathcal{I}_{0}}$, and $E$ which can be restated as saying that they only depend on $\varepsilon$, $C_{\mathcal{I}_{1}}$, and $E$ by our construction. Set $b_{\mathcal{I}_{1}}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}4b_{\mathcal{I}_{0}}>b_{\mathcal{I}_{0}}$ and $\eta_{\mathcal{I}_{0}}>\eta_{\mathcal{I}_{1}}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\tfrac{\eta_{\mathcal{I}_{0}}}{4}>0$. From Section 5 and by the definition of $b_{\mathcal{I}_{1}}$ and $\eta_{\mathcal{I}_{1}}$, we know that there are constants $C$ and $\overline{\sigma}$ depending only on $\varepsilon$, $C_{\mathcal{I}_{1}}$, and $E$ such that the operator $\mathcal{L}\colon W^{2,2}(\Sigma)\to L^{2}(\Sigma)$ is invertible whenever $\Sigma\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ is a surface of constant spacetime mean curvature $\mathcal{H}(\Sigma)\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{1}$ for $\sigma\geq\overline{\sigma}$, and whenever in addition the estimates of Proposition 5.6 and Corollary 5.7 are available on $\Sigma$. Without loss of generality, we may also assume that $C$ and $\overline{\sigma}$ are such that the regularity result in Proposition 4.5 as well as a supplementary result stated in Lemma 6.8 (see Section 6.2 below) apply with $a=\overline{a}=0$, $b=\tfrac{b_{\mathcal{I}_{1}}}{2}$, $\overline{b}=b_{\mathcal{I}_{1}}$, $\eta=2\eta_{\mathcal{I}_{1}}$, $\overline{\eta}=\eta_{\mathcal{I}_{1}}$. We set $\sigma_{\mathcal{I}_{1}}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\max\{\overline{\sigma},4\sigma_{\mathcal{I}_{0}}\}$, and note that by their definition $\sigma_{\mathcal{I}_{1}}$, $b_{\mathcal{I}_{1}}$, and $\eta_{\mathcal{I}_{1}}$ only depend on $\varepsilon$, $C_{\mathcal{I}_{1}}$, and $E$.

Now fix $\sigma_{*}>\sigma_{\mathcal{I}_{1}}$ for the rest of the argument until we start discussing the foliation property when applying Lemma 6.10. Let $Y^{\sigma_{*}}\subseteq[0,1]$ be the maximal subset such that there is a $C^{1}$-map

[TABLE]

with the following properties for every $\tau\in Y^{\sigma_{*}}$:

1. (i)

The surface $\Sigma^{\sigma_{*}}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathcal{F}^{\sigma_{*}}(\tau,\mathbb{S}^{2})$ has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma_{*}}_{\tau})\equiv\nicefrac{{2}}{{\sigma_{*}}}$ with respect to the initial data set $\mathcal{I}_{\tau}$. 2. (ii)

$\partial_{\tau}\mathcal{F}^{\sigma_{*}}(\tau,q)$ is orthogonal to $\Sigma^{\sigma_{*}}_{\tau}$ for every $q\in\mathbb{S}^{2}$. 3. (iii)

The surface $\Sigma^{\sigma_{*}}_{\tau}$ is asymptotically centered in the sense that $\Sigma^{\sigma_{*}}_{\tau}\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$.

Maximality of $Y^{\sigma_{*}}$ is understood here as follows: if the above conditions are satisfied for some $\widetilde{Y}^{\sigma_{*}}\subseteq[0,1]$ and a map $\widetilde{\mathcal{F}}^{\sigma_{*}}\colon\widetilde{Y}^{\sigma_{*}}\times\mathbb{S}^{2}\to M^{3}$, then $\widetilde{Y}^{\sigma_{*}}\subseteq Y^{\sigma_{*}}$ as well as $\mathcal{F}^{\sigma_{*}}|_{\widetilde{Y}^{\sigma^{*}}}\equiv\widetilde{\mathcal{F}}^{\sigma_{*}}$.

Note that for $\tau=0$, Condition (i) is ensured by the assumptions in Theorem 6.1. The same is true for Condition (iii) once one takes into account that $\mathcal{A}(0,b_{\mathcal{I}_{0}},\eta_{\mathcal{I}_{0}})\subseteq\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$. However, Condition (ii) is not automatically satisfied for $\tau=0$ as we do not even know whether the map $\mathcal{F}^{\sigma_{*}}$ exists. The following lemma ensures that $Y^{\sigma_{*}}$ contains an interval $[0,\tau_{0})$ for some $\tau_{0}>0$. In particular, Condition (ii) is satisfied a posteriori for $\tau=0$. More generally, this result shows that $Y^{\sigma_{*}}$ is open around any $\tau_{*}\in Y^{\sigma_{*}}$ such that $\Sigma^{\sigma_{*}}_{\tau_{*}}\in\mathcal{A}(0,b,\eta)$ for $0\leq b<b_{\mathcal{I}_{1}}$ and $\eta_{\mathcal{I}_{1}}<\eta\leq 1$.

Lemma 6.4**.**

Let $0\leq b<\overline{b}\leq b_{\mathcal{I}_{1}}$ and $\eta_{\mathcal{I}_{1}}\leq\overline{\eta}<\eta\leq 1$. For any $\tau_{*}\in[0,1]$ for which there exists a smooth surface $\Sigma^{\sigma_{*}}_{\tau_{*}}\in\mathcal{A}(0,b,\eta)$ satisfying $\mathcal{H}(\Sigma^{\sigma_{*}}_{\tau_{*}})\equiv\nicefrac{{2}}{{\sigma_{*}}}$, there exists an open, connected neighborhood $U_{\tau_{*}}$ of $\tau_{*}$ inside $[0,1]$ and a unique $C^{1}$-map $\mathcal{F}^{\sigma_{*}}\colon U_{\tau_{*}}\times\mathbb{S}^{2}\to M^{3}$ with $\Sigma^{\sigma_{*}}_{\tau_{*}}=\mathcal{F}^{\sigma_{*}}(\tau_{*},\cdot)$ such that (i) and (ii) are satisfied for $\Sigma^{\sigma_{*}}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathcal{F}^{\sigma_{*}}(\tau,\cdot)$, and such that $\Sigma^{\sigma_{*}}_{\tau}\in\mathcal{A}(0,\overline{b},\overline{\eta})$ for all $\tau\in U_{\tau_{*}}$.

Proof.

In order to prove this lemma, suppose that $\tau_{*}\in[0,1]$, and that $b$ and $\eta$ are such as in the statement. As discussed in Section 5.1, in a neighborhood of each $\Sigma^{\sigma_{*}}_{\tau_{*}}$, we may introduce normal geodesic coordinates $y\colon\Sigma^{\sigma_{*}}_{\tau_{*}}\times(-\xi,\xi)\to M^{3}$ near $\Sigma^{\sigma_{*}}_{\tau_{*}}$. Now let $U^{2,2}_{\xi}(\Sigma^{\sigma_{*}}_{\tau_{*}})\subseteq W^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ be an open neighborhood of $0\in W^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ such that $f\in U^{2,2}_{\xi}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ implies $|f|<\xi$; such a neighborhood exists by Lemma A.2.

Next, we consider the graphical spacetime mean curvature map

[TABLE]

which assigns, to every $\tau\in[0,1]$ and every $f\in U_{\xi}^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$, the spacetime mean curvature $\mathcal{H}(\operatorname{graph}f)$ of the geodesic graph, $\operatorname{graph}f=\{y(q,f(q)):q\in\Sigma^{\sigma_{*}}_{\tau_{*}}\}$, with respect to the initial data set $\mathcal{I}_{\tau}$. The Fréchet derivative of the map $h^{\sigma_{*}}$ with respect to the second argument $f$ at the point $(\tau_{*},0)$ is the operator $L^{\mathcal{H}}\colon W^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})\to L^{2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ given by Lemma 5.1, where all geometric quantities are computed with respect to the initial data set $\mathcal{I}_{\tau_{*}}$. As shown in Section 5, the linearized operator $L^{\mathcal{H}}\colon W^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})\to L^{2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ is continuously invertible, because $\sigma_{*}>\sigma_{\mathcal{I}_{1}}$.

By the Implicit Function Theorem, there thus exists a relatively open neighborhood $\widetilde{U}_{\tau_{*}}\subseteq[0,1]$ of $\tau_{*}$ and a unique $C^{1}$-map $\gamma^{\sigma_{*}}\colon\widetilde{U}_{\tau_{*}}\to U_{\xi}^{2,2}(\Sigma^{\sigma_{*}}_{\tau_{*}})$ such that $\gamma^{\sigma_{*}}(\tau_{*})=0$ and

[TABLE]

Thus, by construction, for every $\tau\in\widetilde{U}_{\tau_{*}}$, the surface $\Sigma^{\sigma_{*}}_{\tau}=\operatorname{graph}\gamma^{\sigma_{*}}(\tau)$ has constant spacetime mean curvature $\nicefrac{{2}}{{{\sigma_{*}}}}$ with respect to the initial data set $\mathcal{I}_{\tau}$.

Recall that the surface $\Sigma^{\sigma_{*}}_{\tau_{*}}$ is a graph over some round sphere by our assumptions and by Proposition 4.5, recalling again the a priori bounds on $\sigma_{*}$, $b$, and $\eta$. As $\Sigma^{\sigma_{*}}_{\tau}$ was defined as a graph over $\Sigma^{\sigma_{*}}_{\tau_{*}}$ for every $\tau\in\widetilde{U}_{\tau_{*}}$, composition of these two graphical representations gives us that $\Sigma^{\sigma_{*}}_{\tau}$ is parametrized over a round sphere.

Thus, we may now define the map $\mathcal{F}^{\sigma_{*}}\colon\widetilde{U}_{\tau_{*}}\times\mathbb{S}^{2}\to M^{3}$ by setting $\mathcal{F}^{\sigma_{*}}(\tau,\mathbb{S}^{2})\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Sigma^{\sigma_{*}}_{\tau}$, and defining the parametrization of $\Sigma^{\sigma_{*}}_{\tau}$ by requesting that $\partial_{\tau}\mathcal{F}^{\sigma_{*}}$ be orthogonal to $\Sigma^{\tau}_{\sigma_{*}}$ for all $\tau\in\widetilde{U}_{\tau_{*}}$.

By continuity of $\mathcal{F}^{\sigma_{*}}$ and because $0\leq b<\overline{b}$, $0<\overline{\eta}<\eta\leq 1$, and $\Sigma^{\sigma_{*}}_{\tau_{*}}\in\mathcal{A}(0,b,\eta)$, there exists an open neighborhood $U_{\tau_{*}}\subseteq\widetilde{U}_{\tau_{*}}$ of $\tau_{*}$ such that

[TABLE]

holds for all $\tau\in U_{\tau_{*}}$ as desired. This proves Lemma 6.4. ∎

Choosing $b=b_{\mathcal{I}_{0}}$, $\overline{b}=2b_{\mathcal{I}_{0}}=\tfrac{b_{\mathcal{I}_{1}}}{2}$, and $\eta=\eta_{\mathcal{I}_{0}}$, $\overline{\eta}=\tfrac{\eta_{\mathcal{I}_{0}}}{2}=2\eta_{\mathcal{I}_{1}}$, Lemma 6.4 shows directly via Theorem 6.1 that $0\in Y^{\sigma_{*}}$ and that $Y^{\sigma_{*}}$ is relatively open near [math]. Now we let $X^{\sigma_{*}}$ be the maximal connected subinterval of $Y^{\sigma_{*}}$ containing $\tau=0$. As we have just seen by Lemma 6.4, $X^{\sigma_{*}}$ is relatively open near $\tau=0$. Set $\tau^{*}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\sup X^{\sigma_{*}}$. In Lemma 6.5 below we will show that $\tau^{*}\in X^{\sigma_{*}}$, so that $X^{\sigma_{*}}=[0,\tau^{*}]$ is closed, where $0<\tau^{*}\leq 1$.

Lemma 6.5**.**

The interval $Y^{\sigma_{*}}\subseteq[0,1]$ is closed.

Proof.

Closedness of $Y^{\sigma_{*}}$ can be addressed by following the arguments given in [36, Lemma 5.6] and [37, Lemma 3.14], as the necessary preliminaries are available in the form of Lemma 6.6 and Lemma 6.7 below. Alternatively, one may rely on a more standard method used in [32, Proof of Proposition 6.1], which we describe below. The Sobolev spaces we use throughout the paper are weighted, however, for a given closed, oriented $2$-surface, the weighted Sobolev norms are equivalent to the traditional unweighted ones; we will thus switch to the usual unweighted ones for this proof in order to allow us to use standard results on Sobolev spaces on $2$-surfaces.

Let $\{\tau_{n}\}_{n=1}^{\infty}\subset Y^{\sigma_{*}}$ be a sequence such that $\lim_{n\to\infty}\tau_{n}=\tau\in[0,1]$ and let $\Sigma^{\sigma_{*}}_{\tau_{n}}\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ be a surface with constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma_{*}}_{\tau_{n}})\equiv\nicefrac{{2}}{{\sigma_{*}}}$ with respect to the initial data set $\mathcal{I}_{\tau_{n}}$. By Proposition 4.5 we know that there are functions $f_{n}\colon\mathbb{S}_{r_{n}}(\vec{z}_{n})\to\mathbb{R}$ such that $\Sigma^{\sigma_{*}}_{\tau_{n}}=\operatorname{graph}f_{n}$ where $r_{n}$ and $\vec{z}_{n}$ are the area radius and the coordinate center of $\Sigma^{\sigma_{*}}_{\tau_{n}}$. By the first inequality of (15) and by (21), we know that the sequences $\{r_{n}\}_{n=1}^{\infty}$ and $\{\vec{z}_{n}\}_{n=1}^{\infty}$ are uniformly bounded, so we may assume (up to passing to a subsequence) that $\lim_{n\to\infty}r_{n}=r$ and $\lim_{n\to\infty}\vec{z}_{n}=\vec{z}$. Consequently, in view of (17), we may assume that there is a sequence $\{\tilde{f}_{n}\}_{n=1}^{\infty}$, such that $\tilde{f}_{n}\colon\mathbb{S}_{r}(\vec{z}\,)\to\mathbb{R}$ and $\Sigma^{\sigma_{*}}_{\tau_{n}}=\operatorname{graph}\tilde{f}_{n}$. Again in the view of (17), we may assume that this sequence is uniformly bounded in $W^{2,\infty}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ and hence in $C^{1,\beta}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ for any $0<\beta<1$. Recalling that $\Sigma^{\sigma_{*}}_{\tau_{n}}$ are surfaces of constant spacetime mean curvature, we see that the functions $\tilde{f}_{n}$ satisfy a linear elliptic PDE of the form

[TABLE]

with uniformly bounded coefficients $a_{n}^{\beta\gamma},b_{n}^{\beta},F_{n}\in C^{0,\beta}(\mathbb{S}^{2}_{r}(\vec{z}\,))$, see Appendix B for details. A standard argument using Schauder estimates (see e.g. [22, Theorem 9.19] and [27, Theorem 10.2.1]) allows us to conclude that the functions $\tilde{f}_{n}\in C^{2,\beta}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ are uniformly bounded in $C^{2,\beta}(\mathbb{S}^{2}_{r}(\vec{z}\,))$, and consequently, up to passing to a subsequence, we may assume that $\{\tilde{f}_{n}\}_{n=1}^{\infty}$ converges in $C^{2,\alpha}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ to a limit $f\in C^{2,\alpha}(\mathbb{S}^{2}_{r}(\vec{z}\,))$ for some fixed $0<\alpha<1$. As a consequence of (64) and $C^{2,\alpha}$-convergence, we see that $\Sigma^{\sigma_{*}}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\operatorname{graph}f$ has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma_{*}}_{\tau})\equiv\nicefrac{{2}}{{\sigma^{*}}}$.

Finally, we confirm that $\Sigma^{\sigma_{*}}_{\tau}=\operatorname{graph}f\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ by passing to the limit in the respective inequalities of (15) for $\Sigma^{\sigma_{*}}_{\tau_{n}}=\operatorname{graph}\tilde{f}_{n}\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$. Again, this is possible in the view of the $C^{2,\alpha}$-convergence of the graph functions. ∎

Continuing the proof of Theorem 6.2, we will now use Lemma 6.8 below to show that $\tau^{*}=1$, arguing by contradiction. Suppose instead that $\tau^{*}<1$. Then $\Sigma^{\sigma_{*}}_{\tau}\in\mathcal{A}(0,b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ for all $\tau\in[0,\tau^{*}]$, whereas $\Sigma^{\sigma_{*}}_{0}\in\mathcal{A}(0,b_{\mathcal{I}_{0}},\eta_{\mathcal{I}_{0}})$. Applying Lemma 6.8, we see that in fact $\Sigma^{\sigma_{*}}_{\tau}\in\mathcal{A}(0,\tfrac{b_{\mathcal{I}_{1}}}{2},2\eta_{\mathcal{I}_{1}})$. As a consequence, we may apply Lemma 6.4 with $b=\tfrac{b_{\mathcal{I}_{1}}}{2}$, $\overline{b}=b_{\mathcal{I}_{1}}$, $\eta=2\eta_{\mathcal{I}_{1}}$, and $\overline{\eta}=\eta_{\mathcal{I}_{1}}$ to show that $[0,\tau^{*}+\rho)\subseteq Y^{\sigma_{*}}$ for some $\rho>0$. This contradicts the maximality of the intervals $X^{\sigma_{*}}$, hence $\tau^{*}=1$ and $Y^{\sigma_{*}}=[0,1]$.

Summing up, we have shown that for each $\sigma>\sigma_{\mathcal{I}_{1}}$ there is a surface $\Sigma^{\sigma}_{1}=\mathcal{F}^{\sigma}(1,\mathbb{S}^{2})$ such that its spacetime mean curvature in the initial data set $\mathcal{I}_{1}$ is $\mathcal{H}(\Sigma^{\sigma}_{1})\equiv\nicefrac{{2}}{{\sigma}}$. We may now define $\Psi_{1}\colon(\sigma_{\mathcal{I}_{1}},\infty)\times\mathbb{S}^{2}\to M^{3}$ by setting

[TABLE]

The only remaining thing to check is that the family $\{\Sigma^{\sigma}_{1}\}_{\sigma>\sigma_{\mathcal{I}_{1}}}$ is a foliation, which will be the case if $\Psi_{1}$ is a bijective $C^{1}$-map onto the exterior region $M^{3}\setminus\mathcal{K}_{1}$ of a suitably large compact set $\mathcal{B}\subseteq\mathcal{K}_{1}\subset M^{3}$. This is proven in Lemma 6.10. Note that in this step, we may need to increase $\sigma_{\mathcal{I}_{1}}$, albeit without introducing new dependencies. ∎

6.2. Supplementary lemmas

We will now prove the supplementary lemmas that were used in the proof of Theorem 6.2 above.

Lemma 6.6**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$, with $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ denoting the asymptotic coordinate chart. Assume in addition that $K$ satisfies the potentially stronger decay assumptions $|K|\leq C_{\mathcal{I}}\,|\vec{x}\,|^{-\delta-\varepsilon}$ for some $\delta\geq\tfrac{3}{2}$ and all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$.

Let $\emptyset\neq U\subseteq[0,1]$ be an open subset of $[0,1]$ and define $\mathcal{I}_{\tau}$ as in the proof of Theorem 6.2 for each $\tau\in U$. Let $a\in[0,1)$, $b\geq 0$, $\eta\in(0,1]$ be fixed. Then there exist constants $\overline{\sigma}>0$ and $C>0$, depending only on $\varepsilon$, $\delta$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$, and $E$ such that the following holds for any $\sigma>\overline{\sigma}$: Assume there exists a $C^{1}$-map $\mathcal{F}^{\sigma}\colon U\times\mathbb{S}^{2}\to M^{3}$ such that for every $\tau\in U$ the surface $\Sigma^{\sigma}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathcal{F}^{\sigma}(\tau,\mathbb{S}^{2})$ is in the a priori class $\mathcal{A}(a,b,\eta)$ and has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{\tau}$. Assume further that $\mathcal{F}^{\sigma}$ is a normal variation map in the sense that there exists a continuous lapse function $u=u_{\tau}^{\sigma}\colon\Sigma^{\sigma}_{\tau}\to\mathbb{R}$ such that $\partial_{\tau}\mathcal{F}^{\sigma}=u\,\nu$, where $\nu=\nu^{\sigma}_{\tau}$ is the unit normal to $\Sigma^{\sigma}_{\tau}$ in $(M^{3},g)$. Then we have

[TABLE]

and $\mathcal{L}u=O(\sigma^{1-2\delta-2\varepsilon})$.

Proof.

In this proof, $C>0$ and $\overline{\sigma}>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $\delta$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$, and $E$. The surfaces $\Sigma_{\tau}^{\sigma}$ have constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ in the initial data set $\mathcal{I}_{\tau}$. For clarity, we will write this constant spacetime mean curvature with an explicit reference to the initial data set $\mathcal{I}_{\tau}$ as $\mathcal{H}(\Sigma^{\sigma}_{\tau},\mathcal{I}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ for all $\tau\in U$. Hence

[TABLE]

which gives us the following linear elliptic PDE on the closed surface $\Sigma^{\sigma}_{\tau}$ for the a priori only continuous lapse function $u=u^{\sigma}_{\tau}\colon\Sigma^{\sigma}_{\tau}\to\mathbb{R}$:

[TABLE]

where the elliptic operator $\mathcal{L}$ is (up to a certain factor) the linearization of the spacetime mean curvature operator for the surface $\Sigma^{\sigma}_{\tau}$ in the initial data set $\mathcal{I}_{\tau}$, defined in (28). Then Proposition 5.6 implies that $u\in W^{2,2}(\Sigma^{\sigma}_{\tau})$, and that such a $u$ is unique. Together with (26) and $P=O(\sigma^{-\delta-\varepsilon})$, (67) implies that

[TABLE]

As a consequence, by Corollary 5.7 and (50), we get

[TABLE]

In order to estimate $\|u^{t}\|_{L^{2}(\Sigma_{\tau}^{\sigma})}$, note that by (51) we have for $i=1,2,3$

[TABLE]

where we again use $m_{\mathcal{H}}=m_{\mathcal{H}}(\Sigma^{\sigma}_{\tau})$ as an abbreviation for the Hawking mass and the fact that $m_{\mathcal{H}}\neq 0$ as $E\neq 0$ and $|m_{\mathcal{H}}-E|\leq C\sigma^{-\varepsilon}$ by Proposition 5.6. By (52), we have

[TABLE]

Using the Cauchy–Schwarz Inequality, integration by parts, and (35), together with $f_{i}\approx\sqrt{\tfrac{3}{4\pi r^{4}}}\,x^{i}$, (21), and (68) we obtain

[TABLE]

Combining the last three estimates with a Triangle Inequality, it follows that

[TABLE]

so that

[TABLE]

Recalling (69), we now get a $W^{2,2}$-estimate for $u^{d}$, namely

[TABLE]

From this, as a consequence of (68) and Corollary 5.7, we also have

[TABLE]

∎

Lemma 6.6 enables us to prove the following result.

Lemma 6.7**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$, with $\vec{x}$ denoting the asymptotic coordinate chart. Let $\emptyset\neq U\subseteq[0,1]$ be an open subset of $[0,1]$ and define $\mathcal{I}_{\tau}$ as in the proof of Theorem 6.2 for each $\tau\in U$. Let $a\in[0,1)$, $b\geq 0$, $\eta\in(0,1]$ be fixed. Then there exist constants $\overline{\sigma}>0$ and $C>0$, depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}}$, and $E$ such that the following holds for any $\sigma>\overline{\sigma}$: Assume there exists a $C^{1}$-map $\mathcal{F}^{\sigma}\colon U\times\mathbb{S}^{2}\to M^{3}$ such that for every $\tau\in U$ the surface $\Sigma^{\sigma}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathcal{F}^{\sigma}(\tau,\mathbb{S}^{2})$ is in the a priori class $\mathcal{A}(a,b,\eta)$ and has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{\tau}$. Assume further that $\mathcal{F}^{\sigma}$ is a normal variation map in the sense explained in Lemma 6.6. Then

[TABLE]

Proof.

In this proof, $C>0$ and $\overline{\sigma}>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $a$, $b$, $\eta$, and $C_{\mathcal{I}}$, and $E$. Let $u\colon\Sigma^{\sigma}_{\tau}\to\mathbb{R}$ denote the lapse function as in Lemma 6.6. Then Lemma 6.6 applied with $\delta=\tfrac{3}{2}$ and the Sobolev Embedding Theorem in the form of Lemma A.2 imply that $|\partial_{\tau}\mathcal{F}^{\sigma}|=|u|\leq C\sigma^{1-\varepsilon}$. Then the elementary estimate

[TABLE]

proves (70).

In order to prove (71), we first recall that the mean curvature of $\Sigma^{\sigma}_{\tau}$ satisfies $H=\frac{2}{\sigma}+O(\sigma^{-2-\varepsilon})$, see (26). The first variation of area formula, the fact that the eigenfunctions used to span $L^{2}(\Sigma^{\sigma}_{\tau})$ are $L^{2}(\Sigma^{\sigma}_{\tau})$-orthogonal so that in particular $\int_{\Sigma^{\sigma}_{\tau}}u^{t}\,d\mu=0$, combined with Lemma 6.6 for $\delta=\tfrac{3}{2}$ lead to

[TABLE]

where we also used the Cauchy–Schwarz Inequality.

A very similar analysis demonstrates $\left|\partial_{\tau}|\Sigma^{\sigma}_{\tau}|_{\delta}\right|\leq C\sigma^{\frac{3}{2}-\varepsilon}$. Finally, we prove (72): By definition,

[TABLE]

Using the variation of area formula, the Cauchy–Schwarz Inequality, (20), (21), and Lemma 6.6 with $\delta=\tfrac{3}{2}$ we compute

[TABLE]

This proves (72). ∎

Lemma 6.8**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Let $\emptyset\neq U\subseteq[0,1]$ be an open, connected subset of $[0,1]$ and define $\mathcal{I}_{\tau}$ as in the proof of Theorem 6.2 for each $\tau\in U$.

Let $a,\overline{a}\in[0,1)$, $b,\overline{b}\in[0,\infty)$, $\eta\in(0,2\varepsilon)$, and $\overline{\eta}\in(0,1]$. Then there exists a constant $\overline{\sigma}>0$, depending only on $\varepsilon$, $a$, $\overline{a}$, $b$, $\overline{b}$, $\eta$, $\overline{\eta}$, $C_{\mathcal{I}}$, and $E$ such that the following holds for any $\sigma>\overline{\sigma}$: Assume there exists a $C^{1}$-map $\mathcal{F}^{\sigma}\colon U\times\mathbb{S}^{2}\to M^{3}$ such that for every $\tau\in U$ the surface $\Sigma^{\sigma}_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\mathcal{F}^{\sigma}(\tau,\mathbb{S}^{2})$ is in the a priori class $\mathcal{A}(\overline{a},\overline{b},\overline{\eta})$ and has constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{\tau}$. Assume further that $\mathcal{F}^{\sigma}$ is a normal variation map in the sense explained in Lemma 6.6. Now suppose in addition that $\Sigma^{\sigma}_{\tau_{0}}\in\mathcal{A}(a,b,\eta)$ for some $\tau_{0}\in U$. Then in fact $\Sigma^{\sigma}_{\tau}\in\mathcal{A}(a,b_{\tau},\eta_{\tau})$ with $b_{\tau}=b+O(\sigma^{-\min\{2\varepsilon-\eta,\varepsilon\}})$ and $\eta_{\tau}=\eta+O(\sigma^{-\varepsilon})$ for any $\tau\in U$. Here, the constants in the $O$-notation depend only on $\varepsilon$, $a$, $\overline{a}$, $b$, $\overline{b}$, $\eta$, $\overline{\eta}$, $C_{\mathcal{I}}$, and $E$.

Remark 6.9*.*

Note that the assumption $\eta\in(0,2\varepsilon)$ of the lemma is not restrictive as the inclusion $\mathcal{A}(a,b,\eta_{1})\subseteq\mathcal{A}(a,b,\eta_{2})$ for $0<\eta_{2}<\eta_{1}\leq 1$ implies that we may without loss of generality decrease $\eta\in(0,1]$ to achieve $\eta\in(0,2\varepsilon)$.

Proof.

We drop the explicit reference to $\sigma$ for notational convenience, as $\sigma$ will not be modified in this proof. Let $r_{\tau}$ and $\vec{z}_{\tau}$ denote the area radius and the coordinate center of $\Sigma_{\tau}$, respectively, and let (slightly abusing notation) $\vec{x}_{\tau}$ denote the restriction of the coordinate vector $\vec{x}$ to $\Sigma_{\tau}$, where $\vec{x}$ denotes the asymptotic coordinate chart. The mean curvature of $\Sigma_{\tau}$ is denoted by $H_{\tau}$. In this proof $\overline{\sigma}>0$, $C>0$ and constants involved in the $O$-notation may vary from line to line but depend only on $\varepsilon$, $a$, $\overline{a}$, $b$, $\overline{b}$, $\eta$, $\overline{\eta}$, $C_{\mathcal{I}}$, and $E$.

We first show that there exists $\eta_{\tau}=\eta+O(\sigma^{-\varepsilon})$ such that the second inequality describing the fact that $\Sigma_{\tau}\in\mathcal{A}(a,b_{\tau},\eta_{\tau})$ in (15) holds, namely

[TABLE]

Since $\Sigma_{\tau}\in\mathcal{A}(\overline{a},\overline{b},\overline{\eta})$ for all $\tau\in U$, by the Mean Value Theorem combined with (71) and (21) we have

[TABLE]

hence

[TABLE]

Similarly, combining (70) with (20) and (21) we conclude that

[TABLE]

Since $\Sigma_{\tau_{0}}\in\mathcal{A}(a,b,\eta)$ we have

[TABLE]

which in the view of (74) and (75) can be written as

[TABLE]

Consequently, we have

[TABLE]

Choosing

[TABLE]

(73) follows. Note that by (20) we have

[TABLE]

We will now apply a similar method and adjust the value of the constant $b_{\tau}$ so that the first and the third inequality in (15) describing the fact that $\Sigma_{\tau}\in\mathcal{A}(a,b_{\tau},\eta_{\tau})$, that is

[TABLE]

and

[TABLE]

hold with $\eta_{\tau}$ as defined in (76).

First, we deal with (77). Since $\Sigma_{\tau_{0}}\in\mathcal{A}(a,b,\eta)$ we have

[TABLE]

Combining this with (74) and (72) we obtain

[TABLE]

which in the view of (21) may be further rewritten as

[TABLE]

Since $\varepsilon\leq\tfrac{1}{2}$, we conclude that

[TABLE]

Recall that by our definition (76) of $\eta_{\tau}$ we have $\eta=\eta_{\tau}-\log_{r_{\tau}}(1-C\sigma^{-\varepsilon})$. Hence

[TABLE]

Consequently, we have

[TABLE]

Next, we address (78). We recall that $H_{\tau}^{2}=\mathcal{H}^{2}-\tau^{2}P^{2}=\frac{4}{\sigma^{2}}-\tau^{2}P^{2}$ which implies that $\partial_{\tau}H_{\tau}^{2}=O(\sigma^{-3-2\varepsilon})$. Consequently, we may compute using the variation of area formula and (66) with $\delta=\tfrac{3}{2}$ that

[TABLE]

Again, since $\Sigma_{\tau_{0}}\in\mathcal{A}(a,b,\eta)$ we have

[TABLE]

As before, we use the Mean Value Theorem, (80), and (74) to rewrite this as

[TABLE]

which in the view of (21) may be further rewritten as

[TABLE]

Substituting $\eta=\eta_{\tau}-\log_{r_{\tau}}(1-C\sigma^{-\varepsilon})$ in the view of $\varepsilon\leq\tfrac{1}{2}$ gives

[TABLE]

Consequently, we have

[TABLE]

Together, the inequalities (79) and (81) imply the existence of the constant $b_{\tau}=b+O(\sigma^{-\min\{2\varepsilon-\eta,\varepsilon\}})$ such that (77) and (78) hold. This concludes the proof as we can choose $a_{\tau}=a$ as the above computations show. ∎

Lemma 6.10**.**

Under the assumptions of Theorem 6.2, there exists a constant $\overline{\sigma}>0$, depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}_{1}}$, and $E$, and a compact set $\mathcal{K}\subset M^{3}$ such that the map $\Psi_{1}\colon(\overline{\sigma},\infty)\times\mathbb{S}^{2}\to M^{3}\setminus\mathcal{K}$ defined by (65) is a bijective $C^{1}$-map.

Proof.

To prove the claim, we need to show that $\Psi_{1}$ is $C^{1}$, injective, and surjective onto a suitably chosen exterior region of $M^{3}$. We already proved in Theorem 6.2 that $\mathcal{F}^{\sigma}$ and thus $\Psi_{1}$ is $C^{1}$ with respect to the $\mathbb{S}^{2}$-component. The differentiability with respect to $\sigma$ can be proven following the Implicit Function Theorem argument of Lemma 6.4, where the graphical spacetime mean curvature map is to be interpreted as a function of $\sigma\in(\overline{\sigma},\infty)$ instead of as a function of $\tau\in[0,1]$. This is to be viewed in light of the uniqueness results in Section 6.3.

Injectivity.

In order to show injectivity of $\Psi_{1}$, we need to assert that $\Sigma^{\sigma_{1}}_{1}\cap\Sigma^{\sigma_{2}}_{1}=\emptyset$ for any choice of $\sigma_{1}\neq\sigma_{2}$, $\sigma_{1},\sigma_{2}>\overline{\sigma}$. This can be done by analyzing the lapse function of the variation $\Psi_{1}\colon(\overline{\sigma},\infty)\times\mathbb{S}^{2}\to M^{3}$ with respect to $\sigma$, namely $u=u^{\sigma}_{1}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}g(\partial_{\sigma}\Psi_{1},\nu)$, where $\nu=\nu^{\sigma}_{1}$ denotes the outward unit normal of $\Sigma^{\sigma}_{1}\hookrightarrow(M^{3},g)$. If $u>0$ on $(\overline{\sigma},\infty)\times\mathbb{S}^{2}$, we can conclude that $\Psi_{1}$ is injective.

We will in fact show that $u=1+O(\sigma^{-\frac{1}{2}-\varepsilon})$. Again, $C>0$ and $\overline{\sigma}>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}_{1}}$, and $E$. First, we note that since the spacetime mean curvature of $\Sigma^{\sigma}_{1}$ in the initial data set $\mathcal{I}_{1}$ is constant, $\mathcal{H}(\Sigma^{\sigma}_{1})\equiv\nicefrac{{2}}{{\sigma}}$, we have

[TABLE]

which uniquely determines $u\in W^{2,2}(\Sigma^{\sigma}_{1})$ by Proposition 5.6. Furthermore, by (82), (26), (16), and the asymptotic decay assumptions on $\mathcal{I}_{1}$, we have

[TABLE]

where $H=H(\Sigma^{\sigma}_{1})$, $P=P(\Sigma^{\sigma}_{1})$. This shows that $\mathcal{L}(u-1)=O(\sigma^{-\frac{5}{2}-\varepsilon})$ which is not sufficient for concluding that $u=1+O(\sigma^{-\varepsilon})$ and thus $u>0$ via Corollary 5.7 and Lemma A.2: such an argument would require $\mathcal{L}(u-1)=O(\sigma^{-3-\varepsilon})$ which we obviously do not have.

Instead, we argue as in the proof of Lemma 6.6. For $v\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}u-1$, the above computation shows that $\mathcal{L}v=O(\sigma^{-\frac{5}{2}-\varepsilon})$, which in combination with Corollary 5.7 and (50) gives

[TABLE]

In addition, for $i=1,2,3$, by adding a rich zero and using the orthogonality of $v^{d}$ and $f_{i}$, we have

[TABLE]

where

[TABLE]

by (51), and

[TABLE]

by (52). Using the fact that $\mathcal{L}v=\mathrm{Ric}(\nu,\nu)+O(\sigma^{-3-\varepsilon})$ by (83), Lemma A.3 from [36] (a result showing that $\int_{\Sigma^{\sigma}_{1}}\mathrm{Ric}(\nu,\nu)x^{i}d\mu=O(\sigma^{-\varepsilon})$ which is fully Riemannian and thus directly carries over to our spacetime context), and integration by parts, we obtain as in Lemma 6.6 that

[TABLE]

Combining these estimates, we get, again grouping terms as in Lemma 6.6, that

[TABLE]

which, together with (84), gives

[TABLE]

Finally, Corollary 5.7 gives us

[TABLE]

Then $\|v\|_{W^{2,2}(\Sigma^{\sigma})}\leq C\sigma^{1-\varepsilon}$, so by Lemma A.2 we get that $u=1+v=1+O(\sigma^{-\varepsilon})$ is strictly positive for all $\sigma>\overline{\sigma}$. This shows that $\Psi_{1}$ is indeed injective.

Surjectivity.

By construction, the STCMC-surfaces $\Sigma^{\sigma}=\Psi_{1}(\sigma,\mathbb{S}^{2})$ for $\sigma>\overline{\sigma}$ are in the class of asymptotically centered surfaces $\mathcal{A}(0,b,\eta)$ for some $b>0$ and $\eta\in(0,1]$. In particular, recalling Proposition 4.5, each $\Sigma^{\sigma}$ can be written as a graph over a sphere enclosing the interior region of $M^{3}$. Suppose $p\in M^{3}$ is in the exterior region of some $\Sigma^{\sigma}$ with $\sigma>\overline{\sigma}$. By comparability of the coordinate and mean curvature radii for surfaces in the class $\mathcal{A}(0,b,\eta)$ (see (20) and (21)), we can find $\widetilde{\sigma}>\sigma$ such that $p$ lies in the region enclosed by $\Sigma^{\widetilde{\sigma}}$, and hence in the annulus $A_{\sigma,\widetilde{\sigma}}$ between $\Sigma^{\sigma}=\Psi_{1}(\sigma,\mathbb{S}^{2})$ and $\Sigma^{\widetilde{\sigma}}=\Psi_{1}(\widetilde{\sigma},\mathbb{S}^{2})$. Since $\Psi_{1}\colon[\sigma,\infty)\times\mathbb{S}^{2}\to M^{3}$ is continuous it follows that $A_{\sigma,\widetilde{\sigma}}=\Psi_{1}([\sigma,\widetilde{\sigma}]\times\mathbb{S}^{2})$ hence $p=\Psi_{1}(\hat{\sigma},q)$ for some $\hat{\sigma}\in[\sigma,\widetilde{\sigma}]$ and $q\in\mathbb{S}^{2}$. As $\sigma>\overline{\sigma}$ was arbitrary, this proves surjectivity. ∎

6.3. Uniqueness of the STCMC surfaces.

We close this section by proving that the constant spacetime mean curvature surfaces are unique in the a priori class of asymptotically centered surfaces $\mathcal{A}(a,b,\eta)$. As Brendle and Eichmair [8] constructed examples of asymptotically Euclidean Riemannian manifolds with “off-center” (i.e. not included in the a priori class) CMC-surfaces provided $\mathrm{Scal}\geq 0$ is violated, we will restrict our uniqueness statements to the a priori class — at least when not assuming the dominant energy condition $\mu\geq|J|_{g}$. To the best knowledge of the authors, it is not known whether such examples can also be constructed if the dominant energy condition or its Riemannian analog $\mathrm{Scal}\geq 0$ are satisfied.

The uniqueness result is proven in a similar way as the existence result of Section 6.1, namely by the method of continuity. The “starting point” of the method of continuity is the following result from [36], again adapted to our notation.

Theorem 6.11** (Nerz 2015).**

Let $a\in[0,1)$, $b\geq 0$, and $\eta\in(0,1]$ be constants and let $(M^{3},g)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean manifold viewed as a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}_{0}=(M^{3},g,K\equiv 0,\mu=\tfrac{1}{2}\mathrm{Scal},J\equiv 0)$ with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}_{0}}$ depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}_{0}}$, and $E$, such that for all $\sigma>\sigma_{\mathcal{I}_{0}}$, there is a unique surface $\Sigma^{\sigma}_{0}\in\mathcal{A}(a,b,\eta)$ with constant mean curvature $H(\Sigma^{\sigma}_{0})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to $\mathcal{I}_{0}$.

Our uniqueness result generalizes this to STCMC-surfaces in the spacetime context.

Theorem 6.12** (Uniqueness of STCMC-foliation).**

Let $a\in[0,1)$, $b\geq 0$, and $\eta\in(0,1]$ be constants and let $\mathcal{I}_{1}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Then there is a constant $\sigma_{\mathcal{I}_{1}}$ depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}_{1}}$, and $E$, such that for all $\sigma>\sigma_{\mathcal{I}_{1}}$, there is a unique surface $\Sigma^{\sigma}_{1}\in\mathcal{A}(a,b,\eta)$ with constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{1})\equiv\nicefrac{{2}}{{\sigma_{1}}}$ with respect to $\mathcal{I}_{1}$.

Proof.

We rely on the same type of argument as in the proof of Theorem 6.2. Fix a surface $\Sigma^{\sigma}_{1}$ as in the assumptions, with $\sigma>\sigma_{\mathcal{I}_{1}}$. We now drop the explicit reference to $\sigma$ for notational convenience, as $\sigma$ will not be modified in this proof. Let $Z\subseteq[0,1]$ be a maximal subset such that there is a $C^{1}$-map $\Phi\colon Z\times\mathbb{S}^{2}\to M^{3}$ with the following properties for all $\tau\in Z$:

1. (i)

$\Phi(1,\mathbb{S}^{2})=\Sigma_{1}$, 2. (ii)

$\Sigma_{\tau}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Phi(\tau,\mathbb{S}^{2})$ has constant spacetime mean curvature $\mathcal{H}(\Sigma_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to the initial data set $\mathcal{I}_{\tau}$, where $\mathcal{I}_{\tau}$ is defined as in the proof of Theorem 6.2, 3. (iii)

$\partial_{\tau}\Phi$ is orthogonal to $\Sigma_{\tau}$.

Maximality is understood as in the proof of Theorem 6.2. Arguing as in the proof of Theorem 6.2, we conclude that $Z=[0,1]$ and that there are constants $a_{\mathcal{I}_{1}}\in[0,1)$, $b_{\mathcal{I}_{1}}\geq 0$ and $\eta_{\mathcal{I}_{1}}\in(0,1]$ such that $\Sigma_{\tau}\in\mathcal{A}(a_{\mathcal{I}_{1}},b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ for every $\tau\in[0,1]$, if $\sigma_{\mathcal{I}_{1}}$ suitably large, depending only on $\varepsilon$, $a$, $b$, $\eta$, $C_{\mathcal{I}_{1}}$, and $E$. In particular, we see that $\Sigma_{0}\in\mathcal{A}(a_{\mathcal{I}_{1}},b_{\mathcal{I}_{1}},\eta_{\mathcal{I}_{1}})$ is a surface with constant mean curvature $H(\Sigma_{0})\equiv\nicefrac{{2}}{{\sigma}}$ with respect to $\mathcal{I}_{0}$444Note that although it was assumed throughout the proof of Theorem 6.2 that $a=0$, this proof extends straightforwardly to deal with the general case $a\in[0,1)$. In fact, in the view of Lemma 6.8, we may set $a_{\mathcal{I}_{1}}=a$.. By Theorem 6.11, such a surface is unique in this class. By the method of continuity approach and the local uniqueness in the Implicit Function Theorem, the map $\Phi$ is uniquely determined also by its start value $\Phi(0,\mathbb{S}^{2})=\Sigma_{0}$. It follows directly that $\Sigma_{1}=\Phi(1,\mathbb{S}^{2})$ is uniquely determined by its spacetime mean curvature in $\mathcal{I}_{1}$. ∎

7. The coordinate center of the STCMC-foliation

Let $\{\Sigma^{\sigma}\}_{\sigma>\overline{\sigma}}$ be a foliation of a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ for which $\Sigma^{\sigma}$ grows to the round sphere at infinity as $\sigma\to\infty$. Then we may define the coordinate center of this foliation as the limit $\lim_{\sigma\to\infty}\vec{z}^{\,\sigma}$, where $\vec{z}^{\,\sigma}=\vec{z}\,(\Sigma^{\sigma})$ is the coordinate center of $\Sigma^{\sigma}$ as defined in Definition 4.1, provided that this limit exists (in $\mathbb{R}^{3}$). We would like to draw the attention of the reader to the fact that, while the foliations considered here do not depend on the choice of asymptotic coordinates $\vec{x}$, the coordinate centers $\vec{z}^{\,\sigma}$ and as a consequence also their limit, do depend on $\vec{x}$. We will discuss the subtle consequences of this within this section, too.

Let us first consider the case of the CMC-foliation: In this case, $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{\mathcal{I}_{0}}}$ is the unique foliation of a given $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean manifold $(M^{3},g)$ or initial data set $\mathcal{I}_{0}=(M^{3},g,K\equiv 0,\mathrm{Scal},J\equiv 0)$ by surfaces of constant mean curvature constructed in [26, 32, 36] and discussed in Section 6 above. Under the additional assumption that $(M^{3},g)$ satisfies the Riemannian $C^{2}_{1+\varepsilon}$-Regge–Teitelboim condition (see Definition 2.5), the coordinate center of this foliation is well-defined if and only if the Beig–Ó Murchadha center of mass $\vec{C}_{\text{B\'{O}M}}$ given by (6) is well-defined as was shown by [25, 36]. They also show that in this case, one has

[TABLE]

Now suppose that $\{\Sigma^{\sigma}\}_{\sigma>\sigma_{\mathcal{I}_{1}}}$ is the unique foliation of a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ by surfaces of constant spacetime mean curvature as constructed in Theorem 6.2. One cannot in general expect that (85) also holds for the STCMC-foliation because the foliation is defined in terms of $K$, whereas the Beig–Ó Murchadha center of mass is a purely Riemannian quantity, i.e. independent of $K$. One can only expect that (85) will hold if $K$ falls off very fast, in particular faster than the optimal decay assumed in this paper. In this section we will confirm that this is indeed the case.

7.1. A variational formula for STCMC-surfaces

The following proposition generalizes [36, Proposition 6.5] to the spacetime case.

Proposition 7.1**.**

For $b\geq 0$ and $\eta\in(0,1]$, let $\Sigma\in\mathcal{A}(0,b,\eta)$ be a closed, oriented constant spacetime mean curvature $2$-surface with spacetime mean curvature $\mathcal{H}\equiv\nicefrac{{2}}{{\sigma}}$, with outer unit normal denoted by $\nu$, in a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$. Consider a $C^{1}$-map $\mathcal{F}\colon(-s_{0},s_{0})\times\Sigma\rightarrow M^{3}$ such that $\mathcal{F}(0,\cdot)=\operatorname{Id}_{\Sigma}$. For each $s\in(-s_{0},s_{0})$, we let $\vec{z}_{s}=(z_{s}^{1},z_{s}^{2},z_{s}^{3})$ denote the coordinate center of the surface $\Sigma_{s}=\mathcal{F}(s,\Sigma)$ and consider the lapse $u\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}g((\partial_{s}\mathcal{F})\left|\right._{s=0},\nu)$ of the foliation.

Then there are constants $C>0$ and $\overline{\sigma}>0$ depending only on $\varepsilon$, $b$, $\eta$, $C_{\mathcal{I}}$ such that

[TABLE]

provided that $\sigma>\overline{\sigma}$.

Proof.

Since the coordinate center of a surface is invariant under tangential diffeomorphisms (along $\Sigma$), we may without loss of generality assume that $\mathcal{F}$ is a normal variation of $\Sigma$, such that in particular $\left.(\partial_{s}\mathcal{F})\right|_{s=0}=u\nu$ holds on $\Sigma$. By definition,

[TABLE]

Using the variation of area formula and adding rich zeroes in the third and the forth lines, we compute, dropping the explicit reference to $\delta$ in the denominator,

[TABLE]

Now subtract and add the component $(\nu^{\delta})^{i}$ of the $\delta$-outward unit normal $\nu^{\delta}$ to $\Sigma$ in the bracket of the second term and recall the fact that $\Sigma$ can be written as a graph over $\mathbb{S}^{2}_{r}(\vec{z}\,)$ with graph function $f$ satisfying $\|f\|_{W^{2,\infty}}=O(r^{\frac{1}{2}-\varepsilon})$ by Proposition 4.5. Then, by comparability of $|\vec{x}\,|$ and $r$ as established in Proposition 4.5, we find, recalling $a=0$ in our case,

[TABLE]

Applying the Cauchy–Schwarz Inequality to the above identity for $(\partial_{s}z^{i}_{s})\left|\right._{s=0}$ and using (87), Lemma A.1, and (26) to estimate the individual terms, respectively, we obtain (86). ∎

7.2. STCMC-center of mass

In Section 6, we constructed the unique STCMC-foliation of a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$, i.e. the unique foliation by surfaces $\{\Sigma_{1}^{\sigma}\}_{\sigma>\sigma_{1}}$ of constant spacetime mean curvature $\mathcal{H}(\Sigma_{1}^{\sigma})\equiv\nicefrac{{2}}{{\sigma}}$, provided it has non-vanishing energy $E\neq 0$. This was achieved by deforming the constant mean curvature foliation $\{\Sigma_{0}^{\sigma}\}_{\sigma>\sigma_{0}}$ of the $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean manifold $(M^{3},g)$ from Theorem 6.1 along the curve of initial data sets $\{\mathcal{I}_{\tau}\}_{\tau\in[0,1]}$, where $\mathcal{I}_{\tau}=(M^{3},g,\tau K,\mu_{\tau},\tau J)$ is as described in Section 3.1. We will now apply Proposition 7.1 to find how the coordinate center of a leaf changes under this particular deformation. As a result, we prove Lemma 7.2 relating the respective coordinate centers $\vec{z}^{\,\sigma}_{0}$ and $\vec{z}^{\,\sigma}_{1}$ of the surfaces $\Sigma^{\sigma}_{0}$ and $\Sigma^{\sigma}_{1}$. Note that for the proof of this result, it is necessary to assume that the fall-off rate of $K$ is $K=O(|\vec{x}\,|^{-2})$, which is faster than we originally assumed in Definition 2.1 and in particular faster than one needs for existence and uniqueness of the foliation. See also Conjecture 7.9 below.

Lemma 7.2**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be an STCMC-foliated $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with non-vanishing energy $E\neq 0$. Assume in addition that

[TABLE]

for $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$, with $\vec{x}$ the asymptotic chart. Then there exist constants $C>0$, $\overline{\sigma}>0$ depending only on $\varepsilon$, $C_{\mathcal{I}}$, and $E$ such that for $i=1,2,3$ we have for all $\sigma>\overline{\sigma}$

[TABLE]

where $(z_{0}^{\sigma})^{i}$ and $(z_{1}^{\sigma})^{i}$ denote the components of the coordinate centers $\vec{z}^{\,\sigma}_{0}$ and $\vec{z}^{\,\sigma}_{1}$ of the STCMC-surfaces $\Sigma^{\sigma}_{0}$ and $\Sigma^{\sigma}_{1}$ with respect to $\mathcal{I}_{0}$ and $\mathcal{I}_{1}$ as defined above, respectively.

Remark 7.3*.*

Instead of assuming $K=O(|\vec{x}\,|^{-2})$, one could also assume Regge–Teitelboim conditions on $K$ by carefully tracking all even and odd parts or identify other sufficient decay conditions, as necessary for what one wants to do. For our purposes, it is enough to assume $K=O(|\vec{x}\,|^{-2})$.

Remark 7.4*.*

Assumption (88), Equation (94) in the proof below, and the Mean Value Theorem imply that there exists a constant $C>0$ depending only on $\varepsilon$, $C_{\mathcal{I}}$, and $E$ such that

[TABLE]

However, either or both of the limits

[TABLE]

may fail to exist. On the other hand, $\vec{C}_{\text{STCMC}}=\lim_{\sigma\to\infty}\vec{z}^{\,\sigma}_{1}$ converges if and only if $\lim_{\sigma\to\infty}\left\{\vec{z}^{\,\sigma}_{0}-\frac{1}{32\pi E}\int_{\mathbb{S}^{2}_{\sigma}}\frac{\left(\sum_{k,l}\pi_{kl}x^{k}x^{l}\right)^{2}\vec{x}}{\sigma^{3}}\,d\mu^{\delta}\right\}$ converges. This in particular shows that $K$ can in a sense “compensate” for the diverging coordinate center of the CMC-foliation. See Section 9 for more details on this.

Proof.

First, note that the constants $C_{\mathcal{I}_{\tau}}$ are uniformly bounded by the constant $C_{\mathcal{I}}$. Second, pick constants $b\geq 0$ and $\eta\in(0,1]$ that will remain fixed in this argument and always use the class $\mathcal{A}(0,b,\eta)$ in what follows. Also, $C>0$ and $\overline{\sigma}>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $C_{\mathcal{I}}$, and $|E|$ (as well as on our global choice of $b$ and $\eta$). From now on we assume that $\sigma>\overline{\sigma}$ is fixed. For our choice of $\sigma$ and for $\tau\in[0,1]$, we let $z_{\tau}^{i}=(z_{\tau}^{\sigma})^{i}$, $i=1,2,3$, denote the components of the coordinate center $\vec{z}\,(\Sigma^{\sigma}_{\tau})$ of the unique surface $\Sigma^{\sigma}_{\tau}$ of constant spacetime mean curvature $\mathcal{H}(\Sigma^{\sigma}_{\tau})\equiv\nicefrac{{2}}{{\sigma}}$ in the initial data set $\mathcal{I}_{\tau}$ (see Section 6 for details). Since the index $\sigma$ is assumed to be fixed, it will be suppressed in the remainder of this proof.

According to Proposition 7.1, the variation of the coordinate center with respect to $\tau$ is given by the formula

[TABLE]

where $u_{\tau}$ is the respective lapse function for an arbitrary $\tau\in[0,1]$ and $\nu_{\tau}$ is the outward pointing unit normal to $\Sigma_{\tau}\hookrightarrow(M^{3},g)$. In order to pass from (90) to (89), we will apply Lemma 6.6 with $\delta=2-\varepsilon\geq\tfrac{3}{2}$. By this, we have that $\mathcal{L}u_{\tau}=O(\sigma^{-3})$, $\|(u_{\tau})^{t}\|_{W^{2,2}(\Sigma_{\tau})}\leq C\sigma$ and $\|(u_{\tau})^{d}\|_{W^{2,2}(\Sigma_{\tau})}\leq C\sigma^{\frac{1}{2}-\varepsilon}$.

Next, by (33) and Lemma A.1 we have

[TABLE]

where $f_{\tau}^{i}$ denotes the $i$-th eigenfunction of the operator $-\Delta^{\Sigma_{\tau}}$, see Section 5.2. Then we may rewrite (90) by a Cauchy–Schwarz Inequality and Lemma 6.6 with $\delta=2-\varepsilon$ as

[TABLE]

At the same time, Proposition 5.6 implies that

[TABLE]

where $m_{\mathcal{H}}$ is the Hawking mass of $\Sigma_{\tau}$ with respect to $\mathcal{I}_{\tau}$. Recall that $|m_{\mathcal{H}}-E|\leq C\sigma^{-\varepsilon}$ by Proposition 5.6 so that $m_{\mathcal{H}}\neq 0$ follows from $E\neq 0$. Thus,

[TABLE]

Note that a computation in the proof of Lemma 6.6 shows that

[TABLE]

hence $\int_{\Sigma_{\tau}}f_{\tau}^{i}\,\mathcal{L}(u_{\tau})^{t}\,d\mu=\int_{\Sigma_{\tau}}f_{\tau}^{i}\,\mathcal{L}u_{\tau}\,d\mu+O(\sigma^{-2-2\varepsilon})$. Consequently, in view of (67) and (91), and the Cauchy–Schwarz Inequality, (92) is equivalent to

[TABLE]

Since the expression in the integral is of order $O(\sigma^{-3})$ as mentioned above, and since

[TABLE]

using as before the fact that $m_{\mathcal{H}}=E+O(\sigma^{-\varepsilon})$ along the STCMC-foliation, we conclude by Lemma A.1 that

[TABLE]

Integrating this with respect to $\tau$ over $[0,1]$, we obtain (89). ∎

Theorem 7.5** (STCMC-coordinate expression).**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with respect to an asymptotic coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and decay constant $C_{\mathcal{I}}$, with non-vanishing energy $E\neq 0$. Assume in addition that

[TABLE]

for all $\vec{x}\in\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and that $g$ satisfies the Riemannian $C^{2}_{\nicefrac{{3}}{{2}}+\varepsilon}$-Regge–Teitelboim condition. Then the coordinate center $\vec{C}_{\text{STCMC}}$ of the unique foliation by surfaces of constant spacetime mean curvature is well-defined if and only if the correction term

[TABLE]

limits exist for $i=1,2,3$. In this case, we have

[TABLE]

where $\vec{C}_{\text{B\'{O}M}}$ is the Beig–Ó Murchadha center of mass and $\vec{Z}=(Z^{1},Z^{2},Z^{3})$, or equivalently

[TABLE]

Proof.

Since $g$ satisfies the Riemannian $C^{2}_{\nicefrac{{3}}{{2}}+\varepsilon}$-Regge–Teitelboim condition, $\vec{C}_{\text{CMC}}$ is well-defined and equal to $\vec{C}_{\text{B\'{O}M}}$, see [36, Theorem 6.3]. The result is then a direct consequence of Lemma 7.2. ∎

Remark 7.6*.*

We expect that the coordinate center $\vec{C}_{\text{STCMC}}$ of the spacetime mean curvature foliation translates in a certain sense to the center of mass $\vec{C}_{\text{Sz}}$ as defined by Szabados [41]. Similar to (96), the definition of Szabados takes form $\vec{C}_{\text{Sz}}=\vec{C}_{\text{B\'{O}M}}+\vec{S}$, where $\vec{S}$ is characterized by the extrinsic curvature $K$ of the initial data set and the time function $t$ which realizes this initial data set as a slice in an asymptotically Minkowskian spacetime. The relation between $\vec{C}_{\text{STCMC}}$, $\vec{C}_{\text{Sz}}$, and the Chen–Wang–Yau center of mass defined via optimal isometric embeddings in Minkowski spacetime (see [13]) will be studied in detail in our forthcoming work.

Definition 7.7**.**

We suggest to call the expression $\vec{C}_{\text{mB\'{O}M}}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\vec{C}_{\text{B\'{O}M}}+\vec{Z}$ the modified Beig–Ó Murchadha center of mass.

In Section 9, we will give an example that shows that the contribution of the correction term $\vec{Z}$ is indeed relevant and fixes a problem of the CMC-center of mass uncovered in [10].

Remark 7.8*.*

It is not obvious which decay conditions on $\mathcal{I}$ (e.g. versions of Regge–Teitelboim, faster decay assumptions on $K$, etc.) are sufficient to ensure convergence of the correction term $\vec{Z}$ without forcing it to vanish entirely. This will be studied in detail in our forthcoming work. More importantly, sufficient conditions for convergence of $\vec{C}_{\text{STCMC}}$ that do not force vanishing of $\vec{Z}$ in accordance with the example studied in Section 9 will also be studied in our forthcoming work.

We conjecture the following sufficient conditions, in line with Bartnik’s [3] and Chruściel’s [16] corresponding results for convergence of ADM-energy and ADM-linear momentum.

Conjecture 7.9**.**

We conjecture that the coordinate expression we derived will converge for asymptotic coordinates $\vec{x}$ if $\mu x^{i}\in L^{1}(M^{3})$.

8. Time evolution and Poincaré covariance

of the STCMC-center of mass

8.1. Evolution

In this section, we will study the evolution of the coordinate center of the unique foliation by surfaces of constant spacetime mean curvature under the Einstein evolution equations. We will show that the STCMC-center of mass has the same evolution properties as a point particle in special relativity, evolving according to the formula

[TABLE]

Note that the analogous formula is valid for the CMC-center of mass and also for Chen–Wang–Yau’s center of mass, although under stronger decay assumptions, see [34] and [13], respectively.

Theorem 8.1** (Time evolution of STCMC-foliation).**

Let $(\mathbb{R}\times M^{3},\mathfrak{g})$ be a smooth, globally hyperbolic Lorentzian spacetime satisfying the Einstein equations with energy momentum tensor $\mathfrak{T}$. Suppose that, outside a set of the form $\mathbb{R}\times\mathcal{K}$, $\mathcal{K}\subset M^{3}$ compact, there is a diffeomorphism $\operatorname{Id}_{\mathbb{R}}\times\,\vec{x}\colon\mathbb{R}\times(M^{3}\setminus\mathcal{K})\to\mathbb{R}\times(\mathbb{R}^{3}\setminus\overline{B_{R}(0)})$ which gives rise to asymptotic coordinates $(t,\vec{x})$ on $\mathbb{R}\times(M^{3}\setminus\mathcal{K})$.

Assume that $\mathcal{I}_{0}=(\{0\}\times M^{3},g,K,\mu,J)\hookrightarrow(\mathbb{R}\times M^{3},\mathfrak{g})$ is a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with respect to the coordinate chart $\vec{x}$ and with $E\neq 0$, and suppose additionally that $K=O_{1}(|\vec{x}\,|^{-2})$ with constant $C_{\mathcal{I}}$ as $|\vec{x}\,|\to\infty$. Now consider the $C^{1}$-parametrized family of $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data sets

[TABLE]

with respect to $\vec{x}$ which starts from $\mathcal{I}(0)=\mathcal{I}_{0}$, and which exists for all $t\in(-t_{*},t_{*})$ for some $t_{*}>0$. Assume furthermore that the constants $C_{\mathcal{I}(t)}$ are uniformly bounded on $(-t_{*},t_{*})$, without loss of generality such that $C_{\mathcal{I}(t)}\leq C_{\mathcal{I}_{0}}$.

Assume the foliation $\mathcal{I}(t)$ has initial lapse $N=1+O_{2}(|\vec{x}\,|^{-\frac{1}{2}-\varepsilon})$ as $|\vec{x}\,|\to\infty$ with decay measuring constant denoted by $C_{N}$ and initial shift $X=0$, and suppose furthermore that the initial stress tensor $S$ of $\mathcal{I}_{0}$ satisfies $S=O(|\vec{x}\,|^{-\frac{5}{2}-\varepsilon})$ as $|\vec{x}\,|\to\infty$. There is a constant $\overline{t}>0$, depending only on $\varepsilon$, $C_{\mathcal{I}_{0}}$, $C_{N}$, and $E(0)$ such that the following holds: If the initial data set $\mathcal{I}_{0}$ has well-defined STCMC-center of mass $\vec{C}_{\text{STCMC}}\,(0)$ then the STCMC-center of mass $\vec{C}_{\text{STCMC}}\,(t)$ of $\mathcal{I}(t)$ is also well-defined for $|t|<\overline{t}$. Furthermore, the initial velocity at $t=0$ is given by

[TABLE]

Moreover, we have that $\left.\frac{d}{dt}\right|_{t=0}E=0$ and $\left.\frac{d}{dt}\right|_{t=0}\vec{P}=\vec{0}$.

Remark 8.2*.*

In fact, one gets more information about the evolution of the STCMC-center of mass from the proof of Theorem 8.1: Not only the coordinate expression $\vec{C}_{\text{STCMC}}$ evolves according to (97), but also the individual leaves of the ‹STCMC-foliation evolve in a way more and more close to a translation in direction $\nicefrac{{\vec{P}}}{{E}}$ according to formula (106).

Remark 8.3*.*

We expect that the evolution of the leaves of the STCMC-foliation as well as the evolution of $\vec{C}_{\text{STCMC}}$ can actually be understood when replacing the condition $K=O_{1}(|\vec{x}\,|^{-2})$ by more natural conditions related to integrability criteria on the constraints when integrated against $\vec{x}$. We will investigate this in our forthcoming work, see also Remark 7.8 and Conjecture 7.9.

Remark 8.4*.*

It is straightforward to prove a version of this theorem allowing for non-vanishing shift. As this is not of primary interest here and can also be fixed by a suitable gauge, we will not go in this direction.

Proof.

Throughout this proof, dotted quantities like for example $\dot{E}$ will denote time derivatives at $t=0$, e.g. $\dot{E}=\left.\frac{d}{dt}E\right|_{t=0}$. Moreover, $\overline{\sigma}>0$, $\overline{t}>0$, and $C>0$ denote generic constants that may vary from line to line, but depend only on $\varepsilon$, $C_{\mathcal{I}_{0}}$, and $E(0)$ as well as on $C_{N}(0)$, the constant in the $O$-term of $N$.

The Einstein evolution equations with zero shift are given at $t=0$ by

[TABLE]

where $\mathfrak{Ric}$ is the Ricci tensor of the spacetime $(\mathbb{R}\times M^{3},\mathfrak{g})$, which is completely determined by the stress-energy tensor $\mathfrak{T}$ through the Einsteins equations

[TABLE]

Here, we used that $\mathfrak{Scal}=-\operatorname{tr}_{\mathfrak{g}}\mathfrak{T}=O(|\vec{x}\,|^{-\frac{5}{2}-\varepsilon})$ and thus $\mathfrak{Ric}_{ij}=O(|\vec{x}\,|^{-\frac{5}{2}-\varepsilon})$. In particular, we see from the ADM-formulas (3) and (4) respectively that the energy and linear momentum satisfy $\dot{E}=0$ and $\dot{\vec{P}}=\vec{0}$, with respect to this variation. Also, as $E(0)\neq 0$ and $E(t)$ is continuous, the initial data sets $\mathcal{I}(t)$ have a unique foliation by surfaces of constant spacetime mean curvature near infinity for $|t|<\overline{t}$ and mean curvature radii $\sigma>\overline{\sigma}$. This foliation depends in a $C^{1}$-fashion on $t$ for $|t|<\overline{t}$ which can be seen as follows: Perform a method of continuity procedure around $t=0$ as in the proof of Theorems 6.2, 6.12 or in the proof of Lemma 6.10. Note that, as the initial shift $X$ was chosen to vanish, $X=0$, we in fact know that the STCMC-surface variation is normal at $t=0$. This gives you an STCMC-foliation of $\mathcal{I}(t)$ for each $|t|<\overline{t}$ which depends on $t$ in a $C^{1}$-fashion. By Theorem 6.12, this family of foliations must coincide with the one studied here and must thus depend on $t$ in a $C^{1}$-fashion.

Now fix $\sigma>\overline{\sigma}$ and let $\Sigma_{t}$ denote the unique leaf of the STCMC-foliation with constant spacetime mean curvature $\mathcal{H}(\Sigma_{t},\mathcal{I}(t))\equiv\nicefrac{{2}}{{\sigma}}$ in the initial data set $\mathcal{I}(t)$, where $|t|<\overline{t}$. For this, we use the product rule to see that

[TABLE]

The first term on the right shows how the spacetime mean curvature changes if $\Sigma_{t}$ is considered to be a varying surface in the initial data set $\mathcal{I}(0)$. As the initial shift vanishes and we thus have an initially normal variation, this term exactly gives our well-known linearization $L^{\mathcal{H}}u$, where the operator $L^{\mathcal{H}}$ is as defined in Lemma 5.1 and $u$ is the initial lapse function of the normal variation. The second term on the right shows how the spacetime mean curvature changes if $\Sigma$ is considered to be a fixed surface in the varying initial data set $\mathcal{I}(t)$. More precisely, this term is

[TABLE]

where $H=H(0)$ and $P=P(0)$.

In order to compute $\dot{H}=\left.\frac{d}{dt}\right|_{t=0}H(\Sigma_{0},\mathcal{I}(t))$, we introduce geodesic normal coordinates in a neighborhood $U\subset M^{3}$ of $\Sigma_{0}$, with $y^{n}$ such that $\partial_{n}$ is the outer unit normal to the level set $\{y^{n}=\text{const.}\}$, in particular $\partial_{n}=\nu$ on $\Sigma_{0}$, and $y^{\alpha}$, $\alpha=1,2$, are some coordinates on $\Sigma_{0}$ transported to $U$ along the flow generated by $\partial_{n}$. Note that in this case $g_{nn}=1$, $g_{n\alpha}=0$, and

[TABLE]

in $U$, for all $\alpha,\beta,\gamma=1,2$. We will now drop the index on $\Sigma_{0}$ and just write $\Sigma$ instead for notational convenience. We use the standard formula for the variation of the second fundamental form when the ambient metric is changing (see e.g. Section 3 in [29]555note that our sign convention for the second fundamental form is the opposite of [29].) and compute, using first (98) and (99), second the decay properties of $N$ and the decay estimate for the second fundamental form $A=\tfrac{H}{2}g^{\Sigma}+\mathring{A}$, with $g^{\Sigma}$ the metric induced on $\Sigma$ by $\mathcal{I}(t)$, namely

[TABLE]

from Proposition 4.5, third adding some rich zeros, fourth because $J=O(\sigma^{-3-\varepsilon})$ by assumption, fifth by (101), (102), (103), and finally (26), and Proposition 4.5 to obtain

[TABLE]

where $J$ is the momentum density defined on page 2. Further, let $\eta$ denote the timelike future unit normal vector field to $M^{3}\hookrightarrow(\mathbb{R}\times M^{3},\mathfrak{g})$. Then it is straightforward to check that, by (98), (99), the decay assumptions on the initial data set and on $N$, as well as the definition of $\mu$ and $S$ from page 2

[TABLE]

Summing up and multiplying by $\sqrt{1-(\tfrac{P}{H})^{2}}$, it follows from (100) that

[TABLE]

where the operator $\mathcal{L}$ is given by (28). This uniquely defines $u\in W^{2,2}(\Sigma)$ by the invertibility of $\mathcal{L}$, see Proposition 5.6 as the right hand side is bounded and thus in $L^{2}(\Sigma)$.

In order to compute the initial velocity of $\mathcal{I}(0)$, we first need to compute the initial velocity $\dot{\vec{z}}$ of the Euclidean coordinate center

[TABLE]

of $\Sigma_{t}$. We remind the reader that we chose coordinates $\vec{x}$ which do not depend on $t$. Relying on Proposition 7.1, we will now compute the variation of the coordinate center, $\dot{\vec{z}}$, starting from the variation formula

[TABLE]

The idea is to argue as in the proof of Lemma 7.2 and pass from $u$ to $\mathcal{L}u$, and from $\nu^{i}$ to $f_{i}$ in (105). For this we note that since $\mathcal{I}$ is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean, and since $K$ has faster fall-off $K=O_{1}(|\vec{x}\,|^{-2})$, we have $\mathcal{L}u=O(|\vec{x}\,|^{-3})$ by (104). A computation identical to the one in the proof of Lemma 7.2 yields

[TABLE]

Thus, (104) and integration by parts give us

[TABLE]

Recall that $\Sigma$ is approximated by the coordinate sphere $\mathbb{S}_{\sigma}^{2}$, as described in Section 4. In particular, the functions $f_{i}$, $i=1,2,3$, are close to the respective eigenfunctions $f^{\delta}_{i}$ of the Laplacian $-\Delta^{\mathbb{S}_{\sigma}}$, see (33). Thus

[TABLE]

where $\frac{\vec{x}}{\sigma}$ is the unit normal vector to $\mathbb{S}_{\sigma}^{2}\hookrightarrow(\mathbb{R}^{3},\delta)$. Furthermore, a computation shows that

[TABLE]

Since $f^{\delta}_{i}=\tfrac{\sqrt{3}x^{i}}{\sqrt{4\pi}\sigma^{2}}$ and since $g_{ij}=\delta_{ij}+O(\sigma^{-\frac{1}{2}-\varepsilon})$, we finally arrive at

[TABLE]

Passing to the limit when $\sigma\to\infty$ we obtain the result.∎

8.2. Poincaré-covariance and accordance with Special Relativity

As we have seen before, whether or not a given $2$-surface is STCMC is in fact independent of a choice of slice (as well as of a choice of coordinates). In this sense, STCMC-surfaces are covariant in the sense of General Relativity. The role of the initial data set then is to select a unique family of STCMC-surfaces near the asymptotic end of the spacetime, forming its abstract STCMC-center of mass. In this sense, STCMC-foliations and the associated (abstract) center of mass are Poincaré-covariant.

We will now discuss the transformation behavior of the STCMC-coordinate center under the asymptotic Poincaré group of the ambient spacetime — assuming vanishing angular momentum in the boost case. Dealing with angular momentum and treating the boost case more adequately will be left for our future work. Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be an initial data set which is $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean with respect to asymptotic coordinates $\vec{x}$ and has $E\neq 0$.

Euclidean motions.

Consider the coordinates $\vec{y}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}O\vec{x}+\vec{T}$, with $O$ an orthogonal rotation matrix and $\vec{T}\in\mathbb{R}^{3}$ a translation vector. In other words, $\vec{y}$ arises from $\vec{x}$ through a Euclidean motion. Then, for each leaf $\Sigma^{\sigma}$ of the STCMC-foliation constructed in Theorem 6.2, we find that the Euclidean center of $\Sigma^{\sigma}$ with respect to the $\vec{y}$-coordinates is given by

[TABLE]

Thus, the STCMC-coordinate center

[TABLE]

with respect to the coordinates $\vec{y}$ converges if and only the STCMC-coordinate center

[TABLE]

converges with respect to the coordinates $\vec{x}$ converges and if they converge, we find

[TABLE]

as one would expect from Euclidean Geometry, Newtonian Gravity, and from the description of the spacetime position of a point particle in Special Relativity.

Time translations.

The transformation behavior of $\vec{C}_{\text{STCMC}}$ under asymptotic time translation corresponds to its evolution behavior under the Einstein equations. In other words, Theorem 8.1 tells us under the additional assumption $K=O_{1}(|\vec{x}\,|^{-2})$ that

[TABLE]

which corresponds precisely to the instantaneous law of motion of a point particle in Special Relativity.

Boosts.

The last constituent of the asymptotic Poincaré group of the spacetime are of course the asymptotic boosts. In a given (asymptotic region of a) spacetime

[TABLE]

with asymptotic coordinates $x^{\alpha}=(x^{0},x^{i})$ and suitably decaying lapse $N$, shift $X$, and tensor $h$, a boosted initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)\hookrightarrow(\mathbb{R}\times M^{3},\mathfrak{g})$ is any spacelike hypersurface arising as the set $\{y^{0}=0\}$ with respect to a boosted coordinate system $y^{\alpha}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\Lambda^{\alpha}_{\beta}x^{\beta}$, $y^{\alpha}=(y^{0},\vec{y}\,)$, meaning that the matrix $\Lambda$ is a boost. If the lapse $N$, the shift $X$, and the tensor $h$ decay suitably fast in space and time coordinate directions, the boosted initial data set $\{y^{0}=0\}=\mathcal{I}$ is in fact $C^{1}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean with respect to $\vec{y}$. It is thus reasonable to ask how the STCMC-coordinate centers of the initial data sets $\{x^{0}=0\}$ and $\{y^{0}=0\}$ are related (if they converge). The corresponding question was addressed by Szabados [41] for the BÓM-center of mass although from a slightly different perspective. Of course, we expect that the STCMC-coordinate center boosts like the position of a point particle in Special Relativity, namely

[TABLE]

at least in the absence of angular momentum. This can easily be verified for example for a boosted slice (over the canonical slice) in the Schwarzschild spacetime where in fact the centers both coincide with the center of symmetry $\vec{0}$. Similarly, if one first spatially translates the coordinates on the Schwarzschild spacetime and then considers a boosted slice, the transformation law will be as in (107). In both of these examples, the deviation $\vec{Z}$ introduced in Theorem 7.5 in fact vanishes, so that the transformation law (107) already holds for the CMC-BÓM-center of mass. In view of Section 9 below, it is possible to construct examples of boosted slices in the Schwarzschild spacetime by boosting the example discussed below. One can then see that (107) in fact also applies in this case, but only for the STCMC-center of mass, and not for the CMC-BÓM-center of mass. However, the computation is so tedious that we prefer not to show it here as it is not particularly enlightening.

A proof of a generalized version of (107) incorporating the angular momentum will be given elsewhere, see also Remarks 7.8, 8.3.

9. A concrete graphical example in the Schwarzschild spacetime

As briefly sketched in Sections 2, 3 and analyzed in more detail in Section 7, determining the coordinate center of an asymptotic foliation is tricky and depends on choosing suitable coordinates (see also Conjecture 7.9). In [10, Section 6], this was illustrated by explicitly computing the coordinate center of the CMC-foliation of an asymptotically Euclidean “graphical” time-slice in the Schwarzschild spacetime of mass $m\neq 0$. This example, to be described in more detail below, satisfies all assumptions in [26], in particular those of Theorem 4.2, but yet its CMC-coordinate center does not converge. Equivalently, its BÓM-center also does not converge. After a brief introduction to the graphical example discussed in [10], we will compute that the STCMC-coordinate center (96) does in fact converge in this example and moreover converges to the origin $\vec{0}$, i.e. to the center of symmetry of the spherically symmetric spacetime as one would expect.

We consider the Schwarzschild spacetime $(\mathbb{R}\times M^{3},\mathfrak{g})$ of mass $m\neq 0$ in Schwarzschild coordinates, meaning that

[TABLE]

where $d\Omega^{2}$ denotes the canonical metric on $\mathbb{S}^{2}$. We will freely switch between polar coordinates $(r,\eta)$ and the naturally corresponding Cartesian coordinates $\vec{x}$ defined on $M^{3}$. A graphical time-slice in the (automatically vacuum) Schwarzschild spacetime is an initial data set $(M^{3}_{T},g_{T},K_{T},\mu_{T}\equiv 0,J_{T}\equiv 0)$ arising as the graph of a smooth function $T\colon M^{3}\to\mathbb{R}$ “over” the canonical time-slice $\{{t=0\}}$ (in time-direction), meaning that

[TABLE]

while $g_{T}$ is the Riemannian metric induced on $M_{T}\hookrightarrow(\mathbb{R}\times M^{3},\mathfrak{g})$ and $K_{T}$ is the second fundamental form induced by this embedding with respect to the future pointing unit normal.

Computing the CMC-coordinate center of mass (via the BÓM-center of mass).

Clearly, the center of mass of the canonical time-slice $\{{t=0\}}$ of the Schwarzschild spacetime is the coordinate origin, $\vec{C}_{\text{CMC}}=\vec{C}_{\text{B\'{O}M}}=\vec{0}$. We will now compute this vector for graphical time-slices with the asymptotic decay conditions on $T$ chosen such that $(M^{3}_{T},g_{T},K_{T},\mu_{T}\equiv 0,J_{T}\equiv 0)$ is $C^{2}_{1}$-asymptotically Euclidean with respect to the coordinates $\vec{x}$. To most easily comply with the asymptotic decay conditions specified in Section 2, we will assume that $T=O_{k}(r^{0})$ as $r\to\infty$, with $k\geq 3$.

Now let $\vec{y}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\vec{x}|_{M_{T}}$ denote the induced coordinates on $M_{T}$. As computed in [10, Section 6], the metric $g_{T}$ and second fundamental form $K_{T}$666The corresponding formula for the second fundamental form in [10] has a typo which we corrected here. We thank Axel Fehrenbach for pointing this out to us. are given by

[TABLE]

in the coordinates $\vec{y}$. A straightforward computation shows that the graphical initial data set $(M^{3}_{T},g_{T},K_{T},\mu_{T}\equiv 0,J_{T}\equiv 0)$ is indeed $C^{2}_{1}$-asymptotically Euclidean and in fact has $E=m\neq 0$.

When evaluating the BÓM-center of mass surface integral on a finite coordinate sphere with respect to the $\vec{y}$-coordinates in $M_{T}$, using $s\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}|\vec{y}\,|$ and $\vec{\eta}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}\tfrac{\vec{y}}{s}$, we find

[TABLE]

As in [10, Section 6]777In fact, we are using Schwarzschild coordinates, here, while in [10], isotropic coordinates are used. This allows us to treat the case of $m<0$ as well and does not affect the asymptotic computations. , we pick a fixed vector $\vec{0}\neq\vec{u}\in\mathbb{R}^{3}$ and set

[TABLE]

We point out that this choice of $T$ ensures that $(M^{3},g_{T},\mu_{T}\equiv 0)$ satisfies the Riemannian $C^{2}_{\frac{1}{2}+\varepsilon}$-Regge–Teitelboim conditions so that [34, Cor. 4.2] or [36, Theorem 6.3] apply and ensure that $\vec{C}_{\text{CMC}}=\vec{C}_{\text{B\'{O}M}}$ or that both diverge. One directly computes from the above expression for $\vec{C}_{\text{B\'{O}M}}\left(\mathbb{S}^{2}_{s}(\vec{0}\,)\right)$ that

[TABLE]

which diverges as $s\to\infty$. Hence, the BÓM- and thus also the CMC-coordinate center diverge in this example.

Computing the STCMC-coordinate center of mass (via Formula (96)).

In order to check whether the STCMC-coordinate center of the $C^{2}_{1}$-asymptotically Euclidean initial data set $(M^{3}_{T},g_{T},K_{T},\mu_{T}\equiv 0,J_{T}\equiv 0)$ converges, one needs to compute the STCMC-leaves $\Sigma_{\sigma}$ and the coordinate averages $\vec{z}\,(\Sigma_{\sigma})$ and check whether they converge as $\sigma\to 0$. However, the proof of Theorem 7.5 asserts that $\vec{C}_{\text{STCMC}}$ converges if and only if the coordinate expression given in (96) converges, or in other words if and only if

[TABLE]

converges as $s\to\infty$, where we recall that, using $E=m$ and $(\pi_{T})_{kl}=-(K_{T})_{ij}+\operatorname{tr}_{g_{T}}K_{T}(g_{T})_{ij}$, we know that

[TABLE]

Our (diverging) spacetime correction term $\vec{Z}$ thus precisely compensates for the divergence occurring in $\vec{C}_{\text{CMC}}=\vec{C}_{\text{B\'{O}M}}$. Hence the STCMC-coordinate center of the considered graphical slice converges to $\vec{0}$ as desired.

Appendix A Round surfaces in asymptotically Euclidean manifolds

In this appendix we collect some standard results about closed surfaces in a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set $\mathcal{I}=(M^{3},g,K,\mu,J)$ (see Definition 2.1 and Remark 2.3) that are repeatedly used in this paper.

Lemma A.1**.**

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set with asymptotic coordinate chart $\vec{x}\colon M^{3}\setminus\mathcal{B}\to\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ and let $\Sigma\hookrightarrow M^{3}\setminus\mathcal{B}$ be a closed, oriented $2$-surface. Using the chart $\vec{x}$, we may also view $\Sigma$ as a surface in $\mathbb{R}^{3}\setminus\overline{B_{R}(0)}$ equipped with the Euclidean metric $\delta$, with the induced metric denoted by $\delta^{\Sigma}$. Then there exist positive constants $c$ and $C$ depending only on $\varepsilon$ and $C_{\mathcal{I}}$ such that the following holds, provided that the Euclidean distance to the coordinate origin $|\vec{x}\,|$ on $\Sigma$ satisfies $|\vec{x}\,|\geq c$:

• •

The normals $\nu$ and $\nu^{\delta}$ of $\Sigma$ in the metrics $g$ and $\delta$ satisfy

**

[TABLE]

• •

The volume elements $d\mu$ and $d\mu^{\delta}$ satisfy

**

[TABLE]

• •

The respective second fundamental forms $A$ and $A^{\delta}$ satisfy

**

[TABLE]

and the respective mean curvatures $H$ and $H^{\delta}$ are related via

[TABLE]

Further, if $\|H\|_{L^{2}(\Sigma)}$ is a priori bounded, then the respective trace-free parts of the second fundamental forms satisfy

[TABLE]

where $C$ also depends on the bound on $\|H\|_{L^{2}(\Sigma)}$.

Proof.

See [32, Section 2.4] or [31, Section 1.5], where similar estimates are proven. ∎

The following result is a Sobolev Embedding Theorem which holds for a very general class of $2$-surfaces. In Section 4 this result is applied to $\Sigma$ being a large coordinate sphere $\mathbb{S}^{2}_{r}(\vec{z})\hookrightarrow M^{3}\setminus\mathcal{B}$ in the asymptotic end of an asymptotically Euclidean initial data set. Note that in this case one can without loss of generality replace the area radius $\sqrt{\nicefrac{{|\Sigma|}}{{4\pi}}}$ in the formulation of Lemma A.2 by the coordinate sphere’s radius $r$ as these two radii are uniformly equivalent. In the subsequent sections, this result is applied to $\Sigma\hookrightarrow M^{3}\setminus\mathcal{B}$ being an asymptotically centered closed $2$-surface with constant spacetime mean curvature. Here we are using the fact that (108) is available for large asymptotically centered surfaces in $M^{3}\setminus\mathcal{B}$ in the form of [26, Proposition 5.4], which applies to asymptotically Euclidean initial data sets with general asymptotics as described in Section 2.

Lemma A.2**.**

Let $(\Sigma,g^{\Sigma})$ be a closed, oriented 2-surface with area radius $r=\sqrt{\nicefrac{{|\Sigma|}}{{4\pi}}}$. If there is a constant $C_{S}$ such that for any Lipschitz continuous function $f$ on $\Sigma$, the so-called first Sobolev Inequality

[TABLE]

holds, then we also have the Sobolev Inequality888Note that the constant in (109) is not necessarily optimal.

[TABLE]

for any $f\in W^{2,2}(\Sigma)$.

Proof.

With the first Sobolev Inequality (108) at hand, the following Sobolev Inequalities can be derived for $p>2$:

[TABLE]

see e.g. [35, Proposition II.1.3] for details. As a consequence, for any $f\in W^{2,2}(\Sigma)$ we have

[TABLE]

for any $p>2$. In particular, for $p=4$ we have (109). ∎

The following result is well-known, see e.g. [40, Corollary 2.10] (adapted from [14, Chapter 2]).

Lemma A.3**.**

Let $(\Sigma,g^{\Sigma})$ be a $2$-surface of spherical topology with Gaussian curvature $K$ satisfying

[TABLE]

where $r=\sqrt{\nicefrac{{|\Sigma|}}{{4\pi}}}$ is the area radius of $\Sigma$. Then there is a universal constant $C$ such that for any $f\in W^{2,2}(\Sigma)$ we have

[TABLE]

where $\overline{f}$ denotes the mean value of $f$ on $\Sigma$.

Appendix B The STCMC-condition in normal geodesic coordinates

Let $\mathcal{I}=(M^{3},g,K,\mu,J)$ be a $C^{2}_{\nicefrac{{1}}{{2}}+\varepsilon}$-asymptotically Euclidean initial data set and let $\Sigma\hookrightarrow M^{3}$ be a closed, oriented $2$-surface. Let $(u^{\alpha})$ be coordinates on $\Sigma$ and let $\partial_{\alpha}$ denote the respective tangent vectors to $\Sigma$, for $\alpha\in\{{1,2\}}$. Here and in the rest of this appendix, we use the convention that Greek indices $\alpha,\beta,\gamma\in\{1,2\}$ refer to coordinate vector fields tangential to $\Sigma$.

In a neighborhood of $\Sigma$, the normal geodesic coordinates $y\colon\Sigma\times(-\xi,\xi)\to M^{3}$ are defined for some $\xi>0$, see Section 5.1 for details. In this neighborhood we may write $g=dt^{2}+g_{t}$ where $g_{t}$ is the induced metric on $\Sigma_{t}\mathrel{\raise 0.40903pt\hbox{\rm:}\mkern-5.2mu=}y(\Sigma,t)$.

Consider a $2$-surface $S$ given as the graph of a function $f$ with $|f|<\xi$ over $\Sigma$, i.e.

[TABLE]

Since the vector $-\partial_{t}+\nabla^{g_{t}}f$ is normal to $S$ at the point $(q,t)=(q,f(q))$, the vectors $\partial_{\alpha}+(\partial_{\alpha}f)\partial_{t}$ are tangent to $S$ at this point. As a consequence, the induced metric on $S$ has components given by

[TABLE]

with the components of the inverse given by

[TABLE]

Here and in what follows, $2$-dimensional indices $\alpha$, $\beta$, $\dots$ are raised with respect to $g_{t}$, and all quantities are computed at the point $(q,t)=(q,f(q))$, unless stated otherwise. A straightforward computation then shows that the mean curvature of $S=\operatorname{graph}f$ is given by

[TABLE]

where $A_{t}$ is the second fundamental form of $\Sigma_{t}$. We also have

[TABLE]

Proposition B.1**.**

If $S=\operatorname{graph}f$ is a surface of constant spacetime mean curvature $\mathcal{H}(S)\equiv\nicefrac{{2}}{{\sigma}}$ then $f$ satisfies the equation

[TABLE]

where, with $P$ given by (111), we use the shorthands

[TABLE]