This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

Rigidity, tensegrity and reconstruction of polytopes under metric constraints

Martin Winter Mathematics Institute, University of Warwick, Coventry CV4 7AL, United Kingdom martin.h.winter@warwick.ac.uk
(Date: October 1, 2025)
Abstract.

We conjecture that a convex polytope is uniquely determined up to isometry by its edge-graph, edge lengths and the collection of distances of its vertices to some arbitrary interior point, across all dimensions and all combinatorial types. We conjecture even stronger that for two polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e} with the same edge-graph it is not possible that QQ has longer edges than PP while also having smaller vertex-point distances.

We develop techniques to attack these questions and we verify them in three relevant special cases: PP and QQ are centrally symmetric, QQ is a slight perturbation of PP, and PP and QQ are combinatorially equivalent. In the first two cases the statements stay true if we replace QQ by some graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} of the edge-graph GPG_{P} of PP, which can be interpreted as local resp. universal rigidity of certain tensegrity frameworks. We also establish that a polytope is uniquely determined up to affine equivalence by its edge-graph, edge lengths and the Wachspress coordinates of an arbitrary interior point.

We close with a broad overview of related and subsequent questions.

Key words and phrases:
convex polytopes, reconstruction from the edge-graph, Wachspress coordinates, rigidity and tensegrity of frameworks
2010 Mathematics Subject Classification:
51M20, 52B11, 52C25

1. Introduction

In how far can a convex polytope be reconstructed from partial combinatorial and geometric data, such as its edge-graph, edge lengths, dihedral angles, etc., optionally up to combinatorial type, affine equivalence, or even isometry? Questions of this nature have a long history and are intimately linked to the various notions of rigidity.

In this article we address the reconstruction from the edge-graph and some “graph-compatible” distance data, such as edge lengths. It is well-understood that the edge-graph alone carries very little information about the polytope’s full combinatorics, and trying to fix this by supplementing additional metric data reveals two opposing effects at play.

First and foremost, we need to reconstruct the full combinatorics. As a general rule of thumb, reconstruction from the edge-graph appears more tractable for polytope that have relatively few edges (such as simple polytopes as proven by Blind & Mani [2] and later by Kalai [19])111Though “few edges” is not the best way to capture this in general, see [8] or [17].. At the same time however such polytopes often have too few edges to encode sufficient metric data for reconstructing the geometry. This is most evident for polygons, but happens non-trivially in higher dimensions and with non-simple polytopes as well (see Figure˜1).

Refer to caption
Figure 1. Non-isometric realizations with the same edge lenghts.

In contrast, simplicial polytopes have many edges and it follows from Cauchy’s rigidity theorem that such are determined up to isometry from their edge lengths; if we assume knowledge of the full combinatorics. For simplicial polytopes however, the edge-graph alone is usually not enough to reconstruct the combinatorics in the first place (as evidenced by the abundance of neighborly polytopes).

This leads to the following question: how much and what kind of data do we need to supplement to the edge-graph to permit

  1. ()

    unique reconstruction of the combinatorics, also for polytopes with many edges (such as simplicial polytopes), and at the same time,

  2. ()

    unique reconstruction of the geometry, also for polytopes with few edges (such as simple polytopes).

Also, ideally the supplemented data fits into the structural framework provided by the edge-graph, that is, contains on the order of #edges+#vertices\#\text{edges}+\#\text{vertices} datums.

We propose the following: besides the edge-graph and edge lengths, we also fix a point in the interior of the polytope PP, and we record its distance to each vertex of PP (cf. Figure˜2). We believe that this is sufficient data to reconstruct the polytope up to isometry across all dimensions and all combinatorial types.

Refer to caption
Figure 2. A “pointed polytope”, i.e., a polytope PdP\subset\mathbb{R}^{d} with a point xint(P)x\in\operatorname{int}(P). In addition to the edge lengths we also record the lengths of the gray bars – the “vertex-point distances”.

Here and in the following we can assume that the polytopes are suitably translated so that the chosen point is the origin 0int(P)0\in\operatorname{int}(P).

Conjecture 1.1.

Given polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e} with the origin in their respective interiors, and so that PP and QQ have isomorphic edge-graphs, corresponding edges are of the same length, and corresponding vertices have the same distance to the origin. Then PQP\simeq Q (i.e., PP and QQ are isometric via an orthogonal transformation).

Requiring the origin to lie in the interior is necessary to prevent counterexamples such as the one shown in Figure˜3. This conjecture vastly generalizes several known reconstruction results, such as for matroid base polytopes or highly symmetric polytopes (see Section˜5.1).

Refer to caption
Figure 3. Two non-isometric realizations of a pentagon with the same edge lengths and vertex-point distances. This is possible because the point is not in the interior.

We also make the following stronger conjecture:

Conjecture 1.2.

Given two polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e} with isomorphic edge-graphs, and so that

  1. ()

    0int(Q)0\in\operatorname{int}(Q),

  2. ()

    edges in QQ are at most as long as their counterparts in PP, and

  3. ()

    vertex-origin distances in QQ are at least as large as their counterparts in PP,

then PQP\simeq Q (PP and QQ are isometric via an orthogonal transformations).

Intuitively, ˜1.2 states that a polytope cannot become larger (or “more expanded” as measured in vertex-origin distances) while its edges are getting shorter. It is clear that ˜1.1 is a consequence of ˜1.2, and we shall call the former the “unique reconstruction version” of the latter. Here, the necessity of the precondition 0int(Q)0\in\operatorname{int}(Q) can be seen even quicker: vertex-origin distances can be increased arbitrarily by translating the polytope just far enough away from the origin (see also Figure˜4).

Refer to caption
Figure 4. If xint(Q)x\not\in\operatorname{int}(Q), then it is possible for QQ to have shorter edges than PP, while simultaneously also all vertices farther away from xx.

In this article we develop techniques that we feel confident point us the right way towards a resolution of the conjectures. We then verify the conjectures in the following three relevant special cases:

The last special case clarifies, in particular, the case of 3-dimensional polytopes. Also, our eventual formulations of the first two special cases will in fact be more general, replacing QQ by some embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} of the edge-graph GPG_{P}, where qq is no longer assumed to be the skeleton of any polytope. These results can then also be interpreted as claiming rigidity, local or universal, of certain bar-joint or tensegrity frameworks.

1.1. Notation and terminology

Throughout the article, all polytopes are convex and bounded, in particular, can be written as the convex hull of their vertices:

P=conv{p1,,pn}:={iαipiαΔn},P=\operatorname{conv}\{p_{1},...,p_{n}\}:=\Big{\{}\smash{\sum_{i}\alpha_{i}p_{i}\mid\alpha\in\Delta_{n}}\Big{\}},

where Δn:={x0nx1++xn=1}\Delta_{n}:=\{x\in\mathbb{R}^{n}_{\geq 0}\mid x_{1}+\cdots+x_{n}=1\} denotes the set of convex coefficients.

If not stated otherwise, PdP\subset\mathbb{R}^{d} will denote a polytope in dd-dimensional space for d2d\geq 2, though its affine hull aff(P)\operatorname{aff}(P) might be a proper subspace of d\mathbb{R}^{d}. If dimaff(P)=d\dim\operatorname{aff}(P)=d we say that PP is full-dimensional. Our polytopes are often pointed, that is, they come with a special point xint(P)x\in\operatorname{int}(P) (sometimes also on P\partial P or outside); but we usually translate PP so that xx is the origin. So, instead of distances from the vertices to xx, we just speak of vertex-origin distances.

By (P)\mathcal{F}(P) we denote the face-lattice of PP, and by δ(P)\mathcal{F}_{\delta}(P) the subset of δ\delta-dimensional faces. We shall assume a fixed enumeration 0(P)={p1,,pn}\mathcal{F}_{0}(P)=\{p_{1},...,p_{n}\} of the polytope’s vertices (i.e., our polytopes are labelled), in particular, the number of vertices will be denoted by nn. We also often use a polytope QeQ\subset\mathbb{R}^{e} whose vertices are denoted 0(Q)={q1,,qn}\mathcal{F}_{0}(Q)=\{q_{1},...,q_{n}\}.

The edge-graph of PP is the finite simple graph GP=(V,E)G_{P}=(V,E), where V={1,,n}V=\{1,...,n\} is compatible with the vertex labelling, that is, iVi\in V corresponds to pi0(P)p_{i}\in\mathcal{F}_{0}(P) and ijEij\in E if and only if conv{pi,pj}1(P)\operatorname{conv}\{p_{i},p_{j}\}\in\mathcal{F}_{1}(P). The graph embedding given by ipii\mapsto p_{i} (with edges embedded as line segments) is called (1-)skeleton sk(P)\operatorname{sk}(P) of PP.

When speaking of combinatorially equivalent polytopes PP and QQ, we shall implicitly fix a face-lattice isomorphism ϕ:(P)(Q)\phi\colon\mathcal{F}(P)\xrightarrow{\sim}\mathcal{F}(Q) compatible with the vertex labels, i.e., ϕ(pi)=qi\phi(p_{i})=q_{i}. This also allows us to implicitly associate faces of PP to faces of QQ, for example, for a face σ(P)\sigma\in\mathcal{F}(P) we can write σQ\sigma_{Q} for the corresponding face in (Q)\mathcal{F}(Q). Likewise, if PP and QQ are said to have isomorphic edge-graphs, we implicitly assume a graph isomorphism GPGQG_{P}\xrightarrow{\sim}G_{Q} sending pip_{i} onto qiq_{i}. We will then often say that PP and QQ have a common edge-graph, say, GPG_{P}.

We write PQP\simeq Q to denote that PP and QQ are isometric. Since our polytopes are usually suitably translated, if not stated otherwise, this isometry can be assumed as realized by an orthogonal transformation.

Let us repeat ˜1.2 using our terminology:

Conjecture 1.2.

Given polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e} with the same edge-graph GP=(V,E)G_{P}=(V,E), so that

  1. ()

    0int(Q)0\in\operatorname{int}(Q),

  2. ()

    edges in QQ as most as long as in PP, i.e.,

    qiqjpipj,for all ijE,\|q_{i}-q_{j}\|\leq\|p_{i}-p_{j}\|,\quad\text{for all $ij\in E$},
  3. ()

    vertex-origin distances in QQ are at least as larger as in PP, i.e.,

    qipi,for all iV,\|q_{i}\|\geq\|p_{i}\|,\quad\text{for all $i\in V$},

then PQP\simeq Q.

1.2. Structure of the article

In Section˜2 we prove the instructive special case of ˜1.2 where both PP and QQ are simplices. While comparatively straightforward, the proof helps us to identify a quantity – we call it the expansion of a polytope – that is at the core of a more general approach.

The goal of Section˜3 is to show that the “expansion” of a polytope is monotone in its edge lengths, that is, decreases when the edge lengths shrink. In fact, we verify this in the more general context that replaces QQ by a general embedding q:V(GP)dq\colon V(G_{P})\to\mathbb{R}^{d} of PP’s edge-graph. As a main tool we introduce the Wachspress coordinates (a special class of generalized barycentric coordinates) and discuss a theorem of Ivan Izmestiev.

In Section˜4 we apply these results to prove ˜1.2 for the three special cases: centrally symmetric, close-by and combinatorially equivalent polytopes. We also discuss the special case of inscribed polytopes. We elaborate how our tools can potentially be used to attack ˜1.2.

In Section˜5 we conclude our investigation with further thoughts on our results, notes on connections to the literature, as well as many questions and future research directions. Despite being a conclusion section, it is quite rich in content, as we found it more appropriate to gather many notes there rather than to repeatedly interrupt the flow of the main text.

2. Warmup: a proof for simplices

To get acquainted with the task we discuss the instructive special case of ˜1.2 where both PP and QQ are simplices. The proof is reasonably short but contains already central ideas for the general case.

Theorem 2.1.

Let P,QdP,Q\subset\mathbb{R}^{d} be two simplices so that

  1. ()

    0int(Q)0\in\operatorname{int}(Q),

  2. ()

    edges in QQ are at most as long as in PP, and

  3. ()

    vertex-origin distances in QQ are at least as large as in PP,

then PQP\simeq Q.

Proof.

By (i) we can choose barycentric coordinates αintΔn\alpha\in\operatorname{int}\Delta_{n} for the origin in QQ, that is, 0=iαiqi0=\sum_{i}\alpha_{i}q_{i}. Consider the following system of equalities and inequalities:

iαipi2\displaystyle\sum_{i}\alpha_{i}\|p_{i}\|^{2} =iαipi2+12i,jαiαjpipj2\displaystyle=\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}\!+\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|p_{i}-p_{j}\|^{2}
(2.1)
iαiqi2\displaystyle\sum_{i}\alpha_{i}\|q_{i}\|^{2} =iαiqi2+12i,jαiαjqiqj2\displaystyle=\Big{\|}\sum_{i}\alpha_{i}q_{i}\Big{\|}^{2}\!+\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|q_{i}-q_{j}\|^{2}

The equalities of the first and second row can be verified by rewriting the norms as inner products followed by a straightforward computation. The vertical inequalities follow, from left to right, using (iii), the definition of α\alpha, and (ii) respectively.

But considering this system of (in)equalities, we must conclude that all inequalities are actually satisfied with equality. In particular, equality in the right-most terms yields pipj=qiqj\|p_{i}-p_{j}\|=\|q_{i}-q_{j}\| for all i,jV(GP)i,j\in V(G_{P}) (here we are using αi>0\alpha_{i}>0). But sets of points with pairwise identical distances are isometric. ∎

