Kissing number

In geometry, the kissing number of a mathematical space is defined as the greatest number of non-overlapping unit spheres that can be arranged in that space such that they each touch a common unit sphere. For a given sphere packing (arrangement of spheres) in a given space, a kissing number can also be defined for each individual sphere as the number of spheres it touches. For a lattice packing the kissing number is the same for every sphere, but for an arbitrary sphere packing the kissing number may vary from one sphere to another.

Other names for kissing number that have been used are Newton number (after the originator of the problem), and contact number.

In general, the kissing number problem seeks the maximum possible kissing number for n-dimensional spheres in (n + 1)-dimensional Euclidean space. Ordinary spheres correspond to two-dimensional closed surfaces in three-dimensional space.

Finding the kissing number when centers of spheres are confined to a line (the one-dimensional case) or a plane (two-dimensional case) is trivial. Proving a solution to the three-dimensional case, despite being easy to conceptualise and model in the physical world, eluded mathematicians until the mid-20th century. ^{[1] [2]} Solutions in higher dimensions are considerably more challenging, and only a handful of cases have been solved exactly. For others, investigations have determined upper and lower bounds, but not exact solutions. ^[3]

Known greatest kissing numbers

One dimension

In one dimension, ^[4] the kissing number is 2:

Two dimensions

In two dimensions, the kissing number is 6:

Proof: Consider a circle with center C that is touched by circles with centers C₁, C₂, .... Consider the rays CC_i. These rays all emanate from the same center C, so the sum of angles between adjacent rays is 360°.

Assume by contradiction that there are more than six touching circles. Then at least two adjacent rays, say CC₁ and CC₂, are separated by an angle of less than 60°. The segments C C_i have the same length – 2r – for all i. Therefore, the triangle CC₁C₂ is isosceles, and its third side – C₁C₂ – has a side length of less than 2r. Therefore, the circles 1 and 2 intersect – a contradiction. ^[5]

A highly symmetrical realization of the kissing number 12 in three dimensions is by aligning the centers of outer spheres with vertices of a regular icosahedron. This leaves slightly more than 0.1 of the radius between two nearby spheres.

Three dimensions

In three dimensions, the kissing number is 12, but the correct value was much more difficult to establish than in dimensions one and two. It is easy to arrange 12 spheres so that each touches a central sphere, with a lot of space left over, and it is not obvious that there is no way to pack in a 13th sphere. (In fact, there is so much extra space that any two of the 12 outer spheres can exchange places through a continuous movement without any of the outer spheres losing contact with the center one.) This was the subject of a famous disagreement between mathematicians Isaac Newton and David Gregory. Newton correctly thought that the limit was 12; Gregory thought that a 13th could fit. Some incomplete proofs that Newton was correct were offered in the nineteenth century, most notably one by Reinhold Hoppe, but the first correct proof (according to Brass, Moser, and Pach) did not appear until 1953. ^{[1] [2] [6]}

The twelve neighbors of the central sphere correspond to the maximum bulk coordination number of an atom in a crystal lattice in which all atoms have the same size (as in a chemical element). A coordination number of 12 is found in a cubic close-packed or a hexagonal close-packed structure.

Larger dimensions

In four dimensions, it was known for some time that the answer was either 24 or 25. It is straightforward to produce a packing of 24 spheres around a central sphere (one can place the spheres at the vertices of a suitably scaled 24-cell centered at the origin). As in the three-dimensional case, there is a lot of space left over — even more, in fact, than for n = 3 — so the situation was even less clear. In 2003, Oleg Musin proved the kissing number for n = 4 to be 24. ^{[7] [8]}

The kissing number in n dimensions is unknown for n > 4, except for n = 8 (where the kissing number is 240), and n = 24 (where it is 196,560). ^{[9] [10]} The results in these dimensions stem from the existence of highly symmetrical lattices: the E₈ lattice and the Leech lattice.

If arrangements are restricted to lattice arrangements, in which the centres of the spheres all lie on points in a lattice, then this restricted kissing number is known for n = 1 to 9 and n = 24 dimensions. ^[11] For 5, 6, and 7 dimensions the arrangement with the highest known kissing number found so far is the optimal lattice arrangement, but the existence of a non-lattice arrangement with a higher kissing number has not been excluded.

Some known bounds

The following table lists some known bounds on the kissing number in various dimensions. ^[12] The dimensions in which the kissing number is known are listed in boldface.

Rough volume estimates show that kissing number in n dimensions grows exponentially in n. The base of exponential growth is not known. The grey area in the above plot represents the possible values between known upper and lower bounds. Circles represent values that are known exactly.

