Upper bound of density for packing of congruent hyperballs in hyperbolic 3space^{1}^{1}1Mathematics Subject Classification 2010: 52C17, 52C22, 52B15.
Key words and phrases: Hyperbolic geometry, hyperball packings, packing density.
Abstract
In [26] we proved that to each saturated congruent hyperball packing exists a decomposition of 3dimensional hyperbolic space \mathbb{H}^{3} into truncated tetrahedra. Therefore, in order to get a density upper bound for hyperball packings, it is sufficient to determine the density upper bound of hyperball packings in truncated simplices. In this paper we prove, using the above results and the results of papers [16] and [23], that the density upper bound of the saturated congruent hyperball (hypersphere) packings related to the corresponding truncated tetrahedron cells is realized in a regular truncated tetrahedra with density \approx 0.86338. Furthermore, we prove that the density of locally optimal congruent hyperball arrangement in regular truncated tetrahedron is not a monotonically increasing function of the height (radius) of corresponding optimal hyperball, contrary to the ball (sphere) and horoball (horosphere) packings.
1 Preliminary results
Let X denote a space of constant curvature, either the ndimensional sphere \mathbb{S}^{n}, Euclidean space \mathbb{E}^{n}, or hyperbolic space \mathbb{H}^{n} with n\geq 2. An important question of discrete geometry is to find the highest possible packing density in X by congruent nonoverlapping balls of a given radius [1], [5].
Euclidean cases are the best explored. One major recent development has been the settling of the longstanding Kepler conjecture, part of Hilbert’s 18th problem, by Thomas Hales at the turn of the 21st century. Hales’ computerassisted proof was largely based on a program set forth by L. Fejes Tóth in the 1950’s [8].
In ndimensional hyperbolic geometry occur several new questions concerning the packing and covering problems, e.g. in \mathbb{H}^{n} there are 3 kinds of “generalized balls (spheres)”: the usual balls (spheres), horoballs (horospheres) and hyperballs (hyperspheres). Moreover, the definition of packing density is crucial in hyperbolic spaces as shown by Böröczky [3], for standard examples also see [5], [20]. The most widely accepted notion of packing density considers the local densities of balls with respect to their Dirichlet–Voronoi cells (cf. [3] and [11]). In order to consider ball packings in \overline{\mathbb{H}}^{n}, we use an extended notion of such local density.
In space X^{n} let d_{n}(r) be the density of n+1 mutually touching spheres or horospheres of radius r (in case of horosphere r=\infty) with respect to the simplex spanned by their centres. L. Fejes Tóth and H. S. M. Coxeter conjectured that the packing density of balls of radius r in X^{n} cannot exceed d_{n}(r). This conjecture has been proved by C. A. Rogers for Euclidean space \mathbb{E}^{n}. The 2dimensional spherical case was settled by L.Fejes Tóth in [7].
Ball (sphere) and horoball (horosphere) packings:
In [3] and [4] K. Böröczky proved the following theorem for ball and horoball packings for any n (2\leq n\in\mathbb{N}):
Theorem 1.1 (K. Böröczky)
In an ndimensional space of constant curvature consider a packing of spheres of radius r. In spherical space suppose that r<\frac{\pi}{4}. Then the density of each sphere in its DirichletVoronoi cell cannot exceed the density of n+1 spheres of radius r mutually touching one another with respect to the simplex spanned by their centers.
The above greatest density in \mathbb{H}^{3} is \approx 0.85328 which is not realized by packing with equal balls. However, it is attained by the horoball packing (in this case r=\infty) of \overline{\mathbb{H}}^{3} where the ideal centers of horoballs lie on the absolute figure of \overline{\mathbb{H}}^{3}. This ideal regular tetrahedron tiling is given with CoxeterSchläfli symbol \{3,3,6\}. Ball packings of hyperbolic nspace and of other Thurston geometries are extensively discussed in the literature see e.g. [1], [3], [6], [18] and [34], where the reader finds further references as well.
In a previous paper [12] we proved that the above known optimal horoball packing arrangement in \mathbb{H}^{3} is not unique using the notions of horoballs in same and different types. Two horoballs in a horoball packing are of the ”same type” iff the local densities of the horoballs to the corresponding cell (e.g. DV cell or ideal simplex) are equal, (see [29]). We gave several new examples of horoball packing arrangements based on totally asymptotic Coxeter tilings that yield the above Böröczky–Florian packing density upper bound (see [4])
We have also found that the BöröczkyFlorian type density upper bound for horoball packings of different types is no longer valid for fully asymptotic simplices in higher dimensions n>3 (see [28]). For example in \mathbb{H}^{4}, the density of such optimal, locally densest horoball packing is \approx 0.77038 larger than the analogous BöröczkyFlorian type density upper bound of \approx 0.73046. However, these horoball packing configurations are only locally optimal and cannot be extended to the whole hyperbolic space \mathbb{H}^{4}.
In papers [13] and [14] we have continued our previous investigation in \mathbb{H}^{n} (n\in\{4,5\}) allowing horoballs of different types. We gave several new examples of horoball packing configurations that yield high densities (\approx 0.71645 in \mathbb{H}^{4} and \approx 0.59421 in \mathbb{H}^{5}) where horoballs are centered at ideal vertices of certain Coxeter simplices, and are invariant under the actions of their respective Coxeter groups.
Hyperball (hypersphere) packings:
A hypersphere is the set of all points in \mathbb{H}^{n}, lying at a certain distance, called its height, from a hyperplane, on both sides of the hyperplane (cf. [38] for the planar case).
In hyperbolic plane \mathbb{H}^{2} the universal upper bound of the hypercycle packing density is \frac{3}{\pi}, and the universal lower bound of hypercycle covering density is \frac{\sqrt{12}}{\pi}, proved by I. Vermes in [37, 38, 39]. We note here that independently from him in [15] T. H. Marshall and G. J. Martin obtained similar results to the hypercycle packings.
In [30] and [31] we analysed the regular prism tilings (simply truncated Coxeter orthoscheme tilings) and the corresponding optimal hyperball packings in \mathbb{H}^{n} (n=3,4) and we extended the method developed of paper [31] to 5dimensional hyperbolic space (see [32]). In paper [33] we studied the ndimensional hyperbolic regular prism honeycombs and the corresponding coverings by congruent hyperballs and we determined their least dense covering densities. Furthermore, we formulated conjectures for candidates of the least dense hyperball covering by congruent hyperballs in 3 and 5dimensional hyperbolic spaces.
In [25] we discussed congruent and noncongruent hyperball packings of the truncated regular tetrahedron tilings. These are derived from the Coxeter simplex tilings \{p,3,3\} (7\leq p\in\mathbb{N}) and \{5,3,3,3,3\} in 3 and 5dimensional hyperbolic space. We determined the densest hyperball packing arrangement and its density with congruent hyperballs in \mathbb{H}^{5} and determined the smallest density upper bounds of noncongruent hyperball packings generated by the above tilings in \mathbb{H}^{n},\leavevmode\nobreak\ (n=3,5).
In [24] we deal with the packings derived by horo and hyperballs (briefly hyphor packings) in ndimensional hyperbolic spaces \mathbb{H}^{n} (n=2,3) which form a new class of the classical packing problems. We constructed in the 2 and 3dimensional hyperbolic spaces hyphor packings that are generated by complete Coxeter tilings of degree 1 and we determined their densest packing configurations and their densities. We proved using also numerical approximation methods that in the hyperbolic plane (n=2) the density of the above hyphor packings arbitrarily approximate the universal upper bound of the hypercycle or horocycle packing density \frac{3}{\pi} and in \mathbb{H}^{3} the optimal configuration belongs to the \{7,3,6\} Coxeter tiling with density \approx 0.83267. Furthermore, we analyzed the hyphor packings in truncated orthoschemes \{p,3,6\} (6<p<7,\leavevmode\nobreak\ p\in\mathbb{R}) whose density function is attained its maximum for a parameter which lies in the interval [6.05,6.06] and the densities for parameters lying in this interval are larger that \approx 0.85397.
In [23] we proved that if the truncated tetrahedron is regular, then the density of the densest packing is \approx 0.86338. This is larger than the BöröczkyFlorian density upper bound but our locally optimal hyperball packing configuration cannot be extended to the entirety of \mathbb{H}^{3}. However, we described a hyperball packing construction, by the regular truncated tetrahedron tiling under the extended Coxeter group \{3,3,7\} with maximal density \approx 0.82251.
Recently, (to the best of author’s knowledge) the candidates for the densest hyperball (hypersphere) packings in the 3,4 and 5dimensional hyperbolic space \mathbb{H}^{n} are derived by the regular prism tilings which have been studied in papers [30], [31] and [32].
In [26] we considered hyperball packings in 3dimensional hyperbolic space and developed a decomposition algorithm that for each saturated hyperball packing provides a decomposition of \mathbb{H}^{3} into truncated tetrahedra. Therefore, in order to get a density upper bound for hyperball packings, it is sufficient to determine the density upper bound of hyperball packings in truncated simplices.
In [35] we studied hyperball packings related to truncated regular octahedron and cube tilings that are derived from the Coxeter simplex tilings \{p,3,4\} (7\leq p\in\mathbb{N}) and \{p,4,3\} (5\leq p\in\mathbb{N}) in 3dimensional hyperbolic space \mathbb{H}^{3}. We determined the densest hyperball packing arrangement and its density with congruent and noncongruent hyperballs related to the above tilings. Moreover, we prove that the locally densest congruent or noncongruent hyperball configuration belongs to the regular truncated cube with density \approx 0.86145. This is larger than the BöröczkyFlorian density upper bound for balls and horoballs. We described a noncongruent hyperball packing construction, by the regular cube tiling under the extended Coxeter group \{4,3,7\} with maximal density \approx 0.84931.
In [36] we examined congruent and noncongruent hyperball packings generated by doubly truncated Coxeter orthoscheme tilings in the 3dimensional hyperbolic space. We proved that the densest congruent hyperball packing belongs to the Coxeter orthoscheme tiling of parameter \{7,3,7\} with density \approx 0.81335. This density is equal – by our conjecture – with the upper bound density of the corresponding noncongruent hyperball arrangements.
Remark 1.2
If we try to define the density of system of sets in hyperbolic space as we did in Euclidean space, i.e. by the limiting value of the density with respect to a sphere C(r) of radius r with a fixed centre O. But since for a fixed value of h the volume of spherical shell C(r+h)C(r) is the same order of magnitude as the volume of C(r), the argument used in Euclidean space to prove that the limiting value is independent of the choice of O is does not work in hyperbolic space. Therefore the definition of packing density is crucial in hyperbolic spaces \mathbb{H}^{n} as shown by K. Böröczky [3], for nice examples also see [5], [20]. The most widely accepted notion of packing density considers the local densities of balls with respect to their Dirichlet–Voronoi cells (cf. [3] and [11]), but in our cases these cells are infinite hyperbolic polyhedra. The other possibility: the packing density \delta can be defined (see [38], [39], [30], [32]) as the reciprocal of the ratio of the volume of a fundamental domain for the symmetry group of a tiling to the volume of the ball pieces contained in the fundamental domain (\delta<1). Similarly is defined the covering density \Delta>1. In the present paper our aim is to determine a density upper bound for saturated, congruent hyperball packings in \mathbb{H}^{3} therefore we use an extended notion of such local density.
2 Saturated hyperball packings in \mathbb{H}^{3} and their density upper bound
We use for \mathbb{H}^{3} (and analogously for \mathbb{H}^{n}, n\geq 3) the projective model in the Lorentz space \mathbb{E}^{1,3} that denotes the real vector space \mathbf{V}^{4} equipped with the bilinear form of signature (1,3), \langle\mathbf{x},\leavevmode\nobreak\ \mathbf{y}\rangle=x^{0}y^{0}+x^{1}y^{1% }+x^{2}y^{2}+x^{3}y^{3}, where the nonzero vectors \mathbf{x}=(x^{0},x^{1},x^{2},x^{3})\in\mathbf{V}^{4}\ \ \text{and}\ \ \mathbf% {y}=(y^{0},y^{1},y^{2},y^{3})\in\mathbf{V}^{4}, are determined up to real factors, for representing points of \mathcal{P}^{n}(\mathbb{R}). Then \mathbb{H}^{3} can be interpreted as the interior of the quadric Q=\{(\mathbf{x})\in\mathcal{P}^{3}\langle\mathbf{x},\leavevmode\nobreak\ % \mathbf{x}\rangle=0\}=:\partial\mathbb{H}^{3} in the real projective space \mathcal{P}^{3}(\mathbf{V}^{4},\mbox{\boldmath$V$}\!_{4}) (here \mbox{\boldmath$V$}\!_{4} is the dual space of \mathbf{V}^{4}). Namely, for an interior point \mathbf{y} holds \langle\mathbf{y},\leavevmode\nobreak\ \mathbf{y}\rangle<0.
Points of the boundary \partial\mathbb{H}^{3} in \mathcal{P}^{3} are called points at infinity, or at the absolute of \mathbb{H}^{3}. Points lying outside \partial\mathbb{H}^{3} are said to be outer points of \mathbb{H}^{3} relative to Q. Let (\mathbf{x})\in\mathcal{P}^{3}, a point (\mathbf{y})\in\mathcal{P}^{3} is said to be conjugate to (\mathbf{x}) relative to Q if \langle\mathbf{x},\leavevmode\nobreak\ \mathbf{y}\rangle=0 holds. The set of all points which are conjugate to (\mathbf{x}) form a projective (polar) hyperplane pol(\mathbf{x}):=\{(\mathbf{y})\in\mathcal{P}^{3}\langle\mathbf{x},% \leavevmode\nobreak\ \mathbf{y}\rangle=0\}. Thus the quadric Q induces a bijection (linear polarity \mathbf{V}^{4}\rightarrow\mbox{\boldmath$V$}\!_{4}) from the points of \mathcal{P}^{3} onto their polar hyperplanes.
Point X(\mathbb{x}) and hyperplane \alpha(\mbox{\boldmath$a$}) are incident if \mathbb{x}\mbox{\boldmath$a$}=0 (\mathbb{x}\in\mathbb{V}^{4}\setminus\{\mathbf{0}\},\ \mbox{\boldmath$a$}\in% \mbox{\boldmath$V$}_{4}\setminus\{\mbox{\boldmath$0$}\}).
The hypersphere (or equidistance surface) is a quadratic surface at a constant distance from a plane (base plane) in both halfspaces. The infinite body of the hypersphere, containing the base plane, is called hyperball.
The half hyperball with distance h to a base plane \beta is denoted by \mathcal{H}^{h}_{+}. The volume of a bounded hyperball piece \mathcal{H}^{h}_{+}(\mathcal{A}), delimited by a 2polygon \mathcal{A}\subset\beta, and its prism orthogonal to \beta, can be determined by the classical formula (2.1) of J. Bolyai [2].
\mathrm{Vol}(\mathcal{H}^{h}_{+}(\mathcal{A}))=\frac{1}{4}\mathrm{Area}(% \mathcal{A})\left[k\sinh\frac{2h}{k}+2h\right],  (2.1) 
The constant k=\sqrt{\frac{1}{K}} is the natural length unit in \mathbb{H}^{3}, where K denotes the constant negative sectional curvature. In the following we may assume that k=1.
Let \mathcal{B}^{h} be a hyperball packing in \mathbb{H}^{3} with congruent hyperballs of height h.
The notion of saturated packing follows from that fact that the density of any packing can be improved by adding further packing elements as long as there is sufficient room to do so. However, we usually apply this notion for packings with congruent elements.
In [26] we modified the classical definition of saturated packing for noncompact ball packings with generalized balls (horoballs, hyperballs) in ndimensional hyperbolic space \mathbb{H}^{n} (n\geq 2 integer parameter):
Definition 2.1
A ball packing with noncompact generalized balls (horoballs or/and hyperballs) in \mathbb{H}^{n} is saturated if no new noncompact generalized ball can be added to it.
To obtain hyperball (hypersphere) packing upper bound it obviously suffices to study saturated hyperball packings (using the above definition) and in what follows we assume that all packings are saturated unless otherwise stated.
We take the set of hyperballs \{\mathcal{H}^{h}_{i}\} of a saturated hyperball packing \mathcal{B}^{h} (see Definition 2.1). Their base planes are denoted by \beta_{i}. Thus in a saturated hyperball packing the distance between two ultraparallel base planes d(\beta_{i},\beta_{j}) is at least 2h (where for the natural indices holds i<j and d is the hyperbolic distance function).
In [26] we described a procedure to get a decomposition of 3dimensional hyperbolic space \mathbb{H}^{3} into truncated tetrahedra corresponding to a given saturated hyperball packing whose main steps were the following:

Using the radical planes of the hyperballs \mathcal{H}^{h}_{i}, similarly to the Euclidean space, can be constructed the unique DirichletVoronoi (in short DV) decomposition of \mathbb{H}^{3} to the given hyperball packing \mathcal{B}^{h}.

We consider an arbitrary proper vertex P\in\mathbb{H}^{3} of the above DV decomposition and the hyperballs \mathcal{H}^{h}_{i}(P) whose DV cells meet at P. The base planes of the hyperballs \mathcal{H}^{h}_{i}(P) are denoted by \beta_{i}(P), and these planes determine a noncompact polyhedron \mathcal{D}^{i}(P) with the intersection of their halfspaces containing the vertex P. Moreover, denote A_{1},A_{2},A_{3},\dots the outer vertices of \mathcal{D}^{i}(P) and cut off \mathcal{D}^{i}(P) with the polar planes \alpha_{j}(P) of its outer vertices A_{j}. Thus, we obtain a convex compact polyhedron \mathcal{D}(P). This is bounded by the base planes \beta_{i}(P) and ”polar planes” \alpha_{j}(P). Applying this procedure for all vertices of the above DirichletVoronoi decomposition, we obtain an other decomposition of \mathbb{H}^{3} into convex polyhedra.

We consider \mathcal{D}(P) as a tile of the above decomposition. The planes from the finite set of base planes \{\beta_{i}(P)\} are called adjacent if there is a vertex A_{s} of \mathcal{D}^{i}(P) that lies on each of the above plane. We consider nonadjacent planes \beta_{k_{1}}(P),\beta_{k_{2}}(P),\beta_{k_{3}}(P),\dots\beta_{k_{m}}(P)\in\{% \beta_{i}(P)\} (k_{l}\in\mathbb{N}^{+},\leavevmode\nobreak\ l=1,2,3,\dots m) that have an outer point of intersection denoted by A_{k_{1}\dots k_{m}}. Let N_{\mathcal{D}(P)}\in\mathbb{N} denote the finite number of the outer points A_{k_{1}\dots k_{m}} related to \mathcal{D}(P). It is clear, that its minimum is 0 if \mathcal{D}^{i}(P) is tetrahedron. The polar plane \alpha_{k_{1}\dots k_{m}} of A_{k_{1}\dots k_{m}} is orthogonal to planes \beta_{k_{1}}(P),\beta_{k_{2}}(P),\dots\beta_{k_{m}}(P) (thus it contain their poles B_{k_{1}}, B_{k_{2}},…B_{k_{m}}) and divides \mathcal{D}(P) into two convex polyhedra \mathcal{D}_{1}(P) and \mathcal{D}_{2}(P).