Why can’t we apply this proof to general polytopes? The right-most sum in \tagform@2.1 iterates over all vertex pairs and measures, if you will, a weighted average of pairwise vertex distances in PP. In simplices each vertex pair forms an edge, and hence, if all edges decrease in length, this average decreases as well. In general polytopes however, when edge become shorter, some “non-edge vertex distances” might still increase, and so the right-most inequality cannot be obtained in the same term-wise fashion. In fact, there is no reason to expect that the inequality holds at all.

It should then be surprising to learn that it actually does hold, at least in some controllable circumstances that we explore in the next section. This will allow us to generalize Theorem˜2.1 beyond simplices.

3. α\alpha-expansion, Wachspress coordinates and the Izmestiev matrix

Motivated by the proof of the simplex case (Theorem˜2.1) we define the following measure of size for a polytope (or graph embedding p:V(GP)dp\colon V(G_{P})\to\mathbb{R}^{d}):

Definition 3.1.

For αΔn\alpha\in\Delta_{n} the α\alpha-expansion of PP is

Pα2:=12i,jαiαjpipj2.\|P\|_{\alpha}^{2}:=\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|p_{i}-p_{j}\|^{2}.

The sum in the definition iterates over all pairs of vertices and so the α\alpha-expansion measures a weighted average of vertex distances, in particular, Pα\|P\|_{\alpha} is a translation invariant measure. If all pairwise distances between vertices decrease, so does the α\alpha-expansion.

The surprising fact, and main result of this section (Theorem˜3.2), is that for a carefully chosen αΔn\alpha\in\Delta_{n} the α\alpha-expansion decreases already if only the edge lengths decrease, independent of what happens to other vertex distances.

In fact, this statement holds true in much greater generality and we state it already here (it mentions Wachspress coordinates which we define in the next section; one should read this as “there exist αΔn\alpha\in\Delta_{n} so that …”):

Theorem 3.2.

Let PdP\subset\mathbb{R}^{d} be a polytope with edge-graph GP=(V,E)G_{P}=(V,E) and let αΔn\alpha\in\Delta_{n} be the Wachspress coordinates of some interior point xint(P)x\in\operatorname{int}(P). If q:Vq\colon V e\to\mathbb{R}^{e} is some embedding of GPG_{P} whose edges are at most as long as in PP, then

Pαqα,\|P\|_{\alpha}\geq\|q\|_{\alpha},

with equality if and only if qq is an affine transformation of the skeleton sk(P)\operatorname{sk}(P), all edges of which are of the same length as in PP.

Indeed, Theorem˜3.2 is not so much about comparing PP with another polytope, but actually about comparing the skeleton sk(P)\operatorname{sk}(P) with some other graph embedding qq that might not be the skeleton of any polytope and might even be embedded in a lower- or higher-dimensional Euclidean space. Morally, Theorem˜3.2 says: polytope skeleta are maximally expanded for their edge lengths, where “expansion” here measures an average of vertex distances with carefully chosen weights.

The result clearly hinges on the existence of these so-called Wachspress coordinates, which we introduce now.

3.1. Wachspress coordinates

In a simplex PdP\subset\mathbb{R}^{d} each point xPx\in P can be expressed as a convex combination of the simplex’s vertices in a unique way:

x=iαipi.x=\sum_{i}\alpha_{i}p_{i}.

The coefficients αΔn\alpha\in\Delta_{n} are called the barycentric coordinates of xx in PP.

In a general polytope PdP\subset\mathbb{R}^{d} there are usually many ways to express a point xPx\in P as a convex combination of the polytope’s vertices. In many applications however it is desirable to have a canonical choice, so to say “generalized barycentric coordinates”. Various such coordinates have been defined (see [9] for an overview), one of them being the so-called Wachspress coordinates. Those were initially defined by Wachspress for polygons [28], and later generalized to general polytopes by Warren et al. [30, 32]. A construction, with a geometric interpretation due to [18], is given in Section˜3.3 below.

The relevance of the Wachspress coordinates for our purpose is however not so much in their precise definition, but rather in their relation to a polytope invariant of “higher rank” that we introduced next.

3.2. The Izmestiev matrix

At the core of our proof of Theorem˜3.2 is the observation that the Wachspress coordinates are merely a shadow of a “higher rank” object that we call the Izmestiev matrix of PP; an (n×n)(n\times n)-matrix associated to an nn-vertex polytope with 0int(P)0\in\operatorname{int}(P), whose existence and properties in connection with graph skeleta were established by Lovász in dimension three [21], and by Izmestiev in general dimension [15]. We summarize the findings:

Theorem 3.3.

Given a polytope PdP\subset\mathbb{R}^{d} with 0int(P)0\in\operatorname{int}(P) and edge-graph GP=(V,G_{P}=(V, E)E), there exists a symmetric matrix Mn×nM\in\mathbb{R}^{n\times n} (the Izmestiev matrix of PP) with the following properties:

  1. ()

    Mij>0M_{ij}>0 if ijEij\in E,

  2. ()

    Mij=0M_{ij}=0 if iji\not=j and ijEij\not\in E,

  3. ()

    dimkerM=d\dim\ker M=d,

  4. ()

    MXP=0MX_{P}=0, where XP=(p1,,pn)d×nX_{P}^{\top}=(p_{1},...,p_{n})\in\mathbb{R}^{d\times n}, and

  5. ()

    MM has a unique positive eigenvalue (of multiplicity one).

Izmestiev provided an explicit construction of this matrix that we discuss in Section˜3.3 below. Another concise proof of the spectral properties of the Izmestiev matrix can be found in the appendix of [22].

Observation 3.4.

Each of the properties (i) to (v) of the Izmestiev matrix will be crucial for proving Theorem˜3.2 and we shall elaborate on each point below:

  1. ()

    Theorem˜3.3 (i) and (ii) state that MM is some form of generalized adjacency matrix, having non-zero off-diagonal entries if and only if the polytope has an edge between the corresponding vertices. Note however that the theorem tells nothing directly about the diagonal entries of MM.

  2. ()

    Theorem˜3.3 (iii) and (iv) tell us precisely how the kernel of MM looks like, namely, kerM=spanXP\ker M=\operatorname{span}X_{P}. The inclusion kerMspanXP\ker M\supseteq\operatorname{span}X_{P} follows directly from (iv). But since PP has at least one interior point (the origin) it must be a full-dimensional polytope, meaning that rankXP=d\operatorname{rank}X_{P}=d. Comparison of dimensions (via (iii)) yields the claimed equality.

  3. ()

    let {θ1>θ2>>θm}\{\theta_{1}>\theta_{2}>\cdots>\theta_{m}\} be the spectrum of MM. Theorem˜3.3 (v) then tells us that θ1>0\theta_{1}>0, θ2=0\theta_{2}=0 and θk<0\theta_{k}<0 for all k3k\geq 3.

  4. ()

    M:=M+γIdM^{\prime}:=M+\gamma\operatorname{Id} is a non-negative matrix if γ>0\gamma>0 is sufficiently large, and is then subject to the Perron-Frobenius theorem (see Theorem˜A.1). Since the edge-graph GPG_{P} is connected, the matrix MM^{\prime} is irreducible. The crucial information provided by the Perron-Frobenius theorem is that MM^{\prime} has an eigenvector znz\in\mathbb{R}^{n} to its largest eigenvalue (that is, θ1+γ\theta_{1}+\gamma), all entries of which are positive. By an appropriate scaling we can assume zint(Δn)z\in\operatorname{int}(\Delta_{n}), which is a θ1\theta_{1}-eigenvector to the Izmestiev matrix MM, and in fact, spans its θ1\theta_{1}-eigenspace.

Note that the properties (i) to (v) in Theorem˜3.3 are invariant under scaling of MM by a positive factor. As we verify in Section˜3.3 below, i,jMij>0\smash{\sum_{i,j}M_{ij}>0}, and so we can fix the convenient normalization i,jMij=1\smash{\sum_{i,j}M_{ij}=1}. In fact, with this normalization in place we can now reveal that the Wachspress coordinates emerge simply as the row sums of MM:

(3.1) αi:=jMij,for all i{1,,n}.\alpha_{i}:=\sum_{j}M_{ij},\quad\text{for all $i\in\{1,...,n\}$}.

This connection has previously been observed in [18, Section 4.2] for 3-dimensional polytopes, and we shall verify the general case in the next section (Corollary˜3.6).

3.3. The relation between Wachspress and Izmestiev

The Wachspress coordinates and the Izmestiev matrix can be defined simultaneously in a rather elegant fashion: given a polytope PdP\subset\mathbb{R}^{d} with d2d\geq 2 and 0int(P)0\in\operatorname{int}(P), as well as a vector 𝐜=\mathbf{c}= (c1,,cn)n(c_{1},...,c_{n})\in\mathbb{R}^{n}, consider the generalized polar dual

P(𝐜):={xdx,pici for all iV(GP)}.P^{\circ}(\mathbf{c}):=\big{\{}x\in\mathbb{R}^{d}\mid\langle x,p_{i}\rangle\leq c_{i}\text{ for all $i\in V(G_{P})$}\big{\}}.

We have that P(𝟏)P^{\circ}(\mathbf{1}) with 𝟏=(1,,1)\mathbf{1}=(1,...,1) is the usual polar dual. The (unnormalized) Wachspress coordinates α~n\tilde{\alpha}\in\mathbb{R}^{n} of the origin and the (unnormalized) Izmestiev matrix M~\tilde{M} \in n×n\mathbb{R}^{n\times n} emerge as the coefficients in the Taylor expansion of the volume of P(𝐜)P^{\circ}(\mathbf{c}) at 𝐜=𝟏\mathbf{c}=\mathbf{1}:

(3.2) vol(P(𝐜))=vol(P)+𝐜𝟏,α~+12(𝐜𝟏)M~(𝐜𝟏)+.\operatorname{vol}\!\big{(}P^{\circ}(\mathbf{c})\big{)}=\operatorname{vol}(P^{\circ})+\langle\mathbf{c}-\mathbf{1},\tilde{\alpha}\rangle+\tfrac{1}{2}(\mathbf{c}-\mathbf{1})^{\top}\!\tilde{M}(\mathbf{c}-\mathbf{1})+\cdots.

In other words,

(3.3) α~i:=vol(P(𝐜))ci|𝐜=𝟏andM~ij:=2vol(P(𝐜))cicj|𝐜=𝟏.\tilde{\alpha}_{i}:=\frac{\partial\operatorname{vol}(P^{\circ}(\mathbf{c}))}{\partial c_{i}}\Big{|}_{\mathbf{c}=\mathbf{1}}\quad\text{and}\quad\tilde{M}_{ij}:=\frac{\partial^{2}\operatorname{vol}(P^{\circ}(\mathbf{c}))}{\partial c_{i}\partial c_{j}}\Big{|}_{\mathbf{c}=\mathbf{1}}.

In this form, one might recognize the (unnormalized) Izmestiev matrix of PP as the Alexandrov matrix of the polar dual PP^{\circ}.

Geometric interpretations for \tagform@3.3 were given in [18, Section 3.3] and [15, proof of Lemma 2.3]: for a vertex pi0(P)p_{i}\in\mathcal{F}_{0}(P) let Fid1(P)F_{i}\in\mathcal{F}_{d-1}(P^{\circ}) be the corresponding dual facet. Likewise, for an edge eij1(P)e_{ij}\in\mathcal{F}_{1}(P) let σijd2(P)\sigma_{ij}\in\mathcal{F}_{d-2}(P^{\circ}) be the corresponding dual face of codimension 2. Then

(3.4) α~i=vol(Fi)piandM~ij=vol(σij)pipjsin(pi,pj),\tilde{\alpha}_{i}=\frac{\operatorname{vol}(F_{i})}{\|p_{i}\|}\quad\text{and}\quad\tilde{M}_{ij}=\frac{\operatorname{vol}(\sigma_{ij})}{\|p_{i}\|\|p_{j}\|\sin\sphericalangle(p_{i},p_{j})},

where vol(Fi)\operatorname{vol}(F_{i}) and vol(σij)\operatorname{vol}(\sigma_{ij}) are to be understood as relative volume. The expression for α~\tilde{\alpha} is (up to a constant factor) the cone volume of FiF_{i} in PP^{\circ}​, i.e., the volume of the cone with base face FiF_{i} and apex at the origin. As such it is positive, which confirms again that we can normalize to αΔn\alpha\in\Delta_{n}, and we see that αi\alpha_{i} measures the fraction of the cone volume at FiF_{i} in the total volume of PP^{\circ}. That M~\smash{\tilde{M}} can be normalized follows from the next statement, which is a precursor to \tagform@3.1:

Proposition 3.5.

jM~ij=(d1)α~i\sum_{j}\tilde{M}_{ij}=(d-1)\tilde{\alpha}_{i}.

Proof.

Observe first that vol(P(𝐜))\operatorname{vol}(P^{\circ}(\mathbf{c})) is a homogeneous function of degree dd, i.e.,

vol(P(t𝐜))=vol(tP(𝐜))=tdvol(P(𝐜))\operatorname{vol}(P^{\circ}(t\mathbf{c}))=\operatorname{vol}(tP^{\circ}(\mathbf{c}))=t^{d}\operatorname{vol}(P^{\circ}(\mathbf{c}))

for all t0t\geq 0. Each derivative vol(P(𝐜))/ci\partial\operatorname{vol}(P^{\circ}(\mathbf{c}))/\partial c_{i} is then homogeneous of degree d1d-1. Euler’s homogeneous function theorem (Theorem˜E.1) yields

