arxiv:math/ v3 [math.mg] 1 Mar 2006

Transcription

1 arxiv:math/ v3 [math.mg] 1 Mar 2006 THE KISSING PROBLEM IN THREE DIMENSIONS Oleg R. Musin Abstract The kissing number k(3) is the maximal number of equal size nonoverlapping spheres in three dimensions that can touch another sphere of the same size. This number was the subject of a famous discussion between Isaac Newton and David Gregory in The first proof that k(3) = 12 was given by Schütte and van der Waerden only in In this paper we present a new solution of the Newton-Gregory problem that uses our extension of the Delsarte method. This proof relies on basic calculus and simple spherical geometry. Keywords: Kissing numbers, thirteen spheres problem, Newton-Gregory problem, Legendre polynomials, Delsarte s method 1 Introduction The kissing number k(d) is the highest number of equal nonoverlapping spheres in R d that can touch another sphere of the same size. In three dimensions the kissing number problem is asking how many white billiard balls can kiss (touch) a black ball. The most symmetrical configuration, 12 billiard balls around another, is achieved if the 12 balls are placed at positions corresponding to the vertices of a regular icosahedron concentric with the central ball. However, these 12 outer balls do not kiss each other and may all be moved freely. So perhaps if you moved all of them to one side a 13th ball would possibly fit in? This problem was the subject of a famous discussion between Isaac Newton and David Gregory in 1694 (May 4, 1694; see interesting article [21] for details of this discussion). Most reports say that Newton believed the answer was 12 balls, while Gregory thought that 13 might be possible. However, Casselman [5] found some puzzling features in this story. This problem is often called the thirteen spheres problem. Hoppe [9] thought he had solved the problem in But there was a mistake - an analysis of this mistake was published by Hales in 1994 [8] (see also [20]). Finally this problem Institute for Math. Study of Complex Systems, Moscow State University, Moscow, Russia omusin@mail.ru 1

2 was solved by Schütte and van der Waerden in 1953 [19]. A subsequent twopage sketch of an elegant proof was given by Leech [11] in Most people agree that Leech s proof is correct, but there are gaps in his exposition, many involving sophisticated spherical trigonometry. (Leech s proof was presented in the first edition of the well-known book by Aigner and Ziegler [1], the authors removed this chapter from the second edition because a complete proof would have had to include so much spherical trigonometry.) The thirteen spheres problem continues to be of interest, and new proofs have been published in the last few years by Hsiang [10], Maehara [13], Böröczky [3] and Anstreicher [2]. The main progress in the kissing number problem in high dimensions was at the end of the 1970s. Levenshtein [12], and independently, Odlyzko and Sloane [16] (= [6, Chap.13]) using Delsarte s method in 1979 proved that k(8) = 240 and k(24) = This proof is surprisingly short, clean, and technically easier than all proofs in three dimensions. However, d = 8, 24 are the only dimensions in which this method gives a precise result. For other cases (for instance, d = 3, 4) the upper bounds exceed the lower. We found an extension of the Delsarte method in 2003 [14](see details in [15]) that allowed us to prove the bound k(4) < 25, i.e. k(4) = 24. This extension also yields a proof k(3) < 13. The first version of these proofs was relatively short, but used a numerical solution of some nonconvex constrained optimization problems. Later [15] these calculations were reduced to calculations of roots of polynomials in one variable. In this paper we present a new proof of the Newton-Gregory problem. This proof needs just basic calculus and simple spherical geometry. 2 k(3) = 12 Let us recall the definition of Legendre polynomials P k (t) by the recurrence formula: P 0 = 1, P 1 = t, P 2 = 3 2 t2 1 2,..., P k = 2k 1 k or equivalently P k (t) = 1 d k 2 k k! dt k (t2 1) k t P k 1 k 1 k (Rodrigues formula). P k 2 Lemma 1. Let X = {x 1, x 2,..., x n } be any finite subset of the unit sphere S 2 in R 3. By φ i,j = dist(x i, x j ) we denote the spherical (angular) distance between x i and x j. Then P k (cos(φ i,j )) 0. i=1 j=1 2

