REGULAR POLYTOPES, SPHERE PACKINGS AND APOLLONIAN SECTIONS

IVAN´ RASSKIN

Abstract. In this paper we explore the connection between polytopes and sphere packings. After introducing the notion of Apollonian section, we prove our main result, that every integral tetrahedral, octahedral or cubical Apollonian circle packing is contained as an Apollonian section of an integral or- thoplicial Apollonian section. This result implies that the set of curvatures in every integral tetrahedral, octahedral or cubical Apollonian circle packing is contained in the set of curvatures of an integral ortho- plicial Apollonian sphere packing. We also prove that cross polytopes in dimension ≥ 3 and the 24-cell are M¨obiusunique, we give a new set of non edge-scribable 4-polytopes and we obtain a generalization of the octahedral Descartes’ Theorem of Guettler and Mallows in every dimension.

1. Introduction Sphere packings appear diversely in chemistry, physics and mathematics in crystal structures, granu- lar systems or number theory [UN18; Gaj+21; FSZ19]. There is a long trajectory on the numerical and analytic approaches, but many questions about the algebraic and combinatorial aspects of non-congruent sphere packings are still open. This work is a continuation of [RR21a] where Ram´ırezAlfons´ınand the author introduced a polytopal approach to investigate d-dimensional sphere packings. This yielded, in particular, to a number of results on generalizations of Apollonian circle packings and Descartes’ The- orem based on the Platonic solids. Here, we focus our attention to the third dimensional analogue by investigating the connection between sphere packings and the 4-dimensional regular polytopes.

A particular case of sphere packing based on a regular 4-polytope was given by Nakamura and Dias in two independent works [Nak14; Dia14]. Nakamura called the construction the orthoplicial Apollonian sphere packing and is obtained as follows. Consider a packing of 8 spheres whose tangency graph is the 1-skeleton of an orthoplex. For each subset of four pairwise tangent spheres, there is a unique dual sphere orthogonal to the four. Invert the whole configuration trough every dual sphere, and then repeating this process infinitely (see Figure 1). Nakamura and Dias used the orthoplicial Apollonian sphere packing to prove the corresponding analogue of the local-to-global principle of tetrahedral Apollonian circle packings, conjectured by Graham, Lagarias, Mallow, Wilks and Yan in [Gra+03]. The original conjecture states that for every integral tetrahedral Apollonian packing, every integer large enough satisfying a certain modulo conditions, appears as a curvature. Nowadays, this conjecture is still open [FSZ19]. arXiv:2109.00655v1 [math.CO] 2 Sep 2021

Figure 1. The orthoplicial Apollonian sphere packing at 0, 1 and 2 iterations. Each color represents an orbit.

2010 Mathematics Subject Classification. 51M20, 05B40, 52C17. Key words and phrases. Sphere packings, Apollonian packings, Descartes’ Theorem, Polytopes, Number theory. 1 Cross sections are a natural way to study Apollonian configurations in high dimensions [Bar18]. In order to recognize planar structures in sphere packings, not only geometrically but also arithmetically, we shall present an algebraic generalization of cross sections of Apollonian clusters that we call Apollonian sections. Through this notion we shall prove the following theorem, which implies that the curvatures in every integral tetrahedral, octahedral or cubical Apollonian packing are contained in an integral orthoplicial Apollonian packing (see Figures 2, 15 and 16).

Theorem 1. Let Ω(BP ) be either a tetrahedral, an octahedral or a cubical Apollonian packing. There is an orthoplicial Apollonian packing Ω(BO4 ) containing an Apollonian section arithmetically equivalent to Ω(BP ). Moreover, Ω(BP ) is integral if and only if Ω(BO4 ) is integral.

-3

Figure 2. (Left) An integral tetrahedral Apollonian circle packing Ω(BT 3 ), (center) an integral orthoplicial Apollonian sphere packing Ω(BO4 ), (right) an Apollonian section of Ω(BO4 ) arithmetically equivalent to Ω(BT 3 ).

1.1. Main results. In addition to Theorem 1, our main contributions are the following: (i) New representations of the regular 4-polytopes called Centered Ball Packing projections (Figures 3 and 4). (ii) We prove that cross polytopes in dimension ≥ 3 as well as the 24-cell are M¨obiusunique (Theorem 2 and Proposition 2.2). (iii) We present a new set of non edge-scribable 4-polytopes (Proposition 2.3). (iv) We obtain a generalization of the octahedral Descartes’ Theorem of Guettler and Mallows in every higher dimension (Corollary 2.3). (v) We relate integral polytopal d-ball packings based on the (d + 1)-simplex, (d + 1)-cross polytope and (d + 1)-cube to integers solutions of three diaphontine equations (Corollary 2.5).

1.2. Organization of the paper. In Section 2, after recalling the main tools introduced in [RR21a], we show the Centered Ball Packing projections of the regular 4-polytopes. Then, we extend the results on M¨obiusunicity for the family of cross polytopes in dimension 3 and above and the 24-cell. We briefly discuss the M¨obiusspectra of regular polytopes and use the M¨obiusunicity of the orthoplex, cube and 24-cell to obtain a class of 4-polytopes which are not edge-scribable. After that, we restate the polytopal generalization of Descartes’ Theorem given in [RR21a] in terms of quadratic forms. Through this new approach, we obtain analogues of the equation of Descartes’ Theorem for the family of the simplex, cross polytopes and hypercubes in every dimension, that turns out to be useful to find primitive solutions to three different diophantine equations.

In Section 3, we show that the Nakamura’s definition of orthoplicial sphere packings are in fact polytopal. We then discuss some aspect concerning the duality of orthoplicial sphere packings and we introduce the notion of orthoplicial trinity. Then, we compute the symmetrized Apollonian group and use it to give an alternative matrix representation of the orthoplicial Apollonian group.

Finally, in Section 4, we focus our attention to integral packings. After introducing the notion of Apollonian sections of Apollonian clusters we prove the main theorem. We end with some conjectures and concluding remarks. 2 2. Preliminaries on polytopal d-ball packings. The construction given in [RR21a] to obtain d-ball packings from (d + 1)-polytopes is described as follows. Take a (d+1)-polytope P whose vertices lie outside the unit sphere of the Euclidean space Ed+1. Then, put at each vertex of P a “light source” illuminating the sphere. By projecting stereoreographi- cally the illuminated areas, one obtains a d-ball arrangement in Rcd called the ball-arrangement projection of P, denoted by B(P). If P is edge-scribed, that is, with every edge tangent to the sphere, then B(P) is d-ball packing. In this case, every d-ball packing BP equivalent to B(P) under M¨obiustransformations is said to be polytopal and P is called the tangency polytope of BP . The tangency graph of BP corresponds to the 1-skeleton of P [Che16]. We notice that there are d-ball packings which are not polytopal [RR21a].

One of the main features of polytopal d-ball packings is that they admit a consistent definition of duality in every dimension, contrarily to d-ball packings where the combinatorics is given by a graph or a simplicial complex. If µ is the Mˆobiustransformation sending B(P) to BP , then the dual of BP , ∗ ∗ ∗ denoted by BP , is given by µ(B(P )) where P denotes the polar of P. Duality is then needed to general- ize the Apollonian group and its variants to every polytopal d-ball packing (see [RR21a] for more details).

We say that a polytopal d-ball packing is regular if its tangency polytope is regular. Let us recall the list of regular polytopes in every dimension. All 2-polytopes (polygons) admits a regular realization. The Platonic solids, namely, the tetrahedron T 3, the octahedron O3, the cube C3, the icosahedron I3 and the dodecahedron D3 are the five regular 3-polytopes (polyhedra). The 4-simplex T 4, the orthoplex O4, the hypercube C4, the 600-cell I4 and the 120-cell D4 are five regular 4-polytopes which can be thought as a 4-dimensional analogue of the Platonic solids. The remaining regular 4-polytope, the 24-cell R4 (the notation is not standard), completes the list of regular 4-polytopes. For every d ≥ 2, we shall denote by T d+1, Od+1 and Cd+1 the (d + 1)-dimensional simplex, octahedron and cube. It is well-known that in dimension 5 or above, these three families are the only regular polytopes [Cox73]. The following table contains all the information needed for the paper about the regular (d + 1)-polytopes for d ≥ 1.

d Name Notation Schl¨aflisymbol Eigenvalues Multiplicities Midsphere ratio d = 1 p-gone (p ≥ 3) - {p} Not M¨obius unique for p > 3 tan(π/p) {3,..., 3, 3} −d 1 q d ≥ 1 (d + 1)-Simplex T d+1 d+2 | {z } d d−1 2 d + 1 −2(d + 1) 1 d+1 {3,..., 3, 4} (d + 1)-Cross polytope O | {z } d + 1 14 d−1 0 d d ≥ 2 −2d+1d 1 d+1 {4, 3,..., 3} d+1 −1/2 (d + 1)-Cube C | {z } 2 d + 1 d d−1 0 2d+1 − d − 2 −12ϕ2 1 Icosahedron I3 {3, 5} 4(ϕ2 + 1) 3 ϕ−1 0 8 d = 2 −20ϕ4 1 Dodecahedron D3 {5, 3} 20ϕ2 3 ϕ−2 0 16 −72 1 24-cell R4 {3, 4, 3} 24 4 3−1/2 0 19 d = 3 600-cell I4 {3, 3, 5} M¨obiusunique? 5−1/4ϕ−3/2 120-cell D4 {5, 3, 3} M¨obiusunique? 3−1/2ϕ−3

Table 1. Notations, M¨obiusspectra and midsphere ratio of regular (d + 1)-polytopes for d ≥ 1.

For the latter computations on the space of d-balls, we shall use the inversive coordinates and the inversive product of Wilker [Wil81]. We refer the reader to [RR21a, Section 2] and the reference therein for more details about how to compute center, curvatures and inversions on d-balls in inversive coordinates. We shall refer polytopal 2 (resp. 3)-ball packings as polytopal disk (resp. sphere) packings. 3 2.1. CBP projections of 4-polytopes. For every d ≥ 1, a Centered Ball Packing projection of a regular (d + 1)-polytope P is a polytopal d-ball packing obtained by the ball-arrangement projection of an edge-scribed regular realization of P containing a i-face fi centered in the ray going from the origin to the North Pole of Sd (∞ in Rcd). Such a packing will be called a i-CBP projection of P. The CBP projections of the Platonic solids were presented in [RR21a]. We show in Figures 3 and 4 the CBP projections of the six regular 4-polytopes.

The curvatures in each layer in a i-CBP projection can be expressed as a linear combination of two numbers κP and hi, where κP is the mean of all the curvatures, and hi is the minimum of the positive −1 heights (in the direction of the vector (0, 0, 0, 1)), among all the vertices of P. We notice that κP = `P where `P is the midsphere ratio of P (see [RR21a, Lemma 3]).