jcj2vol(P(𝐜))cicj=(d1)vol(P(𝐜))ci.\sum_{j}c_{j}\frac{\partial^{2}\operatorname{vol}(P^{\circ}(\mathbf{c}))}{\partial c_{i}\partial c_{j}}=(d-1)\frac{\partial\operatorname{vol}(P^{\circ}(\mathbf{c}))}{\partial c_{i}}.

Evaluating at 𝐜=1\mathbf{c}=1 and using \tagform@3.3 yields the claim. ∎

We immediately see that i,jM~ij>0\sum_{i,j}\tilde{M}_{ij}>0 and that we can normalize to i,jMij=1\sum_{i,j}M_{ij}=1. For the normalized quantities then indeed holds \tagform@3.1:

Corollary 3.6.

iMij=αj\sum_{i}M_{ij}=\alpha_{j} for all j{1,,n}j\in\{1,...,n\}.

Lastly, the following properties of the Wachspress coordinates and the Izmestiev matrix will be relevant and can be inferred from the above.

Remark 3.7.
  1. ()

    The Wachspress coordinates of the origin and the Izmestiev matrix depend continuously on the translation of PP, and their normalized variants can be continuously extended to 0P0\in\partial P. If the origin lies in the relative interior of a face σ(P)\sigma\in\mathcal{F}(P), then αi>0\alpha_{i}>0 if and only if piσp_{i}\in\sigma. In particular, if 00\in int(P)\operatorname{int}(P), then αintΔn\alpha\in\operatorname{int}\Delta_{n}.

  2. ()

    The Wachspress coordinates of the origin and the Izmestiev matrix are invariant under linear transformation of PP. This can be inferred from \tagform@3.4 via an elementary computation, as was done for the Izmestiev matrix in [35, Proposition 4.6.].

3.4. Proof of Theorem˜3.2

Recall the main theorem.

Theorem 3.2.

Let PdP\subset\mathbb{R}^{d} be a polytope with edge-graph GP=(V,E)G_{P}=(V,E) and let αΔn\alpha\in\Delta_{n} be the Wachspress coordinates of some interior point xint(P)x\in\operatorname{int}(P). If q:Vq\colon V e\to\mathbb{R}^{e} is some embedding of GPG_{P} whose edges are at most as long as in PP, then

Pαqα,\|P\|_{\alpha}\geq\|q\|_{\alpha},

with equality if and only if qq is an affine transformation of the skeleton sk(P)\operatorname{sk}(P), all edges of which are of the same length as in PP.

The proof presented below is completely elementary, using little more than linear algebra. In Appendix˜F the reader can find a second shorter proof based on the duality theory of semi-definite programming.

Proof.

At the core of this proof is rewriting the α\alpha-expansions Pα\|P\|_{\alpha} and qα\|q\|_{\alpha} as a sum of terms, each of which is non-increasing when transitioning from PP to qq:

(3.5) Pα2\displaystyle\|P\|_{\alpha}^{2}\, =ijEMijpipj2iαipi2+tr(MXPXP)\displaystyle=\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2}-\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}+\operatorname{tr}(MX_{P}X_{P}^{\top})
ijEMijqiqj2iαiqi2+tr(MXqXq)=qα2.\displaystyle\;\,\phantom{\geq}\sum_{ij\in E}M_{ij}\|q_{i}-q_{j}\|^{2}\,-\Big{\|}\sum_{i}\alpha_{i}q_{i}\Big{\|}^{2}+\operatorname{tr}(MX_{q}X_{q}^{\top})=\,\|q\|_{\alpha}^{2}.

Of course, neither the decomposition nor the monotonicity of the terms is obvious; yet their proofs use little more than linear algebra. We elaborate on this now.

For the setup, we recall that the α\alpha-expansion is a translation invariant measure of size. We can therefore translate PP and qq to suit our needs:

  1. ()

    translate PP so that x=0x=0, that is, iαipi=0\sum_{i}\alpha_{i}p_{i}=0.

  2. ()

    since then 0int(P)0\in\operatorname{int}(P), Theorem˜3.3 ensures the existence of the Izmestiev matrix Mn×nM\in\mathbb{R}^{n\times n}.

  3. ()

    Let θ1>θ2>>θm\theta_{1}>\theta_{2}>\cdots>\theta_{m} be the eigenvalues of MM, where θ1>0\theta_{1}>0 and θ2=0\theta_{2}=0. By Observation˜3.4 (iv) there exists a unique θ1\theta_{1}-eigenvector zint(Δn)z\in\operatorname{int}(\Delta_{n}).

  4. ()

    translate qq so that iziqi=0\sum_{i}z_{i}q_{i}=0.

We are ready to derive the decompositions shown in \tagform@3.5: the following equality can be verified straightforwardly by rewriting the square norms as inner products:

12i,jMijpipj2\displaystyle\tfrac{1}{2}\sum_{i,j}M_{ij}\|p_{i}-p_{j}\|^{2} =i(jMij)pi2i,jMijpi,pj,\displaystyle=\sum_{i}\Big{(}\sum_{j}M_{ij}\Big{)}\|p_{i}\|^{2}-\sum_{i,j}M_{ij}\langle p_{i},p_{j}\rangle,

We continue rewriting each of the three terms:

  • on the left: Mijpipj2M_{ij}\|p_{i}-p_{j}\|^{2} is only non-zero for ijEij\in E (using Theorem˜3.3 (ii)). The sum can therefore be rewritten to iterate over the edges of GPG_{P} (where we consider ij,jiEij,ji\in E the same and so we can drop the factor 1/2\nicefrac{{1}}{{2}})

  • in the middle: the row sums of the Izmestiev matrix are exactly the Wachspress coordinates of the origin, that is, jMij=αi\smash{\sum_{j}M_{ij}=\alpha_{i}}.

  • on the right: recall the matrix XPd×nX_{P}\in\mathbb{R}^{d\times n} whose rows are the vertex coordinates of PP. The corresponding Gram matrix XPXPX_{P}X_{P}^{\top} has entries (XPXP)ij=pi,pj(X_{P}X_{P}^{\top})_{ij}=\langle p_{i},p_{j}\rangle.

By this we reach the following equivalent identity:

ijEMijpipj2\displaystyle\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2} =iαipi2i,jMij(XPXP)ij,\displaystyle=\sum_{i}\alpha_{i}\|p_{i}\|^{2}-\sum_{i,j}M_{ij}(X_{P}X_{P}^{\top})_{ij},

We continue rewriting the terms on the right side of the equation:

  • in the middle: the following transformation was previously used in the simplex case (Theorem˜2.1) and can be verified by straightforward expansion of the squared norms:

    iαipi2=12i,jαiαjpipj2+iαipi2.\sum_{i}\alpha_{i}\|p_{i}\|^{2}=\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|p_{i}-p_{j}\|^{2}+\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}.

    Note that the middle term is just the α\alpha-expansion Pα2\|P\|_{\alpha}^{2}.

  • on the right: the sum iterates over entry-wise products of the two matrices MM and XPXPX_{P}X_{P}^{\top}, which can be rewritten as tr(MXPXP)\operatorname{tr}(MX_{P}X_{P}^{\top}).

Thus, we arrive at

ijEMijpipj2\displaystyle\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2} =Pα2+iαipi2tr(MXPXP).\displaystyle=\|P\|_{\alpha}^{2}+\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}\!-\operatorname{tr}(MX_{P}X_{P}^{\top}).

This clearly rearranges to the first line of \tagform@3.5. An analogous sequence of transformations works for qq (we replace pip_{i} by qiq_{i} and XPX_{P} by XqX_{q}, but we keep the Izmestiev matrix of PP). This yields the second line of \tagform@3.5. It remains to verify the term-wise inequalities.

For the first term we have

ijEMijpipj2ijEMijqiqj2\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2}\geq\sum_{ij\in E}M_{ij}\|q_{i}-q_{j}\|^{2}

by term-wise comparison: we use that the sum is only over edges, that Mij>0M_{ij}>0 for ijEij\in E (by Theorem˜3.3 (i)), and that edges in qq are not longer than in PP.

Next, by the wisely chosen translation in setup (i) we have iαipi=0\sum_{i}\alpha_{i}p_{i}=0, thus

iαipi2=0iαiqi2.-\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}\!=0\geq-\Big{\|}\sum_{i}\alpha_{i}q_{i}\Big{\|}^{2}.

The final term requires the most elaboration. By Theorem˜3.3 (iv) the Izmestiev matrix satisfies MXP=0MX_{P}=0. So it suffices to show that tr(MXqXq)\operatorname{tr}(MX_{q}X_{q}^{\top}) is non-positive, as then already follows

(3.6) tr(MXPXP)=0?tr(MXqXq).\operatorname{tr}(MX_{P}X_{P}^{\top})=0\overset{?}{\geq}\operatorname{tr}(MX_{q}X_{q}^{\top}).

To prove tr(MXqXq)0\smash{\operatorname{tr}(MX_{q}X_{q}^{\top})\leq 0} consider the decomposition Xq=Xq1+\smash{X_{q}=X_{q}^{1}+\cdots} +Xqm\smash{+\,X_{q}^{m}} where MXqk=θkXqkMX_{q}^{k}=\theta_{k}X_{q}^{k} (since MM is symmetric, its eigenspaces are orthogonal and XqkX_{q}^{k} is the column-wise orthogonal projecting of XqX_{q} onto the θk\theta_{k}-eigenspace). We compute

tr(MXqXq)\displaystyle\operatorname{tr}(MX_{q}X_{q}^{\top}) =k,tr(MXqk(Xq))\displaystyle=\sum_{k,\ell}\operatorname{tr}(MX_{q}^{k}(X_{q}^{\ell})^{\top})
=k,θktr(Xqk(Xq))\displaystyle=\sum_{k,\ell}\theta_{k}\operatorname{tr}(X_{q}^{k}(X_{q}^{\ell})^{\top}) |by MXqk=θkXqk\displaystyle|\;\text{by $MX_{q}^{k}=\theta_{k}X_{q}^{k}$}
=k,θktr((Xq)Xqk)\displaystyle=\sum_{k,\ell}\theta_{k}\operatorname{tr}((X_{q}^{\ell})^{\top}X_{q}^{k}) |by tr(AB)=tr(BA)\displaystyle|\;\text{by $\operatorname{tr}(AB)=\operatorname{tr}(BA)$}
=kθktr((Xqk)Xqk).\displaystyle=\sum_{k}\theta_{k}\operatorname{tr}((X_{q}^{k})^{\top}X_{q}^{k}). |since (Xq)Xqk=0 when k\displaystyle|\;\text{since $(X_{q}^{\ell})^{\top}X_{q}^{k}=0$ when $k\not=\ell$}

Again, we have been wise in our choice of translation of qq in setup (iv): iziqi=0\sum_{i}z_{i}q_{i}=0 can be written as zXq=0z^{\top}X_{q}=0. Since zz spans the θ1\theta_{1}-eigenspace, the columns of XqX_{q} are therefore orthogonal to the θ1\theta_{1}-eigenspace, hence Xq1=0X_{q}^{1}=0. We conclude

(3.7) tr(MXqXq)\displaystyle\operatorname{tr}(MX_{q}X_{q}^{\top}) =k2θktr((Xqk)Xqk)0,\displaystyle=\sum_{k\geq 2}\theta_{k}\operatorname{tr}((X_{q}^{k})^{\top}X_{q}^{k})\leq 0,

where the final inequality follows from two observations: first, the Izmestiev matrix MM has a unique positive eigenvalue θ1\theta_{1}, thus θk0\theta_{k}\leq 0 for all k2k\geq 2 (Theorem˜3.3 (v)); second (Xqk)Xqk(X_{q}^{k})^{\top}X_{q}^{k} is a Gram matrix, hence is positive semi-definite and has a non-negative trace.

This finalizes the term-wise comparison and established the inequality \tagform@3.5. It remains to discuss the equality case. By now we see that the equality Pα=qα\|P\|_{\alpha}=\|q\|_{\alpha} is equivalent to term-wise equality in \tagform@3.5; and so we proceed term-wise.

To enforce equality in the first term

ijEMijpipj2=!ijEMijqiqj2\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2}\overset{!}{=}\sum_{ij\in E}M_{ij}\|q_{i}-q_{j}\|^{2}

we recall again that Mij>0M_{ij}>0 whenever ijEij\in E by Theorem˜3.3 (i). Thus, we require equality pipj=qiqj\|p_{i}-p_{j}\|=\|q_{i}-q_{j}\| for all edges ijEij\in E. And so edges in qq must be of the same length as in PP.

We skip the second term for now and enforce equality in the last term:

0=tr(MXPXP)=!tr(MXqXq)=\tagform@3.7k2θktr((Xqk)Xqk).0=\operatorname{tr}(MX_{P}X_{P}^{\top})\overset{!}{=}\operatorname{tr}(MX_{q}X_{q}^{\top})\overset{\hyperref@@ii[eq:tr_final_step]{\textup{\tagform@{\ref*{eq:tr_final_step}}}}}{=}\sum_{k\geq 2}\theta_{k}\operatorname{tr}((X_{q}^{k})^{\top}X_{q}^{k}).

Since θk<0\theta_{k}<0 for all k3k\geq 3 (cf. Observation˜3.4 (iii)), for the sum on the right to vanish we necessarily have

tr((Xqk)Xqk)=0 for all k3Xqk=0 for all k3.\operatorname{tr}((X_{q}^{k})^{\top}X_{q}^{k})=0\;\text{ for all $k\geq 3$}\quad\implies\quad X_{q}^{k}=0\;\text{ for all $k\geq 3$}.