If N_{\mathcal{D}_{1}(P)}\neq 0 and N_{\mathcal{D}_{2}(P)}\neq 0 then N_{\mathcal{D}_{1}(P)}<N_{\mathcal{D}(P)} and N_{\mathcal{D}_{2}(P)}<N_{\mathcal{D}(P)} then we apply the point 3 for polyhedra \mathcal{D}_{i}(P),\leavevmode\nobreak\ i\in\{1,2\}.

If N_{\mathcal{D}_{i}(P)}\neq 0 or N_{\mathcal{D}_{j}(P)}=0 (i\neq j,\leavevmode\nobreak\ i,j\in\{1,2\}) then we consider the polyhedron \mathcal{D}_{i}(P) where N_{\mathcal{D}_{i}(P)}=N_{\mathcal{D}(P)}1 because the vertex A_{k_{1}\dots k_{m}} is left out and apply the point 3.

If N_{\mathcal{D}_{1}(P)}=0 and N_{\mathcal{D}_{2}(P)}=0 then the procedure is over for \mathcal{D}(P). We continue the procedure with the next cell.

We have seen in steps 3, 4, 5 and 6 that the number of the outer vertices A_{k_{1}\dots k_{m}} of any polyhedron obtained after the cutting process is less than the original one, and we have proven in step 7 that the original hyperballs form packings in the new polyhedra \mathcal{D}_{1}(P) and \mathcal{D}_{2}(P), as well. We continue the cutting procedure described in step 3 for both polyhedra \mathcal{D}_{1}(P) and \mathcal{D}_{2}(P). If a derived polyhedron is a truncated tetrahedron then the cutting procedure does not give new polyhedra, thus the procedure will not be continued. Finally, after a finite number of cuttings we get a decomposition of \mathcal{D}(P) into truncated tetrahedra, and in any truncated tetrahedron the corresponding congruent hyperballs from \{\mathcal{H}^{h}_{i}\} form a packing. Moreover, we apply the above method for the further cells.
From the above algorithm we obtained the following
Theorem 2.2 (J. Sz. [26])
The in [26] described algorithm provides for each congruent saturated hyperball packing a decomposition of \mathbb{H}^{3} into truncated tetrahedra. \square
Remark 2.3
Przeworski, A. proved a similar theorem in [19] but it was true only for cases if the base planes of hyperspheres form “symmetric cocompacts arrangements” in \mathbb{H}^{n}.
In [16] Y. Miyamoto proved the analogue theorem of K. Böröczky’s theorem (Theorem 1.1):
Theorem 2.4 (Y. Miyamoto, [16])
If a region in \mathbb{H}^{n} bounded by hyperplanes has a hyperball (hypersphere) packing of height (radius) r about its boundary, then in some sense, the ratio of its volume to the volume of its boundary is at least that of a regular truncated simplex of (inner) edgelength 2r.
Remark 2.5
Independently from the above paper A. Przeworski proved a similar theorem with other methods in [19].
Therefore, in order to get density upper bound related to the saturated hyperball packings it is sufficient to determine the density upper bound of hyperball packings in truncated regular simplices (see Fig. 1).
Thus, in the following we assume that the ultraparallel base planes \beta_{i} of \mathcal{H}^{h(p)}_{i} (i=1,2,3,4, and 6<p\in\mathbb{R}) generate a “regular truncated tetrahedron” \mathcal{S}(p) with outer vertices B_{i} (see Fig. 1) whose nonorthogonal dihedral angles are equal to \frac{2\pi}{p}, and the distances between two base planes d(\beta_{i},\beta_{j})=:e_{ij} (i<j\in\{1,2,3,4\}) are equal to 2h(p) depending on the angle \frac{\pi}{p}.
The truncated regular tetrahedron \mathcal{S}(p) can be decomposed into 24 congruent simply truncated orthoschemes; one of them \mathcal{O}=Q_{0}Q_{1}Q_{2}P_{0}P_{1}P_{2} is illustrated in Fig. 1 where P_{0} is the centre of the “regular tetrahedron” \mathcal{S}(p), P_{1} is the centre of a hexagonal face of \mathcal{S}(p), P_{2} is the midpoint of a “common perpendicular” edge of this face, Q_{0} is the centre of an adjacent regular triangle face of \mathcal{S}(p), Q_{1} is the midpoint of an appropriate edge of this face and one of its endpoints is Q_{2}.
In our case the essential dihedral angles of orthoschemes \mathcal{O} are the following: \alpha_{01}=\frac{\pi}{p},\ \ \alpha_{12}=\frac{\pi}{3},\ \ \alpha_{23}=\frac{% \pi}{3}. Therefore, the volume \mathrm{Vol}(\mathcal{O}) of the orthoscheme \mathcal{O} and the volume \mathrm{Vol}(\mathcal{S}(p))=24\cdot\mathrm{Vol}(\mathcal{O}) can be computed for any given parameter p (6<p\in\mathbb{R}) by Theorem 2.6.
Theorem 2.6 (R. Kellerhals, [10], Theorem II.)
The volume of a threedimensional hyperbolic complete orthoscheme (except Lambert cube cases, i.e. complete orthoschemes of degree m=2 with outer edge) \mathcal{O}\subset\mathbb{H}^{3} is expressed with the essential angles \alpha_{01},\alpha_{12},\alpha_{23},\ (0\leq\alpha_{ij}\leq\frac{\pi}{2}) in the following form:
\displaystyle\mathrm{Vol}(\mathcal{O})=\frac{1}{4}\{\mathcal{L}(\alpha_{01}+% \theta)\mathcal{L}(\alpha_{01}\theta)+\mathcal{L}(\frac{\pi}{2}+\alpha_{12}% \theta)+  
\displaystyle+\mathcal{L}(\frac{\pi}{2}\alpha_{12}\theta)+\mathcal{L}(\alpha% _{23}+\theta)\mathcal{L}(\alpha_{23}\theta)+2\mathcal{L}(\frac{\pi}{2}% \theta)\}, 
where \theta\in[0,\frac{\pi}{2}) is defined by:
\tan(\theta)=\frac{\sqrt{\cos^{2}{\alpha_{12}}\sin^{2}{\alpha_{01}}\sin^{2}{% \alpha_{23}}}}{\cos{\alpha_{01}}\cos{\alpha_{23}}}, 
and where \mathcal{L}(x):=\int\limits_{0}^{x}\log{2\sin{t}}dt denotes the Lobachevsky function.
In this case for a given parameter p the length of the common perpendiculars h(p)=\frac{1}{2}e_{ij} (i<j, i,j\in\{1,2,3,4\}) can be determined by the machinery of projective metric geometry. (In the following \mathbf{x}\sim c\cdot\mathbf{x} with c\in\mathbb{R}\setminus\{\mathbf{0}\} represent the same point X=(\mathbf{x}\sim c\cdot\mathbf{x}) of \mathcal{P}^{3}.)
The points P_{2}({\mathbf{p}}_{2}) and Q_{2}({\mathbf{q}}_{2}) are proper points of hyperbolic 3space and Q_{2} lies on the polar hyperplane pol(B_{1})(\mbox{\boldmath$b$}^{1}) of the outer point B_{1} thus
The hyperbolic distance h(p) can be calculated by the following formula (see[23]):
\begin{gathered}\displaystyle\cosh{h(p)}=\cosh{P_{2}Q_{2}}=\frac{\langle{% \mathbf{q}}_{2},{\mathbf{p}}_{2}\rangle}{\sqrt{\langle{\mathbf{q}}_{2},{% \mathbf{q}}_{2}\rangle\langle{\mathbf{p}}_{2},{\mathbf{p}}_{2}\rangle}}=\\ \displaystyle=\frac{h_{23}^{2}h_{22}h_{33}}{\sqrt{h_{22}\langle\mathbf{q}_{2}% ,\mathbf{q}_{2}\rangle}}=\sqrt{\frac{h_{22}\leavevmode\nobreak\ h_{33}h_{23}^% {2}}{h_{22}\leavevmode\nobreak\ h_{33}}},\end{gathered} 
where h_{ij} is the inverse of the CoxeterSchläfli matrix
(c^{ij}):=\begin{pmatrix}1&\cos{\frac{\pi}{p}}&0&0\\ \cos{\frac{\pi}{p}}&1&\cos{\frac{\pi}{3}}&0\\ 0&\cos{\frac{\pi}{3}}&1&\cos{\frac{\pi}{3}}\\ 0&0&\cos{\frac{\pi}{3}}&1\\ \end{pmatrix} 
of the orthoscheme \mathcal{O}. We get that the volume \mathrm{Vol}(\mathcal{S}(p)), the maximal height h(p) of the congruent hyperballs lying in \mathcal{S}(p) and \mathrm{Vol}(\mathcal{H}^{h(p)}\cap\mathcal{S}(p)) all depend only on the parameter p of the truncated regular tetrahedron \mathcal{S}(p).
Therefore, the locally optimal density of congruent hyperball packing related to the regular truncated tetrahedron of parameter p is
\delta(\mathcal{S}(p)):=\frac{4\cdot\mathrm{Vol}(\mathcal{H}^{h(p)}\cap% \mathcal{S}(p))}{\mathrm{Vol}({\mathcal{S}(p)})}, 
and \delta(\mathcal{S}(p)) depends only on p (6<p\in\mathbb{R}). Moreover, the total volume of the parts of the four hyperballs lying in \mathcal{S}(p) can be computed by formula (2.1), and the volume of \mathcal{S}(p) can be determined by Theorem 2.6.
Finally, we obtain the plot after careful analysis of the smooth density function (cf. Fig. 2) and we obtain the following
Theorem 2.7 (J. Sz. [23])
The density function \delta(\mathcal{S}(p)), p\in(6,\infty) attains its maximum at p^{opt}\approx 6.13499, and \delta(\mathcal{S}(p)) is strictly increasing in the interval (6,p^{opt}), and strictly decreasing in (p^{opt},\infty). Moreover, the optimal density \delta^{opt}(\mathcal{S}(p^{opt}))\approx 0.86338 (see Fig. 2).
Remark 2.8

