Rigidity results for initial data sets satisfying the dominant energy condition

Basic Information

• Paper ID: 2304.04145
• Title: Rigidity results for initial data sets satisfying the dominant energy condition
• Authors: Christian Bär, Simon Brendle, Tsz-Kiu Aaron Chow, Bernhard Hanke
• Classification: math.DG (Differential Geometry), math-ph (Mathematical Physics), math.AP (Analysis and PDEs), math.MP (Mathematical Physics)
• Publication Date: arXiv v4, November 10, 2025
• Paper Link: https://arxiv.org/abs/2304.04145

Abstract

This paper establishes rigidity theorems for initial data sets satisfying the dominant energy condition. The objects of study include compact smooth spin manifolds with boundary and compact convex polytopes. For smooth boundary manifolds, the proof is based on boundary value problems for the Dirac operator; for convex polytopes, an approximation method using smooth boundary manifolds is employed.

Research Background and Motivation

Research Problem

This paper investigates rigidity problems for initial data sets in general relativity. An initial data set is a triple $(M,g,q)$, where $M$ is a manifold, $g$ is a Riemannian metric, and $q$ is a symmetric $(0,2)$-tensor. The central question is: under the dominant energy condition $\mu \geq |J|$, which initial data sets can be isometrically embedded into Minkowski spacetime $\mathbb{R}^{n,1}$?

Importance of the Problem

1. Physical Significance: The dominant energy condition is a fundamental physical assumption in general relativity, ensuring that the energy density $\mu$ dominates the momentum density $|J|$
2. Geometric Significance: Rigidity results provide complete classifications of geometric objects, revealing deep connections between curvature conditions and global geometric structures
3. Theoretical Value: Connects spin geometry, index theory, and general relativity, generalizing rigidity statements of the positive mass theorem

Limitations of Existing Approaches

• Schoen-Yau (1981) and Witten (1982) established the spacetime positive mass theorem, but primarily for non-compact asymptotically flat manifolds
• Eichmair et al. (2016) and Huang-Lee (2020) proved related results under dimensional restrictions
• Rigidity results for compact manifolds with boundary and polytopes remain incomplete
• Scalar curvature theory on low-regularity spaces (such as polytopes) is still under development

Research Motivation

The authors shift the focus from non-compact asymptotically flat manifolds to compact manifolds with boundary, replacing asymptotic flatness with appropriate boundary conditions. This approach not only generalizes classical results but also provides new tools for studying low-regularity spaces such as polytopes.

Core Contributions

1. Theorem I (Smooth Boundary Case): For a spin manifold $(M,g,q)$ satisfying the dominant energy condition, the boundary mean curvature condition $H \geq \|dN\|_{tr} + |tr(q)-q(\nu,\nu)|$, and with the Gauss map $N:\Sigma \to S^{n-1}$ having positive degree, the manifold must be isometric to the manifold in Example 1.2 (i.e., embeddable in Minkowski spacetime as a graph hypersurface)
2. Corollary 1.6 (Case $q=0$): When $q=0$, Theorem I simplifies to: a compact Riemannian spin manifold satisfying $R \geq 0$ and $H \geq \|dN\|_{tr}$ can be isometrically embedded in $\mathbb{R}^n$ as a convex domain
3. Theorem II (Polytope Case): For a compact convex polytope $\Omega \subset \mathbb{R}^n$, under the dominant energy condition, boundary face mean convexity $H \geq |tr(q)-q(\nu,\nu)|$, and the matching angle hypothesis, the initial data set $(\Omega,g,q)$ can be isometrically immersed into Minkowski spacetime, with each boundary face being totally geodesic as a codimension-2 submanifold
4. Technical Innovations:
   • Development of boundary value problem theory for the Dirac operator (Section 2)
   • Establishment of holomorphic index theorem applications in odd and even dimensions (Proposition 2.4)
   • Construction of modified Weitzenböck formula (Proposition 3.3)
   • Development of refined approximation techniques for smooth manifolds approximating polytopes (Section 5)

Detailed Methodology

Task Definition

Input: Initial data set $(M,g,q)$, where

• $M$ is an $n$-dimensional compact connected spin manifold with non-empty boundary $\partial M = \Sigma$
• $g$ is a Riemannian metric
• $q$ is a symmetric $(0,2)$-tensor
• Gauss map $N:\Sigma \to S^{n-1}$