Since we also already know Xq1=0X_{q}^{1}=0, we are left with Xq=Xq2\smash{X_{q}=X_{q}^{2}}, that is, the columns of XqX_{q} are in the θ2\theta_{2}-eigenspace (aka. the kernel) of MM. In particular, spanXqkerM\operatorname{span}X_{q}\subseteq\ker M =spanXP=\operatorname{span}X_{P}, where the last equality follows by Observation˜3.4 (ii). It is well-known that if two matrices satisfy spanXqspanXP\operatorname{span}X_{q}\subseteq\operatorname{span}X_{P}, then the rows of XqX_{q} are linear transformations of the rows of XPX_{P}, that is, TXP=XqTX_{P}^{\top}=X_{q}^{\top} for some linear map T:deT\colon\mathbb{R}^{d}\to\mathbb{R}^{e}, or equivalently, qi=Tpiq_{i}=Tp_{i} for all iVi\in V (see Theorem˜C.1 in the appendix for a short reminder of the proof). Therefore, qq (considered with its original translation prior to the setup) must have been an affine transformation of sk(P)\operatorname{sk}(P).

Lastly we note that equality in the middle term of \tagform@3.5 yields no new constraints. In fact, by spanXqkerM\operatorname{span}X_{q}\subseteq\ker M we have MXq=0MX_{q}=0 and

iαiqi=i(jMij)qi=j(iMijqi)=0=iαipi.\sum_{i}\alpha_{i}q_{i}=\sum_{i}\Big{(}\sum_{j}M_{ij}\Big{)}q_{i}=\sum_{j}\Big{(}\sum_{i}M_{ij}q_{i}\Big{)}=0=\sum_{i}\alpha_{i}p_{i}.

Thus, identity in the middle term follows already from identity in the last term.

For the other direction of the identity case assume that qq is an affine transformation of sk(P)\operatorname{sk}(P) with the same edge lengths. Instead of setup (iv) assume a translation of qq for which it is a linear transformation of sk(P)\operatorname{sk}(P), i.e., Xq=TXP\smash{X_{q}^{\top}=TX_{P}^{\top}} for some linear map T:deT\colon\mathbb{R}^{d}\to\mathbb{R}^{e}. Hence iαipi=iαiqi=0\smash{\sum_{i}\alpha_{i}p_{i}=\sum_{i}\alpha_{i}q_{i}=0} and MXP=MXq=0MX_{P}=MX_{q}=0, and \tagform@3.5 reduces to

Pα2=ijEMijpipj20+0=ijEMijqiqj20+0=qα2.\|P\|_{\alpha}^{2}=\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2}-0+0=\sum_{ij\in E}M_{ij}\|q_{i}-q_{j}\|^{2}-0+0=\|q\|_{\alpha}^{2}.

As an immediate consequence we have the following:

Corollary 3.8.

A polytope is uniquely determined (up to affine equivalence) by its edge-graph, its edge lengths and the Wachspress coordinates of some interior point.

Proof.

Given polytopes P1P_{1} and P2P_{2} with the same edge-graph and edge lengths as well as points xiint(Pi)x_{i}\in\operatorname{int}(P_{i}) with the same Wachspress coordinates αΔn\alpha\in\Delta_{n}. By Theorem˜3.2 we have P1αP2αP1α\|P_{1}\|_{\alpha}\geq\|P_{2}\|_{\alpha}\geq\|P_{1}\|_{\alpha}, thus P1α=P2α\|P_{1}\|_{\alpha}=\|P_{2}\|_{\alpha}. Then P1P_{1} and P2P_{2} are affinely equivalent by the equality case of Theorem˜3.2. ∎

Remarkably, this reconstruction works across all combinatorial types and dimensions. That the reconstruction is only up to affine equivalence rather than isometry is due to examples such as rhombi and the zonotope in Figure˜5. In general, such flexibility via an affine transformation happens exactly if “the edge directions lie on a conic at infinity” (see [3, Proposition 4.2] or [5, Proposition 1.4]).

Refer to caption
Figure 5. Two affinely equivalent (but not isometric) polytopes with the same edge length. The edge directions trace a circle on the “plane at infinity”.

Lastly, the reconstruction permitted by Corollary˜3.8 is feasible in practice. This follows from a reformulation of Theorem˜3.2 as a semi-definite program, which can be solved in polynomial time. This is elaborated on in the alternative proof given in Appendix˜F.

4. Rigidity, tensegrity and reconstruction

Our reason for pursuing Theorem˜3.2 in Section˜3 was to transfer the proof of the simplex case (Theorem˜2.1) to general polytopes with the eventual goal of verifying the main conjecture and its corresponding “unique reconstruction version”:

Conjecture 1.2.

Given polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e} with common edge-graph. If

  1. ()

    0int(Q)0\in\operatorname{int}(Q),

  2. ()

    edges in QQ are at most as long as in PP, and

  3. ()

    vertex-origin distances in QQ are at least as large as in PP,

then PQP\simeq Q.

Conjecture 1.1.

A polytope PP with 0int(P)0\in\operatorname{int}(P) is uniquely determined (up to isometry) by its edge-graph, edge lengths and vertex-origin distances.

In contrast to our formulation of Theorem˜3.2, both of the above conjectures are false when stated for a general graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} instead of QQ, even if we require 0intconv(q)0\in\operatorname{int}\operatorname{conv}(q). The following counterexample was provided by Joseph Doolittle [7]:

Example 4.1.

The 3-cube P:=[1,1]33P:=[-1,1]^{3}\subset\mathbb{R}^{3} is inscribed in a sphere of radius 3\smash{\sqrt{3}}. Figure˜6 shows an inscribed embedding q:V(GP)2q\colon V(G_{P})\to\mathbb{R}^{2} with the same circumradius and edge lengths, collapsing GPG_{P} onto a path. In the circumcircle each edge spans an arc of length

2sin1(1/3)70.5287>60.2\sin^{-1}(1/\sqrt{3})\approx 70.5287^{\circ}>60^{\circ}.

The three edges therefore suffice to reach more than half around the circle. In other words, the convex hull of qq contains the circumcenter in its interior.

A full-dimensional counterexample in 3\mathbb{R}^{3} can be obtained by interpreting qq as embedded in 2×{0}\mathbb{R}^{2}\times\{0\} follows by a slight perturbation.

Refer to caption
Figure 6. A 2-dimensional embedding of the edge-graph of the 3-dimensional cube with the same circumradius and edge lengths as the unit cube, also containing the origin in its convex hull.

Potential fixes to the “graph embedding versions” of ˜1.1 and ˜1.2 are discussed in Section˜5.2.

While the general ˜1.2 will stay open, we are confident that our methods point the right way and highlight the essential difficulties. We overcome them in three relevant special cases, for some of which we actually can replace QQ with a graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e}. Those are

  1. ()

    PP and qq are centrally symmetric (Theorem˜4.4).

  2. ()

    qq is a sufficiently small perturbation of sk(P)\operatorname{sk}(P) (Theorem˜4.5).

  3. ()

    PP and QQ are combinatorially equivalent (Theorem˜4.7).

4.1. The remaining difficulty

On trying to generalize the proof of Theorem˜2.1 beyond simplices using Theorem˜3.2 we face the following difficulty: Theorem˜3.2 requires the αΔn\alpha\in\Delta_{n} to be Wachspress coordinates of an interior point xint(P)x\in\operatorname{int}(P), whereas in the proof of Theorem˜2.1 we use that α\alpha is a set of barycentric coordinates of 00\in int(Q)\operatorname{int}(Q). While we have some freedom in choosing xint(P)x\in\operatorname{int}(P), and thereby αΔn\alpha\in\Delta_{n}, it is not clear that any such choice yields barycentric coordinates for 0int(Q)0\in\operatorname{int}(Q). In fact, this is the only obstacle preventing us from proving ˜1.2 right away. For convenience we introduce the following map:

Definition 4.2.

Given polytopes PdP\subset\mathbb{R}^{d} and QeQ\subset\mathbb{R}^{e}, the Wachspress map ϕ:PQ\phi\colon P\to Q is defined as follows: for xPx\in P with Wachspress coordinates α(x)Δn\alpha(x)\in\Delta_{n} set

ϕ(x):=iαi(x)qi.\phi(x):=\sum_{i}\alpha_{i}(x)q_{i}.

In cases where we are working with a graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} instead of QQ we have an analogous map ϕ:Pconv(q)\phi\colon P\to\operatorname{conv}(q).

Our previous discussion amounts to checking whether the origin is in the image of int(P)\operatorname{int}(P) w.r.t. ϕ\phi. While this could be reasonably true if dimPdimQ\dim P\geq\dim Q, it is certainly too much to hope for if dimP<dimQ\dim P<\dim Q: the image of ϕ(int(P))Q\phi(\operatorname{int}(P))\subset Q is of a smaller dimension than QQ and easily “misses” the origin. Fortunately, we can ask for less, which we now formalize in (i) of the following lemma:

Lemma 4.3.

Let PdP\subset\mathbb{R}^{d} be a polytope and q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} some embedding. If

  1. ()

    there exists an xint(P)x\in\operatorname{int}(P) with ϕ(x)x\|\phi(x)\|\leq\|x\| (e.g. ϕ(x)=0\phi(x)=0),

  2. ()

    edges in qq are at most as long as in PP, and

  3. ()

    vertex-origin distances in qq are at least as large as in PP,

then qsk(P)q\simeq\operatorname{sk}(P) (via an orthogonal transformation).

Proof.

Choose xintPx\in\operatorname{int}P with ϕ(x)x\|\phi(x)\|\leq\|x\|, and note that its Wachspress coordinates αintΔn\alpha\in\operatorname{int}\Delta_{n} are strictly positive. In the remainder we follow closely the proof of Theorem˜2.1: consider the system of (in)equalities:

iαipi2\displaystyle\sum_{i}\alpha_{i}\|p_{i}\|^{2} =iαipi2+12i,jαiαjpipj2=x2+Pα2\displaystyle=\Big{\|}\sum_{i}\alpha_{i}p_{i}\Big{\|}^{2}\!+\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|p_{i}-p_{j}\|^{2}=\;\;\|x\|^{2}\;\;\;\,+\|P\|_{\alpha}^{2}
iαiqi2\displaystyle\sum_{i}\alpha_{i}\|q_{i}\|^{2} =iαiqi2+12i,jαiαjqiqj2=ϕ(x)2+qα2,\displaystyle=\Big{\|}\sum_{i}\alpha_{i}q_{i}\Big{\|}^{2}\!+\tfrac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|q_{i}-q_{j}\|^{2}\,=\|\phi(x)\|^{2}\,+\,\|q\|_{\alpha}^{2},

where the two rows hold by simple computation and the vertical inequalities follow (from left to right) by (iii), (i), and (ii) + Theorem˜3.2 respectively. It follows that all inequalities are actually equalities. In particular, since αi>0\alpha_{i}>0 we find both pi\|p_{i}\| == qi\|q_{i}\| for all iVi\in V and pipj=qiqj\|p_{i}-p_{j}\|=\|q_{i}-q_{j}\| for all i,jV(GP)i,j\in V(G_{P}), establishing that qq and sk(P)\operatorname{sk}(P) are indeed isometric via an orthogonal transformation. ∎

The only way for Lemma˜4.3 (i) to fail is if ϕ(x)>x\|\phi(x)\|>\|x\| for all xint(P)x\in\operatorname{int}(P). By (ii) and (iii) we have ϕ(x)=x\|\phi(x)\|=\|x\| whenever xx is a vertex or in an edge of PP, which makes it plausible that (i) should not fail, yet it seems hard to exclude in general.

In each of the three special cases of ˜1.2 discussed below we actually managed to verify 0ϕ(int(P))0\in\phi(\operatorname{int}(P)) in order to apply Lemma˜4.3.

4.2. Central symmetry

Let PdP\subset\mathbb{R}^{d} be centrally symmetric (more precisely, origin symmetric), that is, P=PP=-P. This induces an involution ι:V(GP)V(GP)\iota\colon V(G_{P})\to V(G_{P}) with pι(i)=pip_{\iota(i)}=-p_{i} for all iV(GP)i\in V(G_{P}). We say that an embedding q:V(GP)dq\colon V(G_{P})\to\mathbb{R}^{d} of the edge-graph is centrally symmetric if qι(i)=q_{\iota(i)}= qi-q_{i} for all iV(GP)i\in V(G_{P}).

Theorem 4.4 (centrally symmetric version).

Given a centrally symmetric polytope PdP\subset\mathbb{R}^{d} and a centrally symmetric graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e}, so that

  1. ()

    edges in qq are at most as long as in PP, and

  2. ()

    vertex-origin distances in qq are at least as large as in PP,

then qsk(P)q\simeq\operatorname{sk}(P).

Proof.

Since PP is centrally symmetric, we have 0relint(P)0\in\operatorname{relint}(P) and we can find Wachspress coordinates αint(Δn)\alpha\in\operatorname{int}(\Delta_{n}) of the origin in PP. Since the Wachspress coordinates are invariant under a linear transformation (as noted in Remark˜3.7 (ii)), it holds αi=αι(i)\alpha_{i}=\alpha_{\iota(i)}. For the Wachspress map ϕ\phi follow

ϕ(0)=12iVαiqi+12iVαι(i)qι(i)=12iVαiqi12iVαiqi=0,\phi(0)=\tfrac{1}{2}\sum_{i\in V}\alpha_{i}q_{i}+\tfrac{1}{2}\sum_{i\in V}\alpha_{\iota(i)}q_{\iota(i)}=\tfrac{1}{2}\sum_{i\in V}\alpha_{i}q_{i}-\tfrac{1}{2}\sum_{i\in V}\alpha_{i}q_{i}=0,