3 This lemma easily follows from Schoenberg s theorem [18] for Gegenbauer (ultraspherical) polynomials G (d) k. (Note that P k = G (3) k.) For completeness we give a proof of Lemma 1 in the Appendix. Let f(t) = t t t t t t2 + t Remark. This polynomial of degree 9 is satisfying the assumptions of the extended Delsarte s method [14, 15]. An algorithm for calculating suitable polynomials is presented in the Appendix of [15]. Lemma 2. Suppose X = {x 1, x 2,..., x n } S 2. Then S(X) := i=1 j=1 Proof. The expansion of f in terms of P k is f = 9 k=0 f(cos(φ i,j )) n 2. c k P k = P P P P P P P 9. We have c 0 = 1, c k 0, k = 1, 2,..., 9. Using Lemma 1 we get S(X) = 9 n c k k=0 i=1 j=1 P k (cos(φ i,j )) i=1 j=1 c 0 P 0 = n 2. If n unit spheres kiss the unit sphere in R 3, then the set of kissing points is an arrangement on the central sphere such that the (Euclidean) distance between any two points is at least 1. So the kissing number problem can be stated in another way: How many points can be placed on the surface of S 2 so that the angular separation between any two points is at least 60? Lemma 3. Suppose X = {x 1, x 2,..., x n } is a subset of S 2 such that the angular separation φ i,j between any two distinct points x i, x j is at least 60. Then S(X) = i=1 j=1 f(cos(φ i,j )) < 13n. We give a proof of Lemma 3 in the next section. Theorem. k(3) = 12. Proof. Suppose X is a kissing arrangement on S 2 with n = k(3). Then X satisfies the assumptions in Lemmas 2 and 3. Therefore, n 2 S(X) < 13n. From this n < 13 follows, i.e. n 12. From the other side we have k(3) 12, showing that n = k(3) = 12. 3

4 3 Proof of Lemma 3. We need one fact from spherical trigonometry, namely the law of cosines: cosφ = cosθ 1 cosθ 2 + sin θ 1 sin θ 2 cosϕ, for a spherical triangle ABC with sides of angular lengths θ 1, θ 2, φ and BAC = ϕ (Fig. 1). For ϕ = 90, this reduces to the spherical Pythagorean theorem: cosφ = cosθ 1 cosθ 2. Proof. 1. The polynomial f(t) satisfies the following properties (see Fig.2): (i) f(t) is a monotone decreasing function on the interval [ 1, t 0 ]; (ii) f(t) < 0 for t ( t 0, 1/2]; where f( t 0 ) = 0, t These properties hold because f(t) has only one root t 0 on [ 1, 1/2], and there are no zeros of the derivative f (t) (eighth degree polynomial) on [ 1, t 0 ]. B A ϕ θ 1 θ 2 φ C Fig Fig. 2. The graph of the function f(t) Let S i (X) := f(cos(φ i,j )), then S(X) = j=1 S i (X). From this it follows that if S i (X) < 13 for i = 1, 2,...,n, then S(X) < 13n. We obviously have φ i,i = 0, so f(cosφ i,i ) = f(1). Note that our assumption on X (φ i,j 60, i j) yields cosφ i,j 1/2. Therefore, cosφ i,j lies in the interval [-1,1/2]. By (ii) we have f(cosφ i,j ) 0 whenever cosφ i,j [ t 0, 1/2]. Let J(i) := {j : cosφ i,j [ 1, t 0 )}. We obtain i=1 S i (X) T i (X) := f(1) + j J(i) f(cosφ i,j ). (1) Let θ 0 = arccost Then j J(i) iff φ i,j > 180 θ 0, i.e. θ j < θ 0, where θ j = 180 φ i,j. In other words all x i,j, j J(i), lie inside the spherical cap of center e 0 and radius θ 0, where e 0 = x i is the antipodal point to x i. 4