Dimension	Lower bound	Upper bound
1	2
2	6
3	12
4	24 [7]
5	40	44
6	72	78
7	126	134
8	240
9	306	363
10	510	553
11	592	868
12	840	1,355
13	1,154	2,064
14	1,932	3,174
15	2,564	4,853
16	4,320	7,320
17	5,346	10,978
18	7,398	16,406
19	10,668	24,417
20	17,400	36,195
21	27,720	53,524
22	49,896	80,810
23	93,150	122,351
24	196,560

Generalization

The kissing number problem can be generalized to the problem of finding the maximum number of non-overlapping congruent copies of any convex body that touch a given copy of the body. There are different versions of the problem depending on whether the copies are only required to be congruent to the original body, translates of the original body, or translated by a lattice. For the regular tetrahedron, for example, it is known that both the lattice kissing number and the translative kissing number are equal to 18, whereas the congruent kissing number is at least 56. ^[13]

Algorithms

There are several approximation algorithms on intersection graphs where the approximation ratio depends on the kissing number. ^[14] For example, there is a polynomial-time 10-approximation algorithm to find a maximum non-intersecting subset of a set of rotated unit squares.

Mathematical statement

The kissing number problem can be stated as the existence of a solution to a set of inequalities. Let ${\displaystyle x_{n}}$ be a set of ND-dimensional position vectors of the centres of the spheres. The condition that this set of spheres can lie round the centre sphere without overlapping is: ^[15]

$${\displaystyle \exists x\ \left\{\forall _{n}\{x_{n}^{\textsf {T}}x_{n}=1\}\land \forall _{m,n:m\neq n}\{(x_{n}-x_{m})^{\textsf {T}}(x_{n}-x_{m})\geq 1\}\right\}}$$

Thus the problem for each dimension can be expressed in the existential theory of the reals. However, general methods of solving problems in this form take at least exponential time which is why this problem has only been solved up to four dimensions. By adding additional variables, ${\displaystyle y_{nm}}$ this can be converted to a single quartic equation in N(N − 1)/2 + DN variables: ^[16]

$${\displaystyle \exists xy\ \left\{\sum _{n}\left(x_{n}^{\textsf {T}}x_{n}-1\right)^{2}+\sum _{m,n:m<n}{\Big (}(x_{n}-x_{m})^{\textsf {T}}(x_{n}-x_{m})-1-(y_{nm})^{2}{\Big )}^{2}=0\right\}}$$

Therefore, to solve the case in D = 5 dimensions and N =  40  + 1 vectors would be equivalent to determining the existence of real solutions to a quartic polynomial in 1025 variables. For the D = 24 dimensions and N =  196560  + 1, the quartic would have 19,322,732,544 variables. An alternative statement in terms of distance geometry is given by the distances squared ${\displaystyle R_{mn}}$ between the m^th and n^th sphere:

$${\displaystyle \exists R\ \{\forall _{n}\{R_{0n}=1\}\land \forall _{m,n:m<n}\{R_{mn}\geq 1\}\}}$$

This must be supplemented with the condition that the Cayley–Menger determinant is zero for any set of points which forms a (D + 1) simplex in D dimensions, since that volume must be zero. Setting ${\displaystyle R_{mn}=1+{y_{mn}}^{2}}$ gives a set of simultaneous polynomial equations in just y which must be solved for real values only. The two methods, being entirely equivalent, have various different uses. For example, in the second case one can randomly alter the values of the y by small amounts to try to minimise the polynomial in terms of the y.

See also

• Equilateral dimension
• Spherical code
• Soddy's hexlet
• Cylinder sphere packing

Notes