The claim then follows via Lemma˜4.3

It is clear that Theorem˜4.4 can be adapted to work with other types of symmetry that uniquely fix the origin, such as irreducible symmetry groups.

Theorem˜4.4 has a natural interpretation in the language of classical rigidity theory, where it asserts the universal rigidity of a certain tensegrity framework. In this form it was proven up to dimension three by Connelly [4, Theorem 5.1]. We elaborate further on this in Section˜5.3.

It is now tempting to conclude the unique reconstruction version of Theorem˜4.4, answering ˜1.1 for centrally symmetric polytope. There is however a subtlety: our notion of “central symmetry” for the graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} as used in Theorem˜4.4 has been relative to PP, in that it forces qq to have the same pairs of antipodal vertices as PP. It is however not true that any two centrally symmetric polytopes with the same edge-graph have this relation. David E. Speyer [25] constructed a 12-vertex 4-polytope whose edge-graph has an automorphism that does not preserve antipodality of vertices.

4.3. Local uniqueness

Given a polytope PdP\subset\mathbb{R}^{d} consider the space

P:={q:V(GP)d}\mathcal{E}_{P}:=\{q:V(G_{P})\to\mathbb{R}^{d}\}

of dd-dimensional embeddings of GPG_{P}. We then have sk(P)P\operatorname{sk}(P)\in\mathcal{E}_{P}. Since Pn×d\mathcal{E}_{P}\cong\leavevmode\nobreak\ \mathbb{R}^{n\times d} (in some reasonable sense) we can pull back a metric μ:P×P+\mu\colon\mathcal{E}_{P}\times\mathcal{E}_{P}\to\mathbb{R}_{+}.

Theorem 4.5 (local version).

Given a polytope PdP\subset\mathbb{R}^{d} with 0int(P)0\in\operatorname{int}(P), there exists an ϵ>0\epsilon>0 with the following property: if q:V(GP)dq\colon V(G_{P})\to\mathbb{R}^{d} is an embedding with

  1. ()

    qq is ϵ\epsilon-close to sk(P)\operatorname{sk}(P), i.e., μ(q,sk(P))<ϵ\mu(q,\operatorname{sk}(P))<\epsilon,

  2. ()

    edges in qq are at most as long as in PP, and

  3. ()

    vertex-origin distances in qq are at least as large as in PP,

then qsk(P)q\simeq\operatorname{sk}(P).

Theorem˜4.5 as well can be naturally interpreted in the language of rigidity theory (see Section˜5.3). The proof below makes no use of this language.

In order to prove Theorem˜4.5 we again show 0ϕ(int(P))0\in\phi(\operatorname{int}(P)), which requries more work this time: since 0int(P)0\in\operatorname{int}(P), there exists an ϵ\epsilon-neighborhood Bϵ(0)PB_{\epsilon}(0)\subset P of the origin. The hope is that for a sufficiently small perturbation of the vertices of PP the image of Bϵ(0)B_{\epsilon}(0) under ϕ\phi is still a neighborhood of the origin.

This hope is formalized and verified in the following lemma, which we separated from the proof of Theorem˜4.5 to reuse it in Section˜4.4. Its proof is standard and is included in Appendix˜D:

Lemma D.1.

Let KdK\subset\mathbb{R}^{d} be a compact convex set with 0int(K)0\in\operatorname{int}(K) and f:K×[0,1]f\colon K\times[0,1] d\to\mathbb{R}^{d} a homotopy with f(,0)=idKf(\kern 0.70007pt\cdot\kern 0.70007pt,0)=\operatorname{id}_{K}. If the restriction f|K:K×[0,1]df|_{\partial K}\colon\partial K\times[0,1]\to\mathbb{R}^{d} yields a homotopy of K\partial K in d{0}\mathbb{R}^{d}\setminus\{0\}, then 0intf(K,1)0\in\operatorname{int}f(K,1).

In other words: if a “well-behaved” set KK contains the origin in its interior, and it is deformed so that its boundary never crosses the origin, then the origin stays inside.

Proof of Theorem˜4.5.

Since 0int(P)0\in\operatorname{int}(P) there exists a δ>0\delta>0 with Bδ(0)PB_{\delta}(0)\subset P.

Fix some compact neighborhood NPN\subset\mathcal{E}_{P} of sk(P)\operatorname{sk}(P). Then N×PN\times P is compact in P×d\mathcal{E}_{P}\times\mathbb{R}^{d} and the map

N×Pd,(q,x)ϕq(x):=iαi(x)qiN\times P\to\mathbb{R}^{d},\;\;(q,x)\mapsto\phi_{q}(x):=\sum_{i}\alpha_{i}(x)q_{i}

is uniformly continuous: there exists an ϵ>0\epsilon>0 so that whenever μ(q,q)+xx<ϵ\mu(q,q^{\prime})+\|x-x^{\prime}\|<\epsilon, we have ϕq(x)ϕq(x)<δ/2\|\phi_{q}(x)-\phi_{q^{\prime}}(x^{\prime})\|<\delta/2. We can assume that ϵ\epsilon is sufficiently small, so that Bϵ(sk(P))NB_{\epsilon}(\operatorname{sk}(P))\subset N. We show that this ϵ\epsilon satisfies the statement of the theorem.

Suppose that qq is ϵ\epsilon-close to sk(P)\operatorname{sk}(P), then

μ(sk(P),q)+xx<ϵxϕq(x)=ϕsk(P)(x)ϕq(x)<δ/2.\mu(\operatorname{sk}(P),q)+\|x-x\|<\epsilon\;\,\implies\;\,\|x-\phi_{q}(x)\|=\|\phi_{\operatorname{sk}(P)}(x)-\phi_{q}(x)\|<\delta/2.

The same is true when replacing qq by any linear interpolation q(t):=(1t)sk(P)+tqq(t):=(1-t)\operatorname{sk}(P)+tq with t[0,1]t\in[0,1]. Define the following homotopy:

f:Bδ(0)×[0,1]d,(x,t)ϕq(t)(x).f\colon B_{\delta}(0)\times[0,1]\to\mathbb{R}^{d},\quad(x,t)\mapsto\phi_{q(t)}(x).

We have f(,0)=idf(\kern 0.70007pt\cdot\kern 0.70007pt,0)=\operatorname{id}. That is, if xBδ(0)x\in\partial B_{\delta}(0) then f(x,0)=δ\|f(x,0)\|=\delta, as well as

f(x,t)f(x,0)f(x,0)f(x,t)δδ/2=δ/2,\|f(x,t)\|\geq\|f(x,0)\|-\|f(x,0)-f(x,t)\|\geq\delta-\delta/2=\delta/2,

and f(x,t)0f(x,t)\not=0 for all t[0,1]t\in[0,1]. Since Bδ(0)B_{\delta}(0) is compact and convex, the homotopy ff satisfies the conditions of Lemma˜D.1 and we conclude 0f(Bδ(0),1)=ϕq(Bδ(0))0\in f(B_{\delta}(0),1)=\phi_{q}(B_{\delta}(0)) ϕq(int(P))\subseteq\phi_{q}(\operatorname{int}(P)). Then qsk(P)q\simeq\operatorname{sk}(P) follows via Lemma˜4.3. ∎

Refer to caption
Figure 7. If qq in Theorem˜4.5 is not restricted to aff(P)\operatorname{aff}(P), the vertex-origin distances can be increased by moving out of the affine hull.

The polytope PP in Theorem˜4.5 is assumed to be full-dimensional. This is necessary, since allowing sk(P)\operatorname{sk}(P) to deform beyond its initial affine hull already permits counterexamples such as shown in Figure˜7. Even restricting to deformations with 0intconv(q)0\in\operatorname{int}\operatorname{conv}(q) is not sufficient, as shown in the next example:

Example 4.6.

Consider the 3-cube as embedded in 3×{0}3×25\mathbb{R}^{3}\times\{0\}\subset\mathbb{R}^{3}\times\mathbb{R}^{2}\cong\mathbb{R}^{5}. Let p1,,p83p_{1},...,p_{8}\in\mathbb{R}^{3} ×\times {0}\{0\} be its vertices, and let q1,,q8{0}×2q_{1},...,q_{8}\in\{0\}\times\mathbb{R}^{2} be the vertices as embedded in Figure˜6 (on the right). Since both share the same edge lengths and vertex-origin distances, so does the embedding tp+sqtp+sq whenever t2+s2=1t^{2}+s^{2}=1. Observe further that for both pip_{i} and qiq_{i} the origin can be written as a convex combination using the same coefficients αΔn\alpha\in\Delta_{n} (use an α\alpha with αi=αj\alpha_{i}=\alpha_{j} whenever pi=pjp_{i}=-p_{j}). It follows 00\in intconv(tp+sq)\operatorname{int}\operatorname{conv}(tp+sq).

4.4. Combinatorially equivalent polytopes

In this section we consider combinatorially equivalent polytopes P,QdP,Q\subset\mathbb{R}^{d} and prove the following:

Theorem 4.7 (combinatorial equivalent version).

Let P,QdP,Q\subset\mathbb{R}^{d} be combinatorially equivalent polytopes so that

  1. ()

    0int(Q)0\in\operatorname{int}(Q)

  2. ()

    edges in QQ are at most as long as in PP, and

  3. ()

    vertex-origin distances in QQ are at least as large as in PP.

Then PQP\cong Q.

In particular, since the combinatorics of polytopes up to dimension three is determined by the edge-graph, this proves ˜1.2 for d,e3d,e\leq 3.

Once again the proof uses Lemma˜4.3. Since 0int(Q)0\in\operatorname{int}(Q), we can verify 0ϕ(int(P))0\in\phi(\operatorname{int}(P)) by showing that the Wachspress map ϕ:PQ\phi\colon P\to Q is surjective. This statement is of independent interest, since the question whether ϕ\phi is bijective is a well-known open problem for d3d\geq 3 (see ˜5.11). Our proof of surjectivity uses Lemma˜D.1 and the following:

Lemma 4.8.

Given a face σ(P)\sigma\in\mathcal{F}(P), the Wachspress map ϕ\phi sends σ\sigma onto the corresponding face σQ(Q)\sigma_{Q}\in\mathcal{F}(Q). In particular, ϕ\phi sends P\partial P onto Q\partial Q.

Proof.

Given a point xrelint(σ)x\in\operatorname{relint}(\sigma) with Wachspress coordinates αΔn\alpha\in\Delta_{n}, the coefficient αi\alpha_{i} is non-zero if and only if the vertex pip_{i} is in σ\sigma (Remark˜3.7 (i)). The claim ϕ(x)σQ\phi(x)\in\sigma_{Q} follows immediately. ∎

Lemma 4.9.

The Wachspress map ϕ:PQ\phi\colon P\to Q is surjective.

Proof.

We proceed by induction on the dimension dd of PP. For d=1d=1 the Wachspress map is linear and the claim is trivially true. For d>1d>1 recall that ϕ\phi sends P\partial P to Q\partial Q (by Lemma˜4.8). By induction hypothesis, ϕ\phi is surjective on each proper face, thus surjective on all of P\partial P.

To show surjectivity in the interior, we fix xint(Q)x\in\operatorname{int}(Q); we show ximϕx\in\operatorname{im}\phi. Let ψ:QP\psi\colon Q\to P be the Wachspress map in the other direction (which is usually not the inverse of ϕ\phi) and define ρ:=ϕψ:QQ\rho:=\phi\circ\psi\colon Q\to Q. Note that by Lemma˜4.8 ρ\rho sends each face of QQ to itself and is therefore homotopic to the identity on QQ via the following linear homotopy:

f:Q×[0,1]Q,(y,t)(1t)x+tρ(x).f\colon Q\times[0,1]\to Q,\;(y,t)\mapsto(1-t)x+t\rho(x).

Since faces of QQ are closed under convex combination, f(,t)f(\kern 0.70007pt\cdot\kern 0.70007pt,t) sends Q\partial Q to itself for all t[0,1]t\in[0,1]. Thus, ff satisfies the assumptions of Lemma˜D.1 (with KK chosen as QQ), and therefore xf(Q,1)=imρimϕx\in f(Q,1)=\operatorname{im}\rho\subset\operatorname{im}\phi. ∎

The proof of Theorem˜4.7 follows immediately:

Proof of Theorem˜4.7.

Since 0int(Q)0\in\operatorname{int}(Q) and the Wachspress map ϕ:PQ\phi\colon P\to Q is surjective (by Lemma˜4.9), there exists xPx\in P with ϕ(x)=0\phi(x)=0. Since ϕ(P)=Q\phi(\partial P)=\partial Q (by Lemma˜4.8), we must have xint(P)x\in\operatorname{int}(P). PQP\simeq Q then follows via Lemma˜4.3. ∎

Corollary 4.10.

A polytope with the origin in its interior is uniquely determined by its face-lattice, its edge lengths and its vertex-origin distances.

If the origin lies not in PP then a unique reconstruction is not guaranteed (recall Figure˜3). However, if 0P0\in\partial P then we can say more. Recall the tangent cone of PP at a face σ(P)\sigma\in\mathcal{F}(P):

TP(σ):=cone{xyxP,yσ}.T_{P}(\sigma):=\operatorname{cone}\{x-y\mid x\in P,\,y\in\sigma\}.
Theorem 4.11.