Conditions:

1. Dominant energy condition: $\mu = \frac{1}{2}R + \frac{1}{2}tr(q)^2 - \frac{1}{2}|q|^2 \geq |J| = |div(q) - \nabla tr(q)|$
2. Boundary condition: $H \geq \|dN\|_{tr} + |tr(q)-q(\nu,\nu)|$
3. Topological condition: $deg(N) > 0$

Output: Prove that $(M,g,q)$ is isometric to the standard model in Example 1.2

Model Architecture

Step One: Construction of Modified Spin Connection (Section 3)

Define modified connections $\nabla^{S,q}$ and $\nabla^{E,q}$: $\nabla^{S,q}_X := \nabla^S_X + \frac{i}{2}\gamma(q(X))$$\nabla^{E,q}_X s := \nabla^E_X s + \frac{i}{2}(Qs)(X)$ where $Q$ is the homomorphism field $(Qs)(X) := (\gamma(q(X))\otimes id)s$

Step Two: Modified Weitzenböck Formula (Proposition 3.3)

Establish the key identity: $(D + \frac{i}{2}tr(q))(D - \frac{i}{2}tr(q)) = (\nabla^{E,q})^*\nabla^{E,q} + \frac{1}{2}\mu + \frac{i}{2}\gamma(J)\otimes id$

This formula connects the Dirac operator, modified connection, and dominant energy condition.

Step Three: Boundary Value Problems and Index Calculation (Section 2)

Define the boundary operator: $\chi := -\gamma(\nu) \otimes \gamma^*_0(N)$ satisfying $\chi^2 = id$ and $\chi^* = \chi$, decomposing $E|_\Sigma$ into $\pm 1$ eigensubbundles $F_\pm$.

Key Result (Proposition 2.4): In the odd-dimensional case ($n \geq 3$ odd), the operator $D: H^1(M,E,F_+) \to L^2(M,E)$ is a Fredholm operator with index equal to $deg(N)$.

Step Four: Existence Proof

• Odd-Dimensional Case (Proposition 3.7): Using positive index to guarantee existence of non-trivial section $s \in C^\infty(M,E)$ satisfying $\int_M |\nabla^{E,q}s|^2 \leq \frac{1}{2}\int_M (|J|-\mu)|s|^2 + \frac{1}{2}\int_\Sigma (|tr(q)-q(\nu,\nu)| + \|dN\|_{tr} - H)|s|^2$ with $\chi s = s$
• Even-Dimensional Case (Corollary 3.11): Obtained through the $M \times S^1$ technique and compactness arguments

Step Five: Rigidity Analysis (Section 4)

Prove that the existing $\nabla^{E,q}$-parallel section $s$ is invertible at each point (Lemma 4.5, 4.12), and using boundary condition analysis:

1. Eigenvalues of the second fundamental form $h$ equal singular values of $dN$ (Proposition 4.1(b))
2. $q_{tan} = 0$ (Proposition 4.1(c))

Step Six: Isometric Embedding Construction (Section 4.3)

Construct bundle $T = TM \oplus \mathbb{R}$ equipped with signature $(n,1)$ metric, using existence of $\nabla^{E,q}$-parallel sections to construct $\nabla^{T,q}$-parallel frame, yielding local isometric immersion $F: M \to \mathbb{R}^{n,1}$.

Technical Innovations

1. Application of Holomorphic Index Theorem: First application of the Bär-Brendle-Hanke-Wang holomorphic index theorem to initial data set problems, computing index through boundary operator $A = D_\Sigma + \frac{1}{2}\chi \circ (\sum \gamma(e_j)\otimes \gamma^*_0(dN(e_j)))$
2. Modified Connection Technique: Introduction of $q$-dependent modified connection $\nabla^{E,q}$, allowing the Weitzenböck formula to directly reflect the dominant energy condition
3. Even-Dimensional Case Treatment: Through construction of $M \times S^1$ and carefully chosen map $\phi: S^{n-1} \times S^1 \to S^n$ (Lemma 3.8), converting even-dimensional problems to odd-dimensional ones
4. Polytope Approximation Theory: For convex polytope $\Omega = \bigcap_{i \in I}\{u_i \leq 0\}$, construct smooth approximations $\Omega_\lambda = \{\sum_{i \in I} e^{\lambda u_i} \leq 1\}$ establishing refined $L^\sigma$ estimates ($\sigma \in [1,\frac{3}{2})$), proving $\sup_{p,r} \left(r^{\sigma+1-n}\int_{\Sigma_\lambda \cap B_r(p)} \max\{-W_\lambda, 0\}^\sigma\right)^{1/\sigma} \to 0$ where $W_\lambda = V_\lambda - |tr(q)-q(\nu,\nu)|$