In our case \lim_{p\rightarrow 6}(\delta(\mathcal{S}(p))) is equal to the BöröczkyFlorian upper bound of the ball and horoball packings in \mathbb{H}^{3} [4] (observe that the dihedral angles of \mathcal{S}(p) for the case of the horoball equal 2\pi/6).

\delta^{opt}(\mathcal{S}(p^{opt}))\approx 0.86338 is larger than the BöröczkyFlorian upper bound \delta_{BF}\approx 0.85328; but these hyperball packing configurations are only locally optimal and cannot be extended to the entire hyperbolic space \mathbb{H}^{3}.
We obtain the next theorem as the direct consequence of the previous statements:
Theorem 2.9
The density upper bound of the saturated congruent hyperball packings related to the corresponding truncated tetrahedron cells is realized in a regular truncated tetrahedra belonging to parameter p^{opt}\approx 6.13499 with density \approx 0.86338.
We get from the above theorem directly the denial of the A. Przeworski’s conjecture [19]:
Corollary 2.10
The density function \delta(\mathcal{S}(p)), is not an increasing function of h(p) (the height of hyperballs).
Remark 2.11
The hyperball packings in the regular truncated tetrahedra under the extended reflection groups with CoxeterSchläfli symbol \{3,3,p\}, investigated in paper [23], can be extended to the entire hyperbolic space if p is an integer parameter bigger than 6. They coincide with the hyperball packings given by the regular pgonal prism tilings in \mathbb{H}^{3} with extended CoxeterSchläfli symbols \{p,3,3\}, see in [30]. As we know, \{3,3,p\} and \{p,3,3\} are dually isomorphic extended reflection groups, just with the above frustum of orthoscheme as fundamental domain (Fig. 1, matrix (c^{ij}) in formula (2.2)).
In [23] we studied these tilings and the corresponding hyperball packings. Moreover, we computed their metric data for some integer parameters p (6<p\in\mathbb{N}), where \mathcal{A} is a trigonal face of the regular truncated tetrahedron, cf. Fig. 1. In the Table 1 we recalled from [23] important metric data of some “realizable hyperball packings”.
Table 1, p h(p) \mathrm{Vol}(\mathcal{O}) \mathrm{Vol}(\mathcal{H}^{h}_{+}(\mathcal{A})) \delta(\mathcal{S}(p)) 7 0.78871 0.08856 0.07284 0.82251 8 0.56419 0.10721 0.08220 0.76673 9 0.45320 0.11825 0.08474 0.71663 \vdots \vdots \vdots \vdots \vdots 20 0.16397 0.14636 0.06064 0.41431 \vdots \vdots \vdots \vdots \vdots 50 0.06325 0.15167 0.02918 0.19240 \vdots \vdots \vdots \vdots \vdots 100 0.03147 0.15241 0.01549 0.10165 p\to\infty 0 0.15266 0 0
In hyperbolic spaces \mathbb{H}^{n} (n\geq 3) the problems of the densest ball, horoball and hyperball packings have not been settled yet, in general (see e.g. [13], [28], [29]). Moreover, the optimal sphere packing problem can be extended to the other homogeneous Thurston geometries, e.g. \mathbf{Nil}, \mathbf{Sol}, \widetilde{\mathbf{S}\mathbf{L}_{2}\mathbf{R}}. For these nonEuclidean geometries only very few results are known (e.g. [34] and the references given there).
By the above these we can say that the revisited Kepler problem keeps yet several interesting open questions.
I thank Prof. Emil Molnár for his helpful comments and suggestions to this paper.
References
 [1] Bezdek, K. Sphere Packings Revisited, Eur. J. Combin., 27/6 (2006), 864–883.
 [2] Bolyai, J. Appendix, Scientiam spatii absolute veram exhibens, Marosvásárhely, (1831).
 [3] Böröczky, K. Packing of spheres in spaces of constant curvature, Acta Math. Acad. Sci. Hungar., 32 (1978), 243–261.
 [4] Böröczky, K. — Florian, A. Über die dichteste Kugelpackung im hyperbolischen Raum, Acta Math. Acad. Sci. Hungar., 15 (1964), 237–245.
 [5] Fejes Tóth, G.  Kuperberg, W. Packing and Covering with Convex Sets, Handbook of Convex Geometry Volume B, eds. Gruber, P.M., Willis J.M., pp. 799860, NorthHolland, (1983).
 [6] Fejes Tóth, G. — Kuperberg, G. — Kuperberg, W. Highly Saturated Packings and Reduced Coverings, Monatsh. Math., 125/2 (1998), 127–145.
 [7] Fejes Tóth, L. Regular Figures, Macmillan (New York), 1964.
 [8] Hales, T. C. Historical Overview of the Kepler Conjecture, Discrete and Computational Geometry, 35 (2006), 5–20.
 [9] Im Hof, H.C. Napier cycles and hyperbolic Coxeter groups, Bull. Soc. Math. Belgique, 42 (1990), 523–545.
 [10] Kellerhals, R. On the volume of hyperbolic polyhedra, Math. Ann., 285 (1989), 541–569.
 [11] Kellerhals, R. Ball packings in spaces of constant curvature and the simplicial density function, J. Reine Angew. Math., 494 (1998), , 189–203.
 [12] R.T. Kozma  J. Szirmai, Optimally dense packings for fully asymptotic Coxeter tilings by horoballs of different types, Monatsh. Math., 168/1 (2012), 27–47.
 [13] Kozma, R.T., Szirmai, J. New Lower Bound for the Optimal Ball Packing Density of Hyperbolic 4space, Discrete Comput. Geom., 53/1 (2015), 182198, DOI: 10.1007/s0045401496341.
 [14] Kozma, R.T., Szirmai, J. New Lower Bound for the Optimal Ball Packing Density of Hyperbolic 4space, Submitted Manuscript (2018).
 [15] Marshall, T. H., Martin, G. J. Packing strips in the hyperbolic plane, Proc. Edinb. Math. Soc. (2), 46/1 (2003), 67–73.
 [16] Miyamoto, Y. Volumes of hyperbolic manifolds with geodesic boundary, Topology, 33/4 (1994), 613–629.
 [17] Molnár, E. The Projective Interpretation of the eight 3dimensional homogeneous geometries, Beitr. Algebra Geom., 38/2 (1997), 261–288.
 [18] Molnár, E. — Szirmai, J. Top dense hyperbolic ball packings and coverings for complete Coxeter orthoscheme groups, Publications de l’Institut Mathématique, 103(117) (2018), 129–146, DOI: 10.2298/PIM1817129M, arXiv: 1612.04541.
 [19] Przeworski, A. An upper bound on density for packings of collars about hyperplanes in \mathbb{H}^{n}, Geometriae Dedicata, 163/1 (2013), 193–213, DOI: 10.1007/s107110129745x.
 [20] Radin, C. The symmetry of optimally dense packings, NonEucledian Geometries, eds. A. Prékopa, E. Molnár, pp. 197207, Springer Verlag, (2006).
 [21] Stojanović, M. Coxeter Groups as Automorphism Groups of Solid Transitive 3simplex Tilings, Filomat, 28/3 (2014), 557–577, DOI: 10.2298/FIL1403557S.
 [22] Stojanović, M. Hyperbolic space groups and their supergroups for fundamental simplex tilings, Acta Math. Hungar., 153/2 (2017), 276–288, DOI: 10.1007/s104740170761z.
 [23] Szirmai, J. Hyperball packings in hyperbolic 3space, Mat. Vesn., 70/3 (2018), 211–221.
 [24] Szirmai, J. Packings with horo and hyperballs generated by simple frustum orthoschemes, Acta Math. Hungar., 152/2 (2017), 365–382, DOI:10.1007/s1047401707280.
 [25] Szirmai, J. Density upper bound of congruent and noncongruent hyperball packings generated by truncated regular simplex tilings, Rendiconti del Circolo Matematico di Palermo Series 2, 67 (2018), 307–322, DOI: 10.1007/s1221501703168, arXiv:1510.03208.
 [26] Szirmai, J. Decomposition method related to saturated hyperball packings, Ars Math. Contemp., (to appear) (2018).
 [27] Szirmai, J. The optimal ball and horoball packings to the Coxeter honeycombs in the hyperbolic dspace, Beitr. Algebra Geom., 48/1 (2007), 35–47.
 [28] Szirmai, J. Horoball packings to the totally asymptotic regular simplex in the hyperbolic nspace, Aequat. Math., 85 (2013), 471482, DOI: 10.1007/s0001001201586.
 [29] Szirmai, J. Horoball packings and their densities by generalized simplicial density function in the hyperbolic space, Acta Math. Hungar., 136/12 (2012), 39–55, DOI: 10.1007/s1047401202058.
 [30] Szirmai, J. The pgonal prism tilings and their optimal hypersphere packings in the hyperbolic 3space, Acta Math. Hungar., 111 (12) (2006), 65–76.
 [31] Szirmai, J. The regular prism tilings and their optimal hyperball packings in the hyperbolic nspace, Publ. Math. Debrecen, 69 (12) (2006), 195–207.
 [32] Szirmai, J. The optimal hyperball packings related to the smallest compact arithmetic 5orbifolds, Kragujevac J. Math. 40(2) (2016), 260270, DOI:10.5937/KgJMath1602260S.
 [33] Szirmai, J. The least dense hyperball covering to the regular prism tilings in the hyperbolic nspace, Ann. Mat. Pur. Appl. 195/1 (2016), 235–248, DOI: 10.1007/s1023101404600.
 [34] Szirmai, J. A candidate for the densest packing with equal balls in Thurston geometries, Beitr. Algebra Geom., 55/2 (2014), 441 452, DOI: 10.1007/s1336601301582.
 [35] Szirmai, J. Hyperball packings related to octahedron and cube tilings in hyperbolic space, Submitted manuscript, (2018).
 [36] Szirmai, J. Congruent and noncongruent hyperball packings related to doubly truncated Coxeter orthoschemes in hyperbolic 3space, Submitted manuscript, (2018).
 [37] Vermes, I. Über die Parkettierungsmöglichkeit des dreidimensionalen hyperbolischen Raumes durch kongruente Polyeder, Studia Sci. Math. Hungar., 7 (1972), 267–278.
 [38] Vermes, I. Ausfüllungen der hyperbolischen Ebene durch kongruente Hyperzykelbereiche, Period. Math. Hungar., 10/4 (1979), 217–229.
 [39] Vermes, I. Über reguläre Überdeckungen der BolyaiLobatschewskischen Ebene durch kongruente Hyperzykelbereiche, Period. Math. Hungar., 25/3 (1981), 249–261.
Budapest University of Technology and Economics, Institute of Mathematics,\hfil\leavevmode\nobreak\ \@open@row