  Our work proves a rigidity theorem for initial data sets associated with
convex polytopes, subject to the dominant energy condition. The theorem is
established by utilizing an approach that involves approximating the polytope
of interest with smooth convex domains and solving a boundary value problem for
Dirac operators on these domains.