1. Conway, John H.; Neil J.A. Sloane (1999). Sphere Packings, Lattices and Groups (3rd ed.). New York: Springer-Verlag. p.  21. ISBN   0-387-98585-9.
2. Brass, Peter; Moser, W. O. J.; Pach, János (2005). Research problems in discrete geometry. Springer. p.  93. ISBN   978-0-387-23815-9.
3. Mittelmann, Hans D.; Vallentin, Frank (2010). "High accuracy semidefinite programming bounds for kissing numbers". Experimental Mathematics . 19 (2): 174–178. arXiv: 0902.1105 . Bibcode:2009arXiv0902.1105M. doi:10.1080/10586458.2010.10129070. S2CID   218279.
4. Note that in one dimension, "spheres" are just pairs of points separated by the unit distance. (The vertical dimension of one-dimensional illustration is merely evocative.) Unlike in higher dimensions, it is necessary to specify that the interior of the spheres (the unit-length open intervals) do not overlap in order for there to be a finite packing density.
5. See also Lemma 3.1 in Marathe, M. V.; Breu, H.; Hunt, H. B.; Ravi, S. S.; Rosenkrantz, D. J. (1995). "Simple heuristics for unit disk graphs". Networks. 25 (2): 59. arXiv: math/9409226 . doi:10.1002/net.3230250205.
6. Zong, Chuanming (2008). "The kissing number, blocking number and covering number of a convex body". In Goodman, Jacob E.; Pach, J├ínos; Pollack, Richard (eds.). Surveys on Discrete and Computational Geometry: Twenty Years Later (AMS-IMS-SIAM Joint Summer Research Conference, June 18ÔÇô22, 2006, Snowbird, Utah). Contemporary Mathematics. Vol. 453. Providence, RI: American Mathematical Society. pp. 529–548. doi:10.1090/conm/453/08812. ISBN   9780821842393. MR   2405694..
7. O. R. Musin (2003). "The problem of the twenty-five spheres". Russ. Math. Surv. 58 (4): 794–795. Bibcode:2003RuMaS..58..794M. doi:10.1070/RM2003v058n04ABEH000651. S2CID   250839515.
8. Pfender, Florian; Ziegler, Günter M. (September 2004). "Kissing numbers, sphere packings, and some unexpected proofs" (PDF). Notices of the American Mathematical Society: 873–883..
9. Levenshtein, Vladimir I. (1979). "О границах для упаковок в n-мерном евклидовом пространстве" [On bounds for packings in n-dimensional Euclidean space]. Doklady Akademii Nauk SSSR (in Russian). 245 (6): 1299–1303.
10. Odlyzko, A. M.; Sloane, N. J. A. (1979). "New bounds on the number of unit spheres that can touch a unit sphere in n dimensions". Journal of Combinatorial Theory . Series A. 26 (2): 210–214. doi: 10.1016/0097-3165(79)90074-8 .
11. Weisstein, Eric W. "Kissing Number". MathWorld .
12. Machado, Fabrício C.; Oliveira, Fernando M. (2018). "Improving the Semidefinite Programming Bound for the Kissing Number by Exploiting Polynomial Symmetry". Experimental Mathematics . 27 (3): 362–369. arXiv: 1609.05167 . doi:10.1080/10586458.2017.1286273. S2CID   52903026.
13. Lagarias, Jeffrey C.; Zong, Chuanming (December 2012). "Mysteries in packing regular tetrahedra" (PDF). Notices of the American Mathematical Society: 1540–1549.
14. Kammer, Frank; Tholey, Torsten (July 2012). "Approximation Algorithms for Intersection Graphs". Algorithmica. 68 (2): 312–336. doi:10.1007/s00453-012-9671-1. S2CID   3065780.
15. Numbers m and n run from 1 to N. ${\displaystyle x=(x_{n})_{N}}$ is the sequence of the N positional vectors. As the condition behind the second universal quantifier (${\displaystyle \forall }$) does not change if m and n are exchanged, it is sufficient to let this quantor extend just over ${\displaystyle m,n:m<n}$. For simplification the sphere radiuses are assumed to be 1/2.
16. Concerning the matrix ${\displaystyle y=(y_{mn})_{N\times {N}}}$ only the entries having m < n are needed. Or, equivalent, the matrix can be assumed to be antisymmetric. Anyway the matrix has justN(N − 1)/2 free scalar variables. In addition, there are ND-dimensional vectors x_n, which correspondent to a matrix ${\displaystyle x=(x_{nd})_{N\times D}}$ of N column vectors.

References

• T. Aste and D. Weaire The Pursuit of Perfect Packing (Institute Of Physics Publishing London 2000) ISBN   0-7503-0648-3
• Table of the Highest Kissing Numbers Presently Known maintained by Gabriele Nebe and Neil Sloane (lower bounds)
• Bachoc, Christine; Vallentin, Frank (2008). "New upper bounds for kissing numbers from semidefinite programming". Journal of the American Mathematical Society. 21 (3): 909–924. arXiv: math.MG/0608426 . Bibcode:2008JAMS...21..909B. doi:10.1090/S0894-0347-07-00589-9. MR   2393433. S2CID   204096.

External links

• Grime, James. "Kissing Numbers" (video). youtube. Brady Haran. Archived from the original on 2021-12-12. Retrieved 11 October 2018.