Let P,QdP,Q\subset\mathbb{R}^{d} be combinatorially equivalent polytopes with the following properties:

  1. ()

    0relint(σQ)0\in\operatorname{relint}(\sigma_{Q}) for some face σQ(Q)\sigma_{Q}\in\mathcal{F}(Q), σP\sigma_{P} is the corresponding face in PP, and PP and QQ have isometric tangent cones at σP\sigma_{P} and σQ\sigma_{Q}.

  2. ()

    edges in QQ are at most as long as in PP.

  3. ()

    vertex-origin distances in QQ are at least as large as in PP.

Then PQP\cong Q.

Property (i) is always satisfied if, for example, σQ\sigma_{Q} is a facet of QQ, or if σQ\sigma_{Q} is a face of codimension two at which PP and QQ agree in the dihedral angle.

Proof.

The proof is by induction on the dimension dd of the polytopes. The induction base d=1d=1 is clearly satisfied. In the following we assume d2d\geq 2.

Note first that we can apply Theorem˜4.7 to σP\sigma_{P} and σQ\sigma_{Q} to obtain σPσQ\sigma_{P}\simeq\sigma_{Q} via an orthogonal transformation T:ddT\colon\mathbb{R}^{d}\to\mathbb{R}^{d}, in particular 0relint(σP)P0\in\operatorname{relint}(\sigma_{P})\subseteq P. By (i) this transformation extends to the tangent cones at these faces. Let F1,Q,,Fm,Qd1(Q)F_{1,Q},...,F_{m,Q}\in\mathcal{F}_{d-1}(Q) be the facets of QQ that contain σQ\sigma_{Q}, and let F1,P,,Fm,Pd1(P)F_{1,P},...,F_{m,P}\in\mathcal{F}_{d-1}(P) be the corresponding facets in PP. Then Fi,PF_{i,P} and Fi,QF_{i,Q} too have isometric tangent cones at σP\sigma_{P} resp. σQ\sigma_{Q}, and Fi,PFi,QF_{i,P}\simeq F_{i,Q} follows by induction hypothesis.

Now choose a point xQdx_{Q}\in\mathbb{R}^{d} beyond the face σQ\sigma_{Q} (i.e., above all facet-defining hyperplanes that contain σQ\sigma_{Q}, and below the others) so that xP:=TxQx_{P}:=Tx_{Q} is beyond the face σP\sigma_{P}. Consider the polytopes Q:=conv(Q{xQ})Q^{\prime}:=\operatorname{conv}(Q\cup\{x_{Q}\}) and P:=conv(P{xP})P^{\prime}:=\operatorname{conv}(P\cup\{x_{P}\}). Since xQx_{Q} lies beyond σQ\sigma_{Q}, each edge of QQ^{\prime} is either an edge of QQ, or is an edge between xQx_{Q} and a vertex of some facet Fi,QF_{i,Q}; and analogously for PP^{\prime}. The lengths of edges incident to xQx_{Q} depend only on the shape of the tangent cone and the shapes of the facets Fi,QF_{i,Q}, hence are the same for corresponding edges in PP^{\prime}. Thus, PP^{\prime} and QQ^{\prime} satisfy the preconditions of Theorem˜4.7, and we have PQP^{\prime}\simeq Q^{\prime}.

Finally, as xQ0x_{Q}\to 0, we have QQQ^{\prime}\to Q and PPP^{\prime}\to P (in the Hausdorff metric), which shows that PQP\simeq Q. ∎

Thus, if the origin lies in the interior of PP or a facet of PP then Theorem˜4.7 applies. If the origin lies in a face of codimension three, then counterexamples exist.

Example 4.12.

Consider the pentagons from Figure˜3 as lying in the plane 2×{1}\mathbb{R}^{2}\times\{1\}, with their former origins now at (0,0,1)(0,0,1). Consider the pyramids

P:=conv(P{0})andQ:=conv(Q{0}).P^{*}:=\operatorname{conv}(P\cup\{0\})\quad\text{and}\quad Q^{*}:=\operatorname{conv}(Q\cup\{0\}).

These polytopes have the origin in a vertex (a face of codimension three), have the same edge-graphs, edge lengths and vertex-origin distances, yet are not isometric. Examples with the origin in a high-dimensional face of codimension three can be constructed by considering prisms over PP^{*} resp. QQ^{*}.

We do not know whether Theorem˜4.7 holds if the origin lies in a face of codimension two (see ˜5.13).

4.5. Inscribed polytopes

It is worthwhile to formulate versions of Theorem˜4.7 for inscribed polytopes, that is, polytopes where all vertices lie on a common sphere – the circumsphere. For inscribed polytopes we can write down a direct monotone relation between edge lengths and the circumradius.

Corollary 4.13 (inscribed version).

Given two combinatorially equivalent polytopes P,QdP,Q\subset\mathbb{R}^{d} so that

  1. ()

    PP and QQ are inscribed in spheres of radii rPr_{P} and rQr_{Q} respectively,

  2. ()

    QQ contains its circumcenter in the interior, and

  3. ()

    edges in QQ are at most as long as in PP,

Then rPrQr_{P}\geq r_{Q}, with equality if and only if PQP\simeq Q.

Proof.

Translate PP and QQ so that both circumcenters lie at the origin. Suppose that rPrQr_{P}\leq r_{Q}. Then all preconditions of Theorem˜4.7 are satisfied, which yields PQP\simeq Q, hence rP=rQr_{P}=r_{Q}. ∎

This variant in particular has already found an application in proving the finitude of so-called “compact sphere packings” with spheres of only finitely many different radii [20].

Interestingly, the corresponding “unique reconstruction version” does not require any assumptions about the location of the origin or an explicit value for the circumradius. In fact, we do not even need to apply our results, as it already follows from Cauchy’s rigidity theorem (Theorem˜B.1).

Corollary 4.14.

An inscribed polytope of a fixed combinatorial type is uniquely determined, up to isometry, by its edge lengths.

Proof.

The case d=2d=2 is straightforward: given any circle, there is only a single way (up to isometry) to place edges of prescribed lengths. Also, there is only a single radius for the circle for which the edges reach around the circle exactly once and close up perfectly. This proves uniqueness for polygons.

If PP is of higher dimension then its 2-dimensional faces are still inscribed, have prescribed edge lengths, and by the 2-dimensional case above, corresponding 2-faces in PP and QQ are therefore isometric. Then PQP\simeq Q follows from Cauchy’s rigidity theorem (Theorem˜B.1). ∎

5. Conclusion, further notes and many open questions

We conjectured that a convex polytope is uniquely determined up to isometry by its edge-graph, edge lengths and the collection of distances between its vertices and some interior point, across all dimensions and combinatorial types (˜1.1). We also posed a more general conjecture expressing the idea that polytope skeleta, given their edge lengths, are maximally expanded (˜1.2). We developed techniques based on Wachspress coordinates and the so-called Izmestiev matrix that led to us to resolve three relevant special cases: centrally symmetric polytopes (Theorem˜4.4), small perturbations (Theorem˜4.5), and combinatorially equivalent polytopes (Theorem˜4.7). We feel confident that our approach already highlights the essential difficulties in verifying the general case.

In this section we collected further thoughts on our results, notes on connections to the literature, as well as many questions and future research directions.

5.1. Consequences of the conjectures

˜1.1 vastly generalizes several known “reconstruction from the edge-graph” results. The following is a special case of ˜1.1: an inscribed polytopes with all edges of the same length would be uniquely determined by its edge-graph. This includes the following special cases:

  • The reconstruction of matroids from their base exchange graph: a matroid can be identified with its matroid base polytopes, which is a 01-polytopes (hence inscribed) and has all edges of length 2\sqrt{2}. This reconstruction has been initially proven in [14] and recently rediscovered in [24].

  • The reconstruction of simultaneously vertex- and edge-transitive polytopes from their edge-graph: this was proven in [34, 36], essentially using the tools of this article.

It would imply an analogous reconstruction from the edge-graph for classes of polytopes such as the uniform polytopes or higher-dimensional inscribed Johnson solids [16].

Secondly, a positive answer to ˜1.1 would also resolve Question 6.6 in [35] on whether the metric coloring can capture the Euclidean symmetries of a polytope.

5.2. ˜1.2 for graph embeddings

In Example˜4.1 we show that ˜1.2 does not hold when replacing QQ by some more general graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} of GPG_{P}, even if 0intconv(q)0\in\operatorname{int}\operatorname{conv}(q).

Our intuition for why this fails, and also what distinguishes it from the setting of our conjectures and the verified special cases is, that the embedding of Example˜4.1 does not “wrap around the origin” properly. It is not quite clear what this means for an embedding of a graph, except that it feels right to assign this quality to polytope skeleta, to embeddings close to them, and also to centrally symmetric embeddings.

One possible formalization of this idea is expressed in the conjecture below, that is even stronger than ˜1.2 (the idea is due to Joseph Doolittle):

Conjecture 5.1.

Given a polytope PdP\subset\mathbb{R}^{d} and a graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e} of its edge-graph GPG_{P}, so that

  1. ()

    for each vertex iV(GP)i\in V(G_{P}) the cone

    Ci:=qi+cone{qjqiijE(GP)}C_{i}:=q_{i}+\operatorname{cone}\{q_{j}-q_{i}\mid ij\in E(G_{P})\}

    contains the origin in its interior,

  2. ()

    edges in qq are at most as long as in PP, and

  3. ()

    vertex-origin distances in qq are at least as large as in PP,

then sk(P)q\operatorname{sk}(P)\simeq q.

Note that since iCiconv(q)\bigcap_{i}C_{i}\subseteq\operatorname{conv}(q), (i) already implies 0intconv(q)0\in\operatorname{int}\operatorname{conv}(q).

5.3. Classical rigidity of frameworks

We previously commented on a natural interpretation of Theorem˜4.4 and Theorem˜4.5 in the language of classical rigidity theory (we refer to [6] for any rigidity specific terminology used below).

Consider the edges of PP as cables that can contract but not expand, and connect all vertices of PP to the origin using struts that can expand but not contract. This is known as a tensegrity framework, and we shall call it the tensegrity of PP. Theorem˜4.5 then asserts that these tensegrities are always (locally) rigid.

Using the language of rigidity, a number of natural follow up questions arise. So it turns out that swapping cables and struts does not necessarily preserve rigidity; see Figure˜8 for an example. As a consequence, the tensegrity of a polytope is not necessarily infinitesimally rigid, because infinitesimally rigid frameworks stay rigid under swapping cables and struts.

Refer to caption
Figure 8. Already the skeleton of the cube is not rigid if considered as a tensegrity framework with struts for edges and central cables. Twisting the top and bottom face lengthens the edge struts but keeps the central cables of a fixed length. This corresponds to the infinitesimal flex shown on the left.

Lacking first-order rigidity, we might ask for higher-order rigidity instead:

Question 5.2.

Is the tensegrity of a polytope always second-order rigid, or perhaps even prestress stable?

For an interpretation of Theorem˜4.4 as a tensegrity framework, consider a cable at each edge as before, but each central strut now connects a vertex pip_{i} to its antipodal counterpart pi-p_{i}, and is fixed in its center to the origin. Theorem˜4.4 then asserts that this tensegrity framework is universally rigid, i.e., it has a unique realization across all dimensions.

Here too, swapping cables and struts does not preserve universal or even global rigidity (see Figure˜9). It does not preserve local rigidity either (see Example˜5.3).

Refer to caption
Figure 9. An octagon and an embedding of its edge-graph with longer edges but equally long central cables, showing that the respective tensegrity framework is not globally rigid under forced central symmetry.
Example 5.3.

Consider the 4-cube with its “top” and “bottom” facets (which are 3-cubes) embedded in the hyperplanes 3×{±1}\mathbb{R}^{3}\times\{\pm 1\} respectively. We flex the skeleton as follows: deform the top facet as shown in Figure˜8, and the bottom facet so as to keep the framework centrally symmetric, while keeping both inside their respective hyperplanes. The edge struts inside the facets become longer, and the edge struts between the facets have previously been of minimal length between the hyperplanes, can therefore also only increase in length. The lengths of the central cables stay the same.

As a consequence, the centrally symmetric tensegrity frameworks too are not necessarily infinitesimally rigid.

5.4. Schlegel diagrams

Yet another interpretation of the frameworks discussed in Section˜5.3 is as skeleta of special Schlegel diagrams, namely, of pyramids whose base facet is the polytope PP. It is then natural to ask whether a general Schlegel diagram is rigid as well (this was brought up by Raman Sanyal).

Refer to caption
Figure 10. Two Schlegel diagrams of 4-polytopes: of the pyramid with the 3-cube as base facet (left) and of the 4-cube (right).

The question of rigidity for Schlegel diagrams is already interesting in dimension two, that is, for Schlegel diagrams of 3-polytopes. The edge-graphs of many 3-polytopes are too sparse to be generically rigid in 2\mathbb{R}^{2}, and so one might expect that most of their Schlegel diagrams are flexible. Indeed, flexible Schlegel diagrams exist (see Figure˜11, left).

Refer to caption
Figure 11. A flexible Schlegel diagram (left), and a rigid Schlegel diagram (right).

Surprisingly however, this seems to be the exception rather than the rule. For example, we believe that Schlegel diagrams of (2n+1)(2n+1)-gonal prisms are always rigid (see Figure˜11, right). Since Schlegel diagrams are very special realizations (they are projections of convex objects), the generic ones among them might very well be rigid. This is not clear so far.

Question 5.4.

Is a generic Schlegel diagram rigid?

Above we considered Schlegel diagrams as bar-joint frameworks. If we consider them as tensegrity frameworks then it is easy to find generically flexible examples (see Figure˜12). Schlegel diagrams are also not necessarily globally rigid (see Figure˜13).