5 2. Let us consider on S 2 points e 0, y 1,...,y m such that φ i,j := dist(y i, y j ) 60, i j, θ i := dist(e 0, y i ) < θ 0 for 1 i m. (2) Denote by µ the highest value of m such that the constraints in (2) allow a nonempty set of points y 1,...,y m. Suppose that 0 m µ and Y = {y 1,...,y m } satisfies (2). Let H(Y ) = H(y 1,...,y m ) := f(1) + f( cosθ 1 ) f( cosθ m ), h m := sup {H(Y )}, h max := max {h 0, h 1,..., h µ }. Y It is clear that T i (X) h m, where m = J(i). From (1) it follows that S i (X) h m. Thus, if we prove that h max < 13, then we prove Lemma Now we prove that µ 4. Suppose Y = {y 1,..., y m } S 2 satisfies (2). By symmetry we may assume that e 0 is the North pole and y i has polar coordinates (θ i, ϕ i ). Then from the law of cosines we have: cosφ i,j = cosθ i cosθ j + sin θ i sin θ j cos(ϕ i ϕ j ). Note that θ i > 0 for m 2. Conversely, y i = e 0, θ j = φ i,j 60 > θ 0, a contradiction. From (2) we have cosφ i,j 1/2, then Let then cos(ϕ i ϕ j ) 1/2 cosθ i cosθ j sin θ i sin θ j. (3) Q(α, β) := 1/2 cosαcosβ, sin α sin β Q Q(α, β) 2 cosβ cosα α (α, β) = = α 2 sin 2 α sin β. From this it follows that if 0 < α, β θ 0, then cosβ > 1/2 (because θ 0 < 60 ); so then Q α (α, β) > 0, and Q(α, β) Q(θ 0, β) = Q(β, θ 0 ) Q(θ 0, θ 0 ). Therefore, 1/2 cosθ i cosθ j 1/2 cos2 θ 0 sin θ i sin θ j sin 2 = 1/2 t2 0 θ 0 1 t 2. 0 Combining this inequality and (3), we get cos(ϕ i ϕ j ) 1/2 t2 0 1 t 2. 0 Note that arccos((1/2 t 2 0 )/(1 t2 0 )) > 72. This implies that m 4 because no more than four points can lie in a circle with minimum angular separation between any two points greater than 72. 5

6 4. Now we have to prove that h max = max {h 0, h 1, h 2, h 3, h 4 } < 13. We obviously have h 0 = f(1) = < 13. From (i) follows that f( cosθ) is a monotone decreasing function in θ on [0, θ 0 ]. Then for m = 1 : H(y 1 ) = f(1) + f( cosθ 1 ) attains its maximum at θ 1 = 0. So then h 1 = f(1) + f( 1) = < Let us consider for m = 2, 3, 4 an arrangement {e 0, y 1,...,y m } in S 2 that gives H(Y ) = h m. Here y i e 0 (see 3). Note that in this arrangement, points y k cannot be shifted towards e 0 because in this case H(Y ) increases. For m = 2 this yields e 0 y 1 y 2, and dist(y 1, y 2 ) = 60. If e 0 / y 1 y 2, then the whole arc y 1 y 2 can be shifted towards e 0. If dist(y 1, y 2 ) > 60, then y 1 (and y 2 ) can be shifted towards e 0. For m = 3 we prove that 3 = y 1 y 2 y 3 is a spherical regular triangle with edge length 60. As above, e 0 3, otherwise the whole triangle can be shifted towards e 0. Suppose dist(y 1, y i ) > 60, i = 2, 3, then dist(y 1, e 0 ) can be decreased. From this follows that for any y i at least one of the distances {dist(y i, y j )} is equal to 60. Therefore, at least two sides of 3 (say y 1 y 2 and y 1 y 3 ) have length 60. Also dist(y 2, y 3 ) = 60, conversely y 3 (or y 2, if e 0 y 1 y 3 ) can be rotated about y 1 by a small angle towards e 0 (Fig.3). For m = 4 we first prove that 4 := conv Y (the spherical convex hull of Y ) is a convex quadrilateral. Conversely, we may assume that y 4 y 1 y 2 y 3. The great circle through y 4 that is orthogonal to the arc e 0 y 4 divides S 2 into two hemispheres: H 1 and H 2. Suppose e 0 H 1, then at least one y i (say y 3 ) belongs to H 2 (Fig.4). So the angle e 0 y 4 y 3 greater than 90, then (again from the law of cosines) dist(y 3, e 0 ) > dist(y 3, y 4 ). Thus, θ 3 = dist(y 3, e 0 ) > dist(y 3, y 4 ) 60 > θ 0 a contradiction. Arguing as for m = 3 it is easy to prove that for any vertex y i there are at least two vertices y j at the distance 60 from y i. Note that the diagonals of 4 cannot be both of lengths 60. Conversely, at least one side of 4 is of length less than 60. Thus, 4 is a spherical equilateral quadrangle (rhomb) with edge length Now we introduce the function F 1 (ψ), 1 where ψ [60, 2θ 0 ]: F 1 (ψ) := max ψ/2 θ θ 0 { F 1 (θ, ψ)}, F1 (θ, ψ) = f( cosθ) + f( cos(ψ θ)). So if dist(y i, y j ) = ψ, then Therefore, f( cosθ i ) + f( cosθ j ) F 1 (ψ). (4) H(y 1, y 2 ) h 2 = f(1) + F 1 (60 ) < For given ψ, the value F 1 (ψ) can be found as the maximum of the 9th degree polynomial Ω(s) = F 1 (θ, ψ), s = cos (θ ψ/2), on the interval [cos(θ 0 ψ/2), 1]. 6

