Geodesic ball packings generated by regular prism tilings in geometry ^{1}^{1}1AMS Classification 2000: 52C17, 52C22, 53A35, 51M20
Abstract
In this paper we study the regular prism tilings and construct ball packings by geodesic balls related to the above tilings in the projective model of geometry. Packings are generated by action of the discrete prism groups . We prove that these groups are realized by prism tilings in space if and determine packing density formulae for geodesic ball packings generated by the above prism groups. Moreover, studying these formulae we determine the conjectured maximal dense packing arrangements and their densities and visualize them in the projective model of geometry. We get a dense (conjectured locally densest) geodesic ball arrangement related to the parameters where the kissing number of the packing is , similarly to the densest latticelike geodesic ball arrangement investigated by the second author in [22].
1 Introduction and previous results
In mathematics sphere packing problems concern the arrangements of nonoverlapping equal spheres which fill a space. Usually the space involved is the threedimensional Euclidean space where the famous Kepler conjecture was proved by T. C. Hales and S. P. Ferguson in [5].
However, ball (sphere) packing problems can be generalized to the other dimensional Thurston geometries.
In an dimensional space of constant curvature , , let be the density of spheres of radius mutually touching one another with respect to the simplex spanned by the centres of the spheres. L. Fejes Tóth and H. S. M. Coxeter conjectured that in an dimensional space of constant curvature the density of packing spheres of radius can not exceed . This conjecture has been proved by C. Roger in the Euclidean space. The 2dimensional case has been solved by L. Fejes Tóth. In an dimensional space of constant curvature the problem has been investigated by Böröczky and Florian in [2] and it has been studied by K. Böröczky in [3] for dimensional space of constant curvature .
In [6], [7], [24] and [27] we have studied some new aspects of the horoball and hyperball packings in and we have observed that the ball, horoball and hyperball packing problems are not settled yet in the dimensional hyperbolic space.
In [28] we generalized the above problem of finding the densest geodesic and translation ball (or sphere) packing to the other dimensional homogeneous geometries (Thurston geometries)
and in the papers [21], [22], [23], [25], [28] we investigated several interesting ball packing and covering problems in the above geometries. We described in geometry (see [28]) a candidate of the densest geodesic and translation ball arrangement whose density is .
In this paper we consider the geometry that can be derived from W. Heisenberg’s famous real matrix group. This group provides a noncommutative translation group of an affine 3space. E. Molnár proved in [10], that the homogeneous 3spaces have a unified interpretation in the projective 3sphere . In this work we will use this projective model of the geometry.
In [22] we investigated the geodesic balls of the space and computed their volume, introduced the notion of the lattice, parallelepiped and the density of the latticelike ball packing. Moreover, we determined the densest latticelike geodesic ball packing. The density of this densest packing is , may be surprising enough in comparison with the Euclidean result . The kissing number of the balls in this packing is .
In [25] we considered the analogue question for translation balls. The notions of translation curve and translation ball were introduced by initiative of E. Molnár (see [14], [21]). We have studied the translation balls of space and computed their volume. Moreover, we have proved that the density of the optimal latticelike translation ball packing for every natural lattice parameter is in interval and if then the optimal density is . Meanwhile we can apply a nice general estimate of L. Fejes Tóth [4] in our proof. The kissing number of the latticelike ball packings is less than or equal to 14 and the optimal ball packing is realizable in case of equality. We formulated a conjecture for , where the density of the conjectural densest packing is for lattice parameter , larger than the Euclidean one (), but less than the density of the densest latticelike geodesic ball packing in space known till now [22]. The kissing number of the translation balls in that packing is 14 as well.
In [28] we studied one type of lattice coverings in the space. We introduced the notion of the density of considered coverings and gave upper and lower estimation to the density of the lowest latticelike geodesic ball covering. Moreover we formulate a conjecture for the ball arrangement of the least dense latticelike geodesic ball covering and give its covering density .
In this paper we study the regular prism tilings and construct ball packings by geodesic balls related to the prism tilings in the projective model of geometry where the packings are generated by action of the discrete prism groups . We obtain density formulae for calculations for geodesic ball packings. Analyzing these density functions we obtain a conjecture for optimal geodesic ball packing configurations and determine their densities related to the above prismatic tessellations.
The results are summarized in Theorems 3.4, 4.3 and Conjecture 4.4, and the optimal ball configurations are visualized in our model with Figures 2, 3.
2 Basic notions of the geometry
In this Section we summarize the significant notions and denotations of the geometry (see [10], [22]).
The geometry is a homogeneous 3space derived from the famous real matrix group discovered by W. Heisenberg. The Lie Theory with the method of the projective geometry makes possible to investigate and to describe this topic.
The left (rowcolumn) multiplication of Heisenberg matrices
(2.1) 
defines ”translations” on the points of the space . These translations are not commutative in general. The matrices of the form
(2.2) 
constitute the one parametric centre, i.e. each of its elements commutes with all elements of . The elements of are called fibre translations. geometry of the Heisenberg group can be projectively (affinely) interpreted by the ”right translations” on points as the matrix formula
(2.3) 
shows, according to (1.1). Here we consider as projective collineation group with right actions in homogeneous coordinates. We will use the Cartesian homogeneous coordinate simplex ,,, with the unit point which is distinguished by an origin and by the ideal points of coordinate axes, respectively. Moreover, with (or defines a point of the projective 3sphere (or that of the projective space where opposite rays and are identified). The dual system describes the simplex planes, especially the plane at infinity , and generally, defines a plane of (or that of ). Thus defines the incidence of point and plane , as also denotes it. Thus Nil can be visualized in the affine 3space (so in ) as well.
The translation group defined by formula (2.3) can be extended to a larger group of collineations, preserving the fibering, that will be equivalent to the (orientation preserving) isometry group of . In [11] E. Molnár has shown that a rotation trough angle about the axis at the origin, as isometry of , keeping invariant the Riemann metric everywhere, will be a quadratic mapping in to image as follows:
(2.4) 
This rotation formula, however, is conjugate by the quadratic mapping
(2.5) 
i.e. to the linear rotation formula. This quadratic conjugacy modifies the translations in (2.3), as well. We shall use the following important classification theorem.
Theorem 2.1 (E. Molnár [11])

Any group of isometries, containing a 3dimensional translation lattice, is conjugate by the quadratic mapping in (2.5) to an affine group of the affine (or Euclidean) space whose projection onto the (x,y) plane is an isometry group of . Such an affine group preserves a plane point polarity of signature .

Of course, the involutive line reflection about the axis
preserving the Riemann metric, and its conjugates by the above isometries in 1 (those of the identity component) are also isometries. There does not exist orientation reversing isometry.
Remark 2.2
We obtain from the above described projective model a new model of geometry derived by the quadratic mapping . This is the linearized model of space (see [1]).
2.1 Geodesic curves and spheres
The geodesic curves of the geometry are generally defined as having locally minimal arc length between their any two (near enough) points. The equation systems of the parametrized geodesic curves in our model can be determined by the general theory of Riemann geometry. We can assume, that the starting point of a geodesic curve is the origin because we can transform a curve into an arbitrary starting point by translation (2.1);
The arc length parameter is introduced by
i.e. unit velocity can be assumed.
Remark 2.3
Thus we have harmonized the scales along the coordinate axes.
The equation systems of a helixlike geodesic curves if :
(2.6) 
In the cases the geodesic curve is the following:
(2.7) 
The cases are trivial: .
Definition 2.4
The distance between the points and is defined by the arc length of geodesic curve from to .
In our work [22] we introduced the following definitions:
Definition 2.5
The geodesic sphere of radius with centre at the point is defined as the set of all points in the space with the condition . Moreover, we require that the geodesic sphere is a simply connected surface without selfintersection in the space.
Remark 2.6
We shall see that this last condition depends on radius .
Definition 2.7
The body of the geodesic sphere of centre and of radius in the space is called geodesic ball, denoted by , i.e. iff .
Remark 2.8
Henceforth, typically we choose the origin as centre of the sphere and its ball, by the homogeneity of .
We have denoted by the body of the sphere , furthermore we have denoted their volumes by .
In [22] we have proved the the following theorem:
Theorem 2.9
The geodesic sphere and ball of radius exists in the space if and only if
We obtain the volume of the geodesic ball of radius by the following integral (see 2.8):
(2.8) 
The parametric equation system of the geodesic sphere in our model (see [22]):
(2.9) 
We have obtained by the derivatives of these parametrically represented functions (by intensive and careful computations with Maple through the second fundamental form) the following theorem (see [22]):
Theorem 2.10
The geodesic ball is convex in affineEuclidean sense in our model if and only if .
3 prisms and prism tilings
The prisms and prismlike tilings have been thoroughly investigated in and spaces in papers [16], [19], [30]. Here we consider the analogous problem in space. We will use the in 2. section described projective model of geometry. In the following the plane of axis are called base plane of the model and if we say plane then it is a plane in Euclidean sense.
Definition 3.1
Let be an infinite solid bounded by planes, that are determined by fibrelines passing through the points of a gon (, integer parameter) lying in the baseplane. The images of by isometries are called infinite sided prisms.
The common part of with the base plane is defined as the base figure of the prism.
Let be the image of the base plane in the space (see Remark 2.2) and let be a fibre translation (2.2).
Definition 3.2
Let be an infinite sided prism, that is trimmed by the surface and its translated copy . The parts of and inside the infinite prism are called cover faces and are denoted by and .
The sided bounded prism is the part of between the cover faces and .
Definition 3.3
A bounded or infinite sided prism is said to be regular if its side surfaces are congruent to each other under rotations with angle (see (2.4) and (2.5)) about the central fibre line of the prism.
3.1 Regular bounded prism tilings
In this section we will investigate the existence of regular bounded prism tilings of space. In this case the prism tiles are regular bounded prisms having gonal base figures . The prism itself is a topological polyhedron with vertices, and having at every vertex one gonal cover face and two quadrangle side faces (traced by fibre lines). We are looking such prism tilings of space where at each side edge of the prism (which are fibre lines going through vertices of the base figure) meet prisms regularly, by rotations with angle (, integer parameter).
We shall see in Theorem 3.4 that the regular prism tiling exists for some parameters . Let one of its tiles with with vertices . We may assume that lies on the axis. It is clear that the side curves are derived from each other by rotation about the axis. The corresponding vertices are generated by a fibre translation with a positive real parameter. The cover faces , and the side surfaces form a sided regular prism in . will be generated by its rotational isometry group (if these tiling there exist see Theorem 3.4) which is given by its fundamental domain , . Here, , , and is a piecewise linear topological polyhedron. The group presentation can be determined by a standard procedure [8], called Poincaré algorithm. The generators will pair the bent (piecewise linear) faces of :
mapping onto its neighbours , , , respectively. E.g. for the face a point (relative freely, e.g. in the segment ) is taken. Then the union of triangles , , , will be the face .
Then the image is taken in for the face , as usual. The relations are induced by the edge equivalence classes ; ; ; . So we get the group
(3.1) 
where is a rotation about the fibre line through the origin (the axis), is a rotation about a side fibre line of (through a vertex of its base figure). Notice, that is a screw motion, and thus is the fibre translation connecting the cover faces.
Our first question is the following: For which is ?
The following Theorem answers it:
Theorem 3.4
In there exist regular gonal nonfacetoface prism tilings with isometry group for integer parameters :
the regular triangular prism tiling with ,
the regular square prism tiling with ,
the regular hexagonal prism tiling with ,
and each group has a free parameter .
Proof:
Let be a ”bottom” vertex of the regular bounded prism. Then the other vertices of the bottom cover face can be generated by the rotation formula (see (2.4), (2.5)):
Then the condition for the existence of the tiling is the following:
(3.2) 
where is a rotation about the fibre line through the origin and is a rotation about the side fibre line of through the vertex . means that the corresponding points lie on the same fibre lines.
where and are positive integers. This equation only has the following integer solutions:
We obtain from the above computations, that the existence of the above regular prism tilings is independent from the parameter , so we have proven the Theorem.
Remembering that is the ”vertical” translation of the group, we can also compute the height of the regular bounded prism corresponding to the group tiling, since: where is the origin. Using this, we can also give a metric representation of the group, allowing the visualization of the corresponding prism and prism tiling (see Fig. 1. and Fig. 2.).
4 The optimal geodesic ball packings under group
The sphere packing problem deals with the arrangements of nonoverlapping equal spheres, or balls, which fill the space. While the usual problem is in the dimensional Euclideanspace ), it can be generalized to the other dimensional Thurston spaces (see [29]). In this paper we investigate the optimal ball packings of generated by the above described group.
Let (where or as stated above and ) be a regular prism tiling, and let be one of its tiles that is centered at the origin, with a base face given by the vertices . The corresponding vertices of the prism are generated by fibre translations .
We can assume by symmetry, that the optimal geodesic ball is centered at the origin. The volume of a geodesic ball with radius can be determined by the formula (2.8).
We study only one case of the multiply transitive geodesic ball packings where the fundamental domains of the space groups are not prisms. Let the fundamental domains be derived by the Dirichlet — Voronoi cells (DV cells) where their centers are images of the origin. The volume of the times fundamental domain and of the DV cell is the same, respectively, as in the prism case (for any above fixed). It is easy to see by the formulas (2.5), using the quadratic mapping , that the volume of the Dirichlet — Voronoi cell (or the coresponding prism) is
(4.1) 
These locally densest geodesic ball packings can be determined for all possible fixed integer parameters . The optimal radius is
(4.2) 
where is the geodesic distance function of geometry (see Definition 2.4).
Since the congruent images of under the discrete group cover the space, therefore for the density of the ball packing it is sufficient to relate the volume of the ball to the volume of the prism:
Definition 4.1
The maximal density of the above multiply transitive ball packing for given parameters ( and ):
(4.3) 
For every parameters the locally densest geodesic ball packing can be determined.
If we fixed the parameters and then the distance function is a continuous fuction. Therefore it is easy prove the following Lemma:
Lemma 4.2
In space for the rotation group there always exist for given parameters where
(4.4) 
The system of equations (4.4) in Lemma 4.2 can be solved by numerical methods and the corresponding ball arrangements are denoted by . We obtain  using the formulas (4.13)  that in space for the rotation group the metric data of the godesic ball arrangements are the following:
Theorem 4.3
If the system of equation (4.4) holds then the maximal radii and densities of the optimal ball packings are the following:

If , then , with ,

If , then , with ,

If , then , with .
If we vary the parameter in the above cases then the corresponding radius and the density also change. The following table shows that probably the ball packings with maximal kissing numbers provide the optimal ball packing densities.
Radius  Prism volume  Density  Kissing number  
(3,6)  0.5876  4.1446  0.2063  2 
0.6392  4.9032  0.2246  2  
0.6929  5.7616  0.2438  2  
0.7389  6.5517  0.2593  8  
0.7787  7.8111  0.2558  6  
0.8132  9.0201  0.2525  6  
0.8481  10.3641  0.2495  6  
(4,4)  0.9927  7.8849  0.5283  2 
1.0644  9.0650  0.5678  2  
1.1386  10.3729  0.6090  2  
1.2154  11.8175  0.6512  10  
1.2594  13.4079  0.6404  8  
1.3036  15.1538  0.6295  8  
1.3480  17.0647  0.6194  8  
(6,3)  1.6934  34.4141  0.6190  2 
1.7801  38.0287  0.6537  2  
1.8690  41.9209  0.6897  2  
1.9601  46.1044  0.7272  14  
2.0087  50.5935  0.7153  12  
2.0573  55.4028  0.7038  12  
2.1059  60.5470  0.6929  12 
Therefore, we can formulate by the above results the following conjecture:
Conjecture 4.4
The ball arrangements provide the densest ball packing arrangements related to isometry group with parameters .
Remark 4.5
The optimal ball packing in the case of has a kissing number of , which is greater than the maximal kissing number in the Euclidean dimensional space. In fact, this is the second ball packing arrangement in that has this high of a kissing number (see [22]).
References
 [1] K. Brodaczewska: Elementargeometrie in Nil, Dissertation Dr. rer. nat., Fakultät Mathematik und Naturwissenschaften der Technischen Universität Dresden (2014).
 [2] Böröczky, K. – Florian, A. Über die dichteste Kugelpackung im hyperbolischen Raum, Acta Math. Hung., (1964) 15 , 237–245.
 [3] Böröczky, K. Packing of spheres in spaces of constant curvature, Acta Math. Acad. Sci. Hungar., 32 (1978), 243–261.
 [4] Fejes Tï¿th L.: Regular Figures, Pergamon Press (1964).
 [5] Hales T. C. – Ferguson, S. P.: The Kepler conjecture. Discrete and Computational Geometry 36(1), (2006) 1–269.
 [6] Kozma, T. R. – Szirmai, J. Optimally dense packings for fully asymptotic Coxeter tilings by horoballs of different types, Monatshefte für Mathematik, 168, (2012), 2747, DOI: 10.1007/s006050120393x.
 [7] Kozma, T. R. – Szirmai, J. New Lower Bound for the Optimal Ball Packing Density of Hyperbolic 4space, Discrete and Computational Geometry, 53, (2015), 182198, DOI: 10.1007/s0045401496341.
 [8] Molnár, E. – Lucic, Z. Fundamental domains for planar discontinuous groups and uniform tilings, Geometriae Dedicata, (1992) No. 40, 125–143.
 [9] Molnár, E.  Szirmai, J.  Vesnin, A. The optimal packings by translation balls in , J. Geometry,105 (2), (2014) 287–306, DOI: 10.1007/s000220130207x.
 [10] Molnár, E. The projective interpretation of the eight 3dimensional homogeneous geometries. Beitr. Algebra Geom., 38 (1997) No. 2, 261–288.
 [11] Molnár, E. On projective models of Thurston geometries, some relevant notes on orbifolds and manifolds. Siberian Electronic Mathematical Reports, 7 (2010), 491–498, http://mi.mathnet.ru/semr267
 [12] Molnár, E. – Szirmai, J. Symmetries in the 8 homogeneous 3geometries. Symmetry Cult. Sci., 21/13 (2010), 87117.
 [13] Molnár, E. – Szirmai, J. Classification of lattices. Geom. Dedicata, 161/1 (2012), 251275, DOI: 10.1007/s1071101297055.
 [14] Molnár, E. – Szirmai, J. On Nil crystallography, Symmetry: Culture and Science, 17/12 (2006), 55–74.
 [15] Molnár, E – Szirmai, J. Volumes and geodesic ball packings to the regular prism tilings in space Publ. Math. Debrecen, 84/12 (2014), 189–203, DOI: 10.5486/PMD.2014.5832.
 [16] Pallagi, J. – Schultz B. – Szirmai, J. On regular square prism tilings in space, KoG (Scientific and professional journal of Croatian Society for Geometry and Graphics) 16, (2012), 3642.
 [17] Pallagi, J. – Schultz B. – Szirmai, J. Equidistant surfaces in space, Studies of the University of Zilina, Mathematical Series, 25 (2011), 31–40.
 [18] Schultz, B. – Szirmai, J. On parallelohedra of space, Pollack Periodica (Mathematics in Architecture and Civil Engineering Design and Education), 7. Supplement 1 (2012): 129136. ISBN 9789637298448.
 [19] Schultz, B. – Szirmai, J., Densest geodesic ball packings to space groups generated by screw motions, Mediterr. J. Math, (2014), 1–14, DOI: 10.1007/ s000090140513z.
 [20] Scott, P. The geometries of 3manifolds. Bull. London Math. Soc., 15 (1983) 401–487.
 [21] Szirmai, J. The densest translation ball packing by fundamental lattices in space. Beitr. Algebra Geom. 51 No. 2, (2010), 353–373.
 [22] Szirmai, J. The densest geodesic ball packing by a type of lattices. Beitr. Algebra Geom., 48(2) (2007) 383–398.
 [23] Szirmai, J. Geodesic ball packing in space for generalized Coxeter space groups. Beitr. Algebra Geom., 52(2) (2011), 413–430.
 [24] Szirmai, J. Geodesic ball packing in space for generalized Coxeter space groups. Math. Commun., 17/1 (2012), 151170.
 [25] Szirmai J. Latticelike translation ball packings in space. Publ. Math. Debrecen, 80/34 (2012), 427–440 DOI: 10.5486/PMD.2012.5117.
 [26] Szirmai, J. Horoball packings and their densities by generalized simplicial density function in the hyperbolic space. Acta Math. Hung.. 136/12, (2012), 39–55, DOI: 10.1007/s1047401202058.
 [27] Szirmai, J. Horoball packings to the totally asymptotic regular simplex in the hyperbolic space. Aequat. Math.. 85, (2013), 471–482, DOI: 10.1007/s0001001201586.
 [28] Szirmai, J. On lattice coverings of the space by congruent geodesic balls, Mediterr. J. Math. 10/2, (2013) 953–970, DOI: 10.1007/s0000901202117.
 [29] 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.
 [30] Szirmai, J. Regular prism tilings in space, Aequat. Math., (2014) 88/12, 6779, DOI: 10.1007/s000100130221y.
 [31] Thurston, W. P. (and Levy, S. editor) ThreeDimensional Geometry and Topology. Princeton University Press, Princeton, New Jersey, Vol.1 (1997).