Vertex centered at ∞ Edge centered at ∞ Ridge centered at ∞ Facet centered at ∞

T 4 4-Simplex {3, 3, 3}

p κP = 3/5 Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures p 3 h0 = 1/10 1 1 κP − 4h0 1 2 κP − 2 h1 1 3 κP − h2 1 4 κP − h3 p 3 h1 = 2/15 2 4 κP + h0 2 3 κP + h1 2 2 κP + 2 h2 2 1 κP + 4h3 p h2 = 2/15 p h3 = 1/10

O4 Orthoplex {3, 3, 4}

κP =√ 1 Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures h0 = 2 1 1 κP − h0 1 2 κP − h1 1 3 κP − h2 1 4 κP − h3 h1 = 1 2 4 κP 2 4 κP 2 2 κP 2 4 κP + h3 p h2 = 2/3 3 1 κP + h0 3 2 κP + h1 3 3 κP + h2 p h3 = 1/2

C4 Hypercube {4, 3, 3}

√ κP = 3 Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures h0 = 1 1 1 κP − 2h0 1 2 κP − 3h1 1 4 κP − h2 1 8 κP − h3 p 2 4 κ − h 2 6 κ − h 2 8 κ 2 8 κ + h h1 = √ 1/3 P 0 P 1 P P 3 h2 = 2 3 6 κP 3 6 κP + h1 3 4 κP + h2 h3 = 1 4 4 κP + h0 4 2 κP + 3h1 5 1 κP + 2h0

Figure 3. CBP projections of the 4-simplex, orthoplex and hypercube.

4 R4 24-cell {3, 4, 3}

√ κP = 3 Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures h0 = 1 1 1 κ − 2h0 1 2 κ − 3h1 1 3 κ − 2h2 1 6 κ − h3 p P P P P h1 = 1/3 2 8 κP − h0 2 3 κP − 2h1 2 6 κP − h2 2 12 κP p 3 6 κP 3 6 κP − h1 3 6 κP 3 6 κP + h3 h2 = √2/3 4 8 κP + h0 4 2 κP 4 6 κP + h2 h3 = 2 5 1 κP + 2h0 5 6 κP + h1 5 3 κP + 2h2 6 3 κP + 2h1 7 2 κP + 3h1

I4 600-cell {3, 3, 5}

√ 3/2 κP = 5ϕ Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures 3 4 4 h0 = 1 1 1 κP − 2ϕh0 1 2 κ − (ϕ + ϕ)h 1 3 κ − (ϕ + ϕ)h 1 4 κ − ϕ h 1 P 1 P 2 P 3 2 − 2 12 κ − ϕ2h 3 4 3 h1 = (ϕ + 1) 2 P 0 2 5 κP − (ϕ + 1)h1 2 2 κP − (ϕ + 1)h2 2 4 κP − (ϕ + ϕ)h3 3 20 κ − ϕh 3 4 3 − 1 P 0 3 10 κ − ϕ h 3 6 κ − ϕ h 3 6 κ − (ϕ + 1)h h = (ϕ4 + 1) 2 4 12 κ − h P 1 P 2 P 3 2 P 0 4 2 κ − (ϕ2 + 1)h 4 6 κ − (ϕ3 + ϕ)h 4 12 κ − ϕ3h − 1 5 30 κ P 1 P 2 P 3 h = (ϕ3 + 1) 2 P 3 3 2 3 6 12 κP + h0 5 5 κP − (ϕ − 1)h1 5 6 κP − (ϕ + 1)ϕh2 5 12 κP − ϕ h3 2 3 7 20 κP + ϕh0 6 10 κP − ϕ h1 6 6 κP − ϕ h2 6 12 κP − ϕh3 2 3 7 4 κP − h3 8 12 κP + ϕ h0 7 10 κP − ϕh1 7 3 κP − (ϕ − 1)h2 9 1 κ + 2ϕh 8 10 κP − h1 2 8 12 κP P 0 8 12 κP − ϕ h2 9 12 κP 9 4 κP + h3 9 6 κP − ϕh2 10 10 κP + h1 . . . 10 6 κP − h2 ...... 11 8 κP ...... 12 6 κ + h 4 . . . P 2 15 4 κP + ϕ h3 17 2 κ + (ϕ3 + ϕ)h . . . P 1 ...... 4 21 3 κP + (ϕ + ϕ)h2

D4 120-cell {5, 3, 3}

3√ κP = ϕ 3 Layer n Curvatures Layer n Curvatures Layer n Curvatures Layer n Curvatures p  5 −1 6 2 5 2 4 h0 = 1/2 1 1 κ − ϕ − ϕ h 1 2 κP − (ϕ + ϕ )h1 1 5 κP − (ϕ + ϕ )h2 1 20 κP − ϕ h3 1 P 0 4 −   2 6 κ − (ϕ6 + ϕ)h 2 10 κ − (ϕ5 + ϕ)h 2 20 κ − (ϕ3 + ϕ)h h1 = (ϕ + 1) 2 2 4 κ − ϕ5 − 1 h P 1 P 2 P 3 1 P 0 3 3 κ − (ϕ6 + 1)h 3 10 κ − (ϕ5 + 1)h 3 30 κ − (ϕ3 + 1)h 2 −   P 1 P 2 P 3 h2 = (ϕ + 1) 2 4 2 6 5 3 3 12 κP − ϕ + ϕ h0 4 12 κP − ϕ h1 4 20 κP − ϕ h2 4 60 κP − ϕ h3 h3 = 1 6 5 2  4  5 6 κP − (ϕ − 1)h1 5 10 κP − (ϕ − 1)h2 5 60 κP − ϕ h3 4 24 κP − ϕ + ϕ h0 6 4 2 6 12 κ − (ϕ − ϕ)h1 6 20 κ − (ϕ + ϕ )h2 6 60 κP − ϕh3  4  P P 5 12 κ − ϕ + 1 h 5 3 4 7 20 κP − h3 P 0 7 18 κP − (ϕ + ϕ )h1 7 20 κP − (ϕ + ϕ)h2   5 2 4 8 60 κP 4 −1 8 12 κ − (ϕ + ϕ )h1 8 10 κ − (ϕ + 1)h2 6 4 κP − ϕ + ϕ h0 P P 9 20 κP + h3 9 14 κ − (ϕ5 + ϕ)h 9 30 κ − ϕ4h 7 24 κ − ϕ4h P 1 P 2 . . . P 0 5 3 . . .  3  10 12 κP − (ϕ + 1)h1 10 20 κP − (ϕ + ϕ)h2 8 24 κP − ϕ + ϕ h0 5 3 . . . 11 18 κP − ϕ h1 11 20 κP − (ϕ + 1)h2 4  3  5 3 17 20 κP + ϕ h3 9 32 κP − ϕ + 1 h0 12 6 κP − (ϕ − 1)h1 12 30 κP − ϕ h2 4 2 3 3 13 24 κP − (ϕ + ϕ )h1 13 5 κP − (ϕ − 1)h2 10 24 κP − ϕ h0 4 2  2  14 15 κP − (ϕ + ϕ)h1 14 30 κP − ϕ h2 11 12 κP − ϕ + 1 h0 4 15 30 κ − ϕh   15 2 κP − (ϕ + 1)h1 P 2 2 −1 4 16 20 κP − h2 12 24 κP − ϕ + ϕ h0 16 24 κP − ϕ h1 2 3 17 20 κP 13 28 κP − ϕ h0 17 18 κP − (ϕ + ϕ)h1 18 20 κ + h 3 P 2 15 24 κP − ϕh0 18 12 κP − (ϕ + 1)h1 16 24 κ − h 3 . . . P 0 19 18 κP − ϕ h1 . . . 17 54 κP 2 . . . 20 24 κP − ϕ h1 5 2 18 24 κP + h0 35 5 κP + (ϕ + ϕ )h2 21 18 κP − ϕh1 . . . 22 12 κ − h . . . P 1 . . . 23 24 κP  5 −1 24 12 κP + h1 33 1 κP + ϕ − ϕ h0 ...... 6 2 47 2 κP + (ϕ + ϕ )h1 Figure 4. CBP projections of the 24-cell, 600-cell and 120-cell. 5 2.2. M¨obiusunicity of regular polytopes. An edge-scribable (d+1)-polytope P is said to be M¨obius unique if the ball-arrangement projections of all its edge-scribed realizations are M¨obiusequivalent, i.e. connected by a M¨obiustransformation. Equivalently, P is M¨obiusunique if all its edge-scribed realizations are connected by a projective transformation which preserves the unit sphere of Ed+1. As noticed in [RR21a], we have that: 1. The 2-simplex is the only M¨obiusunique 2-polytope [RR21a, Corollary 2.3]. 2. All the 3-polytopes are M¨obius unique (by the Midsphere theorem of Brightwell and Scheinerman [BS93]). 3. For every d ≥ 2, the d-simplex and the (d + 1)-cube are M¨obiusunique [RR21a, Corollaries 2.2 and 2.4]. Here we shall show that 24-cell and the family of cross polytopes, in dimension 3 or above, are also M¨obius unique. For the proof of the M¨obiusunicity of the cross polytopes, we need the following proposition given in [RR21a]. We recall that a d-ball packing has maximal rank if the rank of its Gramian is d + 2. In particular, polytopal d-ball packings have always maximal rank. Proposition 2.1 ([RR21a]). Let B and B0 two d-ball packings with maximal rank. Then B is M¨obius equivalent to B0 if and only if Gram(B) = Gram(B0). Theorem 2. For every d ≥ 2, the (d + 1)-octahedron is M¨obiusunique. Proof. We proceed similarly to the proof of the M¨obiusunicity of the (d + 1)-cube in [RR21a]. Let d+1 d ≥ 2 and let BOd+1 be a polytopal d-ball packing where O is an edge-scribed realization of the (d + 1)-octahedron. We shall prove by induction on d, that for every two vertices u, v of Od+1, we have

hbu, bvi = 1 − 2dG(u, v)(1) d+1 where dG(u, v) is the distance between u and v in the 1-skeleton of O . This would imply that for any d+1 other polytopal d-ball packing B d+1 , where O is any other edge-scribed realization of the (d + 1)- O2 2 octahedron, we can find an ordering such that Gram(BOd+1 ) = Gram(B d+1 ). The M¨obiusunicity then O2 follows from Proposition 2.1.

