ON STABILITY OF POLYTOPES
BY Wendy Finbow
SUBSIITTED IX P-iFtïL-U. FCZFILLMEST OF THE REQLXRE-UIEXTS FOR THE DEGREE OF
@ Copyright by Wendy Finbon*. 1998 The author has granted a non- L'auteur a accordé une licence non exclusive licence allowing the exclusive permettant à la Natio~liflLiirary of Canada to Bibliothèque nationale du Canada de repfoduce, 10- distriMe or seil reproduire, prêter, distribuer ou copies of this thesis in microform, vendre des copies de cette thèse sous papex or electronic formats. la forme de microfiche/£ïh, de reproduction sur papier ou sur format électronique.
The author retains ownership of the L'auteur consewe la propriété du copyright in this thesis. Neidier the droit d'auteur qui protège cette thèse. thesis nor substantial extracts from it Ni la thèse ni des extraits substantiels may be printed or othefwise de celle-ci ne doivent être imprimés repfoduced doutthe auihor's ou autrement reprodrrits sans son permission. airtorisation.
Canada Contents
List of Tables
List of Figures vii
Acknowledgements X
Abstract xi
Abbreviations and Symbols Used xii
1 Introduction 1
2 Preliminary Results 2.1 Minkowski's Theorem ...... 2.2 Further Results ......
3 Monostatic Simplexes in 8, 9, and 10 Dimensions 24 3.1 The Search Method ...... 25 3.2 Cornputer Algorithm ...... 26 3.3 Error .4nalysis...... 33 3.4 The Main Results ...... 37
4 Loaded Polytopes 44 4.1 Preliminary Results ...... 44 4.2 Nonloadable Polytopes ...... 47 4.3 Loadable Polytopes ...... 52 4.4 Related Results ...... 89
Bibliograp hy 91 List of Tables
3.1 Maximum contribution of each simplex with five facets ...... 30 3.2 Maximumcontributionofeachsimplexwithsixfacets ...... 31 3.3 Maximum contribution of each simplex with seven facets ...... 32 3.4 Mawimum contribution of each sirnplex with eight facets ...... 33 3.4 Maximum contribution of each simplex with eight facets, continued. . 34 3.5 Maximum contribution of trees with 5, 6. 7. 8. and 9 nodes ...... 38 3.6 Contributions of vectors ...... 39
4.1 Loadable/Nonloadable Polytopes ...... 48 4.1 Loadable/rYonloadable Polytopes. continued ...... 49 List of Figures
2.1 Illustration of the Projection Criterion (tetrahedron) ......
3.1 Example of a straight tree ...... 3.2 Example of a tree with six nodes ......
4.1 Overview of a truncated tetrahedron ...... 4.2 Overlap of opposite triangular facets (icosahedron) ...... 4.3 Overlapping facets of the icosahedron ...... 4.4 Overlap of opposite pentagonal facets (dodecahedron)...... 4.5 Overlapping facets of the dodecahedron ...... 4.6 Overlapping facets of the tmncated icosahedron ...... 4.7 Overlapping facets of the truncated dodecahedron ...... 4.8 Overlapping facets of the cuboctahedron ...... 4.9 Overlap of a pentagonal facet onto a triangular facet of the icosidodec- ahedron ...... 4.10 Overlap of a triangular facet onto a pentagonal facet of the icosidodec- ahedron ...... 1.11 Overlap of square facets of the rhornbicoçidodecahedron ...... 4.12 Overlap of a square facet onto a pentagonal facet of the rhombicosido- decahedron ...... 4.13 Overlap of square facets onto triangular facets of a rhombicosidodeca- hedron ......
vii 4.14 Overlap of triangular facets onto square facets of the rhornbicosidodec- ahedron ...... 4.13 Overlap of pentagonal facets onto square facets of the rhombicosido- decahedron ...... 4.16 Overlapping facets of the rhombicosidodecahedron ...... 4-17 Overlapping square facets of a snub cube ...... 4-18 Overlap of square facets onto triangular facets of the snub cube ... 4.19 Overlap of triangular facets of the snub cube and snub dodecahedron 4.20 Overlap of triangular facets onto square facets of the snub cube ... 4.21 Overlapping facets of the snub cube ...... 4.22 Overlapping pentagonal facets of the snub dodecahedron ...... 4.23 Overlap of the pentagonal facets of the snub dodecahedron ...... 4.24 Overlap of triangular facets onto pentagonal facets of the snub dodec- ahedron ...... 4.25 Overlapping facets of the snub dodecahedron ...... I dedzcate my thesis to my Dad. Acknowledgements
1 should like to thank my supervisor, Dr. Dawson, for his support and encourage- ment over the last two years, and Dr. Thornpson, for his endless supply of corrections.
1 should also like to thank Dr. M. Heukaeufer' a professor in the department of Modem Languages at Saint Mary's University, and Connie hk. for their help in translating the paper by Minkowski, (111, from Gerrnan to English.
Lastly, 1should like to thank my fnends and famil. for ail their love and support. Abstract
A polytope is said to be monostatic if its center of gravity is above only one facet. We show that no monostatic simplex exists in eight or fewer dimensions?and there is no monostatic simplex in ten dimensions that fdls sequentially onto dl its facets.
We cda polytope loadable if there is a point in the interior of the polytope that is above only one facet. We classify some convex polytopes as either loadable or nonloadable. Abbreviat ions and Symbols Used
e,= (0,. . . 0,1,0,. . . ,O)=, where the 1 is in the ith position of the vector.
0 u 5 V. where u = (ul,.. . , and v = (vlo... v,JT implies that ui 5 viVi =
1, ...,72
xii Chapter 1
Introduction
-4closed and bounded region in space is said to be convez if given aqv tn-O points p and q in the region, then the line joining them is aiso completely contained in the region. (In other w-or&. tp + (1- t )q is inside the region for every t E [O. 11). A conves set wïth nonempty interior is called a convez body. Let a E Rn. and let b E R. Then a hyperplane is the set of d points x = (xl. . . . .x,) satisfiing the equation a - x = b. The vector a = (alt... .a,) points in the direction normal to the hyperplane. and if llall = 1. b is the distance fiom the origin 0. to the point in the hyperplane ciosesr to it.
-4 polytope is definecl. see. for example. [4]to be a conves region in Rn bounded by a finite number of hwyperplanes. -1 h-yperplane that bounds a polytope is cailed a supportzng hyperplane, or a support plane of the polytope. If HI.. . . . Hm is a finite set of hwyperplanesthat bound a polytope P. Hi n P. where i = 1.. . . . m is said to be a facet of P provided that it is maximal in the sense that no Hj n P. j = 1. . . . . m. j # i properly contains it. If H, f~ P is a facet? then H* is called a facet hypeîplane. For any facet F* of a polytope, the facet vector vi is the vector with length equal to the (n - 1)-dimensional volume of the facet, directeci outward and normal to the facet. We may describe a polytope P mith m facets in the following -: The facet Fi lies in the hyperplane ai-x = bit V2 = 1,. . .,m, with 11aJ1 = 1. The facet vector vi that describes Fi is a multiple of the vector Vi = 1, . . . , m. Then if R-edefine Yi = ~~/ll~~ll,and let A = (Yl,.. . ,y,)=. b = (bl,.- . .bJT, P is descrîbed as the region defined by -4.x 5 b. ( For v = (vil .. .,un) E Rn.w = (wl,.. ., w,) E Rn-we Say v 5 w ifui 5 wi, Vz = 1,. . . ,n.)
Let a polytope P E Rn be given. Let density of P. p(x), be a non-oegative, real- valued function defined on the points x E P. If the volume of P is V. then define the mas of P to be m = J Jv p(x)dCi. The center of grauity of P is
If for ail x,y E P, p(x) = p(y), then we say the polytope has unzfom density.
We Say that a point is above a facet of a polytope if the line parallel to the facet vector v and through the point meets the interior of the facet. .A polytope is said to be stable on a facet if and only if its center of gravity is above that facet. Therefore the position of the center of gravity of a polytope determines whether or not the polytope is stable on a particular facet. A regular polytope (constructed) of uniform density will be stable on dl of its facets, since its center of gravity will be located at the center of the polytope, above every facet.
Imagine a pardelogram with two 45' angles, such that two sides are tnice as long as the other ones. Such a pardelogram will be stable on any of the long sides, but if set on a short side, it will fa11 over, since the center of gravity of the object is not above that face. A convex polytope of uniform density will fall from its ith facet to another facet if the center of gravity of the polytope is not above the ith facet. The dihedral angle between the ith facet and the facet onto which the polytope falls must be obtuse, and if the polytope is stable on its jth facet, the center of gravity di be above that facet. Because the angle between any two facets in a polytope such that it falls from one facet to the other must be obtuse, the angle between the facet vectors must be acute. Note that when a polytope falls boom one facet to another. the position of the center of gavity of the polytope is lowered. so the polytope cannot continue to fa11 kom one facet to another facet forever. In fact, a polytope will never fall onto a facet from which it has already fallen. This means a polytope must be stable on at least one of its facets, a resting jacet; a polytope which is stable on only one of its facets is said to be monostatic. Given a monostatic polytope P. we define a pendtirnate facet of P to be a facet from which the polytope fails onto its resting facet.
A simplex in Rn is defined to be the convex hull of n + 1 independant points. Equivalently, it can also be defined in Rnto be a convex body which is bounded by exactly n + 1 facet-planes. In a paper published in 1969 by Conway and Guy. [3]. it was proved that no tetrahedron or polygon is monostatic, and it was asked in how few dimensions a monostatic simplex can exîst, if any. Dawson. [5] showed in 1985. that no monostatic simplex exists in six or fewer dimensions, and he found one that does exist in RLo.In 1996, it was shown by Dawson, Finbow, and Mak, [7] that t here are no monostatic simplexes in seven dimensions or fewer, and a monostatic sirnplex exists in RH that fails through each of its facets in turn.
In Chapter 3, Theorem 3, it is shown that no such simplex exists in eight dimen- sions. In addition, we show that in ten dimensions, no monostatic simplex exists which falls sequentially through al1 of its facets. In Chapter 2, proofs are provided of previously known results which will be needed in Chapter 3. These include an English translation of Minkowski's Theorem on equilibrated facet vecton.
In Chapter 4 we look at polytopes of non-uniform density. -4 polytope is loadable if there is a point in the interior of the polytope that is above only one facet. If the polytope were weighted in such a way as to move the center of gravity to this point, then such a weighting would yield a monostatic polytope. We investigate the polytopes in several well known classes to determine which ones are loadable. These classes are: the regdar solids, the zonotopes, and the Archimedean solids.
Findy, in Chapter 4, we study the construction of density functions which may be useful for optimizing the position of the center of gravity. Let K be a bounded convex body in Rn, and for each density function p on K. let c(p) be the correspond- ing center of gavity of K. We consider the set D of density functions p whose range is {dl,d2}, where dl > 0, and d2 < W. For a given linear functional f on K. we constmct (Theorem 5) p* E D such that j(c(p*))= mau{/(c(p)) : p E D}.
The following notation is used throughout this thesis: ei = (0, .. . .O, 1.0, . . . ,O)=. where the 1 is in the ith position of the vector. O = (0, . . . ,o)~. and 1 = (1, . . . .1)*. If v = (u~,.. .,un), w = (utlt.. .,wn)E Rn,, then V = v/llvll. and v 5 w means Chapter 2
Preliminary Result s
In this chapter, we develop definitions and lemmas that will be used in Chapter 3 to prove various results on monostatic simplexes.
2.1 Minkowski's Theorem
-4 set of vectors {m,xi,. . . ,x,,,) in Rnis equzlzbrated [B. 111 if:
ii) {q} spans Rn,and
iii) (Vz # f,i,j =O, .. . ,m,a 3 O) # axj .
Given an rn x n matrix A, each of whose rows are unit vectors. and a vector b E Rm such that b 2 0, we define the polytope Pa4(b)= {x E Rn : .-lx 5 b)' and V,(b) to be the n-dimensional volume of P4(b). When no ambiguity can arise. the subscript A wilI be dropped. If vi is the ith row of A, then the number rl is called the ith tangentid parameter of the polytope if rl 2 0, and the plane v,-x= rf is a supporting hyperplane of the polytope. If wf is the ith tangential parameter of P(w) then P(w) = P(w*),where w*= (w;,. .. ,w;)=. We cal1 the vector w* a tangential vector of P(w).
Let TA = {b E Rm : 3w E Rmsuch that b is the tangential vector of P4(w)). We note that Tq is a cmex cone in Rm. Indeed. if u, w E TA,then for each i there exist y, z E Rn,SUC~ that y E P(u),and z E P(w) SUC~that vi*y= ui and Vi.2 = Wi. Hence vi*(Ay+ (1 - X)Z) = Xvi*y+ (1 - X)vi.z = Aui + (1 - A)wi (2-1) and therefore Xy+(l-X)z E P(Xu+(l -X)w), and we dsosee that Xy+(l-X)z is contained in the hyperplane vi*x= Aui + (1 - A) wi, thus Au+ (1- A) w is a tangential vector. The set TAis also a cone, since P(Aw) = XP(w) for al1 A.
Theorem 1 is a generalization from 3 to n dimensions of a theorem discovered by Minkowski which originally appeared in 1897, [II], in German. For the proof we dl need two technical lemmas.
Lemma 1 (Minkowski, [Il]) Let {vi}zR=,be a set of equzlzbrated vectors in Rn. Then there ezists a bounded polytope with m fucets, whose ichfacet vector is a multiple of vi
Proof. Let Yi = vi/llvill The system .Ix 5 1, where d = (Yl,Y*, . . . . defines the polytope P(1). The m planes described by the equation =tu = 1 are tangent to the unit ball in Rn,thus P(1)contains the ball. Furthemore, the point of tangency of the i" plane with the unit ball is contained in no other hyperplane, and hence each hyperplane intersects P(1) in a facet. Therefore 1 is a tangential vector for P(1),P(l) has rn facets and its ichfacet vector is a multiple of Yi.
Findy we show that P(1)is bounded. Ifq E P(I),then vi*q= II~illYi*q5 II~ill, for i = 1, . . ., m. On the other hand, since {vi) is equilibrated, for ail q E P(1). But then P(l)is a subset of the bounded parallelepiped B obtained as follows: let S c {vi)ZLbe a basis for Rn and define B = nvEs{x: a, 5 v - x 5
Lemma 2 (Brunn, [2])Let A and B be conuez sets in Rn. and let Q4) and V(B)be the volumes of the sets A and B respectivefg. Then for t E (O. 1) we have dv(t4+ (1 - t)B) 2 t + (1 - t) uith equality if and only if -4 is obtained ftom B by dzlation and translation.
Theorem 1 (Minkowski, [Il])Every equilibrated set of uectors in Rn is the set of fucet vectors for a unique conuez polytope up to translation: and the set of facet vectors of any convez polytope is equilibrated.
Proof. Let P be any convex polytope, and let P have rn facets Fi. mhere the m facets have an arbitras. ordering. The following construction and result ing formulas for the derivative will be useful in both halves of the proof.
Consider any q in the interior of P and define the one parameter family of poly- topes P(t)= tP + (1- t)q. Let C(t)be the volume of P(t),and for i = 1.. . . . m. let C&) be the volume of the convex hd, E(t)of {q) U (tFi + (1 - t)q),d,(t) be the (n - 1)-dimensional volume of tF, + (1- t)q, and hi(t)be the length of the perpen- dicular àropped fiom q to the hyperplane containing tFi + (1- t)q. Then it is easy to see that C(t)= CgR=,Ci(t), k(t) = th(l),and &(t) = tn-'.li(l). NOW,since Pi(t) is a pyramid, Ci(t)= iAih(t) = :tn%(l)hi(l),and hence dCi(t) = tn-lAi(l)hi(l)= Ai(t)hi(l). If V denotes the total volume of P, we have both
V = &hi(l) and dV = didhi. i= 1 n i= 1 Now, let {vilbe the set of facet vectors of such a polytope P. we will show that
{v,) is an equilibrated set. Choose q in the interior of P. there exists E > O such that qk = q + cep is in the interior of P for k = 1: 2,. . . , n? and then b- 2.3 we have both
where qki is the length of the perpendicdar dropped from qk to the hyperplane containing Fi,V k = 1, . . ., n. But then
for k = 1,. . . , n and hence m
Therefore, {y) is equilibrated, for, since P has non-empty interior. {vi)spms Rn. and convexity assures that vi # avi, Vi # j, a 2 O.
Conversely, let {vi)g"=,ean equilibrated set of vectors in Rn. Let -4 be the matrix whose ith row is the vector Ci,where Yi = vi/j(vi11, let T = T4be the convex cone of tangential vectors, and for each x E T, let V(x)= I.A(x) be the n-dimensional volume of P(x) = PA(x).Suppose x = (X~~....X,)~,and set r = mâ~~<~~~{~~}. Xow P(l)is bounded, as shown in Lemma 1, and P(;x) C P(1).hence if we define p so that V(1) = pdn we obtain
This implies that P(x)is bounded.
Let A,(x) be the area of the ith facet of P(x). Xote that .-li(x)and V(x) are continuous functions of x; indeed V is differentiable and by Equation 2.3. 61- =
Setting f(x) = (r,x) x E Tl we have by Lemma 2.
for x, u E Tl and hence the set
T T W = {W = (wl1... ,w,+l) : (utl,.. ., wJT E T and O 5 w,+l 5 /(L.I,.. . . w,) ) is convex. We further note that it is a cone with apex at the origin. For x E T. x =
(xll...,x,)~. Set x = (21,.. .,x,+~)*,where X,+I = f(x) and let L: = {x:x E T). Observe t hat and hence we have
i= 1 is the tangent hyperplane to W at the point x E U. It follows that for dl q E 11,: q =
To prove existence of the polytope, consider ÇVl = {w E W : w,+~ = 1). let d = ( IIv~ll,. - - ,IIvm 11, O)=, and define the real-valued function (T on CVl by set ting O(W) = d-W.First, it is established that the minimum value of o is attained on LC;.
We do this by showing that the set Z = {w E 1.h : d-w 5 p XE1IIv~~~) is com- pact and nonempty. Indeed V(p1) = pnV(l) = 1 so r = (p.. . . .p. 1)= E 11'; and d-r = pEhl Ilvill: and therefore we see that Z # 0. On the other hand, Z is com- pact since d=(wl,.. . ,w,+~)~ 5 p llvi1l implies chat mi 5 fi zz"=,Ivi11. where M = m~~~i~n{llvill).
Note that the minimum value of o is positive: for if w E T is such that TV E 1.b;
~d r = ~u~<~<~{w~),- - then by 2.7, p = JI ~v(w)5 r and hence d-W 2 .Ur 2 Mp.
Let p, = (pl,. . ., P,+~)~ = (pl, . . . p,. l)= be the point where the minimum
OCCUIS, let p = (po,.. . ,Pm), and let the minimum va,lue of 0 be a. Then we have d-po = a. If x E W, = -E IVl, so ~-X/X,+~2 a, which implies d-x 2 ~m+l î/&F) ax,+l. This implies that XEl ll~~ll~~= axrn+i supports the cone FI; at p, and hence must be the same hyperplane as 2.10 with x = po. Comparing coefficients and not ing that V(pj = 1, we see that
Now P(p) is a polytope whose ith facet has n - 1dimensional volume nllvill/a. There- fore P = P (($ p) is a polytope with the correct face vecton.
To see that this po-ytope is unique up to translation, suppose the set of face vectors of both P(p) and P(q) is {vl,.. . ,v,). Since ll~ill> 0 W ~e sec that P and q are in T. Observe that since Ai(p) = Ai(q) = Ilvill, we have by 2-10 and 2-11 which implies V(p) 5 V!q). Symmetrically, V(q) 5 C'(p). so Y(p) = Ce(q).
Xote that this implies by equation 2.10 that the hyperplane that supports CV at p and the hyperplane that supports W at q are identical. and the line {w + (1- t)q: t E [O, 11) also lies on this plane. But then
so /(tp + (1 - t)q)= t f (p) + (1- t)f (q). But then by Lemma 2. rve have that P(p) can be obtained from P(q) \la translation and dilation. and since C'(p) = C'(q). the dilation factor is 1.
2.2 Further Results
In this section, we concentrate on results about simplexes. Given a simplex in Rn with facet vectors Q, . . . ,%, we dehe pi(xj) = ~=x,/ll~ll= IIx~~~ COSO~~~ where Qij is the angle between vectors and xj. In the following lemrna we determine what conditions must hold in order for a simplex to fd from a facet i to a facet j. The Projection Criterion, which originally appeared in [5], gives an inequality involving the area of the ith facet of the simplex, and the projection of the jth face vector onto the ith face vector. Lemma 3 (Projection Criterion, Dawson, [5]) A simplez wzth/acet vectors {v*} will tip from the ith facet to the jl facet if and only zf // < pi(xj)-In particular, the ith facet is smaller (that is, the (n- 1)-dimensional volume) then the jth facet.
Proof. Let F, be the facet described by xi, and vi be the vertex opposite Fi in the simplex. The center of gravity of a simplex, c = ,1 (1:z vi). Let L be the (n - 2) - dimensional hyperplane containing Fin F,, and without loss of generality, assume that the origin O is the foot of the perpendicular dropped from v, to L. See Figure 2.1.
Let pLL(x) denote the vector projection of a point x into LI. and define A(& n 4)to be the (n - 2)- dimensional volume of F, n Fj. Since Ilxi(l is the area of
= l*llpLl(vj)ll.-I(Fin-L n F,), and similarily 11~~11= t herefore follows t hat
Let M be the hyperplane containing LL fi Fi,and define pLtI(x)to be the scalar projection of the point x into M. Then
Given that the simplex falls £rom its ith facet onto its jth facet, its center of gravity c must not be above the ith facet, therefore pM(c)< O. Then holds if and onlv if
and the statement of the lemma foliows.
Figure 2.1: Illustration of the Projection Criterion (enahedron) The following lemma domus to project the facet vectors of a simplex into the plane spanned b- 3 and x,,, thereby simplif'ing any calculations to do with the projection of a facet of the simplex onto W.
Lemma 4 (Dawson, Finbow, Mak, [?]) Let Bi, be the angle between cectors x, and xj in Rn. For any Hven value of Bo, E (0. a). cos Bol cos 012- - - COS em-lm is m&màzed for Bi i+l E (O. ~/2)~helien Bol = B12 = = Bm-lm = Bh/m-
Proot We need to show three things in order to prove this lemma. one. the prod- uct cos OO1 cos BI2 COS 0,- 1 m is ma.-edden the vectors {x, } Lie in a plane. two. the above product is ma.cimized when the vector x, lies in the cone generated b-x,-~ ad and three, it is mashized when OO1 = OI2 = . . . = 0rn-l m.
Take any three vectors {q?x, Y xr } in Rn. If al1 t hree do not lie in a plane. then we can project the jfi vector into the plane spanned by x, and xk. and ob- tain a new vector $ to replace xj. If tj is the angle betn-een x, and x,. then COS B;, COS 5, > cos 19, cos Ojk. This is because q, < Bij. and cosine is a decreasing POS- itive function on (O. ;i/2)' w COS Bij < COS q,,. similady n-e see that COS qk> COS ojk. so COSO,~COSO~~~ahere O*jj'Ojk € (0. ~12).is maximized when the three vectors lie in a plane. Since this is the case for any set of three vectors. it is t hen rrue for al1 vectors.
Sel? we must show that the product cos Ool cos e12- cos Om-l, is maximized when x, is a vector lyhg in the cone generated by q-l and the plane. for if this were SO, we can say that the product is maximized when Bi-1 i-1 = i i Bi i-1. This can be proved by contradictiony so suppose xi does not lie betn-een x,-1 and x+ll Say! q+l lies between q-1 and 4, instead. The angle Bi-l > Bi-l i,l. so cos Bi i-l < COS Oi-, i,, . Therefore interchanging xi and X+I ni11 increase the product of the cosines of those angles.
The argument above tells us that the product of the cosines is maximized n-hen ali the vectors lie in a plane' and the sum of all the angles between the vectors is equal to the angle between the first vector m, and the last vector. h.If a?0 E (O, r), and a + ,O is constant, then 1 COS CE cos ,û = - [cos(cr + 8) + cos(a - O)] (2.18) 2 is maximized when a - /3 = O, that is, when O = B. Hence cos eOlcos - COS for Qi-l i E (0, '2) Vz is maximized when al1 the angles between each of the wctors is the same. 1
In light of Lemma 3, for vectors G,xj, we write x, \A xj (read xi falls onto Xj) if and only if llqll < pi(xj) For i = O, 1,. . ., let
Given a simplex in Rn, with facets Fo, . . . , F, which falls from Fito Fo t hrough other facets, the number ai represents a bound on how small the projection of an individual facet vector x, onto xo can be. The following lemma. which first appeared in [5], determines the bound on how smd ai cm be for a given 2.
Lemma 5 (Dawson, [5]) If x, \ q-i \ - - \ q,then and by the Projection Criterion (Lemma 3), kJ.LIlx,ll < COS Oij,
Multiplying both sides of Equation 2.21 by cos O&, we see that
a- L 11411,me, > min ~,co~e~IIcos~~,+~ llxoll j=o
The product in 2.22 is maximized when the angles Oj j+i lie in a plane by Lemma 4, so c~~~B~~+~2 19,. Adding a to both sides and solving, we find (T - Bot) + 6, j+i 2 r. The product of the cosines is maximized when the angles are equal, that is, when a - Bo, = Oo1 = . . . = 1 ~/i+ 1, and we find that
Since cos(n - O~)= - cos go, ,
Finally, since - (cos (&))i+L < O, we have, by 2.22 and 2.23,
PO(%)= (rn" cos 0 ) llqll
Lemrna 5 implies that PO(%) ai = min - llxall It is also clear that a0 = 1. If xl \ qlthen Bol must be acute, so al 2 0. and by the projection criterion, a-1 2 1.
For a simplex in Rn with facet vectors qt...,%, where llxoll 2 3 ll~~ll, optimizing the projection of each individual vector onto one particular vector will give bounds that enable us to show that no monostatic simplex exists in five or fewer dimensions as follows: Let a set of six vectors in R5 be given and suppose, for a contradiction, that they are the facet vectors of a monostatic simplex. There are two cases to consider, if the simplex has 1 penultimate facet, and if the simplex has more than 1 penultimate facet. For the case where the sirnplex has 1 penul- timate facet xl, assume without loss of generality that llxl 11 = 1- The projection of a onto xl is bounded below by a-1, and the projection of a vector onto it- self is a*. .4.s well, the projection of the ilh vector, qtonto XI is bounded belon. by ai-1. The sum of the projections of the vectors xo, . . . , xj onto xl is bounded below by a-l + a. + al + a:! + a3 + a4 zz 1.278443219 > 0: and therefore the vec- tors cannot be equilibrated. In the other case, where there are two penultimate facets, and we project each vector onto q,so without loss of generality, assume 11% 11 = 1. Then the mm of the projection of the vectors must be bounded below b- a0 + 2a1 + a2 + as + a4 x 0.27843219 > 0, and we see that no set of six vectors in R5 cmbe equilibrated, and therefore no monostatic simplex exists in R5.
Optimizing the projection of two vectors simultaneously onto a will give us better bounds, and allow us to show that no monostatic simplexes exist in sLu and seven dimensions.
The parameter bi is defined to be the minimum of the combined projections of xl and onto the vector W. For any bi we have:
Hence if the system tips from xl to Q, we get bi > COS*^^^ - COS^^^. Setting x = cos OO1, we get bi 1 x2 - x. This is minimized when z = 8,therefore,
\Ire now have bounds on vector projections that are tight enough to show that no monostatic simplex exists in R6. (see [SI). Let a set of seven vectors in R6, 3,. . . ,~g, where 11% 11 2 - 3 11% 11 be given, and suppose, for a contradiction, that they are the facet vectors of a monostatic simplex. Again there are two cases to consider, depending on whether the simplex has 1 or more than 1 penultimate facet. If there is only one penultimate facet, the sum of the projections of al1 the vectors onto the vector describing the penultimate facet, xl, is bounded below by a-1 + a. + a* + a3 + a* + b5 = 1.02843219 > O, contradicting the assumption that the vectors are equilibrated. If there is more than one penultimate facet, the sum of the projections of al1 the vectors onto the vector describing the largest facet, xo, is bounded below by a0 +ai + a2 + as + a4 + bs = 0.28432189 > O, and no set of 7 vectors in R6 which obey the various inequalities required by the projection criterion for a monostatic simplex is equilibrated, therefore there is no monostatic simplex in Rb For sorne calculations, we will need to improve the bound on bi. Lemma 6 (Dawson, Finbow, Mak, [7])Suppose x, \ xi-1 \
bi 2 cos2(ir ) + COS(Z~)cosi (fif), where y = (a - Boi)/i
Proof. -4s was seen, it is straightforward to determine a lower bound on bi = min{po(xL+ q)/llm(l : (3 x.,-*, . . . .XL E Rn)xi\ xi-1 \ . \ q). However. -I4 is a somewhat rough estimate, and when the number of vectors lying in the cone generated by Q and q,and hence the angle Ba is small. then bi will be larger than
-I.4 This is because the angle between any two facets in a simplex such that the simplex falls from one facet to the other, must be obtuse. Hence the angle between the facet vectors must be acute.
First, note that
by the Projection Criterion, and similarily,
Note that COS Bli = - COS(* - Oli)- By Lemma 4, the product in 2.30 is maximized when the vectors lie consecutively in a plane, so cor + + Bi-i < Bo*, and the angles Bi-i are dl equal to &,Ji. This implies that the product is maximized when T - Ooi = Oi-1i = Bol, thus Bo, = Since 7 = r - BO1/i,Bol = 7r - i7 Substituting these into Equation 2.31. Lemma 7 follows.
In particular, 1
To see this, suppose that x2 \ xl \ a.Since b2 = minpo(xl + x2)/llxali, Lemma 7
Differentiating the right-hand side of Equation 2.33 with respect to 7 and minimizing, we see that
O = -4 cos(27) sin(2y) - 2 sin(2y) cos2y - 2 cos(27) sin cos
= - sin y cos y[8(2 cos2 7 - 1) + 4 cos2 nf + 2(2 cos2 - l)] = (12~0s~~-5)(-2suiycosy).
Solving for y, nie find the smaller of the two solutions is
Substituting y = arccos(,/5/12) into Equatioo 2.33, we see that 62 2 2 .
In solving for each remaining bi, i = 3,d, 5, we obtain the following numerical solutions:
Lemma 7 (Dawson, Finbow, Mak, [7]) If O 5 i 5 j, then ai 2 aj, and b, 1 bj. Lemma 7 relies on the fact that the projection of a vector v onto another vector u is bigger than the projection of a third vector w onto u- if the angle between u and w is greater than the angle between u and v, and llvll 2 Ilwll To show bj 2 bk, where j is less then or equal to k, it is sufficient to show that bi 5 biVl for al1 i > 0.
A key step to showing that no monostatic simplex elusts in Riand R8 requires that given two penultimate facets Fiand F', such that the simplex never falls onto F., we find a bound for the minimum of pl(x. + Q), where x. and q are the facet vectors corresponding to facets F. and Fo respectively.
Lemma 8 (Dawson, Finbow, Mak, [7])Suppose xt \ % and x. \ Q. Then
To prove Lemma 8, we consider two cases: one when the angle between facet vectors xl and x. is les than or equal to g, and the other case, when BI. > 5. Let a monostatic simplex with two penultimate facets described by the vectors xl and x. which both fall onto a final vector be given. The goal is to find an upper bound on how small pl (x. + q)can be.
If 01, 5 a, then cos&. -> O and so pl(x.) 2 O. This implies that pl(m + x*) = 1 pi (m) +pl (4)2 pi (m) Since $ cos Bol < i, we actually have cos Ool > $ cos Bol - 5. ho,çince = cos solthen Noting that BI. = Bol +Bo., and taking the derivative with respect to 80. to minimize the right hand side of Equation 2.36, we find that
O = - sin 4. cos(80i + 60.) - cos 80. sin(Ooi + Bo.) = - sin(Ool + 2Bo.). (2.37)
Since xl \ ~0 and x. \ a,801 and Bo. must be acute, so Bol + 2eo. = T!and solving for Oo,, we find that
00. = (T - eol)/z. (2.38)
Substituting this into equation 2.36, and using half angle formulas to simpli@. we find that
7r - 001 n + 001 = coseoi +cos ( * ) cos ( )
For the next reçult we need the following definition, which modifies the deletion of ai. We cm now prove that in R7,no monostatic simplex exists. (See [il). Let a set of eight vectors in R~,Q,. .. ,XT, where ll~lj2 - - > Ilx711 be @en, and suppose. (for a contradiction) that they are the facet vectors of a monostatic simplex. The proof has six different cases, depending on how many penultimate Eacets there are. and on how many facets fa11 onto them. Five cases are straightforward. and can be verified in the same way as we showed that no monostatic simplex exists in fewer than 7 dimensions above. The sixth case, in which there is a penultimate facet onto which the simplex never falls is more difficult, and requires the use of Lemma 8. There is a penultimate facet onto which the simplex never falls. Without loss of generalit- label the corresponding facet vector x~.If the other penultimate facet is described by xa, then the sum of the projections of these vectors onto xo is greater than. by Lemma 3,
Note that if the set {q)is equilibrated, the left hand side is equal to zero. and Bo, > arccos 113 = 70.52878". Projecting onto m, by Lemrna 3 we find
PO(^') > - cos Bo, cos' llvoll (" -i"~) . 1 L and a2(Ooa) > -119, a3(OOa)> -0.173211.. . , a4(aoa)> 9 + 6;.i '-: -0.207336 - . .. a5(Boa)> -0.229257. . . , and as(Ooa) > -0.24606. . .. Therefore the surn of the projection of dl the vectors ont0 the largest vector, vo is bounded below by a0 + 2al + ~~(0,~)+ a3(e0,) + a4(OOa)+ a5(Boo)+ a6(Ooa) = 0.03478 > 0 adthe set of vectors is not equilibrated. Hence no monostatic simplex exists in R'. Xote that the sum of the projections of eight vectors is baxely positive. Csing the lemmas and theorems established in this chapter, we are unable to show that monostatic simplexes do not exist in R8. Chapter 3
Monostatic Simplexes in 8, 9, and 10 Dimensions
In Chapter 2 ae shoaed how. optimizing the projection of one vector individually. or rao vectors added together onto a single vector. a-e could prove no monostatic simplex exists in 7 dimensions or fewer. The techniques del-eloped in the previous chapter are not sufficient on their onto show that no monostatic simples esists in R8. and a natural extension of them wodd be to try to optimize the projection of the sum of three or more vecsors onto a single vector. This is veq- dficult analgicallc so in this chapter ae use a cornputer to optimize the projecrion of the sum of three or more vectors onto a single vector.
The aigorithm for the program that fin& the overd maximum of the projection of the sum of 3 or more wcton onto a single vector is discussed in Section 3.1 and appears in Section 3.2. In 3.3. ae provide a suitable upper bound on the error that may occur as a rdtof using a cornputer search ro calculase this maximum. and in Section 34we use our resdts to show no monostatic simples e-dts in 8 dimensions. and various other resdts about simpleves in 9 and 10 dimensions. 3.1 The Search Method
Let a sirnplex in Rn which falls through each of its facets in turn be given. Let vo, . . . , v, be its facet vectors. Then without los of generality. llvo11 2 IIvl II 2 2 Ilv, 11. Suppose vo points dong the negative 2-&S. Then if llvn11 = 1, the projection of v, onto vo, po (vn)= cos a,where a E (O. ?r). We want to find a formula for the projection of each of the facet vectors onto one single facet vector. Define fo(a)= cos (Y. Chooçe rn = 1, .. . : n - 1, and ,/3 E (O. r). We place the rn + l vectors v,-,, . . . ,v, so that vn-, makes an angle of 4 with the x-ais, and has length 1: and the other m vectors are at angles les than 4. We define fm(3)= max C:.,-,+, po(vi).Recall the Projection Critenon. which states that in order for a simplex to fall from its ichfacet to its jch facet, llxill < ll~,llcos0~~.Since the simplex we are looking at falls sequentially through al1 its facets. then IIviii 11 < llvill If v,-,+l has inclination LY and length 1. then by the Projection Criterion, IIv~-~+~II< COS(Q-p). Then f,+,(a) = cosa+rnaiq,g<,cos(a-i3) fm(,3). Therefore, for a simplex in Rn which falls through each of its facets sequentially, the projection of the sum of its facet vectors onto the largest facet vector. vo, ivill be
where fo(a)= cos a.
Given a directed tree with nodes x and y, x is said to be the parent of y if there exists a directed edge from x to y. Similady, y is the child of x if there exists a directed edge from x to y. We define the falling pattern of a monostatic sirnplex to be the directed tree with one root node such that each node in the tree represents a facet of the simplex, and an edge joining two nodes indicates that the simplex falls from the child node onto the parent node. Note that the root node dlrepresent the largest facet of the simplex, since by the projection criterion, a simplex will only fa11 from a smaller facet onto a larger one. Given a simplex with falling pattern r, the root node of T describes the largest facet of the simplex. If the vector representing the facet corresponding to the root of r has length 1 and lies dong the negative x-axis,we dehe f,(~)to be the maximum projection of the facet vectors of the simplex onto the laxgest facet. Then. as in 3.1 and 3.2, we see that