Refer to caption
Figure 12. The skeleton of Schlegel diagram of a triangular prism with cables on the outside and struts on the inside is not rigid. Twisting the inner triangle increases the lengths of struts and fixes all other lengths.
Refer to caption
Figure 13. Folding the Schlegel diagram of the 3-cube along a diagonal preserves all edge lengths.

5.5. Stoker’s conjecture

Stoker’s conjecture asks whether the dihedral angles of a polytope determine its face angles, and thereby its overall shape to some degree. Recall that dihedral angles are the angle at which facets meet in faces of codimension two, whereas face angles are the dihedral angles of the facets. Stoker’s conjecture was asked in 1968 [26], and a proof was claimed recently by Wang and Xie [29]:

Theorem 5.5 (Wang-Xie, 2022).

Let P1P_{1} and P2P_{2} be two combinatorially equivalent polytopes such that corresponding dihedral angles are equal. Then all corresponding face angles are equal as well.

Our results allow us to formulate a semantically similar statement. The following is a direct consequence of Corollary˜4.10 when expressed for the polar dual polytope:

Corollary 5.6.

Let P1P_{1} and P2P_{2} be two combinatorially equivalent polytopes such that corresponding dihedral angles and facet-origin distances are equal. Then P1P2P_{1}\simeq P_{2}.

While the assumptions in Corollary˜5.6 are unlike stronger compared to Stoker’s conjecture (we require facet-origin distances), we also obtain isometry instead of just identical face angles. While related, we are not aware that either of Theorem˜5.5 or Corollary˜5.6 follows from the other one easily.

5.6. Pure edge length constraints

Many polytopes cannot be reconstructed up to isometry from their edge-graph and edge lengths alone (recall Figure˜1). However, for all we know the following is open:

Question 5.7.

Is the combinatorial type of a polytope uniquely determined by its edge-graph and edge lengths?

This alone would already prove ˜1.1 (by Corollary˜4.10). It would also imply a positive answer to the following:

Question 5.8.

Is a polytope uniquely determined up to isometry by its 2-skeleton (i.e., the face-lattice cut off at, but including dimension two) and the shape of each 2-face?

Note that this is a particular strengthening of Cauchy’s rigidity theorem, which requires the face-lattice to be prescribed in its entirety, rather than on some lower levels only.

Let us now fix the combinatorial type. We are aware of three types of polytopes that are not determined (up to isometry) by their face-lattice and edge lengths:

  1. ()

    nn-gons with n4n\geq 4.

  2. ()

    Minkowski sums: if P=Q+RP=Q+R and QQ and RR are generically oriented w.r.t. each other, then a slight reorientation of the summands changes the shape of PP but keeps its edge lengths (see Figure˜14).

  3. ()

    polytopes having all edge directions on a “conic at infinity”: this implies an affine flex (see [5]). This is most easily implemented for zonotopes (recall Figure˜5), but happens for other polytopes as well, such as 3-polytopes with up to five edge directions (see Figure˜15).

Refer to caption
Figure 14. The cuboctahedron can be written as the Minkowski sum of two simplices, and twisting these simplices leads to a flex of the cuboctahedron that preserves edge lengths.
Refer to caption
Figure 15. A cuboid sliced at an angle in an appropriate way has only five edge directions and has an edge length preserving affine flex. The flex deforms the bottom face into a rhombus and keeps the vertical edges vertical.

We are not aware of other examples of polytopes that flex in this way and so we wonder whether this is already a full characterization.

Question 5.9.

If a polytope is not determined up to isometry by its combinatorial type and edge lengths, is it necessarily a polygon, a non-trivial Minkowski sum or has all its edge directions on a conic at infinity? Is this true at least up to dimension three?

In how far a 3-polytope is determined by local metric data at its edges was reportedly discussed in an Oberwolfach question session (as communicated by Ivan Izmestiev on MathOverflow [33]), where the following more general question was asked:

Question 5.10.

Given a simplicial 3-polytope and at each edge we prescribe either the length or the dihedral angle, in how far does this determine the polytope?

Having length constraints at every edge determines a simplicial polytope already up to isometry via Cauchy’s rigidity theorem (Theorem˜B.1). The angles-only version is exactly the 3-dimensional Stoker’s conjecture (Section˜5.5). We are not aware that this question has been addressed in the literature beyond these two extreme cases.

Note also that ˜5.10 is stated for simplicial 3-polytopes, but actually includes general 3-polytopes via a trick: if PP is not simplicial, triangulate every 2-face, and at each new edge created in this way prescribe a dihedral angle of 180180^{\circ} to prevent the faces from folding at it.

5.7. Injectivity of the Wachspress map ϕ\boldsymbol{\phi}

In Lemma˜4.9 we proved that the Wachspress map ϕ:PQ\phi\colon P\to Q (cf. Definition˜4.2) between combinatorially equivalent polytopes is surjective. In contrast, the injectivity of the Wachspress map has been established only in dimension two by Floater and Kosinka [10] and is conjectured for all d3d\geq 3.

Conjecture 5.11.

The Wachspress map ϕ:PQ\phi\colon P\to Q is injective.

If true, the Wachspress map would provide an interesting and somewhat canonical homeomorphism (in fact, a rational map, see [31]) between any two combinatorially equivalent polytopes.

5.8. What if 𝟎𝐢𝐧𝐭(𝑸)\boldsymbol{0\not\in\mathrm{int}(Q)}?

If 0Q0\not\in Q then Figures˜4 and 3 show that our conjectures fail. We do however not know whether in the “unique reconstruction” case the number of solutions would be finite.

Question 5.12.

Given edge-graph, edge lengths and vertex-origin distances, are there only finitely many polytopes with these parameters?

This is in contrast to when we replace QQ with a graph embedding q:V(GP)eq\colon V(G_{P})\to\mathbb{R}^{e}, which can have a continuum of realizations (see Figure˜16).

Refer to caption
Figure 16. The square based pyramid (left) is flexible as a framework (since then the bottom face needs not stay flat). Likewise, the framework of the square based frustum with this particular choice of origin (right) flexes. It is however (locally) rigid as a polytope.

In Section˜4.4 we showed that reconstruction from the face-lattice, edge lengths and vertex-origin distances is possible even if the origin lies only in the inside of a facet of PP, but that it can fail if it lies in a face of codimension three. We do not know what happens for a face of codimension two.

Question 5.13.

Is a polytope uniquely determined by its face-lattice, edge lengths and vertex-origin distances if the origin is allowed to lie in the inside of faces of codimension 0,10,1 and 22?

Appendix A Perron-Frobenius theory

The following fragment of the Perron-Frobenius theorem is relevant to this article. Recall that a matrix is irreducible if no simultaneous row-column permutation brings it in a block-diagonal form with more than one block; or equivalently, if it is not the (weighted) adjacency matrix of a disconnected graph. See also [11].

Theorem A.1 (Perron-Frobenius).

Let Mn×nM\in\mathbb{R}^{n\times n} be a non-negative irreducible symmetric matrix, then

  1. ()

    the largest eigenvalue θ\theta of MM is positive and has multiplicity one.

  2. ()

    there is a θ\theta-eigenvector znz\in\mathbb{R}^{n} with strictly positive entries.

Appendix B Cauchy’s rigidity theorem

Cauchy’s famous rigidity theorem was initially formulated in dimension three and is often quoted briefly as follows:

3-polytopes with isometric faces are themselves isometric.

Generalizations to higher dimensions have been proven by Alexandrov [1] where one assumes isometric facets to conclude global isometry (see also its proof in [23]). We state a rigorous version that only requires isometric 2-faces and that can be easily derived from the facet versions using induction by dimension:

Theorem B.1 (Cauchy’s rigidity theorem, version with 2-faces).

Given two combinatorially equivalent polytopes P,QdP,Q\subset\mathbb{R}^{d} and a face-lattice isomorphism ϕ:(P)(Q)\phi\colon\mathcal{F}(P)\to\mathcal{F}(Q). If ϕ\phi extends to an isometry on every 2-face σ2(P)\sigma\in\mathcal{F}_{2}(P), then ϕ\phi extends to an isometry on all of PP, that is, PQP\simeq Q.

Appendix C Some linear algebra

Theorem C.1.

Given two matrices Ad×nA\in\mathbb{R}^{d\times n} and Bd×mB\in\mathbb{R}^{d\times m} with spanBspanA\operatorname{span}B\subseteq\operatorname{span}A, there exists a linear transformation T:nmT\colon\mathbb{R}^{n}\to\mathbb{R}^{m} with TA=BTA^{\top}=B^{\top}.

Proof.

Set UA:=spanAU_{A}:=\operatorname{span}A and dA:=dimUA=rankAd_{A}:=\dim U_{A}=\operatorname{rank}A. Respectively, set UB:=spanBUAU_{B}:=\operatorname{span}B\subset U_{A} and dB:=dimUB=rankBd_{B}:=\dim U_{B}=\operatorname{rank}B. We can assume that the columns of AA and BB are sorted so that a1,,adAUAa_{1},...,a_{d_{A}}\in U_{A} form a basis, and likewise b1,,bdBUBb_{1},...,b_{d_{B}}\in U_{B} form a basis. Let T~:UAUB\tilde{T}\colon U_{A}\to U_{B} be the uniquely determined linear map that maps T~ai=bi\tilde{T}a_{i}=b_{i} for i{1,,dB}i\in\{1,...,d_{B}\} and T~ai=0\tilde{T}a_{i}=0 for i{dB+1,,dA}i\in\{d_{B}+1,...,d_{A}\}. Then T~A=B\tilde{T}A=B.

The Moore-Penrose pseudo inverse An×dA^{\dagger}\in\mathbb{R}^{n\times d} of AA satisfies AA=πUAAA^{\dagger}=\pi_{U_{A}}, where πUA\pi_{U_{A}} is the orthogonal projection onto UAU_{A}. We set T:=(AT~A)T:=(A^{\dagger}\tilde{T}A)^{\top} and verify

TA=(AT)=(AAT~A)=(πUAT~A)=(πUAB)=B,TA^{\top}=(AT^{\top})^{\top}=(AA^{\dagger}\tilde{T}A)^{\top}=(\pi_{U_{A}}\tilde{T}A)^{\top}=(\pi_{U_{A}}B)^{\top}=B^{\top},

where for the last equality we used that all columns of BB are already in UAU_{A} and the projection acts as identity. ∎

Appendix D A topological argument

Lemma D.1.

Let KdK\subset\mathbb{R}^{d} be a compact convex set, xint(K)x\in\operatorname{int}(K) a point and f:K×[0,1]df\colon K\times[0,1]\to\mathbb{R}^{d} a homotopy with f(,0)=idKf(\kern 0.70007pt\cdot\kern 0.70007pt,0)=\operatorname{id}_{K}. If the restriction f|K:K×[0,1]df|_{\partial K}\colon\partial K\times[0,1]\to\mathbb{R}^{d} yields a homotopy of K\partial K in d{x}\mathbb{R}^{d}\setminus\{x\}, then xintf(K,1)x\in\operatorname{int}f(K,1).

Proof.

Suppose that xintf(K,1)x\not\in\operatorname{int}f(K,1). Since f\partial f is a homotopy in d{x}\mathbb{R}^{d}\setminus\{x\}, we actually have xf(K,1)x\not\in f(K,1). We derive a contradiction.

Construct a map g:KKg\colon K\to\partial K as follows: for yKy\in K consider the unique ray emanating from xx passing through f(y,1)f(y,1). Let g(x)g(x) be the unique intersection of this ray with K\partial K. Likewise, construct the map h:K×[0,1]Kh\colon\partial K\times[0,1]\to\partial K: for yKy\in\partial K and t[0,1]t\in[0,1], let h(y,t)h(y,t) be the intersection of K\partial K with the unique ray emanating from xx and passing through f(y,t)f(y,t). Note that h(,0)=idKh(\kern 0.70007pt\cdot\kern 0.70007pt,0)=\operatorname{id}_{\partial K} and h(,1)=g|Kh(\kern 0.70007pt\cdot\kern 0.70007pt,1)=g|_{\partial K}. In other words, g|Kg|_{\partial K} is homotopic to the identity on K\partial K.

The existence of such a map g:KKg\colon K\to\partial K is a well-known impossibility. This can be quickly shown by considering the following commutative diagram (left) and the diagram induced on the \mathbb{Z}-homology groups (right):

K{K}K{{\partial K}}K{{\partial K}}g|K\scriptstyle{g|_{\partial K}}i\scriptstyle{i}g\scriptstyle{g}   H(K){H_{\bullet}(K)}H(K){{H_{\bullet}(\partial K)}}H(K){{H_{\bullet}(\partial K)}}(g|K)\scriptstyle{(g|_{\partial K})_{*}}i\scriptstyle{i_{*}}g\scriptstyle{g_{*}}

Since g|Kg|_{\partial K} is homotopic to the identity, the arrow (g|K)(g|_{\partial K})_{*} is an isomorphism, and so must be the arrows above it. This is impossible because

Hd1(K)=0=Hd1(K).H_{d-1}(\partial K)=\mathbb{Z}\not=0=H_{d-1}(K).

Appendix E Euler’s homogeneous function theorem

Theorem E.1 (Euler’s homogeneous function theorem).

Let f:nf\colon\mathbb{R}^{n}\to\mathbb{R} be a homogeneous function of degree d1d\geq 1 i.e., f(t𝐱)=tdf(𝐱)f(t\mathbf{x})=t^{d}f(\mathbf{x}) for all t0t\geq 0. Then

ixif(𝐱)xi=df(𝐱).\sum_{i}x_{i}\frac{\partial f(\mathbf{x})}{\partial x_{i}}=d\cdot f(\mathbf{x}).
Proof.