The case d = 2 can be easily checked in a single octahedral disk packing. Since, 3-polytopes are M¨obiusunique, then (1) holds for every edge-scribed realization of O3. Let us now suppose that (1) holds for any edge-scribed realization of the (d+1)-octahedron for some d ≥ 2. Let BOd+2 be a polytopal (d + 1)-ball packing where Od+2 ⊂ Ed+2 is an edge-scribed realization of the (d + 2)-octahedron. We give to the vertices of Od+2 an antipodal labelling d+2 V (O ) = {v1, . . . , vd+2, v−1, . . . , v−(d+2)} d+1 where vi and vj are connected by an edge of O if and only if j 6= −i. For every 1 ≤ i < j ≤ d + 1, we consider the following collection of (d + 1)-balls

Bi,j := {b1, . . . , bd+1, b−i, b−j} ⊂ BOd+2 d+2 d+2 d+2 where bk := bvk . Since BO is polytopal, the tangency graph of BO is the 1-skeleton of O . Therefore, bd+2 and b−(d+2) are tangent to every bk ∈ Bi,j. Translating the tangency conditions with the inversive product (see [RR21b]), we have that bd+2 and b−(d+2) satisfies

(2) hbk, bi = −1 for every bk ∈ Bi,j,

and b ∈ {bd+2, b−(d+2)}. In inversive coordinates, (2) becomes the following linear system

(3) Bi,jQd+3X = −1d+3

where Bi,j is the matrix of the inversive coordinates of Bi,j, Qd+3 is the matrix of the inversive product, X is a (d + 3)-column vector, and 1d+3 is the (d + 3)-column vector of only 1’s. Since bd+2 and b−(d+2) are distinct, (3) has more than one solution. Therefore, Bi,j is singular, which implies that there is a d+2 hyperplane Hi,j of E such that

Vi,j := {v1, . . . , vd+1, v−i, v−j} ⊂ Hi,j. 0 d+2 Moreover, since, for every 1 ≤ j ≤ d + 1, the hyperplanes Hi,j and Hi,j0 share d + 2 points of E , then they must be the same hyperplane. Therefore, there is one hyperplane H containing all the vertices d+2 V (O ) \{vd+2, v−(d+2)}.

6 We thus can find a M¨obiustranformation µ ∈ M¨ob(Sd+1) which the corresponding projective trans- d+2 d+1 formation sends H to the hyperplane {xd+2 = 0} ⊂ E . After identifying µ(H) with E we obtain that µ(H ∩Od+2) becomes an edge-scribed realization of the (d+1)-octahedron Od+1. The identification µ(H) ' Ed+1 preserves the inversive product of the (d + 1)-balls corresponding to the points lying in H. Moreover, the distance between u and v in the graph of Od+2 is equal to the distance in the graph of Od+1. By the invariance of the inversive product under M¨obiustransformations and the induction hypothesis, we have that (1) holds for any two (d + 1)-balls of BOd+2 \{bd+2, b−(d+2)}.

The same arguments work if we exchange b1 by bd+2 in Bi,j, so (1) holds in the remaining cases.  Proposition 2.2. The 24-cell is M¨obiusunique. Proof. It is well-known that the 1-skeleton of a 24-cell admits a 3-coloring such that the vertices of each color span an orthoplex. By taking the vertices of two out of the three colors, we span a hypercube. Moreover, the edges of this hypercube are also edges of the initial 24-cell. Therefore, every polytopal sphere packing whose tangency polytope is the 24-cell contains a hypercubical sphere packing.

Let BR4 be the 0-CBP projection of the 24-cell. The even layers of BR4 gives a 3-CBP projection BC4 of the hypercube (see Figure 3). Let BP be a polytopal sphere packing where P is another edge-scribed realization of the 24-cell. Let Q ⊂ P be one of the hypercubes contained in P. Since the hypercube is M¨obiusunique, there is a M¨obiustransformation µ sending the packing BQ to BC4 . Since M¨obius transformations preserve the tangency relation, every sphere b ∈ BR4 \BC4 must be tangent to 6 spheres 4 of BC4 corresponding to the vertices of a facet of C . This condition forces µ to send the spheres in BP \BQ to BR4 \BC4 implying that BP and BR4 are M¨obiusequivalent.  2.3. M¨obiusspectra. The M¨obiusspectra of an edge-scribable M¨obiusunique (d + 1)-polytope P, introduced in [RR21a], is a spectral invariant defined as the spectrum of the Gramian of the ball- arrangement projection of an edge-scribed realization of P. We believe that M¨obius spectra may shed light on the properties of edge-scribable polytopes. We show in Table 1 the M¨obiusspectra of the regular polytopes that are known to be M¨obiusunique. 2.4. Non edge-scribable 4-polytopes. The first examples of 4-polytopes not admitting an edge- scribed realization were given by Schulte in [Sch87]. In [EKZ02, Corollary 9], Eppstein, Kupperberg and Ziegler extended the list by noticing that every stacked 4-polytope with more than 6 vertices is not edge-scribable. We recall that a stacked polytope is a polytope obtained by applying consecutively connected sums of simplices, where the connected sum of two polytopes P and P0 along a facet f, is 0 0 another polytope P#f P whose faces are the union of the faces of P and P minus f. The main ingredient for the proof of this impossibility result is the following. Lemma 1. The consecutive connected sum of three 4-simplices is not edge-scribable. 0 00 0 00 Proof. We consider the 4-polytope (T #f1 T )#f2 T where T , T and T are three 4-simplices. The 0 00 facets f1 and f2 must intersect in a common ridge r (a triangle) of T , T and T . We label the vertices by V (r) = {1, 2, 3}, V (T ) = {1, 2, 3, 4, 5}, V (T 0) = {1, 2, 3, 5, 6} and V (T 00) = {1, 2, 3, 6, 7}. Let us 0 00 suppose that (T #f1 T )#f2 T admits an edge-scribed realization P. By a applying the ball-arrangement projection to P and then an inversion on a sphere centered at the tangency point of the spheres b1 and b2 we obtain a polytopal sphere packing BP Euclidean congruent to the packing depicted at Figure 5.

b1 b2

b5 b6

b4 b3 b7 bf0

Figure 5. The sphere packing BP obtained by inverting the ball-arrangement projection of a glueing of three 4-simplices.