7 y 3 y e y 1 Fig. 3 y 3 y 4 H 1 e 0 y 1 y 2 Fig. 4 H 2 y θ y 3 c θ y 1 1 e θ 2 0 y 2 Fig Let m = 4, d 1 = dist(y 1, y 3 ), d 2 = dist(y 2, y 4 ), Since 4 = y 1 y 2 y 3 y 4 is a spherical rhomb, we have cos(d 1 /2)cos(d 2 /2) = 1/2 (Pythagorean theorem, the diagonals y 1 y 3, y 2 y 4 of 4 are orthogonal). So if then ρ(s) := 2 arccos 1 2 cos(s/2), ρ(d 1 ) = d 2, ρ(d 2 ) = d 1, ρ(90 ) = 90, ρ(ρ(s)) = s. Suppose d 1 d 2. The inequalities θ i θ 0 yield d 2 2θ 0. Then ρ(2θ 0 ) d 1 90 d 2 2θ 0. Now we consider two cases: 1) ρ(2θ 0 ) d 1 < 77, and 2) 77 d ) Clearly, F 1 (ψ) is a monotone decreasing function in ψ. Then (4) implies f( cosθ 1 ) + f( cosθ 3 ) F 1 (d 1 ) F 1 (ρ(2θ 0 )), so then f( cosθ 2 ) + f( cosθ 4 ) F 1 (d 2 ) = F 1 (ρ(d 1 )) < F 1 (ρ(77 )), H(Y ) < f(1) + F 1 (ρ(2θ 0 )) + F 1 (ρ(77 )) < 13. 2) In this case we have Thus, h 4 < 13. H(Y ) f(1) + F 1 (77 ) + F 1 (90 ) < Our last step is to show that h 3 < Since 3 is a regular triangle, H(Y ) = f(1) + f( cosθ 1 ) + f( cosθ 2 ) + f( cosθ 3 ) is a symmetric function in the θ i, so it is sufficient to consider the case θ 1 θ 2 θ 3 θ 0. In this case R 0 θ 3 θ 0, where R 0 = arccos 2/ is the (spherical) circumradius of 3. 2 A more detailed analysis shows h , h

8 Let y c be the center of 3. We have γ := y 1 y 3 y c = y 2 y 3 y c. Using the low of cosines for the triangle y 1 y 3 y c, we get γ = arccos 2/3, i.e. γ = R 0. Denote the angle e 0 y 3 y c by u. Then (see Fig.5) cosθ 1 = cos60 cosθ 3 + sin 60 sinθ 3 cos(r 0 u), cosθ 2 = cos60 cosθ 3 + sin 60 sinθ 3 cos(r 0 + u), where 0 u u 0 := arccos(cotθ 3 / 3) R 0. Note that if u = u 0, then θ 2 = θ 3 ; u = 0 yields θ 1 = θ 2 ; and if 0 < u < u 0, then θ 1 < θ 2 < θ 3. For fixed θ 3 = ψ, H(y 1, y 2 ) is a polynomial of degree 9 in s = cosu. Denote by F 2 (ψ) the maximum of this polynomial on the interval [cosu 0, 1]. Let {ψ 1,..., ψ 6 } = {R 0, 38, 41, 44, 48, θ 0 }. It is clear that F 2 (ψ) is a monotone increasing function in ψ on [R 0, θ 0 ]. From the other side, f( cosψ) is a monotone decreasing function in ψ. Therefore for θ 3 [ψ i, ψ i+1 ] we have Since, H(Y ) = H(y 1, y 2 ) + f( cosθ 3 ) < w i := F 2 (ψ i+1 ) + f( cosψ i ). {w 1,...,w 5 } { , , , , }, we get h 3 < max{w i } < 13. Thus, h m < 13 for all m as required. Appendix. Proof of Lemma 1. In this proof we are using Schoenberg s original proof [18] which is based on the addition theorem for Gegenbauer polynomials. 3 The addition theorem for Legendre polynomials was discovered by Laplace and Legendre in : where = P k (cosθ 1 )P k (cosθ 2 ) + 2 = P k (cos θ 1 cosθ 2 + sin θ 1 sin θ 2 cosϕ) k m=1 (k m)! (k + m)! P m k (cos θ 1)P m k (cosθ 2) cosmϕ k c m,k Pk m (cos θ 1 )Pk m (cosθ 2 ) cosmϕ, m=0 (See details in [4] and [7].) P m k (t) = (1 t2 ) m 2 d m dt m P k(t). 3 Pfender and Ziegler[17] give a proof as a simple consequence of the addition theorem for spherical harmonics. This theorem is not so elementary. The addition theorem for Legendre polynomials can be proven by elementary algebraic calculations. 8