Experimental Setup

This is a pure mathematics theoretical paper with no numerical experiments. All results are obtained through rigorous mathematical proofs.

Theoretical Verification Framework

The paper verifies theoretical correctness through the following levels:

1. Special Case Verification: Examples 1.1 and 1.2 provide explicit examples satisfying all conditions
2. Dimensional Induction: The $n=1$ case can be directly verified through ODE methods (Remark 1.5)
3. Limit Process: Polytope case verified through limit process of smooth approximations (Section 6)

Experimental Results

Complete Statement of Main Theorems

Theorem I: Let $(M,g,q)$ be an initial data set with $M$ an $n \geq 2$-dimensional compact connected spin manifold with non-empty boundary $\partial M = \Sigma$, and $N:\Sigma \to S^{n-1}$ a smooth map. Assume:

• $(M,g,q)$ satisfies the dominant energy condition
• $H \geq \|dN\|_{tr} + |tr(q)-q(\nu,\nu)|$ holds along $\Sigma$
• $N:\Sigma \to S^{n-1}$ has positive degree

Then $(M,g)$ is isometric to the manifold in Example 1.2, $q$ is the second fundamental form of $M$ in $\mathbb{R}^{n,1}$, and $N$ agrees with the outward normal $\partial K \to S^{n-1}$ up to $O(n)$ elements.

Key Corollaries

Corollary 1.6: When $q=0$, conditions simplify to $R \geq 0$ and $H \geq \|dN\|_{tr}$, with conclusion that $M$ can be isometrically embedded in $\mathbb{R}^n$ as a compact convex domain with smooth boundary.

This generalizes Gromov's question on scalar curvature comparison to the case of manifolds with boundary.

Results for Polytopes

Theorem II: Let $\Omega \subset \mathbb{R}^n$ be a compact convex polytope, $g$ a Riemannian metric defined on an open set containing $\Omega$, and $q$ a symmetric $(0,2)$-tensor. Assume:

• $(\Omega,g,q)$ satisfies the dominant energy condition
• $H \geq |tr(q)-q(\nu,\nu)|$ holds on boundary faces
• The Matching Angle Hypothesis is satisfied

Then $(\Omega,g,q)$ can be isometrically immersed into Minkowski spacetime $\mathbb{R}^{n,1}$, $q$ is the second fundamental form, and each boundary face is totally geodesic as a codimension-2 submanifold.

Technical Results

1. Index Formula (Proposition 2.4): $ind(D: H^1(M,E,F_+) \to L^2(M,E)) = deg(N)$
2. Energy Estimate (Proposition 3.7): Under the dominant energy condition, $\int_M |\nabla^{E,q}s|^2 \leq \frac{1}{2}\int_M (|J|-\mu)|s|^2 + \frac{1}{2}\int_\Sigma (|tr(q)-q(\nu,\nu)| + \|dN\|_{tr} - H)|s|^2$ The right-hand side is non-positive when conditions are satisfied, ensuring $s$ is $\nabla^{E,q}$-parallel.
3. Approximation Estimate (Corollary 5.9): For $\sigma \in [1,\frac{3}{2})$, $\sup_{p \in \mathbb{R}^n} \sup_{r \leq 1} \left(r^{\sigma+1-n}\int_{\Sigma_\lambda \cap B_r(p)} \max\{-W_\lambda,0\}^\sigma\right)^{1/\sigma} \to 0 \text{ as } \lambda \to \infty$

Related Work

History of Positive Mass Theorem