7 In this packing, the dual sphere bf0 , corresponding to the facet f0 of T1 with vertices {1, 2, 3, 4}, cuts orthogonally the sphere b7. This implies that the vertex v7 of T3 lies in the affine hull of f0, so f0 is not a 0 00 face of (T #f1 T )#f2 T . This contradicts the condition on the set of faces in the definition of connected 0 00 sum. Therefore, (T #f1 T )#f2 T is not edge-scribable.  Chen mentioned in [Che16, Section 5] similar arguments as above, by considering a 2-CBP projection of the 4-simplex instead of a 1-CBP projection. A natural generalization of stacked polytopes are the stacked P-polytopes, introduced by Chen and Padrol in [CP17], as polytopes obtained by connected sums of several copies of a given polytope P. With the new results on M¨obiusunicity and the following lemma given in [RR21b], we can generalize the construction of Eppstein, Kuperberg and Ziegler by applying the same arguments to the stacked P-polytopes, where P = O4, C4, R4. We recall that a d-ball packing B is i said to be standard if it contains the half-spaces bi = {xd ≥ 1} and bj = {xd ≤ −1}, denoted by [B]j. Lemma 2. [RR21b, Lemma 2] Let B and B0 be two d-ball packings with same tangency graph G and 0 i 0 i let ij be an edge of G. Then B and B are M¨obiusequivalent if and only if [B]j and [B ]j are Euclidean congruent. Proposition 2.3. The following 4-polytopes are not edge-scribable: • The connected sum of two orthoplexes. • The consecutive connected sum of three hypercubes sharing a ridge. • The consecutive connected sum of three 24-cell sharing a ridge. Proof. We may apply the same arguments as those in the proof of Lemma 1. By combining the M¨obius unicity of the orthoplex, hypercube and 24-cell and Lemma 2, we have that after applying the corre- sponding inversion and a proper rescaling, we must obtain a packing which is Euclidean similar to a glueing by reflections of 1-CBP projections of the corresponding 4-polytope. Then, as above, we would have a dual sphere cutting orthogonally spheres of the different components in the connected sum, and the same contradiction arises.  2.5. The Polytopal Descartes’ Theorem. In [RR21a], a generalization of the Descartes’ theorem for regular polytopal d-ball packings was used to construct integral Apollonian packings of the Platonic solids. In this section, we restate this result with quadratric forms which will be useful to obtain further consequences. Let us first define and recall the terminology needed for the rest of the paper.

We define the polytopal curvatures of a polytopal d-ball packing BP as the real numbers 1 X (4) κ := κ(b ). f |V (f)| v v∈V (f) where f is a face of P, V (f) is the set of vertices of f and κ(bv) is the curvature of the d-ball bv ∈ BP . From the point of view of the polytope, these curvatures correspond to the Lorentzian curvatures of the faces of P, which were defined in [RR21a] by the following Lorentzian product:

(5) κf = −hxN , xf i

d+1,1 where xN is the vector ed+1 + ed+2 ∈ L and xf is the Lorentzian barycenter of f, that is, 1 X (6) x := x f |V (f)| bv v∈V (f)

where xbv is the Lorentzian vector corresponding to bv. The polytopal curvature of a vertex v is exactly the curvature of the corresponding d-ball bv.

The next lemma applies to every centrally symmetric regular (d + 1)-polytopes. These are the p-gons with p even when d = 1 and all the regular polytopes not belonging to the simplex family when d ≥ 2. The case d = 2 was also given in [RR21a].

Lemma 3 (Antipodal relation). Let BP be polytopal d-ball packing where P is a regular edge-scribed (d + 1)-polytope which is centrally symmetric. Then, for any two vertices v, v¯ at maximal distance in the 1-skeleton of P, we have κ + κ κ = v v¯ . P 2 8 1 Proof. Since P is centrally symmetric, then 2 (xbv + xbv¯ ) is the Lorentzian barycenter of P. The Lemma follows from (5) and linearity. 

We define the flag quadratic form of a regular (d + 1)-polytope P as

d 2 2 X (xi − xi+1) x (7) Φ (x , x , . . . , x ) := + d+1 P 0 1 d+1 L (i) − L (i + 1) L (d + 1) i=0 P P P where

 −1 i = 0  LP (i) := 0 i = 1  `−2 if 2 ≤ i ≤ d + 1 fi

where `fi denotes the half edge-length of a regular edge-scribed realization of fi. We call this value the midsphere ratio of a regular polytope since it can be obtained from every regular realization P (not necessarily edge-scribed) by the ratio between the half edge-length of P and the radius of the midsphere (sphere tangent to every edge) of P. The midsphere ratios of every regular polytope are given in the Table 1 and were adapted from [Cox73]. The flag quadratic form of the (d + 1)-simplex, (d + 1)-cross polytope and (d + 1)-cube, for every d ≥ 1, are

d   X i + 2 2 d + 2 2 (8) Φ d+1 (x , . . . , x ) = − (x − x ) + x T 0 d+1 2 i i+1 d d+1 i=0 d−1   X i + 2 2 d + 1 2 2 (9) Φ d+1 (x , . . . , x ) = − (x − x ) − (x − x ) + x O 0 d+1 2 i i+1 2 d d+1 d+1 i=0 d X 2 1 2 (10) Φ d+1 (x , . . . , x ) = − (x − x ) + x . C 0 d+1 i i+1 d d+1 i=0

In terms of the flag quadratic form, the polytopal generalization of the Descartes’ Theorem given in [RR21a, Theorem 4.1] is restated as follows.

Theorem 3 (Polytopal Descartes’ Theorem, [RR21a]). For d ≥ 1, let BP a polytopal d-ball packing where P is a regular edge-scribed (d + 1)-polytope. Then, for any flag (v, . . . , f, P), the polytopal curvatures κv, . . . , κf , κP of BP satisfies

(11) ΦP (κv, . . . , κf , κP ) = 0.

Corollary 2.1. Let BP+ and BP− two regular polytopal d-ball packings where one is obtained from the other by the inversion on a dual d-ball bf . Then,

 2 r `f −2  2 2  (12) κP± = κf ± `P `f − `P Φf (κv, . . . , κf ). `P

Proof. It follows from the definition (7) that

 2 2  2 `f xd − `P xd+1 (13) ΦP (x0, x1, . . . , xd+1) = Φf (x0, x1, . . . , xd) − 2 2 `f − `P

By combining this with (11) and then resolving for κP we obtain (12).  9 Let us explain in detail why Theorem 3 generalizes the Descartes’ Theorem. We define the simplicial, hyperoctahedral and hypercubical quadratic form as  1 1  (14) T (u , . . . , u ) :=Φ d+1 u , (u + u ),..., (u + ... + u ) d+1 1 d+2 T 1 2 1 2 d + 2 1 d+2 d+2 d+2 ! 1 1 X X = ( u )2 − u2 2 d i i i=1 i=1  1 1  (15) O (u , . . . , u ) :=Φ d+1 u , (u + u ),..., (u + ... + u ), u d+1 1 d+2 O 1 2 1 2 d + 1 1 d+1 d+2 d+1 1 X =u2 − (u − u )2 d+2 2 i d+2 i=1  1 1  (16) C (u , . . . , u ) :=Φ d+1 u , (u + u ),..., (u + u ) d+1 1 d+2 C 1 2 1 2 2 1 d+2 d+1 ! 1 1 X = (u + u )2 − (u − u )2 4 d 1 d+2 i i+1 i=1

Corollary 2.2 (Soddy-Gosset Theorem). Let κ1, . . . , κd+2 be the curvatures of a polytopal d-ball packing BT d+1 . Then, d+2 d+2 X 2 X 2 (17) d κi = ( κi) i=1 i=1 Proof. Since T d+1 is M¨obiusunique, the Polytopal Descartes’ Theorem holds for any polytopal d-ball d+1 packing BT d+1 . For every i = 0, . . . , d + 1, let vi be the vertex of T corresponding to the d-ball of d+1 curvature κi. Since T is a (d + 1)-neighborly polytope, that is, every set of vertices span a face, we d+1 can find a flag (f0, . . . , fd, fd+1 = T ) where the vertices of fi are v0, . . . , vi. Therefore, we have 1 (18) κ = (κ + ... + κ ). fi i 1 i Then, by the Polytopal Descartes’ Theorem,

d+1 Td+1(κ1, . . . , κd+2) = ΦT (κf0 , . . . , κfd+1 ) = 0 which is equivalent to (17). 

Remark 1. In the above corollary, “polytopal d-ball packing BT d+1 ” can be replaced by a “d-ball packing made by d+2 pairwise tangent d-balls”, as in the original statement. By Proposition 2.1, both definitions are equivalent.

Corollary 2.3. Let κ1, . . . , κd+1 be the curvatures of d + 1 pairwise tangent d-balls of a polytopal d-ball packing BOd+1 . Then, d+1 X 2 2 (19) (κi − κOd+1 ) = 2κOd+1 i=1 Proof. We can apply the same arguments as in the proof of Corollary 2.2. The vertices corresponding to the curvatures are the vertices of a facet of Od+1, which is a d-simplex. Therefore, we can find a a d+1 flag (f0, . . . , fd, O ) where (18) is satisfied for every i = 0, . . . , d. The Polytopal Descartes’ Theorem combined with hyperoctahedral quadratic form give the result.  Remark 2. Equation (19) for d = 2, 3 is equivalent to the equation given in the generalizations of Descartes’ Theorem of Guettler-Mallows [GM08], Nakamura [Nak14] and Dias [Dia14], respectively. The extra equations needed to obtain the remaining curvatures correspond to the antipodal relation.

Corollary 2.4. For every d ≥ 1, let κ1, . . . , κd+2 be the curvatures of d + 2 consecutive tangent d-balls of a polytopal d-ball packing BCd+1 where κ1 and κd+2 are the curvatures of two d-balls at distance d + 1 in the tangency graph of BCd+1 . Then we have d+1 X 2 2 (20) d (κi − κi+1) = (κ1 + κd+2) . i=1 10 d+1 Proof. By M¨obiusunicity, we can consider that C is regular. Let (v1, . . . , vd+2) a path passing through d+1 the vertices of C corresponding to the curvatures. There is a flag (f0, . . . , fd+1) where fi is the unique d i-face containing the vertices v1, . . . , vi. For every i = 2, . . . , d+1, the intersection fi ∩S gives a polytopal i-ball packing BCi . By applying the antipodal relation to BCi+1 we obtain that, for every i = 1, . . . , d + 2, 1 (21) κ = (κ + κ ). fi 2 1 i and therefore, by the Polytopal Descartes’ Theorem,

d+1 Cd+1(κ1, . . . , κd+2) = ΦC (κf0 , . . . , κfd , κfd+1 ) = 0 which is equivalent to (20).  In the case when all the curvatures in the d-ball packing in Corollaries 2.2, 2.3 and 2.4 are integers, then we obtain a geometric method to find solutions to three Diophantine equations.

Corollary 2.5. Let d ≥ 1. If there is an integral polytopal d-ball packing BT d+1 , BOd+1 or BCd+1 then the following Diophantine equations 2 2 2 (22) d(m1 + ··· + md+2) = n 2 2 2 (23) m1 + ··· + md+1 = 2n 2 2 2 (24) d(m1 + ··· + md+1) = n , admit interger solutions. Proof. Equations (22), (23) and (24) are obtained by additioning the numbers inside the parenthesis in the equalities (17), (19) and (20), respectively. 

10 7 90 89 3 88 12 59 64 15 67 62 22 9 8 15 10

2(22 + 22 + 32 + 152) = 222 392 + 402 + 412 = 2 × 492 2(32 + 52 + 522) = 742

Figure 6. Three primitive solutions to Equations (22), (23), (24) for d = 2, obtained by the relations between the curvatures in a integral tetrahedral, octahedral and cubical disk packing, respectively.

3. The orthoplicial sphere packing. In two independent works, Nakamura [Nak14] and Dias [Dia14] defined the orthoplicial sphere packings as a class of packings whose tangency graph is the 1-skeleton of an orthoplex. These packings were also studied by Sheydvasser in [A19] as a packing arising from a quaternary quadratic field. They can also be obtained as a particular case of Boyd-Maxwell packing [CL15] or a crystallographic sphere packing [KN19]. In this section, we shall study orthoplicial sphere packings as polytopal sphere packings, i.e. as packings M¨obiusequivalent to the ball-arrangement projection of an edge-scribed orthoplex. As we show below, orthoplicial sphere packings, in the sense of Nakamura and Dias, are polytopal. Lemma 4. Orthoplicial sphere packings are polytopal. 11 Proof. For any orthoplicial sphere packing B, the tancency conditions gives all the entries of Gram(B) except for the entries corresponding to the inversive product of disjoint spheres, which must be strictly less than −1. Therefore, the rank of Gram(B) must be equal to 5, implying that the entries of disjoint spheres must be −3. By Equation (1), we can reorder the spheres of B such that Gram(B) = Gram(BO4 ), where BO4 is the ball-arrangement projection of an edge-scribed orthoplex. By Proposition 2.1, B and BO4 are M¨obiusequivalent, and therefore B is polytopal.  ∗ 3.1. Orthoplicial trinities. By polarity, the dual BO4 of an orthoplicial sphere packing is M¨obius ∗ equivalent to the ball-arrangement projection of a ridge-scribed hypercube. Therefore, BO4 is not a packing. However, by alternating the vertices of the hypercube, we can split the dual in two orthoplicial ∗ 0 00 0 00 sphere packings BO4 = BO4 ∪BO4 . Such a collection of three orthoplicial sphere packings {BO4 , BO4 , BO4 } will be called an orthoplicial trinity. 0 00 0 00 Lemma 5. Let {BO4 , BO4 , BO4 } be an orthoplicial trinity. Then, for any B ∈ {BO4 , BO4 , BO4 }, we have ∗ 0 00 that B = {BO4 , BO4 , BO4 }\B. Proof. By the M¨obiusunicity of the orthoplex, combined with the fact that M¨obiustransformations preserves duality, implies that it is enough to prove the result in a particular case. Let us consider the standard orthoplicial sphere packing B0 given in Figure 7.

b1 Spheres Inversive coordinates b1 0 0 1 1 1 b2 0 0 −1 1 1 b−4 b−3 b3 1 1 0 0 1 b−2 b4 −1 1 0 0 1 b−1 b−1 0 0 −1 −1 1 b3 b4 b−2 0 0 1 −1 1 b−3 −1 −1 0 0 1 b−4 1 −1 0 0 1 b2

Figure 7. The standard orthoplicial sphere packing B0.

∗ It can be checked that B0 can be split into two orthoplicial ball packings, both obtained from B0 by a π rotation of 2 degree, one around the x-axis and the other around the y-axis, as shown in Figure 8. 

Figure 8. An orthoplicial trinity containing the standard orthoplicial sphere packing B0.

The three packings in the previous orthoplicial trinity are 1-CBP projections of the orthoplex. A similar case arises when an orthoplicial trinity contains a 0-CBP projection of the orthoplex. In this case the other two packings must be 3-CBP projections as it is shown in Figure 9. We notice that any orthoplicial trinity is M¨obiusequivalent to the ball-arrangement projection of a circunscribed 24-cell. 12 Figure 9. From left to right: a 3-CBP projection BO4 of the orthoplex; a 0-CBP projection ∗ ∗ contained in BO4 ; the other 3-CBP projection contained in BO4 ; the orthoplicial trinity containing BO4 .

3.2. Apollonian groups of orthoplicial sphere packings. In [RR21a], several groups generalizing the notion of the Apollonian, SuperApollonian and symmetric groups of [Gra+05; Gra+06a; Gra+06b] were introduced for polytopal d-ball packings. We recall that for any polytopal d-ball packing BP , the ∗ Apollonian group A(BP ) is the group generated by the inversions on the dual spheres of BP , the sym- metric group of BP is the group generated by the inversion leaving invariant BP , and the symmetrized Apollonian group SA(BP ) is the group generated by the previous two groups. When P is M¨obiusunique, then the three groups can be defined canonically for P.

In this context, the orthoplicial Apollonian group A(BO4 ), introduced by Nakamura in [Nak14] and also by Dias in [Dia14], corresponds to the Apollonian group of the standard orthoplicial sphere packing B0. Lemma 6. The symmetrized orthoplicial Apollonian group is the hyperbolic Coxeter group with Coxeter graph and it is isomorphic to the group generated to the following five matrices: 4 4 −1 0 0 0 0  1 1 −1 −1 0 1 0 0 0 0 1 0 0 0 0 1 0 0 0 0 0 1 0 0 0 1 1 1 1 1 0 0 1 0 0 0 0 1 0 0 0 0 −1 0 −2 2 V =  0 0 1 0 0 E = −1 1 1 −1 0 R = 0 0 −1 0 0 F = 0 0 0 1 0 S = 0 0 1 0 0 .  0 0 0 1 0 2 −1 1 −1 1 0 0 0 0 1 0 0 0 1 0 0 0 −2 0 −1 2 0 0 0 0 1 0 0 0 0 2 0 0 0 0 1 0 0 0 0 1 0 −2 0 −2 3

Proof. Let B0 be the standard orthoplicial ball packing. Let {rv, re, rr, rf } be the fundamental generators 4 of the symmetric group of B0, with respect to the flag (v, e, r, f, O ), where v = 1, e = 12, r = 123 and f = 1234. Then, {rv, re, rr, rf } are represented by the following M¨obiustransformations:

• rv is the reflection on the plane {x = 0}. • re is the inversion on the sphere with center (−1, 1, −1) and radius 2. • rr is the reflection on the plane {z = 0}. √ • rf is the inversion on the sphere with center (0, 0, 1) and radius 2.

Since the orthoplex is regular, then SA(B0) is generated by {rv, re, rr, rf , sf }, where sf is the inversion on the sphere orthogonal to b1, b2, b3 and b4. By using the inversive coordinates, we obtain a linear faithful ↑ representation of SA(B0) as the discrete subgroup of O4,1(Q), where the generators {rv, re, rr, rf , sf } are represented by the matrices {V, E, R, F, S}, respectively. The relations in the Coxeter graph can be checked by straightforward computations on the matrices. By the M¨obiusunicity of the orthoplex, the symmetrized Apollonian group of any other orthoplicial ball packing is isomorphic to SA(B0). 

13 ↑ Corollary 3.1. The orthoplicial Apollonian group is isomorphic to the subgroup of O4,1(Z) generated by the following 16 matrices:

1 0 0 0 0 1 0 0 0 0 0 −1 0 −2 2 0 −1 −2 0 2 S1234 =0 0 1 0 0 = S, S = 0 −2 −1 0 2 = FS1234F, 0 −2 0 −1 2 1234 0 0 0 1 0 0 −2 0 −2 3 0 −2 −2 0 3 1 0 0 0 0  1 0 0 0 0  0 −1 2 0 2 0 −1 0 2 2 S =0 2 −1 0 −2 = RS R, S = 0 0 1 0 0  = FS F, 1234 0 0 0 1 0  1234 1234 0 2 0 −1 −2 1234 0 −2 2 0 3 0 −2 0 2 3 −1 0 0 −2 2 −1 0 −2 0 2 0 1 0 0 0 0 1 0 0 0 S = 0 0 1 0 0 = ES E, S = −2 0 −1 0 2 = FS F, 1234 −2 0 0 −1 2 1234 1234  0 0 0 1 0 1234 −2 0 0 −2 3 −2 0 −2 0 3 −1 0 2 0 2  −1 0 0 2 2  0 1 0 0 0 0 1 0 0 0 S = 2 0 −1 0 −2 = RS R, S =  0 0 1 0 0  = FS F, 1234  0 0 0 1 0  1234 1234  2 0 0 −1 −2 1234 −2 0 2 0 3 −2 0 0 2 3 −1 0 0 2 −2 −1 0 2 0 −2 0 1 0 0 0 0 1 0 0 0 S = 0 0 1 0 0  = VS V, S =  2 0 −1 0 2  = FS F, 1234  2 0 0 −1 2  1234 1234  0 0 0 1 0  1234 2 0 0 −2 3 2 0 −2 0 3 −1 0 −2 0 −2 −1 0 0 −2 −2 0 1 0 0 0 0 1 0 0 0 S =−2 0 −1 0 −2 = RS R, S =  0 0 1 0 0  = FS F, 1234  0 0 0 1 0  1234 1234 −2 0 0 −1 −2 1234 2 0 2 0 3 2 0 0 2 3 1 0 0 0 0  1 0 0 0 0  0 −1 0 2 −2 0 −1 2 0 −2 S =0 0 1 0 0  = ES E, S = 0 2 −1 0 2  = FS F, 1234 0 2 0 −1 2  1234 1234 0 0 0 1 0  1234 0 2 0 −2 3 0 2 −2 0 3 1 0 0 0 0  1 0 0 0 0  0 −1 −2 0 −2 0 −1 0 −2 −2 S =0 −2 −1 0 −2 = RS R, S = 0 0 1 0 0  = FS F. 1234 0 0 0 1 0  1234 1234 0 −2 0 −1 −2 1234 0 2 2 0 3 0 2 0 2 3

↑ The previous representation of the orthoplicial Apollonian group in O4,1(Z) can be obtained by con- ↑ jugating the matrix representation OW (Z) given by Nakamura in [Nak14], where W is the matrix of the inversive product in augmented curvature-center coordinates coordinates (see [LMW01]). Both represen- tations satisfy the following relations: 2 (R1) Sijkl = 1 for every (i, j, k, l) = (±1, ±2, ±3, ±4). 2 0 0 0 0 (R2)( SijklSi0j0k0l0 ) = 1 if the labels ijkl and i j k l differ from only one letter.

4. Apollonian sections of the orthoplicial Apollonian sphere packing.

Let BP a polytopal d-ball packing and let Ω(BP ) be the Apollonian cluster of BP , defined as the orbit of A(BP ) on BP . We say that a subset Σ ⊂ Ω(BP ) is an Apollonian section of Ω(BP ) if there is subgroup 0 0 0 Γ < A(BP ) and a subset X ⊆ BP such that Σ = Γ · X. Let Σ = Γ · X and Σ = Γ · X be two Apollonian 0 sections of two polytopal Apollonian clusters Ω(BP ) and Ω(BP0 ), respectively. We say that Σ and Σ are algebraically equivalent if Γ and Γ0 are isomorphic and there is equivariant bijection between Γ · X and Γ0 · X0 with respect to the actions. More specifically, Σ and Σ0 are algebraically equivalent if there exists two bijections φ :Γ → Γ0 and ψ :Σ → Σ0 such that (i) φ :Γ → Γ0 is a group isomorphism. (ii) For all g ∈ Γ and all b ∈ X, ψ(g · b) = φ(g) · ψ(b). If in addition ψ preserves curvatures, we say that Σ and Σ0 are arithmetically equivalent.

14 In [RR21a], the Platonic Apollonian groups were defined as the Apollonian groups of the polytopal disk packings whose tangency polytope is one of the Platonic solids. We shall say that an Apollonian cluster has a tetrahedral, octahedral, cubical, icosahedral or dodecahedral Apollonian section if it contains an Apollonian section algebraically equivalent to an Apollonian packing of the corresponding Platonic solid. Theorem 4. Every orthoplicial Apollonian packing contains a tetrahedral, octahedral and cubical Apol- lonian section. Proof. Since the orthoplex is M¨obiusunique, it is enough to find the desired Apollonian sections in a particular orthoplicial Apollonian packing. First, we shall construct a tetrahedral and octahedral Apollonian section in Ω(B0) (see Figure 10).

Figure 10. The orthoplicial Apollonian packing Ω(B0) at depth ≤ 3, seen from above.

(Tetrahedral section) Let ΣT = ΓT ·XT where XT = {b1, b2, b3, b4} ⊂ B0,ΓT = hs1234, s1234, s1234, s1234i < ∗ A(B0) and let b1234 ∈ B0 . One can check that ΓT leaves invariant b1234. Let HT be the boundary of 2 b1234. In B0, HT is the plane {y = 1}. After identifying HT with Rc, we obtain that

XT ∩ HT := {b1 ∩ HT , b2 ∩ HT , b3 ∩ HT , b4 ∩ HT } becomes a tetrahedral disk packing BT 3 . Moreover, the restriction of ΓT on HT induces an isomorphism between ΓT and A(BT 3 ) given by

φT :ΓT −→ A(BT 3 )

sb 7−→ sb∩HT where sb denotes the inversion on b. On the other hand, we can construct a bijection

ψT :ΣT −→ Ω(BT 3 )

g · b 7−→ (g · b) ∩ HT By inspection, one can check that for every ijkl ∈ {1234, 1234, 1234, 1234} and for every m = 1, 2, 3, 4, we have

ψT (sijkl · bm) = (sijkl · bm) ∩ HT

= φT (sijkl) · (bm ∩ HT )

= φT (sijkl) · ψT (bm) which gives the equivariance of ψT . Therefore, ΣT is algebraically equivalent to the tetrahedral Apollo- nian packing Ω(BT 3 ).

(Octahedral section) We may apply a similar strategy to construct an octahedral section in Ω(B0) by intersecting with the plane HO := {x − y = 1}, orthogonal to every b ∈ XO := B0 \{b4, b−4}. The intersection HO ∩ XO gives an octahedral disk packing BO3 . We define then ΣO := ΓO · XO ⊂ Ω(B0) where

ΓO = ht123, t123, t123, t123, t123, t123, t123 | tijk := sijk4sijk4i. 15 In this case, the group isomorphism φO :ΓO 7→ A(BO3 ) is given by tijk 7→ sijk, where sijk denotes the inversion on the disk orthogonal to the three disks {bi ∩HO, bj ∩HO, bk ∩HO}. The equivariant bijection ψO :ΣO → Ω(BO3 ) is then given by g(b) 7→ g(b) ∩ HO. We illustrate in Figure 11 the tetrahedral and octahedral sections ΣT and ΣO.

Figure 11. A tetrahedral (left) and octahedral (right) section of the orthoplicial Apollo- nian packing Ω(B0).

(Cubical section) This case is a bit more tricky. We consider the orthoplicial ball packing B1 obtained by a 3-CBP projection of the orthoplex, with the labelling and coordinates given in Figure 12.

Spheres √ Inversive√ coordinates√ √ b1 1/√2 1/√2 1/√2 1/√2 1 b−4 b1 b2 1/√2 −1/√2 1/√2 −1/√2 1 b2 b−3 b3 1/√2 −1/√2 −1/√2 1/√2 1 b4 1/√2 1/√2 −1/√2 −1/√2 1 b −1 b4 b−1 −1/√2 −1/√2 −1/√2 −1/√2 1 b b3 −2 b−2 −1/√2 1/√2 −1/√2 1/√2 1 b−3 −1/√2 1/√2 1/√2 −1/√2 1 b−4 −1/ 2 −1/ 2 1/ 2 1/ 2 1

Figure 12. The orthoplicial ball packing B1.

In this labelling, i is positive if and only if the first coordinate of the center of bi is positive. The orthoplicial Apollonian packing Ω(B1) and a cubical Apollonian section of Ω(B1) are shown in Figure 14. Let us now describe the construction. First, we consider the rigde-scribed hypercube 1 C4 = conv{√ (±1, ±1, ±1, ±1)}. 2 We then split C4 into two edge-scribed orthoplexes induced the classes of a 2-coloring of its vertices. Let O4 be the orthoplex in the class of the vertex √1 (1, 1, 1, 1). We have that B is the ball-arrangement 2 1 4 4 3 4 projection of O . Let π : E → E the orthographic projection onto the hyperplane {x1 = 0} ⊂ E . We 4 3 1 have that π(O ) is the edge-scribed cube C with vertices √ (±1, ±1, ±1). Let B 3 be the cubical disk 2 C 3 4 packing obtained by the ball-arrangement projection of C . By mapping bv → bπ(v), for every v ∈ O , we construct a bijection πe : B1 → BC3 .

∗ 3 Let Bx be the packing made by the six spheres of B1 which are orthogonal to HC := {x = 0} ⊂ Rc (see the blue packing in Figure 9), and let ΓC < A(B1) be the group generated by the inversions on the spheres of Bx. According to the labelling of B1, we have that ΓC corresponds to the parabolic subgroup of A(B1) hs1234, s1234, s1234, s1234, s1234, s1234i. 16 2 ∗ By intersecting with HC, and then identifying HC to Rc, we map Bx to BC3 (see Figure 13).

∗ Figure 13. (Left) B1 in gray with Bx in blue; (right) BC3 in gray with BC3 in blue.

Therefore, we can define a group isomorphism φC :ΓC → A(BC3 ) by mapping the inversion on every ∗ sphere b ∈ Bx to the inversion on the circle (b ∩ HC) ∈ BC3 . Let ΣC := ΓC ·B1. It can be checked that for every g ∈ ΓC and every b ∈ B1,

ψC :ΣC → Ω(BC3 ) g · b 7→ φ(g) · πe(b) defines an equivariant bijection with respect to the action of ΓC on B1. Therefore, ΣC is a cubical Apollonian section of Ω(B1). 

We notice that the bijections ψT and ψO described above preserve both the inversive product which is not the case for ψC. Indeed, spheres which are tangent in the cubical Apollonian section may correspond to disjoint disks in the cubical Apollonian packing. To see this, we might consider the cubical Apollonian 0 packing Ω(BC3 ) given above with a 2-coloring. Then, two disks ψC(b) and ψC(b ) have the same color 0 if and only if the centers of b and b both lie on the same side of the plane HC. When two disks with same color (and therefore disjoint), correspond to two vertices lying in a same square-face, then the corresponding spheres are tangent (see Figure 14).

Figure 14. From left to right: the orthoplicial Apollonian sphere packing Ω(B1); a cubical Apollonian section ΣC of Ω(B1); ΣC zoomed; the corresponding cubical Apollonian disk packing with a minimal coloration.

On the arithmetic side, the bijections ψT and ψ√O preserve the curvatures, contrarily to ψC. However, 2 by composing ψC with a rescaling of Rb of factor 2 we obtain a bijection that does preserve the curva- tures. These properties will be important in the next section in order to construct integral orthoplicial Apollonian packings containing an Apollonian section with prescribed curvatures.

17 4.1. Construction of orthoplicial Apollonian packings containing a given integral section. In this section we shall prove Theorem 1. For the integrality part, we will need the following three propositions.

Proposition 4.1. Let κ1, κ2, κ3, κ4 be the curvatures of four pairwise tangent spheres of an orthoplicial p sphere packing BO4 . Then, Ω(BO4 ) is integral if and only if κ1, κ2, κ3, κ4 and T3(κ1, κ2, κ3, κ4) are integers. Proof. Nakamura gave a proof in [Nak14] for all orthoplicial sphere packings M¨obiusequivalent to the standard B0. Since the orthoplex is M¨obiusunique, Nakamura’s proposition applies to every orthoplicial sphere packing. 

Proposition 4.2. Let κ1, κ2, κ3 be the curvatures of three pairwise tangent disks of on octahedral disk p packing BO3 . Then, Ω(BO3 ) is integral if and only if κ1, κ2, κ3 and 2T2(κ1, κ2, κ3) are integers. Proof. The sufficiency part was treated in [RR21a]. For the necessity, we shall mimic the same method used by Nakamura to show the validity of Proposition 4.1. Let us suppose that Ω(BO3 ) is integral. Then

(25) T2(κ1, κ2, κ3) = κ1κ2 + κ2κ3 + κ1κ3 ∈ Z. By the antipodal relation we have that

2κO3 = κv + κv¯ ∈ Z. On the other hand, by Corollary 2.1, we have p κO3 = κ1 + κ2 + κ3 ± 2T2(κ1, κ2, κ3) Therefore, p 2κO3 = 2(κ1 + κ2 + κ3) ± 2 2T2(κ1, κ2, κ3) ∈ Z p which implies that 2 2T2(κ1, κ2, κ3) is an integer. Let us show that is an even integer. If there is m ∈ Z p such that 2 2T2(κ1, κ2, κ3) = 2m + 1, then we would have 1 1 p2T (κ , κ , κ ) = m + ⇔ 2T (κ , κ , κ ) = m2 + m + 6∈ , 2 1 2 3 2 2 1 2 3 4 Z p p contradicting (25). Hence, 2 2T2(κ1, κ2, κ3) is an even integer, so 2T2(κ1, κ2, κ3) is also an integer. 

Proposition 4.3. Let κ1, κ2, κ3 be the curvatures of three consecutive tangent disks in a cubical disk p packing BC3 . Then Ω(BC3 ) is integral if and only if κ1, κ2, κ3 and 2C2(κ1, κ2, κ3) are integers.

Proof. Idem to Proposition 4.2.  We may now prove Theorem 1. Let us recall the statement.

Theorem 1. Let Ω(BP ) be either a tetrahedral, an octahedral or a cubical Apollonian packing. There is an orthoplicial Apollonian packing Ω(BO4 ) containing an Apollonian section arithmetically equivalent to Ω(BP ). Moreover, Ω(BP ) is integral if and only if Ω(BO4 ) is integral.

Proof. We recall that the curvature of a d-ball b can be obtained from its inversive coordinates xb by T κ(b) = kd+2xb where kd+2 is the (d + 2)-column matrix (0,..., 0, −1, 1), c.f. Equation (5).

3 Let BO4 be the orthoplicial ball packing obtained from B0 by applying a translation of R of vector (−1, 0, 0), so HT becomes the plane given by {x = 0}. Let ΣT = ΓT ·XT be the corresponding tetrahedral section of Ω(BO4 ) and let BT 3 be the tetrahedral disk packing given by ψT (XT ). In inversive coordinates, −1 the bijection ψT : Ω(BT 3 ) → ΣT ⊂ Ω(BO4 ) is given by −1 ψT : Ω(BT 3 ) → ΣT xb 7→ (0, xb)

18 0 Let BT 3 be any tetrahedral disk packing. Since the tetrahedron is M¨obiusunique, there exists µ ∈ 2 0 0 ↑ M¨ob(Rc) such that µ ·BT 3 = BT 3 . Therefore, µ · Ω(BT 3 ) = Ω(BT 3 ). Let M ∈ O3,1(R) be the matrix corresponding to µ and letµ ˜ ∈ M¨ob(Rc3) be the M¨obiustransformation corresponding to the matrix

 T  1 04 ↑ ∈ O4,1(R) 04 M

0 where 04 is the null column-matrix of size 4. We have that ΣT := µe · ΣT is a tetrahedral section of the 0 orthoplicial Apollonian packing Ω(BO4 ) = µe · Ω(BO4 ). We define the bijection

−1 ψeT := µ ◦ ψT ◦ µe

0 0 0 mapping ΣT to Ω(BT 3 ). Let us show that ψeT preserves curvatures. For every disk b ∈ Ω(BT 3 ), we have

−1 −1 −1 κ(ψeT (b)) = κ(µe ◦ ψT ◦ µ (b))  T   T  T 1 04 04 −1 = k5 M xb where I4 is the identity matrix of size 4 04 M I4   T 0 = k5 xb T = k4 xb = κ(b)

0 0 0 and therefore, Ω(BT 3 ) and ΣT ⊂ Ω(BO4 ) are arithmetically equivalent.

We may use the beginning of the proof to show the octahedral case in a similar way as above. In this case, we may consider the orthoplicial sphere packing obtained from B0 after a rotation around the π axis directed by (0, 0, 1) of angle 4 , so the plane HO becomes the plane {x = 0}. Then, the bijection −1 −1 ψO : Ω(BO3 ) → ΣO has the same expression in inversive coordinates as ψT and the rest of computa- tions are identical as above.

The cubical case is more delicate since the corresponding Apollonian section is not obtained by a cross section. Let B1 be the orthoplicial sphere packing of Figure 12 and let ΣC = ΓC ·B1 be the cubical section described in the proof of Theorem 4. Let BC3 be the cubical disk packing given by ψC(B1). We define the function  : BC3 → {1, −1} by (b) = sign(x1) where x1 is the first coordinate of the vertex −1 −1 corresponding to ψC (b). Then, in inversive coordinates, the bijection ψC : BC3 → B1 and the group −1 isomorphism φC : A(BC3 ) → ΓC are given by

−1 −1 0 φ : A(BC3 ) → ΓC ψ : B 3 −→ B C (26) C C 1  1 0T  x 7−→ √1 ((b), x ) 4 b 2 b A 7→ 04 A

−1 The equivariance of ψC allow us to extend ψC : Ω(BC3 ) → ΣC by

−1 −1 −1 ψC (g · b) = φC (g) · ψC (b)(27)

0 for every g ∈ A(BC3 ) and every b ∈ BC3 . Let BC3 be any cubical disk packing. As before, the M¨obius 2 0 unicity of the cube implies that there is µ ∈ M¨ob(Rc) such that µ · Ω(BC3 ) = Ω(BC3 ). We define 3 3 µ ∈ M¨ob(Rc) analogously to the tetrahedral case described above. Let µ = λ √1 ◦ µ ∈ M¨ob(Rc) where e b 2 e 3 √1 0 λ √1 is the rescaling of Rc with factor . We have that ΣC := µ · ΣC is a cubical section of the 2 2 b 0 −1 orthoplicial Apollonian packing Ω(BO4 ) = µb · Ω(B1). We define the bijection ψbC = µ ◦ ψC ◦ µb mapping 0 0 ΣC to Ω(BC3 ). As above, we need to prove that ψbC preserves curvatures. For, we first notice that for 0 0 0 0 0 0 0 any disk b ∈ Ω(BC3 ), there is g ∈ A(BC3 ) and b ∈ BC3 such that b = g · b , by definition. Thus, for any 0 b ∈ Ω(BC3 ), we have 19  −1   −1 0 0  κ ψbC (b) = κ ψbC (g · b ) −1 −1 0 0
= κ µb ◦ ψC ◦ µ (g · b )  −1 −1 0 0  = κ λ √1 ◦ µe ◦ ψC ◦ µ (g · b ) √ 2 −1 −1 −1 0
−1 0 = 2κ µ ◦ ψC ◦ µ (µgµ ) · (b ) where g = µ g µ ∈ A(BC3 ) √ e −1 −1 −1 0
= 2κ µe ◦ φC (g) · ψC ◦ µ (b ) by (27) √  T   T   −1 0  T 1 04 1 04 1 (µ (b )) = 2k5 √ where A is the matrix of g 04 M 04 A 2 xµ−1(b0)  −1 0  T (µ (b )) = k5 MAxµ−1(b0) T = k4 MAxµ−1(b0) T −1 = k4 MAM xb0 = κ(µgµ−1 · b0) = κ(g0 · b0) = κ(b)

0 0 0 Let us now suppose that Ω(BT 3 ), Ω(BO3 ) and Ω(BC3 ) are integral. We shall show that in the three 0 cases we can find four pairwise tangent spheres of BO4 with curvatures κ1, κ2, κ3 and κ4 satisfying that p 0 T3(κ1, κ2, κ3, κ4) ∈ Z. The integrality of Ω(BO4 ) then follows from Proposition 4.1.

0 (Tetrahedral section) Let κ1, κ2, κ3, κ4 be the curvatures of the four disks of BT 3 . By Descartes’ theorem, we have that

T3(κ1, κ2, κ3, κ4) = 0.

−1 Since ψeT preserves curvatures, then κ1, κ2, κ3, κ4 are the curvatures of four pairwise tangent spheres of 0 p BO4 satisfying that T3(κ1, κ2, κ3, κ4) ∈ Z.

0 (Octahedral section) Let κ1, κ2, κ3, κ−1, κ−2, κ−3 be the curvatures of the six disks of BO3 under an antipodal labelling, so κi and κ−i are the curvatures of non-tangent disks. By (19) , we have

2 2 2 2 (κ1 − κO3 ) + (κ2 − κO3 ) + (κ3 − κO3 ) = 2κO3

1 −1 where κO3 = 2 (κi+κ−i) for every i = 1, 2, 3. Since ψeO preserves curvatures, then κ1, κ2, κ3, κ−1, κ−2, κ−3 0 are also the curvatures of six spheres of BO4 . Let κ4 and κ−4 be the curvatures of the remaining spheres 0 0 of BO4 . One can check that the labelling of the curvatures in BO4 is also antipodal. By the antipodal 0 0 relation, which holds for BO3 and BO4 , we have

κ1 + κ−1 (28) κ 3 = = κ 4 . O 2 O

0 We also have that κ1, κ2, κ3, κ4 are the curvatures of four pairwise tangent spheres of BO4 . By combining (28) with Corollary 2.3 we obtain

2 2 2 2 2 (κ1 − κO4 ) + (κ2 − κO4 ) + (κ3 − κO4 ) + (κ4 − κO4 ) = 2κO4 2 2 2 2 2 ⇔ (κ1 − κO3 ) + (κ2 − κO3 ) + (κ3 − κO3 ) + (κ4 − κO3 ) = 2κO3

⇔ κ4 = κO3

20 Therefore,

p p T3(κ1, κ2, κ3, κ4) = T3(κ1, κ2, κ3, κO3 ) q 1 2 2 2 2 2
= (κ + κ + κ + κ 3 ) − 2 κ + κ + κ + κ 2 1 2 3 O 1 2 3 O3 q 1 2 2 2 2 = 2(κ κ + κ κ + κ κ ) + 2κ − (κ − κ 3 ) − (κ − κ 3 ) − (κ − κ 3 ) 2 1 2 2 3 1 3 O3 1 O 2 O 3 O 1 = p2(κ κ + κ κ + κ κ ) 2 1 2 2 3 1 3 1 = p2T (κ , κ , κ ) 2 2 1 2 3

0 p Since Ω(B 3 ) is integral, then, 2T2(κ1, κ2, κ3) is an even integer. Moreover, by Proposition 4.2, 2T2(κ1, κ2, κ3) O p is also an even integer implying thus T3(κ1, κ2, κ3, κ4) is an integer.

(Cubical section) Let κ1, κ2, κ3, κ4 be the curvatures of four consecutive tangent disks b1, b2, b3, b4 of 0 BC3 with b4 tangent to b1. Then κ1, κ2, κ3, κ4, by applying the antipodal relation on a square-face of κC3 we have

κ1 + κ3 = κ2 + κ4.

−1 Since ψbC preserves curvatures, then κ1, κ2, κ3, κ4 are the curvatures of four pairwise tangent spheres of BO4 . Then, we have

p p T3(κ1, κ2, κ3, κ4) = T3(κ1, κ2, κ3, κ1 + κ3 − κ2) q 2 = κ1κ2 + κ2κ3 + κ1κ3 − κ2 p = 2C3(κ1, κ2, κ3)

0 p By the integrality of Ω(BC3 ) and Proposition 4.3, we have that T3(κ1, κ2, κ3, κ4) is an integer. 

The proof of Theorem 1 is constructive giving thus a method to obtain integral orthoplicial Apollonian packings containing a given tetrahedral, octahedral or cubical section. We have applied this method to produce the three orthoplicial Apollonian packings showed in Figures 2, 15 and 16.

293 365 365 293 398 229 229 398 294 173 173 294 206 125 300 300 125 206 134 85 212 212 85 134

356 356 381 140 140 381

309 309 260 53 53 260 340 78 78 340 84 398 398 84 180 180 189 382357 357382 189 220 29 302 302 29 220

358 358 141 364 364 141 373 222 222 373 116 38 262 262 38 116 341 156 237 237 156 341 356 356 238 325 325 238 294 294 300 44 44 300

198189 189198 158 158

340 340 308 308 166 166

221236 236221

389 68 68 389 390 390 268 268 13 374365 365374 13 69 84 348 348 84 69 228 349 278 110 110 278 349 228 244 244 173 126 126 173 380 117 117 380 325 325 388 388 150156 317 317 156150 157 157

277 277 366350 350366 269 269 397 397 388 20 20 388 93 93 133 133 133 133

180 373 373 180 164 124 284 94 94 284 124 164 14 118 334 358 4 358 334 118 14 78 245 245 78 292396 396292 246172 172246 342 342 349 302 302 349 380 380 45 388 388 45 229182 182229 230 220 220 230 268 268 221 52 52 221

317286 286317 270 270 270 270 148 148 36 205 205 36 333 333 206 397 397 206 398 398 382 254 254 382 62 197 389 389 197 62 101108 69 262 262 69 108101 70 70 70 70

181188 262269 269262 188181 204 204 116 116

285292 292285 174 174

244 244

326 326

214 214 278285 12 285278 76 76

292 292 77 221 221 77 70 286 286 70 236 236 237 237 172 172 318 318 55 356 356 270 270 196 196 332 332 325 325 134 382 382 134 220 220 220 220 316 316 213 317 317 213 84 132 102 102 132 84 374 124 124 374

125 342364 125 125 364342 125 286238 388 158166 166158 388 238286 301 381 381 301 46 68 181 181 68 46

293356 356293 308 61 28 28 61 308 173 173 252 229 229 252 230 230 188 188 230 230 189 222 166 166 222 189 389 20 20 389

269 380 132 132 380 269 358 358 109 262 262 109 110 332317 317332 110 388 388 365 365 276 21 6 21 276 76 76 390 390

276 276 206 206 253 68 68 253 156 156 349 53 53 349 334 334 148 101 268 268 101 148 260 165 165 260 245 245 -2 341 341

Figure 15. From left to right: an integral octahedral Apollonian packing with initial curvatures (−2, 4, 5); an integral orthoplicial Apollonian packing with initial curvatures (−2, 4, 5, 5); an octahedral Apollonian section of the packing on the center, which is arith- metically equivalent to packing on the left.

21 462 462 404 404 350 470 470 350 412 412 300 358 358 300 254 308 308 254 212 469 469 212 477262 262477 174389 220 220 389174 397 397 317 182 182 317 489 140 325 325 140 489 253 148 148 253 261 261 377490110 110490377 464 464 197 118 118 197 84 205 205 84 281368 433 433 368281 346433 433346 149 92 92 149 157 157 424 424 489 489 476 62 402337 337402 62 476 201266 266201 422 485 485 422 313 493 493 313 469 248 70 70 248 469 109 109 117 322 322 117 478 478 486 369 369 486 413 257 257 413 421 304 304 421

326373 44 44 373326 466 466

184 349 349 184 137 357 52 52 357 137 217 382390 390382 217 424 424 284 462 462 284 457 170 240 240 170 457 317 77 77 317 485 85 273 193 193 273 85 485 340348 226 226 348340 420 420 453 453 293301 301293

389422 30 38 38 30 422389 381 381 245 245253 348 348 253245 245 212 488 488 212 337 276 276 337 304 268 178 178 268 304 413 145 145 413 122 436 201 201 436 122 458 5 458 145 168 481 318 318 481 168 145 89 61 341 341 61 89 498 53 53 498 406 406 274 465 465 274 112 373 238246 246238 373 112 182 297 297 182 205 205213 213205 205 494 494 309 309 332 453461 461453 332

470 470 493 493 493 493

293270 270293 28 401378 378401 28 136113 113136 413 413 413 413 405 360252 252360 405 430 430 377 445 445 377 181 488 153 269 269 153 488 181 173 130 462 462 130 173 198 497 497 198 412 412 432 190 389 389 190 432 20 449 449 20

413 413 413 413 390 45 430 430 45 390 268 228 228 268 245 405 245 320 320 245 405 245 382 157 277 457 337 222 397 397 222 337 457 277 157 382 494 434 314 314 434 494 337 294 180 97 412 412 97 180 294 337 354 237 237 354 134 217 172 157 114 329 329 114 157 172 217 134 194 80 80 194 493 37 37 493 401 476 112 112 476 401 356 149 433 433 149 356 416 448 97 378 468 341 341 448 448 341 341 468 378 97 448 416 465 333 129 129 333 465 425 276 356 425 425 356 276 425 442 336 22 442 442 22 336 442 333 333 313 350 350 313 57 74 253 253 74 57 386 386 401 401 386 317 317 386 176 436 436 176 326 193 309 141 141 309 193 326 445 445 428 156 37 37 156 428 445 445 133 133 385368 368385 133 133 116 422 20 422 116 148 414 498481 481498 414 148 333 449 449 333 316 464 464 316 330 412 150 150 412 330 214 345 133 133 345 214 197 404 404 197 406293 293406 317 142 398285 285398 142 317 300 125 125 300 442 442 442 442 425 425 425 425 329 329 312 312

438453 29 29 453438 14 205220 220205 14 264 264 249 437 437 249 214 420 420 214 258 197 197 258 365 73 241 241 73 365 394 409 409 394 277 460 460 277 377 58 345 345 58 377 445 182 330 330 182 445 260428 277109153 328 328 153109277 428260 490 313 313 490 165 260 473 473 260 165 494 92 136 56 56 136 92 494 385 385 385 385 341 472 341 341 472 341 457 474 474 457 246290 145 430 430 145 290246 476 41 130 101 101 130 41 476 280369 12 369280 325 253 253 325 462 236 164 164 236 462 86 238 238 86 397245 245397 149 382 382 149 452 344 344 452 230 461 461 230 329 329 -3 446 446

Figure 16. From left to right: an integral cubical Apollonian packing with initial cur- vatures (−3, 5, 12); an integral orthoplicial Apollonian packing with initial curvatures (−3, 5, 12, 20); a cubical Apollonian section of the packing on th center which is arith- metically equivalent to the packing on the left.

5. Open questions and concluding remarks In [RR21a] it was also conjectured that for d ≥ 2, all edge-scribable (d + 1)-polytopes are M¨obius unique. In view of the family of polytopes that are known to be M¨obiusunique, a natural next step would be to show the validity of the conjecture for regular polytopes. Question 1. Are the 600-cell and the 120-cell M¨obiusunique? We believe that the techniques developed in this work can be useful to study open question on the behaviour of the curvatures in integral Apollonian packings such as the local-global principle. Similar techniques can be applied to construct integral Apollonian packings of the hypercube and 24-cell con- taining Platonic Apollonian sections.

As instance, we can consider the polytopal sphere packing BR4 obtained by rescaling a 1-CBP projec- tion of the 24-cell by a factor of √1 (see Figure 17). Numerical experiments suggest that every integer 3 is represented without no modulo restriction as a curvature of Ω(BR4 ). We end with the following two conjectures.

Conjecture 1. The set of curvatures of Ω(BR4 ) is N.

Conjecture 2. There is a sequence of consecutive tangent spheres (b0, b1,...) ⊂ Ω(BR4 ) such that, for every n ∈ N, the curvature of bn is n.

Figure 17. The sphere packing BR4 (left) and its Apollonian packing Ω(BR4 ) (right) with the curvatures.

22 Figures: All the figures in this paper were done by an own software coded with Mathematica 12.2 [Inc].

Acknowledgements: This work is a part of the content of the PhD thesis of the author at the University of Montpellier. The author would like to thank his advisor, Prof. Jorge Ram´ırezAlfons´ınfor several helpful suggestions and many fruitful discussions.

References [A19] Sheydvasser A. “Quaternion orders and sphere packings”. In: Journal of Number Theory 204 (2019), pp. 41–98. issn: 0022-314X. doi: https://doi.org/10.1016/j.jnt.2019.03.014. url: https://www.sciencedirect.com/science/article/pii/S0022314X19301192. [Bar18] Arthur Baragar. “Higher dimensional Apollonian packings, revisited”. In: Geometriae Ded- icata 195.1 (Aug. 2018), pp. 137–161. issn: 1572-9168. doi: 10.1007/s10711-017-0280-7. url: https://doi.org/10.1007/s10711-017-0280-7. [BS93] G. R. Brightwell and E. R. Scheinerman. “Representations of Planar Graphs”. In: SIAM Journal on Discrete Mathematics 6.2 (1993), pp. 214–229. doi: 10.1137/0406017. eprint: https://doi.org/10.1137/0406017. url: https://doi.org/10.1137/0406017. [Che16] H. Chen. “Apollonian Ball Packings and Stacked Polytopes”. In: Discrete & Computational Geometry 55.4 (2016), pp. 801–826. issn: 1432-0444. doi: 10.1007/s00454-016-9777-3. url: https://doi.org/10.1007/s00454-016-9777-3. [CL15] H. Chen and J.P. Labb´e.“Lorentzian Coxeter systems and Boyd–Maxwell ball packings”. In: Geometriae Dedicata 174.1 (Feb. 2015), pp. 43–73. issn: 1572-9168. doi: 10.1007/s10711- 014-0004-1. url: https://doi.org/10.1007/s10711-014-0004-1. [CP17] H. Chen and A. Padrol. “Scribability problems for polytopes”. In: European Journal of Combinatorics 64 (2017), pp. 1–26. issn: 0195-6698. doi: https://doi.org/10.1016/ j.ejc.2017.02.006. url: https://www.sciencedirect.com/science/article/pii/ S0195669817300136. [Cox73] H.S.M. Coxeter. Regular Polytopes. Dover books on advanced mathematics. Dover Publica- tions, 1973. isbn: 9780486614809. [Dia14] D. Dias. “The Local-Global Principle for Integral Generalized Apollonian Sphere Packings”. In: arXiv: Number Theory (2014). [EKZ02] D. Eppstein, G. Kuperberg, and G.M. Ziegler. “Fat 4-polytopes and fatter 3-spheres”. In: arXiv: Combinatorics (2002). [FSZ19] Elena Fuchs, Katherine E. Stange, and Xin Zhang. “Local-global principles in circle pack- ings”. In: Compositio Mathematica 155.6 (2019), pp. 1118–1170. doi: 10.1112/S0010437X19007139. [Gaj+21] Parmesh Gajjar, Chris Johnson, James Carr, Kevin Chrispeels, J. Gray, and Philip Withers. “Size segregation of irregular granular materials captured by time-resolved 3D imaging”. In: Scientific Reports 11 (Apr. 2021), p. 8352. doi: 10.1038/s41598-021-87280-1. [Gra+03] R. Graham, J. C. Lagarias, C. L. Mallows, A. R. Wilks, and C. Yan. “Apollonian circle packings: number theory”. In: Journal of Number Theory 100.1 (2003), pp. 1–45. issn: 0022-314X. doi: https://doi.org/10.1016/S0022- 314X(03)00015- 5. url: https: //www.sciencedirect.com/science/article/pii/S0022314X03000155. [Gra+06a] R. Graham, J. C. Lagarias, C. L. Mallows, A. R. Wilks, and C. Yan. “Apollonian Circle Packings: Geometry and Group Theory II. Super-Apollonian Group and Integral Packings”. In: Discrete & Computational Geometry 35 (Sept. 2006), pp. 1–36. doi: 10.1007/s00454- 005-1195-x. [Gra+06b] R. Graham, J. C. Lagarias, C. L. Mallows, A. R. Wilks, and C. Yan. “Apollonian Circle Packings: Geometry and Group Theory III. Higher Dimensions”. In: Discrete & Computa- tional Geometry 35 (2006). doi: 10.1007/s00454-005-1197-8. [Gra+05] R. Graham, J. C. Lagarias, C. L. Mallows, A. R. Wilks, and C. H. Yan. “Apollonian Circle Packings: Geometry and Group Theory I. The Apollonian Group”. In: Discrete & Com- putational Geometry 34.4 (2005), pp. 547–585. doi: 10.1007/s00454-005-1196-9. url: https://doi.org/10.1007/s00454-005-1196-9. [GM08] G. Guettler and C. L. Mallows. “A generalization of Apollonian packing of circles”. In: Journal of Combinatorics 1 (Jan. 2008). doi: 10.4310/JOC.2010.v1.n1.a1. [Inc] Wolfram Research Inc. Mathematica, Version 12.2. Champaign, IL, 2021. url: https : //www.wolfram.com/mathematica.

23 [KN19] A. Kontorovich and K. Nakamura. “Geometry and arithmetic of crystallographic sphere packings”. In: Proceedings of the National Academy of Sciences 116.2 (2019), pp. 436–441. issn: 0027-8424. doi: 10 . 1073 / pnas . 1721104116. eprint: https : / / www . pnas . org / content/116/2/436.full.pdf. url: https://www.pnas.org/content/116/2/436. [LMW01] J. Lagarias, C. L. Mallows, and A. Wilks. “Beyond the Descartes Circle Theorem”. In: The American Mathematical Monthly 109 (2001). doi: 10.2307/2695498. [Nak14] K. Nakamura. The local-global principle for integral bends in orthoplicial Apollonian sphere packings. 2014. arXiv: 1401.2980 [math.NT]. [RR21a] J. L. Ram´ırezAlfons´ınand I. Rasskin. “A polytopal generalization of Apollonian packings and Descartes’ theorem”. In: (2021). arXiv: 2107.09432 [math.CO]. [RR21b] J. L. Ram´ırezAlfons´ınand I. Rasskin. “Ball packings for links”. In: European Journal of Combinatorics 96 (2021), p. 103351. issn: 0195-6698. doi: https://doi.org/10.1016/ j.ejc.2021.103351. url: https://www.sciencedirect.com/science/article/pii/ S0195669821000433. [Sch87] E. Schulte. “Analogues of Steinitz’s theorem about non-inscribable polytopes”. In: Colloq. Math. Soc. J´anosBolyai 48 (Jan. 1987). [UN18] A. Umayahara and M. Nespolo. “Derivative structures based on the sphere packing:” in: Zeitschrift f¨urKristallographie - Crystalline Materials 233.3-4 (2018), pp. 179–203. doi: doi:10.1515/zkri-2017-2125. url: https://doi.org/10.1515/zkri-2017-2125. [Wil81] J. B. Wilker. “Inversive Geometry”. In: (1981). Ed. by Chandler Davis, Branko Gr¨unbaum, and F. A. Sherk, pp. 379–442.

IMAG, Univ. Montpellier, CNRS, Montpellier, France Email address: [email protected]

24