Differentiate both sides of f(t𝐱)=tdf(𝐱)f(t\mathbf{x})=t^{d}f(\mathbf{x}) w.r.t. tt

dtd1f(𝐱)=t(tdf(𝐱))=tf(t𝐱)=if(t𝐱)(txi)(txi)t=ixif(t𝐱)(txi)dt^{d-1}f(\mathbf{x})=\frac{\partial}{\partial t}(t^{d}f(\mathbf{x}))=\frac{\partial}{\partial t}f(t\mathbf{x})=\sum_{i}\frac{\partial f(t\mathbf{x})}{\partial(tx_{i})}\frac{\partial(tx_{i})}{\partial t}=\sum_{i}x_{i}\frac{\partial f(t\mathbf{x})}{\partial(tx_{i})}

and evaluate at t=1t=1. ∎

Appendix F An alternative proof of Theorem˜3.2 using semi-definite optimization

The following proof of Theorem˜3.2 does not address the equality case.

Proof.

Theorem˜3.2 can be equivalently phrased as the claim that the following program attains its optimum if we choose qiq_{i} to be the skeleton of PP:

maxqαs.t.qiqjpipj,for all ijEq1,,qnn\begin{array}[]{rl}\max&\|q\|_{\alpha}\\ \text{s.t\hbox to0.0pt{.\hss}}&\|q_{i}-q_{j}\|\leq\|p_{i}-p_{j}\|,\quad\text{for all $ij\in E$}\\ &q_{1},...,q_{n}\in\mathbb{R}^{n}\end{array}

Since qα2=12i,jαiαjqiqj2\|q\|_{\alpha}^{2}=\frac{1}{2}\sum_{i,j}\alpha_{i}\alpha_{j}\|q_{i}-q_{j}\|^{2} =iαiqi2iαiqi2=\sum_{i}\alpha_{i}\|q_{i}\|^{2}-\|\!\sum_{i}\alpha_{i}q_{i}\|^{2}, we obtain the following equivalent program:

maxiαiqi2=:e(q)s.t.iαiqi=0qiqjpipj,for all ijEq1,,qnn\begin{array}[]{rl}\max&\sum_{i}\alpha_{i}\|q_{i}\|^{2}=:e(q)\\ \text{s.t\hbox to0.0pt{.\hss}}&\sum_{i}\alpha_{i}q_{i}=0\\ &\|q_{i}-q_{j}\|\leq\|p_{i}-p_{j}\|,\quad\text{for all $ij\in E$}\\ &q_{1},...,q_{n}\in\mathbb{R}^{n}\end{array}

This particular program has been studied extensively (see e.g. [13, 12, 27]). It can be rewritten as a semi-definite program (which we do not repeat here) with the following dual:

minijEwijpipj2=:d(w)s.t.Lwdiag(α)+μαα0w0,μ free\begin{array}[]{rl}\min&\sum_{ij\in E}w_{ij}\|p_{i}-p_{j}\|^{2}=:d(w)\\ \text{s.t\hbox to0.0pt{.\hss}}&L_{w}-\operatorname{diag}(\alpha)+\mu\alpha\alpha^{\top}\succeq 0\\ &w\geq 0,\mu\text{ free}\end{array}

where LwL_{w} is the Laplace matrix of GPG_{P} with edge weights ww (that is Lij=wijL_{ij}=-w_{ij} and Lii=jiwijL_{ii}=\sum_{j\not=i}w_{ij}), diag(α)\operatorname{diag}(\alpha) is the diagonal matrix with α\alpha on its diagonal, and X0X\succeq 0 asserts that XX is a positive semi-definite matrix.

Recall the following property of a dual program: if there are qinq_{i}\in\mathbb{R}^{n}, w0w\geq 0 and μ\mu\in\mathbb{R} so that the primal and the dual program attain the same objective value, then we know that there is no duality gap and we found optimal solutions for both programs. We now claim that such a choice can be made using qi:=piq_{i}:=p_{i}, wij:=Mijw_{ij}:=M_{ij} (where MM is the Izmestiev matrix of PP), and with a value for μ\mu to be determined later. We first verify that the objective values agree:

d(M)\displaystyle d(M) =ijEMijpipj2=12i,jMijpipj2\displaystyle=\sum_{ij\in E}M_{ij}\|p_{i}-p_{j}\|^{2}=\tfrac{1}{2}\sum_{i,j}M_{ij}\|p_{i}-p_{j}\|^{2}
=i(jMij)pi2i,jMijpi,pj\displaystyle=\sum_{i}\Big{(}\sum_{j}M_{ij}\Big{)}\|p_{i}\|^{2}-\sum_{i,j}M_{ij}\langle p_{i},p_{j}\rangle
=iαipi2tr(MXP=0XP)\displaystyle=\sum_{i}\alpha_{i}\|p_{i}\|^{2}-\operatorname{tr}(\,\underbrace{MX_{P}}_{=0}X_{P}^{\top})
=iαipi2=e(p),\displaystyle=\sum_{i}\alpha_{i}\|p_{i}\|^{2}=e(p),

where XP:=(p1,,pn)n×dX_{P}^{\top}:=(p_{1},...,p_{n})\in\mathbb{R}^{n\times d}, MXP=0MX_{P}=0 by Theorem˜3.3 (iv), as well as iMij\sum_{i}M_{ij} =αj=\alpha_{j} by Corollary˜3.6.

It only remains to verify that there exists μ\mu\in\mathbb{R} so that LMdiag(α)+μαα0L_{M}-\operatorname{diag}(\alpha)+\mu\alpha\alpha^{\top}\succeq 0. Set D:=diag(α11/2,,αn1/2)D:=\operatorname{diag}(\alpha_{1}^{\smash{-1/2}},...,\alpha_{n}^{\smash{-1/2}}) and observe that the matrices XX and DXDDXD have the same signature. It therefore suffices to verify

0D(LMdiag(α)+μαα)D=DLMDId+μ(D1𝟏)(D1𝟏).0\preceq D(L_{M}-\operatorname{diag}(\alpha)+\mu\alpha\alpha^{\top})D=DL_{M}D-\operatorname{Id}+\mu(D^{-1}\mathbf{1})(D^{-1}\mathbf{1})^{\top}.

First we claim that LMdiag(α)=ML_{M}-\operatorname{diag}(\alpha)=-M. Since both sides agree on the off-diagonal, it suffices to compare their row sums. And in fact, since LM𝟏=0L_{M}\mathbf{1}=0 we have (LMdiag(α))𝟏=α=M𝟏(L_{M}-\operatorname{diag}(\alpha))\mathbf{1}=-\alpha=-M\mathbf{1}. Hence, DLMDIdDL_{M}D-\operatorname{Id} has the same signature as M-M, i.e., a unique negative eigenvalue, and one can check that the corresponding eigenvector is D1𝟏D^{-1}\mathbf{1}. We see that the term μ(D1𝟏)(D1𝟏)\mu(D^{-1}\mathbf{1})(D^{-1}\mathbf{1})^{\top} just shifts this smallest eigenvalue of DLMDIdDL_{M}D-\operatorname{Id} up or down, while not changing the other eigenvalues, and so we can choose μ\mu large enough to make this eigenvalue positive. ∎

The formulation of Theorem˜3.2 as a semi-definite program allows for a simultaneous reconstruction (cf.  Corollary˜3.8) of both the polytope and its Izmestiev matrix from only the edge-graph, the edge lengths and the Wachspress coordinates of some interior point. Since semi-definite programs can be solved in polynomial time, this approach is actually feasible in practice.


Funding. This work was supported by the Engineering and Physical Sciences Research Council [EP/V009044/1]


Acknowledgements. I thank Raman Sanyal, Joseph Doolittle, Miek Messerschmidt, Bernd Schulze, James Cruickshank, Robert Connelly and Albert Zhang for many fruitful discussions on the topic of this article, many of which lead to completely new perspectives on the results and to numerous subsequent questions.

References

  • [1] A. D. Alexandrov. Convex polyhedra, volume 109. Springer, 2005.
  • [2] R. Blind and P. Mani-Levitska. Puzzles and polytope isomorphisms. Aequationes mathematicae, 34:287–297, 1987.
  • [3] R. Connelly. Generic global rigidity. Discrete & Computational Geometry, 33:549–563, 2005.
  • [4] R. Connelly. Stress matrices and M matrices. Oberwolfach Reports, 3:678–680, 2006.
  • [5] R. Connelly, S. J. Gortler, and L. Theran. Affine rigidity and conics at infinity. International Mathematics Research Notices, 2018(13):4084–4102, 2018.
  • [6] R. Connelly and W. Whiteley. Second-order rigidity and prestress stability for tensegrity frameworks. SIAM Journal on Discrete Mathematics, 9(3):453–491, 1996.
  • [7] J. Doolittle. Answer to “Given the skeleton of an inscribed polytope. If I move the vertices so that no edge increases in length, can the circumradius still get larger?” https://mathoverflow.net/a/419107/108884. version: 2022-03-28.
  • [8] J. Doolittle. Reconstructing nearly simple polytopes from their graph. arXiv preprint arXiv:1701.08334, 2017.
  • [9] M. S. Floater. Generalized barycentric coordinates and applications. Acta Numerica, 24:161–214, 2015.
  • [10] M. S. Floater and J. Kosinka. On the injectivity of wachspress and mean value mappings between convex polygons. Advances in Computational Mathematics, 32(2):163–174, 2010.
  • [11] G. Frobenius, F. G. Frobenius, F. G. Frobenius, F. G. Frobenius, and G. Mathematician. Über matrizen aus nicht negativen elementen. 1912.
  • [12] F. Göring, C. Helmberg, and M. Wappler. Embedded in the shadow of the separator. SIAM Journal on Optimization, 19(1):472–501, 2008.
  • [13] F. Göring, C. Helmberg, and M. Wappler. The rotational dimension of a graph. Journal of Graph Theory, 66(4):283–302, 2011.
  • [14] C. A. Holzmann, P. Norton, and M. Tobey. A graphical representation of matroids. SIAM Journal on Applied Mathematics, 25(4):618–627, 1973.
  • [15] I. Izmestiev. The colin de verdiere number and graphs of polytopes. Israel Journal of Mathematics, 178(1):427–444, 2010.
  • [16] N. W. Johnson. Convex polyhedra with regular faces. Canadian Journal of Mathematics, 18:169–200, 1966.
  • [17] M. Joswig and G. M. Ziegler. Neighborly cubical polytopes. Discrete & Computational Geometry, 24:325–344, 2000.
  • [18] T. Ju, S. Schaefer, J. D. Warren, and M. Desbrun. A geometric construction of coordinates for convex polyhedra using polar duals. In Symposium on Geometry Processing, pages 181–186, 2005.
  • [19] G. Kalai. A simple way to tell a simple polytope from its graph. Journal of combinatorial theory, Series A, 49(2):381–383, 1988.
  • [20] E. Kikianty and M. Messerschmidt. On compact packings of euclidean space with spheres of finitely many sizes. arXiv preprint arXiv:2305.00758, 2023.
  • [21] L. Lovász. Steinitz representations of polyhedra and the colin de verdiere number. Journal of Combinatorial Theory, Series B, 82(2):223–236, 2001.
  • [22] H. Narayanan, R. Shah, and N. Srivastava. A spectral approach to polytope diameter. arXiv preprint arXiv:2101.12198, 2021.
  • [23] I. Pak. Lectures on discrete and polyhedral geometry. Manuscript (http://www. math. ucla. edu/˜ pak/book. htm), 2010.
  • [24] G. Pineda-Villavicencio and B. Schröter. Reconstructibility of matroid polytopes. SIAM Journal on Discrete Mathematics, 36(1):490–508, 2022.
  • [25] D. E. Speyer. Does the edge-graph of a centrally symmetric polytope determine which vertices are antipodal? https://mathoverflow.net/q/440534/108884. version: 2023-02-27.
  • [26] J. J. Stoker. Geometrical problems concerning polyhedra in the large. Communications on pure and applied mathematics, 21(2):119–168, 1968.
  • [27] J. Sun, S. Boyd, L. Xiao, and P. Diaconis. The fastest mixing markov process on a graph and a connection to a maximum variance unfolding problem. SIAM review, 48(4):681–699, 2006.
  • [28] E. L. Wachspress. A rational finite element basis. 1975.
  • [29] J. Wang and Z. Xie. On gromov’s dihedral ridigidity conjecture and stoker’s conjecture. arXiv preprint arXiv:2203.09511, 2022.
  • [30] J. Warren. Barycentric coordinates for convex polytopes. Advances in Computational Mathematics, 6(1):97–108, 1996.
  • [31] J. Warren. On the uniqueness of barycentric coordinates. Contemporary Mathematics, 334:93–100, 2003.
  • [32] J. Warren, S. Schaefer, A. N. Hirani, and M. Desbrun. Barycentric coordinates for convex sets. Advances in computational mathematics, 27(3):319–338, 2007.
  • [33] M. Winter. Is the dodecahedron flexible (as a polytope with fixed edge-lengths)? https://mathoverflow.net/q/434771. version: 2022-11-17.
  • [34] M. Winter. Symmetric and spectral realizations of highly symmetric graphs, 2020.
  • [35] M. Winter. Capturing polytopal symmetries by coloring the edge-graph. arXiv preprint arXiv:2108.13483, 2021.
  • [36] M. Winter. Spectral Realizations of Symmetric Graphs, Spectral Polytopes and Edge-Transitivity. PhD thesis, Technische Universität Chemnitz, 2021.