1. Schoen-Yau (1981): Proved 3-dimensional spacetime positive energy theorem using minimal surface techniques
2. Witten (1982): Proved positive mass theorem in all dimensions using spinor methods
3. Parker-Taubes (1982): Provided rigorous mathematical form of Witten's argument
4. Chruściel-Maerten (2006), Beig-Chruściel (1996): Proved rigidity statements

Research on Manifolds with Boundary

1. Eichmair-Huang-Lee-Schoen (2016): Proved spacetime positive mass theorem in dimensions up to 7
2. Huang-Lee (2020): Proved rigidity statement of this theorem
3. Eichmair-Galloway-Mendes (2021): First obtained rigidity results for boundary case
4. Glöckle (2023): Studied initial data rigidity through Dirac-Witten operator

Scalar Curvature on Polytopes

1. Gromov (2014): Pioneering work on scalar curvature bounds for polytopes
2. Li (2020, 2024): Proved comparison theorem for 3-dimensional polytopes and dihedral angle rigidity for $n$-prisms
3. Brendle (2024): Proved scalar curvature rigidity for convex polytopes
4. Brendle-Wang (2023): Studied Gromov rigidity theorem for acute-angled polytopes
5. Tsang (2021): Studied dihedral angle rigidity for cubic initial data sets

Dirac Operator Boundary Value Problems

1. Hörmander (1994): Early research
2. Bär-Ballmann (2012, 2016): Systematic development of theory
3. Bär-Brendle-Hanke-Wang (2024): Holomorphic index theorem

Advantages of this paper:

• Unified treatment of smooth manifolds and polytope cases
• First application of holomorphic index theorem to initial data sets
• Handles general $(g,q)$ rather than special cases like $q=0$ or $q=g$
• Provides complete rigidity characterization rather than only comparison results

Conclusions and Discussion

Main Conclusions

1. Complete Characterization: Under the dominant energy condition and appropriate boundary conditions, initial data sets on compact spin manifolds are completely characterized by graph hypersurfaces embeddable in Minkowski spacetime
2. Topological Necessity: The positive degree condition $deg(N) > 0$ is essential, guaranteeing existence of critical sections through index theory
3. Geometric Rigidity: The boundary condition $H \geq \|dN\|_{tr} + |tr(q)-q(\nu,\nu)|$ achieves equality if and only if the manifold is the standard model in Example 1.2
4. Polytope Generalization: The matching angle hypothesis allows generalization to non-smooth boundaries, providing tools for low-regularity geometry

Limitations

1. Spin Assumption: Theorem I requires $M$ to be a spin manifold, excluding certain topological types (such as $\mathbb{RP}^{2n}$). This is an essential limitation of the method, as the Dirac operator depends on spin structure
2. Dimensional Restrictions: While theoretically applicable to all $n \geq 2$, the $n=1$ case requires separate treatment (Remark 1.5), and proof techniques differ significantly between odd and even dimensions
3. Boundary Regularity: Theorem II requires the metric $g$ and tensor $q$ to be defined on an open set containing $\Omega$, which may be overly restrictive in some applications
4. Matching Angle Hypothesis: While physically natural, this hypothesis may be difficult to verify in pure geometric settings
5. Non-Compact Case: The paper focuses on compact cases, complementing but not overlapping with the non-compact asymptotically flat setting of classical positive mass theorem

Future Directions

1. Non-Spin Case: Explore use of other elliptic operators (such as Dirac-Witten operator) for non-spin manifolds
2. More General Boundaries: Study cases with corners or lower regularity boundaries
3. Time Evolution: Generalize results to evolution problems in spacetime
4. Numerical Methods: Develop numerical algorithms to verify or construct initial data sets satisfying the conditions
5. Physical Applications: Apply these rigidity results to constrain initial data in numerical relativity

In-Depth Evaluation

Strengths

1. Methodological Innovation:
   • Creative application of holomorphic index theorem to initial data sets, marking its first use in general relativity
   • Elegant construction of modified connection $\nabla^{E,q}$ encoding information of second fundamental form $q$
   • Treatment of even-dimensional case through $M \times S^1$ technique demonstrates profound geometric insight
2. Technical Depth:
   • Proof synthesizes spin geometry, index theory, PDEs, and convex geometry
   • Modified Weitzenböck formula (Proposition 3.3) is the technical core, converting physical conditions to analytical ones
   • $L^\sigma$ estimates for polytope approximation (Section 5) are technically refined, controlling $\lambda$-dependence