9 Proof. Let X = {x 1,...,x n } S 2 and x i has spherical (polar) coordinates (θ i, ϕ i ). Then from the law of cosines we have: which yields cosφ i,j = cosθ i cosθ j + sin θ i sin θ j cosϕ i,j, ϕ i,j := ϕ i ϕ j, P k (cosφ i,j ) = i,j i,j = m k c m,k Pk m (cosθ i )Pk m (cosθ j )cosmϕ i,j m=0 c m,k u m,i u m,j cosmϕ i,j, u m,i = Pk m (cosθ i ). i,j Let us prove that for any real u 1,...,u n u i u j cosmϕ i,j 0. i,j Pick n vectors v 1,..., v n in R 2 with coordinates v i = (cosmϕ i, sinmϕ i ). If v = u 1 v u n v n, then 0 v 2 = v, v = i,j u i u j cosmϕ i,j. This inequality and the inequalities c m,k > 0 complete our proof. References [1] M. Aigner and G.M. Ziegler, Proofs from THE BOOK, Springer, 1998 (first ed.) and 2002 (second ed.) [2] K. Anstreicher, The thirteen spheres: A new proof, Discrete and Computational Geometry, 31(2004), [3] K. Böröczky, The Newton-Gregory problem revisited, Proc. Discrete Geometry, Marcel Dekker, 2003, [4] B.C. Carlson, Special functions of applied mathematics, Academic Press, [5] B. Casselman, The difficulties of kissing in three dimensions, Notices Amer. Math. Soc., 51(2004), [6] J.H. Conway and N.J.A. Sloane, Sphere Packings, Lattices, and Groups, New York, Springer-Verlag, 1999 (Third Edition). [7] A. Erdélyi, editor, Higher Transcendental Function, McGraw-Hill, NY, 3 vols, 1953, Vol. II, Chap. XI. 9

10 [8] T. Hales, The status of the Kepler conjecture, Mathematical Intelligencer 16(1994), [9] R. Hoppe, Bemerkung der Redaction, Archiv Math. Physik (Grunet) 56 (1874), [10] W.-Y. Hsiang, Least action principle of crystal formation of dense packing type and Kepler s conjecture, World Scientific, [11] J. Leech, The problem of the thirteen spheres, Math. Gazette 41 (1956), [12] V.I. Levenshtein, On bounds for packing in n-dimensional Euclidean space, Sov. Math. Dokl. 20(2), 1979, [13] H. Maehara, Isoperimetric theorem for spherical polygons and the problem of 13 spheres, Ryukyu Math. J., 14 (2001), [14] O.R. Musin, The problem of the twenty-five spheres, Russian Math. Surveys, 58(2003), [15] O.R. Musin, The kissing number in four dimensions, preprint, September 2003, math. MG/ [16] A.M. Odlyzko and N.J.A. Sloane, New bounds on the number of unit spheres that can touch a unit sphere in n dimensions, J. of Combinatorial Theory A26(1979), [17] F. Pfender and G.M. Ziegler, Kissing numbers, sphere packings, and some unexpected proofs, Notices Amer. Math. Soc., 51(2004), [18] I.J. Schoenberg, Positive definite functions on spheres, Duke Math. J., 9 (1942), [19] K. Schütte and B.L. van der Waerden, Das Problem der dreizehn Kugeln, Math. Ann. 125 (1953), [20] G.G. Szpiro, Kepler s conjecture, Wiley, [21] G.G. Szpiro, Newton and the kissing problem, 10