3. Result Completeness:
   • Proves not only existence but complete isometric classification
   • Generalization from smooth case to polytopes is logically clear
   • Provides multiple concrete examples (Examples 1.1, 1.2) and boundary case discussions (Remarks 1.3-1.5)
4. Clarity of Exposition:
   • Well-organized structure: introduction → boundary value problems → spin geometry → proof → approximation → polytopes
   • Key steps supported by detailed lemmas
   • Consistent notation system with clear definitions

Weaknesses

1. Accessibility Challenges:
   • Requires deep background in spin geometry and index theory for non-specialists
   • Some technical details (such as construction of map $\phi$ in Lemma 3.8) lack geometric intuition
   • Separate treatment of odd and even dimensions increases reading burden
2. Insufficient Examples:
   • Beyond Example 1.2 of graph hypersurfaces, lacks other explicit examples
   • No counterexamples when boundary conditions fail
   • Polytope case lacks concrete computational examples
3. Technical Assumptions:
   • Geometric meaning of matching angle hypothesis could be explained more thoroughly
   • Necessity of spin assumption deserves more discussion
   • Geometric/physical significance of positive degree condition $deg(N) > 0$ could be explored more deeply
4. Generalization Potential:
   • Method highly dependent on Dirac operator, difficult to directly extend to non-spin case
   • Polytope proof relies on smooth approximation, applicability to more general singular spaces unclear

Impact

1. Contribution to Field:
   • Provides new analytical tools for initial data set theory (modified connections and boundary value problems)
   • Advances scalar curvature theory on low-regularity spaces
   • Connects general relativity and convex geometry
2. Practical Value:
   • Provides theoretical foundation for numerical relativity, useful for verifying physical reasonableness of initial data
   • Polytope results may apply to discretization methods
3. Reproducibility:
   • Proofs are completely self-contained with clear citations of external results
   • Technical approach is clear, available for subsequent research
   • Appendices A and B provide independently valuable auxiliary results
4. Open Problems:
   • Stimulates subsequent research on non-spin manifolds and lower regularity boundaries
   • Provides new perspective on rigidity theory of spacetime positive mass theorem

Applicable Scenarios

1. General Relativity:
   • Verify whether initial data in numerical simulations correspond to submanifolds of Minkowski spacetime
   • Study geometric constraints on gravitational wave initial data
2. Differential Geometry:
   • Serves as paradigm for studying rigidity problems on manifolds with boundary
   • Generalizable to other curvature conditions (such as Ricci curvature)
3. Convex Geometry:
   • Study relationship between intrinsic geometry and extrinsic embedding of convex polytopes
   • Generalize to more general convex bodies
4. Numerical Analysis:
   • Provide theoretical guarantees for finite element methods
   • Design discretization schemes preserving geometric properties

References

The paper cites 32 key references, primarily including:

1. Classical Positive Mass Theorem Literature:
   • 29 Schoen-Yau (1981): Minimal surface method
   • 32 Witten (1982): Spinor method
   • 27 Parker-Taubes (1982): Rigorous formulation of Witten's argument
2. Boundary Value Problem Theory:
   • 2,3 Bär-Ballmann (2012, 2016): Dirac operator boundary value problems
   • 4 Bär-Brendle-Hanke-Wang (2024): Holomorphic index theorem
3. Scalar Curvature on Polytopes:
   • 17,18,19 Gromov (2014, 2023, 2024): Pioneering work
   • 7 Brendle (2024): Convex polytope rigidity
   • 23,24,25 Li (2020, 2022, 2024): Polytope comparison theorems
4. Initial Data Sets:
   • 13 Eichmair-Huang-Lee-Schoen (2016)
   • 21 Huang-Lee (2020)
   • 12 Eichmair-Galloway-Mendes (2021)

Overall Assessment: This is a high-quality mathematics paper deeply combining spin geometry, index theory, and general relativity, solving important problems in initial data set theory. The proof techniques are sophisticated, and results have both geometric and physical significance. While technically demanding, it provides powerful tools and clear directions for subsequent research. The paper's main value lies in methodological innovation (application of holomorphic index theorem, modified connection techniques) and result completeness (unified treatment from smooth manifolds to polytopes).