## ORBITAL COUNTING OF CURVES ON ALGEBRAIC SURFACES AND SPHERE PACKINGS

## IGOR DOLGACHEV

To the memory of Andrey Todorov

Abstract. We realize the Apollonian group associated to an integral Apollonian circle packings, and some of its generalizations, as a group of automorphisms of an algebraic surface. Borrowing some results in the theory of orbit counting, we study the asymptotic of the growth of degrees of elements in the orbit of a curve on an algebraic surface with respect to a geometrically finite group of its automorphisms.

## 1. Introduction

Let Γ be a discrete subgroup of isometries of a hyperbolic space H n and O ( Γ x 0 ) be its orbit. Consider a family of compact subsets B T of H n whose volume tends to infinity as T →∞ . The problem of finding the asymptotic of #O Γ ( x 0 ) ∩B T is a fundamental problem in harmonic analysis, number theory, ergodic theory and geometry. In this paper we discuss an application of this problem to algebraic geometry.

Let X be a smooth projective algebraic surface and let Num( X ) be the group of divisor classes on X modulo numerical equivalence. For any divisor class D we denote by [ D ] its image in Num( X ). The intersection form on divisor classes defines a non-degenerate symmetric bilinear form on Num( X ) of signature (1 , n ). Let Γ ′ be a group of automorphisms of X such that its image Γ in the orthogonal group O(Num( X )) is an infinite group. This implies that X is birationally isomorphic to an abelian surface, a K3 surface, an Enriques surface or the projective plane.

For any divisor class D we denote by O Γ ( D ) the Γ-orbit of [ D ]. Fix an ample divisor class H and an effective divisor class C . For any positive real number T , let

N T ( H,C ) = # [ { C ′ ] ∈ O ( Γ C ) : H C · ′ ≤ T } = # { H ′ ∈ O ( Γ H ) : H ′ · C ≤ T } .

We are interested in the asymptotic of this function when T goes to infinity. As far as I know, this problem was first considered by Arthur Baragar in his two papers [3], [4]. 1

glyph[negationslash]

In the present paper we explain how does this problem relate to a general orbital counting problem in the theory of discrete subgroups of Lie groups. More precisely, we consider the group G = SO(1 , n ) 0 realized as the group of orientation preserving isometries of the hyperbolic space H n ⊂ P ( V ) associated with the real inner product vector space V = Num( X ) R of signature (1 , n ). We will represent points of H n (resp. points in its boundary ∂ H n , resp. the points in P n \ H n ) by vectors v ∈ V with ( v, v ) = 1 , ( v, h 0 ) &gt; 0 (resp. positive rays of vectors with ( v, v ) = 0 , ( v, h 0 ) &gt; 0, resp. vectors with ( v, v ) = -1 ( , v, h 0 ) &gt; 0), where h 0 is a fixed vector with ( h , h 0 0 ) = 1. Let B T ( e ) be the set of points x in H n whose hyperbolic distance from the point e ∈ H n (resp. from a fixed horosphere with center at e if e ∈ ∂ H n , resp. from the orthogonal hyperplane H e with normal vector e if e ∈ H n ) is less than or equal to T . We use the following result (see [28]) which represents the state of the art in the study of orbital counting.

Theorem 1.1. Assume that Γ is a non-elementary geometrically finite discrete subgroup of orientation preserving isometries of H n . Let δ Γ be the Hausdorff dimension of the limit set Λ(Γ) of Γ . If ( e, e ) ≤ 0 , we additionally assume that the orbit O Γ ([ e ]) is a discrete set and δ Γ &gt; 1 if ( e, e ) &lt; 0 . Then there exists a positive constant c Γ ,h 0 ,e depending only on Γ [ , h ] and e such that

$$\lim _ { T \to \infty } \frac { \# O _ { \Gamma } ( h _ { 0 } ) \cap \mathcal { B } _ { T } } { \exp ( T ) ^ { \delta _ { \Gamma } } } = c _ { \Gamma, h _ { 0 }, e }.$$

We apply this theorem to our situation to obtain the following theorem.

Theorem 1.2. Let X be an abelian surface, or a K3 surface, or an Enriques surface, or a rational surface. Let Γ ′ be a group of automorphisms of X such that its image in O Num ( ( X )) is a non-elementary geometrically finite discrete group. Fix an ample numerical divisor class H ∈ Num ( X ) and an effective numerical divisor class C ∈ Num ( X ) . Assume that δ Γ &gt; 1 if C 2 0 &lt; 0 . Then there exists a positive constant c Γ ,H,C depending only on Γ [ , H ] and [ C ] such that

$$\lim _ { T \to \infty } \frac { N _ { T } ( H, C ) } { T ^ { \delta _ { \Gamma } } } = c _ { \Gamma, H, C }.$$

1 I thank Serge Cantat for these references.

As we see, the main ingredient of the asymptotic expression is the Hausdorff dimension δ Γ of the limit set Λ(Γ) ⊂ ∂ H n of Γ. One of the first and the most beautiful example where δ Γ was computed with many decimals is the example where Γ is an Apollonian group in H 3 whose limit set is the closure of a countable union of circles that intersect in at most one point, an Apollonian gasket . It was proved by D. Boyd [7] that the asymptotic of curvatures of circles in an Apollonian gasket is equal to c T ( ) T δ Γ for some function c T ( ) such that lim T →∞ c T ( ) T = 0. He also showed that 1 &lt; δ Γ &lt; 2. It is known now that δ is about 1 305686729 . ... [39]. An Apollonian sphere packing is a special case of a Boyd-Maxwell sphere packing defined by certain Coxeter groups of hyperbolic type introduced by G. Maxwell in [26]. . 2 The Hausdorff dimension of the limit set of the Boyd-Maxwell Coxeter groups coincides with the sphere packing critical exponent equal to

$$( 1. 1 ) \quad \inf \{ s \, \colon \sum _ { S \in \mathfrak { P } } r ( S ) ^ { s } < \infty \} = \sup \{ s \, \colon \sum _ { S \in \mathfrak { P } } r ( S ) ^ { s } = \infty \},$$

where r S ( ) denotes the radius of a sphere S in a sphere packing P . Although the Hausdorff dimension of the limit set of a geometrically finite group is usually difficult to compute, the bounds for sphere packing critical exponent is easier to find.

We give several concrete examples where the Boyd-Maxwell groups, and, in particular the Apollonian groups, are realized as automorphism groups of algebraic surfaces so that our counting problem can be solved in terms of the sphere packing critical exponent. Since this paper addresses a problem in algebraic geometry, we give a short introduction to the hyperbolic geometry and the theory of sphere packings in a language which we think is more suitable for algebraic geometers. Whenever the details or references are omitted, they can be found in [1], [2], [34], or [40].

It is my pleasure to thank Peter Sarnak for his inspiring talks in Ann Arbor, in April 2014, that gave rise to the present paper. The paper could never be written without the assitance of Jeffrey Lagarias who patiently answered my numerous questions on the theory of Apollonian sphere packings, of Boris Apanasov and Serge Cantat who helped me to correct some of my illiteracy in the theory of Kleinian groups, and of Hee Oh for explaining to me the state of the art in the general orbital counting problem.

2 We borrow the terminology from [9]

I am also thankful to the organizers of the conference at Schiermonnikoog2014 for giving me an opportunity to talk about the beautiful subject of Apollonian circle packings.

## 2. Hyperbolic space

Let R n, 1 denote the Minkowski space defined by the quadratic form (the fundamental quadratic form )

$$q = - t _ { 0 } ^ { 2 } + \sum _ { i = 1 } ^ { n } t _ { i } ^ { 2 }$$

of signature ( n, 1). 3 We denote by ( v, w ) the value of the associated symmetric form on two vectors v, w ∈ R n, 1 , so that q v ( ) = ( v, v ). Let

$$\mathbb { H } ^ { n } = \{ v \in \mathbb { R } ^ { n, 1 } \colon q ( v ) < 0 \} / \mathbb { R } ^ { * } \subset \mathbb { P } ^ { n + 1 } ( \mathbb { R } )$$

be the n -dimensional hyperbolic (or Lobachevsky ) space, a Riemannian space of constant negative curvature. Its points are lines [ v ] = R · v spanned by a vector v ∈ R n, 1 with ( v, v ) &lt; 0. The set of real points of the fundamental quadric Q : q = 0 is called the absolute . The union of H n and the absolute is the closure H n of H n in P n ( R ). In affine coordinates y i = x /x i 0 , one represents H n by the interior of the n -dimensional ball

$$\cdot \stackrel { \dots } { \Phi } B _ { n } = \{ ( y _ { 1 }, \dots, y _ { n } ) \colon \sum _ { \substack { i = 1 \\ \Phi } } ^ { n } y _ { i } ^ { 2 } < 1 \}. \\ \cdot \stackrel { \dots } { \Phi } \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Psi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi\Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phis \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \ Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi \Phi$$

This is called the Klein model of H n . In this model the Riemannian metric of constant curvature -1 is given by

$$d s ^ { 2 } = \frac { ( 1 - | y | ^ { 2 } ) \sum _ { i = 1 } ^ { n } \bar { d y } _ { i } ^ { 2 } + ( \sum _ { i = 1 } ^ { n } y _ { i } d y _ { i } ) ^ { 2 } } { ( 1 - | y | ^ { 2 } ) ^ { 2 } },$$

where | y | 2 = ∑ n i =1 y 2 i . One can also represent any point in H n by a vector v = ( x , . . . , x 0 n ) ∈ R n, 1 with ( v, v ) = -1 , x 0 &gt; 0. Thus, we obtain a model of H n as one sheet of a two-sheeted hyperboloid Y n in R n, 1

$$- t _ { 0 } ^ { 2 } + \sum _ { i = 1 } ^ { n + 1 } t _ { i } ^ { 2 } = - 1, \ t _ { 0 } > 0. \\......... \.. \.... \ e \prod _ { i = 1 } ^ { n + 1 } \ e \prod _ { i = 1 } ^ { n + 1 }.$$

This model is called the vector model of H n . The hyperbolic distance d v, w ( ) in this model between two points in H n satisfies

$$\cosh d ( v, w ) = - ( v, w ).$$

3 Note that we use the signature ( n, 1) instead of (1 , n ) to agree with the standard notations in hyperbolic geometry.

Here and later we use that, for any two vectors v, w in the vector model of H n , we have ( v, w ) &lt; 0 (see [34], Theorem 3.1.1). The maps

↦

$$B ^ { n } \rightarrow Y ^ { n }, \ ( y _ { 1 }, \dots, y _ { n } ) \mapsto ( \frac { 1 } { \sqrt { 1 - | y | ^ { 2 } } }, \frac { y _ { 1 } } { \sqrt { 1 - | y | ^ { 2 } } }, \dots, \frac { y _ { n } } { \sqrt { 1 - | y | ^ { 2 } } } )$$

and

↦

$$Y ^ { n } \to B _ { n }, \ ( t _ { 0 }, \dots, t _ { n + 1 } ) \mapsto ( \frac { t _ { 1 } } { t _ { 0 } }, \dots, \frac { t _ { n } } { t _ { 0 } } )$$

are inverse to each other.

The induced metric of R n, 1 ds 2 = -dt 2 0 + ∑ n i =1 dt 2 i , after the scaling transformation y i = t /q i ( ) t 1 2 / , defines a Riemannian metric on H n of constant curvature -1

$$( 2. 2 ) \quad \ \ d s ^ { 2 } = \frac { 1 } { ( 1 - | y | ^ { 2 } ) ^ { 2 } } \left ( ( 1 - | y | ^ { 2 } ) \sum _ { i = 1 } ^ { n } d y _ { i } ^ { 2 } + ( \sum _ { i = 1 } ^ { n } y _ { i } d y _ { i } ) ^ { 2 } \right ),$$

where | y | 2 = y 2 1 + · · · + y 2 n . Consider the section of B n by the the hyperplane y n = 0. It is isomorphic to the n -1-dimensional ball

$$( 2. 3 ) \quad \ \ ^ { \prime } B _ { n - 1 } = \{ ( y _ { 1 }, \dots, y _ { n - 1 } ) \in \mathbb { R } ^ { n - 1 } \colon \sum _ { i = 1 } ^ { n - 1 } y _ { i } ^ { 2 } < 1 \}.$$

The restriction of the metric (2.2) to ′ B n -1 coincides with the metric of H n -1 . Thus we may view ′ B n -1 as a ball model of a geodesic hypersurface of B n = H n .

↦

Let p = [1 0 , , . . . , 0 , -1] ∈ Q ( R ) be the southern pole of the absolute. The projection P n glyph[axisshort]glyph[axisshort]glyph[arrowaxisright] P n -1 from this point is given by the formula [ t 0 , . . . , t n ] → [ t 1 , . . . , t n -1 , t 0 + t n ]. Since the hyperplanes t 0 + t n = 0 , t 0 = 0 do not intersect Q ( R ), Q ( R ) \ { p } maps bijectively onto R n -1 with coordinates u i = y i y n +1 . The rational map defined by the formula (2.4)

↦

Φ : P n -1 → Q , [ x , . . . , x 0 n -1 ] → [ x 2 0 + | x | 2 , 2 x x , . . . , 0 1 2 x x 0 n -1 , x 2 0 -| x | 2 ] , where | x | 2 = x 2 1 + · · · + x 2 n -1 , blows down the hyperplane x 0 = 0 to the point p and equals the inverse of the restriction of the projection Q ( R ) \{ } → p P n -1 ( R ). In affine coordinates u i = x /x i 0 in P n -1 \{ x 0 = 0 } and affine coordinates ( y , . . . , y 1 n ) in Q \{ t 0 = 0 } the map is defined by the formula (2.5)

↦

$$\lambda _ { \lambda } ^ { \dots \cdot \cdot \cdot / } ( u _ { 1 }, \dots, u _ { n - 1 } ) \mapsto ( y _ { 1 }, \dots, y _ { n } ) = ( \frac { 2 u _ { 1 } } { 1 + | u | ^ { 2 } }, \dots, \frac { 2 u _ { n - 1 } } { 1 + | u | ^ { 2 } }, \frac { 1 - | u | ^ { 2 } } { 1 + | u | ^ { 2 } } ),$$

where | u | 2 = u 2 1 + · · · + u 2 n -1 . Let Q + be the part of Q ( R ) where y n &gt; 0 (the northern hemisphere). The preimage of Q + under the map Φ is the subset of R n -1 defined by the inequality | u | 2 &lt; 1. It is a unit ball

B n -1 of dimension n -1. Consider the orthogonal projection α of Q + to the ball B n -1 from (2.3). The composition

↦

$$( 2. 6 ) \ \alpha \circ \Phi \, \colon B _ { n - 1 } \to ^ { \prime } B _ { n - 1 }, \ ( u _ { 1 }, \dots, u _ { n - 1 } ) \mapsto ( \frac { 2 u _ { 1 } } { 1 + | u | ^ { 2 } }, \dots, \frac { 2 u _ { n - 1 } } { 1 + | u | ^ { 2 } } )$$

with the inverse

↦

$$\overset { \prime } { B } _ { n - 1 } \to B _ { n - 1 }, \ ( y _ { 1 }, \dots, y _ { n - 1 } ) \mapsto ( \frac { y _ { 1 } } { 1 + \sqrt { 1 - | y | ^ { 2 } } }, \dots, \frac { y _ { n - 1 } } { 1 - \sqrt { 1 + | y | ^ { 2 } } } )$$

is a bijection which we can use to transport the metric on ′ B n -1 to a metric on B n -1 with curvature equal to -1. In coordinates u , . . . , u 1 n -1 , it is given by the formula

$$d s ^ { 2 } = 4 ( 1 - | u | ^ { 2 } ) ^ { - 2 } \sum _ { i = 1 } ^ { n - 1 } d u _ { i } ^ { 2 }.$$

This is a metric on B n -1 of constant curvature -1, called the Poincar´ e metric . The corresponding model of the hyperbolic space H n -1 is called a conformal or Poincar´ e model of H n -1 .

Yet, another model of H n is realized in the upper half-space

$$\mathcal { H } ^ { n } = \{ ( u _ { 1 }, \dots, u _ { n } ) \in \mathbb { R } ^ { n }, u _ { n } > 0 \}.$$

The isomorphism B n →H n is defined by

$$( u _ { 1 }, \dots, u _ { n } ) = \frac { 1 } { \rho ^ { 2 } } ( 2 y _ { 1 }, \dots, 2 y _ { n - 1 }, 2 ( y _ { n } + 1 ) - \rho ^ { 2 } ),$$

where ρ 2 = y 2 1 + · · · + y 2 n -1 +( y n +1) . The metric on 2 H n is given by

$$d s ^ { 2 } = \frac { 1 } { y _ { n } ^ { 2 } } \sum _ { i = 1 } ^ { n } d u _ { i } ^ { 2 }.$$

A k -plane Π in H n is the non-empty intersection of a k -plane in P n ( R ) with H n . We call it proper if its intersection with H n is nonempty, and hence, it is a geodesic submanifold of H n of dimension k isomorphic to H k .

A geodesic line corresponds to a 2-dimensional subspace U of signature (1 , 1). It has a basis ( f, g ) that consists of isotropic vectors with ( f, g ) = -1 (called a standard basis ). Thus any geodesic line glyph[lscript] intersects the absolute at two different points [ f ] and [ g ]. We can choose a parameterization γ : R → l of glyph[lscript] of the form γ t ( ) = 1 √ 2 ( e f t + e -t g ) so that γ ( -∞ ) = [ f ] and γ (+ ∞ ) = [ ]. g The hyperbolic distance d γ ( (0) , γ ( )) from the point t γ (0) = 1 √ 2 ( f + ) to the point g γ t ( ) satisfies d γ ( (0) , γ ( )) t = -( γ (0) , γ ( )) t = cosh t . Thus t = d γ ( (0) , γ ( )) t is the natural parameter of the geodesic line. The subgroup of SO( U ) 0 ∼ =

SU (1 , 1) 0 is isomorphic to R &gt; 0 and consists of transformation g t given in the standard basis by the matrix ( e t 0 0 e -t ) . It can be extended to a group of isometries of R n, 1 that acts identically on U ⊥ . Comparing with the natural parameterization of l , we see that each transformation g t moves a point x ∈ l to a point g t ( x ) such that d x, g ( t ( x )) = t . For this reason, it is called the hyperbolic translate .

A n -1-plane H is called a hyperplane . We write it in the form

$$H _ { \mathfrak { e } } = \{ v \in \mathbb { R } ^ { n, 1 } \, \sum \{ 0 \} \, \colon ( v, \mathfrak { e } ) = 0, ( v, v ) \leq 0 \} / \mathbb { R } ^ { * }.$$

Since the signature of q is equal to ( n, 1), H n ∩ H e = ∅ if and only if ( e , e ) ≥ 0. If ( e , e ) = 0, the intersection with H n +1 consists of one point [ ]. e We assume that ( e , e ) ≥ 0 and ( e , e ) = 1 if ( e , e ) &gt; 0. These properties define e uniquely, up to multiplication by a constant (equal to ± 1 if ( e , e ) = 1). A choice of one of the two rays of the line R e is called an orientation of the hyperplane. The complement of H e in H n consists of two connected components (half-spaces) defined by the sign of ( v, e ). We denote the closures of these components by H ± e , accordingly. Changing the orientation, interchanges the half-spaces.

Consider a hyperplane H e ′ in ′ B n -1 . It is given by an equation

$$a _ { 0 } t _ { 0 } + \sum _ { i = 1 } ^ { n - 1 } a _ { i } t _ { i } = 0$$

for some ( a , . . . , a 0 n -1 ) ∈ R n -1 1 , with norm 1.

Under the map (2.6), the pre-image of H e ′ is equal to the intersection S ( e ′ ) ∩ ¯ B n -1 , where

$$S ( \mathfrak { e } ^ { \prime } ) \coloneqq - a _ { 0 } \sum _ { i = 1 } ^ { n - 1 } u _ { i } ^ { 2 } + 2 \sum _ { i = 1 } ^ { n - 1 } a _ { i } u _ { i } - a _ { 0 } = 0.$$

glyph[negationslash]

If k = a 0 = 0, we can rewrite this equation in the form

$$k ^ { 2 } \sum _ { i = 1 } ^ { n - 1 } ( u _ { i } - a _ { i } / k ) ^ { 2 } = \sum _ { i = 1 } ^ { n - 1 } a _ { i } ^ { 2 } - a _ { 0 } ^ { 2 } = 1.$$

This is a sphere of radius r = 1 k with the center c = ( a /k, . . . , a 1 n -1 /k ). It is immediately checked that S ( e ′ ) intersects the boundary of B n -1 orthogonally (i.e the radius-vectors are perpendicular at each intersection point in the Euclidean inner product). The intersection of S ( e ′ ) with B n -1 is a hyperbolic codimension 1 subspace of the Poincar´ model of e H n -1 .

Assume k = 0, the equation of S ( e ′ ) becomes ∑ n -1 i =1 a u i i = 0 . The hyperplane S ( e ′ ) should be viewed as a sphere of radius r = ∞ with the center at the point at infinity corresponding to p in the compactification

R n ∪ {∞} ∼ = Q ( R ). For any k ≥ 0, the intersection S ( e ′ ) ∩ B n -1 is a geodesic hypersurface in H n -1 . It contains the center of the ball if and only if a 0 = 0.

Example 2.1. A one-dimensional hyperbolic space can be modeled by an open interval ( -1 1). , The Poincar´ e and the Klein metrics coincide and equal to dy / 2 (1 -y 2 ) 2 . It is isomorphic to the Euclidean onedimensional space in this case.

Assume n = 2. We have two models of a hyperbolic 2-dimensional space H 2 both realized in the unit disk B 2 . The the Klein model is defined by the metric

$$d s ^ { 2 } = \frac { y _ { 1 } d y _ { 2 } ^ { 2 } + y _ { 2 } d y _ { 1 } ^ { 2 } } { ( 1 - y _ { 1 } ^ { 2 } - y _ { 2 } ^ { 2 } ) ^ { 2 } }. \\. \quad. \quad. \quad.$$

The Poincar´ e model is defined by the metric

$$d s ^ { 2 } = \frac { d y _ { 1 } ^ { 2 } + d y _ { 2 } ^ { 2 } } { ( 1 - y _ { 1 } ^ { 2 } - y _ { 2 } ^ { 2 } ) ^ { 2 } }.$$

In the Klein model, a geodesic line in H 2 is a line joining two distinct points on the absolute. In the Poincar´ model, a geodesic line is a circle e arc intersecting the absolute orthogonally at two distinct points or a line passing through the center of the disk.

The group Iso( H n ) of isometries of H n is the index 2 subgroup O( n, 1) ′ of the orthogonal group of R n. 1 that preserves each sheet of the hyperboloid representing the vector model of H n . When n is even, we have O( n, 1) ′ = SO( n, 1). It is also isomorphic to the group PO( n, 1) of projective transformations of P n that leaves invariant the quadric Q = ∂ H n . It consists of two connected components, the connected component of the identity PO( n, 1) 0 is the group of isometries that preserve the orientation. It is isomorphic to the subgroup of SO( n, 1) 0 of elements of spinor norm 1.

Via the rational map (2.4), the group Iso( H n ) becomes isomorphic to the subgroup of the real Cremona group Cr R ( n -1) of P R n -1 . It is generated by affine orthogonal transformations of P R n -1 and the inversion birational quadratic transformation

↦

$$[ x _ { 0 }, \dots, x _ { n - 1 } ] \mapsto [ \sum _ { i = 1 } ^ { n - 1 } x _ { i } ^ { 2 }, x _ { 0 } x _ { 1 }, \dots, x _ { 0 } x _ { n - 1 } ].$$

This group is known classically as the Inversive group in dimension n -1. We can also identify Q with a one-point compactification ˆ E n -1 = R n -1 ∪ {∞} of the Euclidean space of dimension n -1. The point ∞ corresponds to the point on Q from which we project to P n -1 . In this

↦

↦

model the Inversive group is known as the M¨bius group M¨b( o o n -1). If we identify, in the usual way, the 2-dimensional sphere with the complex projective line P 1 ( C ), we obtain that the group M¨b(2) is isomorphic o to PSL ( 2 C ) and consists of M¨bius transformations o z → az + b cz + d . The inversion transformation (2.7) is conjugate to the transformation z → -1 /z .

It is known that the orthogonal group O( n, 1) is generated by reflections

↦

$$s _ { \alpha } \colon v \mapsto v - 2 \frac { ( v, \alpha ) } { ( \alpha, \alpha ) } \alpha,$$

glyph[negationslash]

where ( α, α ) = 0. If ( α, α ) &lt; 0, then s α has only one fixed point in H n , the point [ α ] ∈ H n . If ( α, α ) &gt; 0, then the fixed locus of s α is a hyperplane in H n . The pre-image of its intersection with the absolute ∂ H n is the sphere S α/ α, α ( ( ) 1 2 / ). The M¨ obius transformation corresponding to s α is the inversion transformation that fixes the sphere.

An isometry γ of a hyperbolic space H n is of the following three possible types:

- · hyperbolic if γ has no fixed points in H n but has two fixed points on the boundary;
- · parabolic if γ has no fixed points in H n but has one fixed point on the boundary;
- · elliptic if γ has a fixed point in H n .

An element is elliptic if and only if it is of finite order. In the ˆ E n -1 model of the absolute, a parabolic element corresponds to a translation in the affine orthogonal group AO( n -1).

A discrete subgroup Γ of Iso( H n ) is called a Kleinian group . It acts discontinuously on H n and the extension of this action to the absolute acts discontinuously on the set Ω(Γ) = ∂ H n \ Λ(Γ), where Λ(Γ) is the set of limit points of Γ, i.e. points on the boundary that belongs to the closure of the Γ-orbit of a point x in H n (does not matter which one). The set Ω(Γ) is called the discontinuity set of Γ. The set H n ∪ Ω(Γ) is the largest open subset of H n on which Γ acts discontinuously.

Following [34], § 12.2, we call a Kleinian group to be of the first kind (resp. of the second type ) if its discontinuity set is empty (resp. not empty).

A Kleinian group Γ is called elementary if it has a finite orbit in H n . It is characterized by the property that it contains a free abelian subgroup of finite index.

The limit set Λ(Γ) of an elementary Kleinian group is either empty or consists of at most two points, the fixed points of a parabolic or a

hyperbolic element in Γ. Otherwise Λ(Γ) is a perfect set, i.e. has no isolated points. A fixed point of a hyperbolic or a parabolic transformation is a limit point. The limit set of a non-elementary Kleinian group is equal to the closure of the set of fixed points of parabolic or hyperbolic elements [34], Theorem 12.2.2.

A Kleinian group in dimension 2 is called a Fuchsian group . Usually one assumes that Γ is contained in the connected component of Iso(2), i.e. it preserves the orientation. The discontinuity set Ω(Γ) is either empty, and then Γ is called of the first kind, or consists of the union of open intervals, then γ is of the second kind. The quotient X = H 2 ∪ Ω(Γ) Γ is a Riemann surface. / If Γ is of the first kind, one can compactify X by a finite set of orbits of cuspidal limit points. A Fuchsian group with H 2 / Γ of finite volume is of the first kind. If Γ is of the second kind, then X is a Riemann surface with possibly a non-empty boundary.

## 3. Geometrically finite Kleinian groups

A non-empty subset A of a metric space is called convex if any geodesic connecting two of its points is contained in A . Its closure is called a closed convex subset . A closed convex subset P of H n is equal to the intersection of open half-spaces H ,i e i ∈ I . It is assumed that none of the half-spaces H e i contains the intersection of other halfspaces. A maximal non-empty closed convex subset of the boundary of P is called a side . Each side is equal to the intersection of P with one of the bounding hyperplanes H e i . The boundary of P is equal to the union of sides and two different sides intersect only along their boundaries.

A convex polyhedron in H n is a convex closed subset P with finitely many sides. It is bounded by a finite set of hyperplanes H i .

We choose the vectors e i as in the previous section so that

$$P = \bigcap _ { i = 1 } ^ { N } H _ { \mathfrak { e } _ { i } } ^ { - }.$$

The dihedral angle φ H ,H ( e i e j ) between two proper bounding hyperplanes is defined by the formula

$$\cos \phi ( H _ { \mathfrak { e } _ { i } }, H _ { \mathfrak { e } _ { j } } ) \coloneqq - ( \mathfrak { e } _ { i }, \mathfrak { e } _ { j } ).$$

If | ( e i , e j ) | &gt; 1, the angle is not defined, we say that the hyperplanes are divergent . In this case the distance between H e and H e ′ can be found from the formula

$$& ( 3. 2 ) & & \cosh d ( H _ { \mathfrak { e } }, H _ { \mathfrak { e } ^ { \prime } } ) = | ( \mathfrak { e }, \mathfrak { e } ^ { \prime } ) |.$$

For any subset J of { 1 , . . . , N } of cardinality n -k , such the intersection of H ,i e i ∈ J is a k -plane Π in H n , the intersection P ∩ Π is a polyhedron in Π. It is called a k -face of P . A ( n -1)-face of P is a side of P .

The matrix

$$G ( P ) = ( g _ { i j } ), \ \ g _ { i j } = ( \mathfrak { e } _ { i }, \mathfrak { e } _ { j } ),$$

is called the Gram matrix of P . There is a natural bijection between the set of its k -dimensional proper faces and positive definite principal submatrices of G P ( ) of size n -k . The improper vertices of P correspond to positive semi-definite principal submatrices of size n .

Recall that a closed subset D of a metric space X is called a fundamental domain for a group Γ of isometries of X if

- (i) the interior D o of D is an open set;
- (ii) γ D ( o ) ∩ D o = ∅ , for any γ ∈ Γ \ { 1 } ;
- (iii) the set of subsets of the form γ D ( ) is locally finite ; 4
- (iv) X = ∪ γ ∈ Γ γ D ( );

A group Γ admits a fundamental domain if and only if it is a discrete subgroup of the group of isometries of X . For example, one can choose D to be a Dirichlet domain

$$( 3. 3 ) \quad D ( x _ { 0 } ) = \{ x \in X \, \colon d ( x, x _ { 0 } ) \leq d ( \gamma ( x ), x _ { 0 } ), \text{ for any } \gamma \in \Gamma \},$$

glyph[negationslash]

where x 0 is fixed point in X and d x, y ( ) denotes the distance between two points. Assume X = H n . For any γ ∈ Γ \ Γ x 0 , let H γ be the hyperplane of points x such that d x , x ( 0 ) = ρ x, γ ( ( x 0 )). Then D x ( 0 ) = ∩ γ ∈ \ Γ Γ x 0 H -γ and, for any γ ∈ Γ x 0 , S = γ D x ( ( 0 )) ∩ D x ( 0 ) ⊂ H γ is a side in ∂D x ( 0 ). Each side of D x ( 0 ) is obtained in this way for a unique γ .

A fundamental domain D for a Kleinian group Γ in H n is called polyhedral if its boundary ∂D = D \ D o is contained in the union of a locally finite set of hyperplanes H i and each side S of the boundary is equal to D ∩ γ D ( ) for a unique σ S ∈ Γ. A Dirichlet domain is a convex polyhedral fundamental domain.

A choice of a polyhedral fundamental domain allows one to find a presentation of Γ in terms of generators and relations. The set of generators is the set of elements γ S ∈ Γ, where S is a side of D . A relation γ S t ◦· · ·◦ γ S 1 = 1 between the generators corresponds to a cycle

$$D _ { 0 } = D, D _ { 1 } = \gamma _ { S _ { 1 } } ( D _ { 0 } ), \dots, D _ { t } = \gamma _ { S _ { t } } \circ \dots \circ \gamma _ { S _ { 1 } } ( D _ { 0 } ) = D _ { 0 }.$$

Among various equivalent definitions of a geometrically finite Kleinian group we choose the following one (see [2], Chapter 4, § 1): A Kleinian

4 This means that every point x ∈ H n is contained in a finite set of subsets γ D ( )

group Γ is called geometrically finite if it admits a polyhedral fundamental domain with finitely many sides. 5 It follows from above that such a group is finitely generated and finitely presented. The converse is true only in dimension n ≤ 2. In dimensions n ≤ 3, one can show that Γ is geometrically finite if and only if there exists (equivalently any) a Dirichlet fundamental domain with finitely many sides (loc. cit., Theorem 4.4). On the other hand, for n &gt; 3 there are examples of geometrically finite groups all whose Dirichlet domains have infinitely many sides (loc. cit. Theorem 4.5).

glyph[negationslash]

Example 3.1. A convex polyhedron P in H n is called a Coxeter polyhedron if each pair of its hyperplanes H ,H e i e j is either divergent or the dihedral angle between them is equal to π/m ij , where m ij ∈ Z or ∞ . The Coxeter diagram of a Coxeter polyhedron P is a labeled graph whose pairs of vertices corresponding to divergent hyperplanes are joined by a dotted edge, two vertices corresponding to hyperplanes forming the zero angle are joined by a solid edge or by putting ∞ over the edge, and two vertices corresponding to the hyperplanes with the dihedral angle π/m ij = 0 , π/ 2 are joined by an edge with the label m ij -2, the label is dropped if m ij = 3.

glyph[negationslash]

A Coxeter polyhedron is a fundamental polyhedron of the reflection Kleinian group Γ P generated by the reflections s e i . For each pair of reflections s e i , s e j with m ij = ∞ , the product s e i ◦ s e j is a rotation in angle 2 π/m ij around the codimension 2 face H e i ∩ H e j . It follows that ( s i ◦ s j ) m ij is one of the basic relations. In fact, these relations are basic relations. An abstract Coxeter group is a pair (Γ , S ) consisting of a group Γ and a set S of its generators s i such that ( s s i j ) m ij = 1 are the basic relations, where ( m ij ) is a matrix with diagonal elements equal to 1 and off-diagonal elements equal to integers ≥ 2 or the infinity symbol ∞ if the order of the product is infinite. A realization of a Coxeter group ( G,S ) in the form of Γ P , where P is a Coxeter polyhedron in a hyperbolic space H n is called a hyperbolic Coxeter group. It follows that the cardinality of the set S of generators must be equal to n +1. Compact Coxeter polytopes (resp. of finite volume) are described by Lanner (resp. quasi-Lanner) diagrams which can be found for example in [40], Chapter 5, Table 4.

Let Γ be any discrete group and Γ r be its subgroup generated by reflections. It is a Coxeter group with a fundamental Coxeter polyhedron P . Then Γ r is a normal subgroup of Γ and

$$\Gamma = \Gamma _ { r } \rtimes A ( P ),$$

5 Other equivalent definition is given in terms of the convex core of Γ, the minimal convex subset of H n that contains all geodesics connecting any two points in Λ(Γ).

where A P ( ) is the subgroup of Γ which leaves P invariant.

## 4. Sphere packings

Recall that a real n-sphere S in the projective space P R n over reals is the set of real points Q ( R ) of a quadric Q defined over R that contains a given nonsingular quadric Q 0 with Q 0 ( R ) = ∅ in a fixed hyperplane H ∞ . We fix the projective coordinates to assume that H ∞ : x 0 = 0 and Q 0 : x 0 = x 2 1 + · · · + x 2 n = 0. A 1-sphere is called a circle , it is a conic passing through the points [0 , 1 , ± i ] at infinity. We do not exclude the case n = 1 where H ∞ is the point ∞ = [0 1] and , Q 0 = ∅ . A 0-dimensional sphere is a set of two real points in P R 1 that may coincide.

Let

$$\mathbb { P } _ { \mathbb { R } } ^ { n } \dashrightarrow \mathbb { P } _ { \mathbb { R } } ^ { n + 1 }$$

be the rational map defined by the formula

↦

$$[ x _ { 0 }, \dots, x _ { n } ] & \mapsto [ 2 x _ { 0 } ^ { 2 }, - 2 x _ { 0 } x _ { 1 }, \dots, - 2 x _ { 0 } x _ { n }, \sum _ { i = 1 } ^ { n } x _ { i } ^ { 2 } ]. \\. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

Its image is the quadric Q given by the equation

$$q & = 2 t _ { 0 } t _ { n + 1 } + t _ { 1 } ^ { 2 } + \dots + t _ { n } ^ { 2 } = 0.$$

This map is given by a choice of a basis in the linear system of quadrics containing the quadric Q 0 . Adifferent choice of a basis leads to the map (2.4) whose image is a quadric with equation -t 2 0 + ∑ n +1 i =1 t 2 i = 0. Let e = ( a , . . . , a 0 n +1 ) ∈ P n +1 1 , with ( e , e ) ≥ 0. We assume that ( e , e ) = 1 if ( e , e ) &gt; 0. The pre-image of a hyperplane

$$H _ { \mathfrak { e } } \colon a _ { 0 } t _ { n + 1 } + a _ { n + 1 } t _ { 0 } + \sum _ { i = 1 } ^ { n + 1 } a _ { i } t _ { i } = 0$$

in P R n +1 is a quadric in P R n defined by the equation

$$\mu _ { 0 } \sum _ { i = 1 } ^ { n } x _ { i } ^ { 2 } - 2 a _ { n + 1 } x _ { 0 } ^ { 2 } - 2 x _ { 0 } \sum _ { i = 1 } ^ { n } a _ { i } x _ { i } = 0.$$

glyph[negationslash]

If a 0 = 0, we can rewrite this equation in the form

$$( 4. 2 ) \quad \ a _ { 0 } ^ { 2 } \sum _ { i = 1 } ^ { n } ( \frac { x _ { i } } { x _ { 0 } } - \frac { a _ { i } } { a _ { 0 } } ) ^ { 2 } = 2 a _ { 0 } a _ { n + 1 } + \sum _ { i = 1 } ^ { n } a _ { i } ^ { 2 } = 1.$$

So, we can identify its real points with a n -dimensional sphere S ( e ) in the Euclidean space P n ( R ) \ { x 0 = 0 } of radius square r 2 = 1 /a 2 0 and the center c = [ a 1 a 0 , . . . , a n a 0 ]. It is natural to call | a 0 | the curvature of the sphere. If a 0 = 0, we get the union of two hyperplanes H ∞ and S := ∑ n i =1 a x i i + a n +1 x 0 = 0. In this case we set r := ∞ . Note that if

glyph[negationslash]

we take e = ( a , . . . , a 0 n +1 ) in (4.2) with ( e , e ) = 0 , a 0 = 0, we obtain a quadratic cone with the vertex at c . Its set of real points consists only of the vertex.

One introduces the oriented curvature equal to a 0 . We agree that the positive curvature corresponds to the interior of the sphere S ( e ), i.e. an open ball B ( e ) of radius r . The negative curvature corresponds to the open exterior of the sphere, we also call it the ball corresponding to the oriented sphere S and continue to denote it B ( e ). It can be considered as a ball in the extended Euclidean space ˆ E n .

A sphere packing is an infinite set P = ( S i ) i ∈ I of oriented n -spheres such that any two of them are either disjoint or touch each other. We say that a sphere packing is strict if, additionally, no two open balls B i intersect. An example of a non-strict sphere packing is an infinite set of nested spheres. We assume also that the set P is locally finite in the sense that, for any t &gt; 0, there exists only finitely many spheres of curvature at most t in any fixed bounded region of the space. The condition that two spheres S i and S j are disjoint or touch each other is easily expressed in terms of linear algebra. We have the following.

Lemma 4.1. Let S v ( ) and S w ( ) be two n -dimensional spheres corresponding to hyperplanes H v and H w . Then their interiors do not intersect if and only if

$$( v, w ) \leq - 1.$$

The equality takes place if and only if the spheres are tangent to each other, and hence intersect at one real point.

Proof. Consider the pencil of spheres λS v ( ) + µS w ( ) = H λv ( + µw ) (here, for brevity of notation, we do not normalize the normal vector defining the hyperplane λH v ( )+ µH µ ( )). It is clear that B v ( ) ∩ B w ( ) = ∅ if and only if S v ( ) = S w ( ) = ∅ or the pencil contains a cone with its vertex equal to the tangency points of the spheres. A point of intersection of two spheres is defined by a vector x such that ( x, v ) = ( x, w ) = ( x, x ) = 0. Such x exists if and only if the plane 〈 v, w 〉 spanned by v and w is positive or semi-positive definite. This happens if and only if ( v, v )( w, w ) -( v, w ) 2 ≥ 0. The equality takes place if and only x ∈ 〈 v, w 〉 . This is equivalent to that one of the spheres in the pencil is a degenerate sphere, i.e. a cone with its vertex equal to the tangency point of the two spheres. Thus we proved that | ( v, w ) | ≥ 1. As we noted before, ( v, w ) must be negative. glyph[square]

Let Γ be a geometrically finite Kleinian group in H n +1 that is identified with a subgroup of the M¨bius group M¨b( o o n ). We say that a sphere packing is a sphere Γ -packing if it is invariant with respect to

the action of Γ on spheres and has only finitely many orbits, say N , on the set. We also say that a sphere Γ-packing is clustered if one can choose representatives S , . . . , S 1 N of the orbits such that P is equal to the union of the sets { γ S ( 1 ) , . . . , γ ( S N ) } , γ ∈ Γ. Each such set is called a cluster of the sphere packing. We will see later than Apollonian sphere packings and, their generalizations, Boyd-Maxwell sphere packings, are examples of clustered sphere Γ-packings.

Let Γ i be the stabilizer subgroup of some sphere S i in a sphere Γpacking. The sphere corresponds to a geodesic hyperplane H i in H n +1 , hence the ball B i is the Klein model of H n on which Γ i acts as a geometrically finite Kleinian group. Its limit points are on S i . We know that they belong to the closure of the set of fixed points of nonelliptic elements in Γ , hence they are limit point of the whole Γ. i

Let us introduce the following assumptions.

- (A1) for any sphere S i , Λ(Γ ) = i S . i
- (A2) ⋃ i ∈ I ¯ B i = ˆ E . n

Assumption (A1) implies that

$$\overline { \bigcup _ { i \in I } ^ { S _ { i } } } \subset \Lambda ( \Gamma ).$$

Assumption (A2) means that the packing is maximal , i.e. one cannot find a sphere whose interior is disjoint from the interiors of all other spheres in the packing.

glyph[negationslash]

If the packing is strict, assumptions (A1) and (A1) imply the equality in (4.3). Indeed, suppose x is a limit point of Γ that does not lie in the closure of the union of the spheres. By (A2), it belongs to the closure of some open ball B i . Suppose x belongs to the open ball B i . We mentioned before that the limit set is the closure of the set of fixed points of non-elliptic elements of Γ. Thus we can find some open neighborhood of x inside B i that contains a fixed point of a non-elliptic element γ in Γ. The stabilizer subgroup Γ i of B i in Γ is a Kleinian group acting on B i . It cannot contain a non-elliptic element with a fixed point in the ball. This shows that γ ∈ Γ i and hence γ B ( i ) glyph[subsetnoteql] B i . Thus the set { γ n ( B ,n i ) ∈ Z } consists of nested spheres. This is excluded by the definition of a strict sphere packing.

In the next two sections we give examples of sphere Γ-packings satisfying assumptions (A1) and (A2).

## 5. Boyd-Maxwell sphere packings

Let P be a Coxeter polytope in H n +1 with the Gram matrix G P ( ). Let e 1 , . . . , e n +2 be the normal vectors of its bounding hyperplanes. We

can take the set ( e 1 , . . . , e n +2 ) to be a basis of the space R n +1 1 , with the fundamental quadric (4.1).

Let ω j be a vector in R n, 1 uniquely determined by the condition

$$( 5. 1 ) & & ( \omega _ { j }, \mathfrak { e } _ { i } ) = \delta _ { i j }.$$

We have

$$\sum ^ { n } g ^ { i j } \mathfrak { e } _ { i },$$

$$\omega _ { j } = \sum _ { i = 1 } ^ { \ell \ell } g ^ { i j } \mathfrak { e } _ { i },$$

where

$$G ( P ) ^ { - 1 } = ( g ^ { i j } ) = ( ( \omega _ { i }, \omega _ { j } ) ).$$

We call ω i real if ( ω , ω i i ) = g ii &gt; 0 . In this case we can normalize it to set

$$\bar { \omega } _ { i } \coloneqq \omega _ { i } / \sqrt { g ^ { i i } }.$$

Let J = { j 1 &lt; · · · &lt; j r } be the subset of { 1 , . . . , n +2 } such that ¯ ω j is real if and only if j ∈ J . Consider the union

$$\mathfrak { P } ( P ) = \bigcup _ { k = 1 } ^ { r } O _ { \Gamma _ { P } } ( S ( \bar { \omega } _ { j _ { k } } ) ).$$

A Boyd-Maxwell sphere packing is a sphere packing of the form P ( P ), where P is a Coxeter polytope. By definition, it is clustered with clusters

$$( \gamma ( S ( \bar { \omega } _ { j _ { 1 } } ) ), \dots, \gamma ( S ( \bar { \omega } _ { j _ { r } } ) ) ) = ( S ( \gamma ( \bar { \omega } _ { j _ { 1 } } ) ), \dots, S ( \gamma ( \bar { \omega } _ { j _ { r } } ) ).$$

We call the cluster ( S ω ( j 1 ) , . . . , S ( ω j r )) the initial cluster . It is called non-degenerate if all ω i are real. This is equivalent to that all principal maximal minors of the matrix G P ( ) are negative.

glyph[negationslash]

Following G. Maxwell [26], we say that the Coxeter diagram is of level l if, after deleting any l of its vertices, we obtain a Coxeter diagram of Euclidean or of parabolic type describing a Coxeter polyhedron in an Euclidean or a spherical geometry. They can be characterized by the property that all m ij = ∞ (except in the case n = 1) and the symmetric matrix ( m ij ) is non-negative definite, definite in the Euclidean case, and having a one-dimensional radical in the parabolic case. Coxeter diagrams of level l = 1 correspond to Lanner and quasi-Lanner Coxeter groups. Coxeter diagrams of level 2 have been classified by G. Maxwell (with three graphs omitted, see [9]). They occur only in dimension n ≤ 10. For n ≥ 4 they are obtained from quasi-Lanner diagrams by adding one vertex.

For example, a Coxeter polyhedron in H 10 with Coxeter diagram of level 2

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-afa5b531-0e6b-4cd6-b0f2-dceddaf43db3.jpg)

defines a sphere packing with only one real ω i corresponding to the extreme vertex on the right.

The following theorem is proven, under a certain assumption, in loc.cit., Theorem 3.3. The assumption had been later removed in [27], Theorem 6.1. 6

Theorem 5.1. Let P be a Coxeter polyhedron in H n . Then P ( P ) is a sphere Γ P -packing if and only if the Coxeter diagram of P is of level ≤ 2 .

If P is a quasi-Lanner Coxeter polytope, the Kleinian group Γ P has a fundamental polyhedron of finite volume. It follows that Γ P is a Kleinian group of the first kind ([34], Theorem 12.2.13). Thus any point on the absolute is its limit point. The stabilizer Γ P,i of each S ω (¯ ) i is the reflection group of the Coxeter group defined by the Coxeter diagram of level 1, hence the limit set of Γ P,i is equal to S ω (¯ ). i Thus the assumption (A1) is satisfied.

Let P ( P ) be a non-degenerate Boyd-Maxwell sphere packing and let k = ( k , . . . , k 1 n +2 ) be the vector of the curvatures of the spheres S v ( i ) in its cluster ( S v ( 1 ) , . . . , S ( v n +2 )).

Theorem 5.2. Let G P ( ) be the Gram matrix of P . Then

$$\ ^ { t } { \mathbf k } \cdot G ( P ) \cdot { \mathbf k } = 0.$$

Proof. Let J n +2 be the matrix of the symmetric bilinear form defined by the fundamental quadratic form q . It satisfies J n +2 = J -1 n +2 . Let X be the matrix whose j th column is the vector of coordinates of the vector v j . Recall that the first coordinate of each vector v j is equal to the curvature of the sphere S v ( j ) . Let v i = γ ω (¯ ) for some i γ ∈ Γ . We P have

$$( v _ { i }, v _ { j } ) = ( \gamma ( \bar { \omega } _ { i } ), \gamma ( \bar { \omega } _ { j } ) ) = ( \bar { \omega } _ { i }, \bar { \omega } _ { i } ).$$

Thus, hence

$$^ { t } X ^ { - 1 } \cdot G ( P ) ^ { - 1 } \cdot X ^ { - 1 } = J _ { n + 2 }.$$

6 I am grateful to Hao Chen for the reference.

$$^ { t } X \cdot J _ { n + 2 } \cdot X = G ( P ) ^ { - 1 },$$

Taking the inverse, we obtain

$$X \cdot G ( P ) \cdot ^ { t } X = J _ { n + 2 } ^ { - 1 } = J _ { n + 2 }.$$

The first entry a 11 of the matrix in the right-hand side is equal to zero. Hence

$$\L ^ { t } { \mathbf k } \cdot G ( P ) \cdot { \mathbf k } = 0,$$

It is clear that we have the following possibilities for a cluster in P ( P ).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-47daa083-b814-43ce-896d-29860e97b782.jpg)

Assume n ≥ 2. Suppose P ( P ) is a non-degenerate Boyd-Maxwell sphere packing. Let ¯ ω , . . . , ω 1 ¯ n +2 be defined as (5.2). Then the hyperplanes H ω (¯ ) bound a convex polyhedron with no proper vertices on i the absolute. Let us denote it by P ⊥ . Since P ( P ) is a sphere packing, no two bounding hyperplanes of P ⊥ have a common point in the absolute. By Lemma 4.1, glyph[negationslash]

$$( 5. 5 ) & & ( \bar { \omega } _ { i }, \bar { \omega } _ { j } ) \leq - 1, \ \ i \neq j.$$

Thus any two bounding hyperplanes of P ⊥ either diverge or have the dihedral angle equal to 0. In particular, we obtain that P ⊥ is a Coxeter polyhedron and hence its reflection group Γ ⊥ P := Γ P ⊥ is a Kleinian group. The Coxeter diagram of P ⊥ has either solid or dotted edges. As a non-weighted graph, it is a complete graph K n +2 with n + 2 vertices. It is clear that it is of level l ≤ 2 if and only if n ≤ 3. In this case we can define the dual Boyd-Maxwell sphere packing

$$\mathfrak { P } ( P ) ^ { \perp } \coloneqq \mathfrak { P } ( P ^ { \perp } ).$$

glyph[negationslash]

Note that the property ( e i , ω ¯ ) = 0 for k i = k is interpreted in terms of spheres as the property that the spheres S ( e i ) and S ω (¯ ) are orthogonal j to each other (see, for example, [13], Proposition 7.3). This implies that the dual sphere packing can be defined by a cluster that consists of spheres orthogonal to all spheres of a cluster defining P ( P ) except one.

Let A be the matrix of a transformation from Γ P ⊥ in a basis ( ω , . . . , ω 1 n +1 ). Then

$$^ { t } A \cdot G ( P ) ^ { - 1 } \cdot A = G ( P ).$$

Taking the inverse, we obtain

$$A ^ { - 1 } \cdot G ( P ) \cdot ^ { t } A ^ { - 1 } = G ( P ).$$

This shows that under the homomorphism A → t A -1 , the group Γ ⊥ P becomes isomorphic to a subgroup of the orthogonal group of the quadratic form defined by the matrix G P ( ).

The group Γ P is generated by the reflections s e i that act on the dual basis ( ω , . . . , ω 1 n +2 ) by the formula

$$& ( 5. 7 ) \quad s _ { \mathfrak { e } _ { i } } ( \omega _ { j } ) = \omega _ { j } - 2 ( \mathfrak { e } _ { i }, \omega _ { j } ) \mathfrak { e } _ { i } = \omega _ { j } - 2 \delta _ { i j } \mathfrak { e } _ { i } = \omega _ { j } - 2 \delta _ { i j } \sum _ { k = 1 } ^ { n + 2 } g _ { k i } \omega _ { k }, \\ & \cdot \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad. \quad.$$

where δ ij is the Kronecker symbol. This gives explicitly the action of Γ P on the clusters of the sphere packing and also on the set of their curvature vectors.

The group Γ P ⊥ is generated by the reflections s ω i = s ¯ ω i that act on the basis ( e 1 , . . . , e n +2 ) by the formula

$$s _ { \omega _ { i } } ( \mathfrak { c } _ { j } ) = \mathfrak { c } _ { j } - 2 ( \bar { \omega } _ { i }, \mathfrak { c } _ { j } ) \bar { \omega } _ { i } = \mathfrak { c } _ { j } - \frac { 2 \delta _ { i j } } { g ^ { i i } } \omega _ { i } = \mathfrak { c } _ { j } - \frac { 2 \delta _ { i j } } { g ^ { i i } } \sum _ { k = 1 } ^ { n + 2 } g ^ { k i } e _ { k }.$$

Example 5.3. This is the most notorious and widely discussed beautiful example of a Boyd-Maxwell sphere packing (see, for example, [18], [19], [20],[17],[36]).

We start with a set of pairwise touching non-degenerate n -spheres S v ( 1 ) , . . . , S ( v n +2 ) with positive radii that touch each other. By Lemma 4.1, the Gram matrix G = (( v , v i j )) has 1 at the diagonal and -1 off the diagonal. Thus,

$$G = c i r ( 1, - 1, \dots, - 1 ),$$

where cir( α , . . . , α 1 n +2 ) denotes the circulant matrix with the first row ( α , . . . , α 1 n +2 ). It is easy to compute the inverse of G to obtain

$$G ^ { - 1 } = \frac { 1 } { \varepsilon _ { n } } c i r ( n - 1, - 1, \dots, - 1 ).$$

$$^ {. 1 } = \frac { 1 } { 2 n } \text{cir} ( n - 1, - 1, \dots, - 1 ).$$

Assume n ≥ 2. Then

$$\frac { 2 n } { n - 1 } G ^ { - 1 } = \text{cir} ( 1, \frac { 1 } { 1 - n }, \dots, \frac { 1 } { 1 - n } ). \quad \ _ { - } \.$$

It is equal to the Gram matrix of a Coxeter polytope P if and only if n = 2 3. , The group Γ P is called the Apollonian group and will be denoted by Ap . It is a Kleinian group only for n n = 2 3. We call , P an

Apollonian polyhedron . Its Coxeter diagram is equal to the complete graph K 4 with solid edges if n = 2 and a complete graph K 5 with simple edges if n = 3. In each case the Coxeter diagram is of level 2. Thus, P defines a Boyd-Maxwell sphere packing in dimensions n = 2 or 3. In the case n = 2, we see that all circles in the cluster defined by the vectors ¯ ω , . . . , ω 1 ¯ n +2 mutually touch each other. So, all clusters in this packing define a set of n +2 mutually touching circles. Any two circles in the packing are either disjoint or touch each other.

We have G P ( ) -1 = n -1 2 n G , thus P ⊥ has the Gram matrix equal to G . It follows that P ⊥ is always a Coxeter polyhedron with Coxeter diagram equal to a complete graph K n +2 with solid edges. It defines a sphere packing only if n ≤ 3. Clearly, it is non-degenerate, and hence by Theorem 5.1, it is maximal. If n = 2 3, it is , the dual to the Apollonian packing. Its clusters consist of n +2 mutually tangent n -spheres. If n = 2, the dual of an Apollonian circle packing is an Apollonian circle packing. Its clusters are orthogonal to clusters of the original packing as shown on Figure 1 below. The group Ap ⊥ n is called the dual Apollonian group .

Figure 1. Dual Apollonian circle clusters (from [26])

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-da8e143d-6317-4468-b358-b839741d603c.jpg)

Let H ¯ ω 1 , . . . , H ¯ ω n +2 be the bounding hyperplanes of an Apollonian Coxeter polyhedron and s e i be the corresponding reflections generating the Apollonian group Ap n . Let S ( e 1 ) , . . . , S ( e n +2 ) be the cluster of the dual Apollonian sphere packing. The formula (5.7) specializes to give

s e i (¯ ) = ¯ , and ω j ω j glyph[negationslash]

$$s _ { \mathfrak { e } _ { i } } ( \bar { \omega } _ { i } ) = - \bar { \omega } _ { i } + \frac { 2 } { n - 1 } \sum _ { j \neq i } \bar { w } _ { j }.$$

For example, s e 1 is represented in the basis (¯ ω , . . . , ω 1 ¯ n +2 ) by the matrix

$$\begin{array} { c c c c } \bullet & & \bullet & \bullet _ { \nu } \\ & & A _ { 1 } = \begin{pmatrix} - 1 & 0 & 0 & \dots & 0 \\ \frac { 2 } { n - 1 } & 1 & 0 & \dots & 0 \\ \vdots & \vdots & \vdots & \vdots \\ \frac { 2 } { n - 1 } & 0 & 0 & \dots & 1 \end{pmatrix}. \\ \bullet & & \bullet & \bullet & \bullet & \bullet & \end{array}$$

So, the Apollonian group is generated by n +2 matrices A , . . . , A 1 n +2 of this sort.

glyph[negationslash]

Formula (5.8) shows that the dual Apollonian group Ap ⊥ n acts on the Apollonian packing by acting on spheres S ( e i ) via the reflections s ¯ ω i . We have s ¯ ω i ( e j ) = e j if j = i and glyph[negationslash]

$$s _ { \bar { \omega } _ { i } } ( \mathfrak { e } _ { i } ) = - \mathfrak { e } _ { i } + 2 \sum _ { k = 1, k \neq i } ^ { n + 2 } \mathfrak { e } _ { j }$$

For example, s ¯ ω 1 is represented by the matrix

$$B _ { 1 } = \begin{pmatrix} - 1 & 2 & 0 & \dots & 2 \\ 0 & 1 & 0 & \dots & 0 \\ \vdots & \vdots & \vdots & \vdots \\ \dot { 0 } & \dot { 0 } & \dot { 0 } & \dots & 1 \end{pmatrix}.$$

The equation from Theorem 5.2 becomes

$$( 5. 1 0 ) \quad \ \ n ( k _ { 1 } ^ { 2 } + \dots + k _ { n + 2 } ^ { 2 } ) - ( k _ { 1 } + \dots + k _ { n + 2 } ) ^ { 2 } = 0.$$

It is known as Descartes's equation or Soddy's equation 7

glyph[negationslash]

Using (5.6), we see that Ap ⊥ n acts on the solutions of the equation (5.4) by the contragradient representations. It is generated by the matrices t A -1 = t A , . . . , 1 t A -1 n +2 = t A n +2 . However, if n = 2, it does not leave the Apollonian packing invariant.

Assume n &gt; 1. Consider a solution ( k , . . . , k 1 n +2 ) of the Descartes' equation which we rewrite in the form

$$\sum _ { i = 1 } ^ { n + 2 } k _ { i } ^ { 2 } - \frac { 2 } { n - 1 } \sum _ { 1 \leq i \leq j \leq n + 2 } k _ { i } k _ { j } = 0.$$

7 From a 'poem proof' of the theorem in the case n = 2 in 'Kiss Precise' by Frederick Soddy published in Nature, 1930:

Four circles to the kissing come. / The smaller are the bender. / The bend is just the inverse of / The distance from the center. / Though their intrigue left Euclid dumb / Theres now no need for rule of thumb. / Since zero bends a dead straight line / And concave bends have minus sign, / The sum of the squares of all four bends / Is half the square of their sum.

Thus, k n +2 satisfies the quadratic equation

$$t ^ { 2 } - \frac { 2 t } { n - 1 } \sum _ { i = 1 } ^ { n + 1 } k _ { i } - \frac { 2 } { n - 1 } \sum _ { 1 \leq i < k \leq n + 1 } ^ { n + 2 } k _ { i } = 0.$$

It expresses the well-know fact (the Apollonian Theorem ) that, given n +1 spheres touching each other, there are two more spheres that touch them. Thus, starting from an Apollonian cluster ( S , . . . , S 1 n +2 ), we get a new Apollonian cluster ( S , . . . , S 1 n +1 , S ′ n +2 ) such that the curvatures of S n +2 and S ′ n +2 are solution of the above quadratic equation. The curvatures of S ′ n +2 is equal to

$$k _ { n + 2 } ^ { \prime } = - k _ { n + 2 } + \frac { 2 } { n - 1 } \sum _ { i = 1 } ^ { n + 1 } k _ { i }.$$

So the new cluster ( S , . . . , S 1 n +1 , S ′ n +2 ) is the cluster obtained from the original one by applying the transformation s e n +2 .

Figure 2. Apollonina circle cluster with curvature vector ( -10 18 23 27) , , ,

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-01b6fb98-c661-4167-a309-a66b7f848d0e.jpg)

If we start with a cluster as in the picture, then all circles will be inclosed in a unique circle of largest radius, so our circle packing will look as in following Figure 3.

Another possibility is when we start with a cluster of three circles touching a line. Then, the circle packing will look like in the following Figure 4.

The following example is taken from [6].

Example 5.4. We take P with the Gram matrix

$$G ( P ) = \text{cir} ( 1, - 1, 0, - 1 ).$$

It is a Coxeter polyhedron with the Coxeter diagram

Figure 3. Apollonina circle packing

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-78e00927-9be4-4ff2-8d83-6fa58b85fefc.jpg)

Figure 4. Apollonina circle packing in a band

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-d331fce2-da69-4c50-9349-17fd22e07313.jpg)

$$\begin{smallmatrix} & & \infty & & \infty \\ \dddot { \longrightarrow } & & \infty & & \\ \end{smallmatrix}$$

The dual polyhedron P ⊥ has the Gram matrix G P ( ) -1

$$G ( P ) ^ { - 1 } = \frac { 1 } { 2 } \text{cir} ( 1, - 1, - 2, - 1 ).$$

Example 5.5. Let us consider the exceptional case of Apollonian packings when n = 1 (this case was briefly mentioned in [17]). Here 0dimensional spheres are sets of two or one points in P 1 ( R ) that we can identify with the boundary of the unit disk model of the hyperbolic plane H 2 . We start with three 0-dimensional spheres touching each other. This means that each pair has a common point. Let p 12 , p 13 , p 23 be the common points. We use the Poincar´ model of e H 2 . Then the corresponding hyperplanes H ω ( 1 ) , H ( ω 1 ) , H ( ω 1 ) are three geodesics joining pairs of these points. They form a ideal triangle with vertices on the boundary.

Figure 5. Ideal triangle

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-b3d7ad0f-d39f-4342-997a-65f50ef1d234.jpg)

The Apollonian group is not defined in this case. However, the dual Apollonian group Ap ⊥ 1 is defined. It is a Coxeter group with Coxeter diagram equal to the complete graph K 3 with solid edges. It is a quasiLanner Coxeter group generated by the reflections in sides of the ideal triangle.

The fundamental polyhedron of Ap ⊥ 1 is equal to the closure of the interior of the triangle. Its Gram matrix is equal to cir(1 , -1 , -1). Since it is of finite volume, the group is a Fuchsian group of the first kind with 3 cusps. The limit set of Ap ⊥ 1 is the whole circle. However, the packing is not strict since the interiors of our 0-dimensional spheres are one of the two intervals in P 1 ( R ) with end points equal to the sphere. Obviously we have a nested family of intervals. The set of 0-dimensional spheres in the sphere Ap ⊥ 1 -packing is a proper subset of Λ(Ap ⊥ 1 ). It consists of cusp points and is countable.

The Descartes equation becomes

$$& ( 5. 1 2 ) & & k _ { 1 } k _ { 2 } + k _ { 1 } k _ { 3 } + k _ { 2 } k _ { 3 } = 0.$$

## 6. Integral packings

Let K ⊂ R be a totally real field of algebraic numbers of degree d = [ K : Q ] and o K be its ring of integers. A free o K -submodule L of rank n +2 of R n +1 1 , is called a o K -lattice. Let ( e 1 , . . . , e n +2 ) be a basis of L . Assume that the entries of its Gram matrix G = (( e i , e j )) = ( g ij ) belong to o K Let σ i : K ↪ → R , i = 1 , . . . , d -1 , be the set of nonidentical embeddings of K into R . We assume additionally that the matrices G σ i = ( σ i ( g ij )) are positive definite. Let

$$f = \sum _ { 1 \leq i, j \leq n + 2 } g _ { i j } t _ { i } t _ { j }$$

be the quadratic form defined by the matrix G and O( f, o K ) be the subgroup of GL n +2 ( o K ) of transformations that leave f invariant. The group O( f, o K ) is a discrete subgroup of O( R n +1 1 , ) with a fundamental polyhedron of finite volume, compact if f does not represent zero. By passing to the projective orthogonal group, we will view such groups as Kleinian subgroups of Iso( H n +1 ). A subgroup of Iso( H n +1 ) which is commensurable with a subgroup of the form PO( f, o K ) is called arithmetic . Two groups PO( f , ′ o ′ K ) and PO( f, o K ) are commensurable if and only if K = K ′ and f is equivalent to λf ′ for some positive λ ∈ K . In particular, any subgroup of finite index in PO( f, o K ) is an arithmetic group.

We will be interested in the special case of an integral lattice where K = Q . A Kleinian group Γ in H n +1 is called integral if it is commensurable with a subgroup of PO( f, Z ) for some integral quadratic form f of signature ( n + 1 1). , We will be dealing mostly with nonarithmetic geometrically finite Kleinian groups. However, it follows from [5], Proposition 1, that such a group is always Zariski dense in PSO( n +1 1) or PO( , n +1 1) if it acts irreducibly in , H n +1 . In terminology due to P. Sarnak, Γ is a thin group .

We will use the usual terminology of the theory of integral quadratic forms, or more abstractly, integral quadratic lattices , free abelian groups L of finite rank equipped with a symmetric bilinear form with values in Z . By tensoring with R , we embed L in L R := L ⊗ R . We assume that the quadratic form of L extends to a quadratic form on L R of signature ( n +1 1), this will identify , L with an integral lattice in the previous sense. The set L ∨ of all vectors in L R such that ( x, v ) ∈ Z for all x ∈ L is called the dual lattice . It is a free abelian group spanned by the dual basis of a basis of L in L R . The Gram matrix of this basis is the inverse of the Gram matrix of the basis of L . The quotient group L /L ∨ is called the discriminant group of L . Its order is equal to the minus of the determinant of the Gram matrix of any basis of L . Note that L ∨ ⊂ L R is not, in general, an integral lattice. However, it becomes integral if we multiply its quadratic form by the exponent of the discriminant group of L . We denote by O( L ) the orthogonal group of L , the subgroup of GL( L ) that preserves the bilinear form. A reflection isometry of R N defined by a vector v ∈ L belongs to O( L ) if and only if 2( v,x ) ( v,v ) ∈ Z for all x ∈ L .

Remark 6.1 . Let Γ be an integral Kleinian group of isometries of H n of finite covolume, for example, the orthogonal group of some integral quadratic lattice L . It is obviously geometrically finite, however, for

n &gt; 3 it may contain finitely generated subgroups which are not geometrically finite. In fact, for any lattice L of rank ≥ 5 that contains a primitive sublattice defined by the quadratic form q = -x 2 0 + x 2 1 + x 2 2 + x 2 3 , the orthogonal group O( L ) contains finitely generated but not finitely presented subgroups (see [21]). 8

A sphere packing P = ( S i ) i ∈ I is called integral if the following conditions are satisfied.

- (i) The corresponding norm one vectors v i in R n +1 1 , span an integral lattice L ;
- (ii) There exists a positive integer λ such that λ v , v ( i j ) ∈ Z for all i, j ∈ I . The smallest such c is called the exponent .
- (iii) The curvatures of spheres S i are integers.

After multiplying the quadratic form of L , the restriction of the fundamental quadratic form to L , by the exponent we obtain an integral quadratic lattice int( L ). If Γ is a Kleinian group that leaves an integral sphere packing invariant, then it obviously leaves int( L ) invariant and becomes a subgroup of O(int( L )).

It follows from formula (5.7) that a Boyd-Maxwell sphere packing defined by a Coxeter polytope P with Gram matrix G P ( ) = ( g ij ) is integral if and only if the curvatures of the spheres from the initial cluster are integers and g ij ∈ Q . The quadratic lattice is defined by the matrix eG P ( ), where e is the smallest positive integer such that eg ij ∈ Z .

Example 6.2. Let P be an Apollonian polyhedron in H n +1 , n &gt; 1. It is defined by the matrix cir(1 , 1 1 -n , . . . , 1 1 -n ) whose inverse is equal to n -1 2 n cir(1 , -1 , . . . , -1). Recall that it is a Coxeter polyhedron only for n = 2 3. , The Apollonian packing is an integral Coxeter packing with exponent 1. The Apollonian group Ap n becomes a subgroup of the orthogonal group O( A p n ( )), where A p n ( ) is an integral quadratic lattice of rank n +2 with the quadratic form

$$q _ { n } = \sum _ { i = 1 } ^ { n + 2 } t _ { i } ^ { 2 } - 2 \sum _ { 1 \leq i < j \leq n + 2 } t _ { i } t _ { j }.$$

Recall that an integral quadratic lattice is called even if all its values are even integers, it is called odd otherwise. Obviously, the lattice A p n ( ) is odd. For any lattice L , we denote by L ev its largest even sublattice. Obviously, O( L ) ⊂ O( L ev ). Let us find an orthogonal decomposition of the lattice A p n ( ) ev .

8 I thank Boris Apanasov who informed me about this fact.

The quadratic form (6.1) is given in a basis ¯ ω , . . . , ω 1 ¯ n +2 of A p n ( ). Consider a new basis f , . . . , f 1 n +2 formed by the vectors ¯ ω , ω 1 ¯ +¯ 1 ω , 2 -¯ ω 1 -¯ ω , ω 3 ¯ 1 + ¯ ω 2 + ¯ ω 3 -¯ ω , ω 4 ¯ 4 -¯ ω , . . . , ω 5 ¯ n +1 -¯ ω n +2 . The Gram matrix of this basis is the following matrix.

$$\begin{pmatrix} 1 \, 0 \, 0 \, 0 \, 0 \, 0 \, 0 \, 0 \, \dots \, 0 \, 0 \\ 0 \, 0 \, 2 \, 0 \, 0 \, 0 \, 0 \, 0 \, \dots \, 0 \, 0 \\ 0 \, 2 \, 0 \, 0 \, 0 \, 0 \, 0 \, \dots \, 0 \, 0 \\ 0 \, 0 \, 0 \, 4 \, - 2 \, 0 \, 0 \, \dots \, 0 \, 0 \\ 0 \, 0 \, 0 \, - 2 \, 4 \, - 2 \, 0 \, \dots \, 0 \, 0 \\ \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \\ \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \vdots \end{pmatrix}.$$

This shows that, for any n &gt; 1,

$$\mathcal { A } p ( n ) & \cong \langle 1 \rangle \oplus U ( 2 ) \oplus \mathbf A _ { n - 1 } ( 2 ).$$

Here we use the standard notation from the theory of integral quadratic lattices: U is the integral hyperbolic plane, A n is the even positive definite root lattice defined by the Cartan matrix of type A n , 〈 k 〉 is the lattice of rank 1 spanned by a vector v with ( v, v ) = k , and M t ( ) denotes the quadratic lattice M with its integral quadratic form multiplied by a rational number t .

We have

$$( 6. 3 ) \quad \mathcal { A } p ( n ) ^ { e v } = \mathcal { A } _ { 1 } ( 2 ) \oplus \mathsf U ( 2 ) \oplus \mathsf A _ { n - 1 } ( 2 ) \cong ( \mathsf A _ { 1 } \oplus \mathsf U \oplus \mathsf A _ { n - 1 } ) ( 2 ).$$

The discriminant group of A p n ( ) (resp. A p n ( ) ev ) is isomorphic to ( Z / 2 Z ) n ⊕ Z / 2 n Z (resp. ( Z / 2 Z ) n ⊕ Z / 2 n Z ⊕ Z / 4 Z ).

Next we consider the dual Apollonian sphere packing. Its lattice is the dual lattice of A p n ( ) ∨ defined by the quadratic form from the left-hand side of (5.10) multiplied by n -1 2 n .

$$q _ { n } ^ { \perp } \coloneqq ( n - 1 ) \sum _ { i = 1 } ^ { n + 2 } t _ { i } ^ { 2 } - 2 \sum _ { 1 \leq i < j \leq n + 2 } t _ { i } t _ { j }.$$

Its exponent is equal to 2 n and coincides with the exponent of the discriminant group of A p n ( ). We have

$$\mathcal { A } p ( n ) ^ { \perp } \coloneqq \text{int} ( \mathcal { A } p ( n ) ^ { \vee } ) = \mathcal { A } p ( n ) ^ { \vee } ( 2 n ) \cong ( \langle 1 \rangle \oplus \mathsf U ( 1 / 2 ) \oplus \mathsf A ^ { \vee } _ { n - 1 } ( 1 / 2 ) ) ( 2 n )$$

$$= \langle 2 n \rangle \oplus \mathsf U ( n ) \oplus \mathsf A _ { n - 1 } ^ { \vee } ( n ).$$

The lattice A ∨ n -1 is generated by A n -1 and a vector

$$\beta = \frac { 1 } { n } ( \omega _ { 1 } + 2 \omega _ { 2 } + \cdots + ( n - 1 ) \omega _ { n - 1 } ),$$

where ω , . . . , ω 1 n -1 is the standard basis defined by simple roots. In the basis ( ω , . . . , ω 2 n -1 , β ), the matrix of A ∨ n -1 ( n ) is equal to

$$\begin{pmatrix} \dots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \\ ( 6. 4 ) & \left ( \begin{matrix} 2 n & - n & 0 & 0 & \dots & 0 & 0 & 0 \\ - n & 2 n & - n & 0 & \dots & 0 & 0 & 0 & 0 \\ 0 & - n & 2 & - n & \dots & 0 & 0 & 0 & 0 \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \\ 0 & 0 & 0 & 0 & \dots & - n & 2 n & - n & 0 \\ 0 & 0 & 0 & 0 & \dots & 0 & - n & 2 n & n \\ 0 & 0 & 0 & 0 & \dots & 0 & 0 & n & n - 1 \end{matrix} \right ) \\ \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots & \vdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots. \end{pmatrix}$$

Thus we obtain that Ap ⊥ n is even if n is odd, and, for n even

$$( 6. 5 ) \quad \ \ ( \mathcal { A } p ( n ) ^ { \perp } ) ^ { \text{ev} } \cong \langle 2 n \rangle \oplus \mathbf U ( n ) \oplus ( A _ { n - 1 } ^ { \vee } ( n ) ) ^ { \text{ev} },$$

where ( A ∨ n -1 ( n )) ev is defined by the matrix (6.4) in which the submatrix ( 2 n n n n -1 ) is replaced with ( 2 n 2 n 2 n 4 n -4 ).

Note that the Apollonian packing in dimension n = 1 exists but it is not a Boyd-Maxwell sphere packing. We have

$$\mathcal { A } p _ { 1 } \cong \langle 1 \rangle \oplus U ( 2 ), \quad \mathcal { A } p _ { 1 } ^ { \text{ev} } \cong ( \mathbf A _ { 1 } \oplus U ) ( 2 ).$$

The lattice A p ev 1 (1 / 2) = A U 1 is from Nikulin's list of integral even hyperbolic quadratic lattices of rank 3 and signature (2 , 1) such that the subgroup of the orthogonal group generated by reflections s , v ( v, v ) = 2 is of finite index (a 2-reflexive lattice) [31]. Each such lattice is a sublattice of A 1 ⊕ U . According to Nikulin's classification of even 2reflexive hyperbolic lattices of rank ≥ 5, none of the lattices A p n ( ) ev or A p n ( ) ev is 2-reflexive [29].

If n = 2, the lattices A p 2 and A p ⊥ 2 are isomorphic. According to Vinberg's list of 2-reflexive even hyperbolic lattice of rank 4 [41], the lattice A p ev 2 is not 2-reflexive.

Note that all examples from [6] are examples of integral Boyd-Maxwell sphere packings.

## 7. Orbital counting

Let A be a subset of the Euclidean space R n . Recall that the Hausdorff dimension δ A ( ) of A is defined to be the infimum for all s ≥ 0 for which

$$\mu _ { s } ( A ) = \inf _ { A \subset \cup _ { j } B _ { j } } ( \sum _ { j } r ( B _ { j } ) ^ { s } ) = 0.$$

Here ( B j ) is a countable set of open balls of radii r B ( j ) which cover A . For example, if the Lebesque measure of A is equal to 0, then (7.1) holds for s = n , hence δ A ( ) ≤ n . The Hausdorff dimension coincides with the Lebesque measure if the latter is positive and finite. A countable set has the Hausdorff measure equal to zero. Also it is known that

the topological dimension of A is less than or equal to its Hausdorff dimension.

The Hausdorff dimension is closely related to the fractal dimension of a fractal set A , i.e. a set that can be subdivided in some finite number N λ ( ) of subsets, all congruent (by translation or rotation) to one another and each equal to a scaled copy of A by a linear factor λ . It is defined to be equal to log N s ( ) log(1 /λ ) . For example, the Cantor set consists of two parts A 1 and A 2 (contained in the interval [0 , 1 3] and / [1 / 3 1]), , each rescaled version of the set with the scaling factor 1 / 3. Thus its fractal dimension is equal to log 2 / log 3 &lt; 1. It is clear that the Hausdorff dimension of a bounded fractal set of diameter D = 2 R is less than or equal than the fractal dimension. In fact, such a set can be covered by N λ ( ) balls of radius λR , or by N λ ( ) 2 balls of radius λ R 2 . Since λ &lt; 1, we get µ s ( A ) = lim k →∞ N λ ( ) k ( λ R k ) s which is zero if N λ λ ( ) s &lt; 1 or s &gt; log N/ log(1 /λ ). In fact, the Hausdorff dimension of the Cantor set coincides with its fractal dimension log 2 / log 3.

By a theorem of D. Sullivan [37], for a geometrically finite nonelementary discrete group Γ, the Hausdorff dimension δ Γ of Λ(Γ) is positive and coincides with the critical exponent of Γ equal to

$$\inf \{ s > 0 \colon \sum _ { g \in \Gamma } e ^ { - s d ( x _ { 0 }, g ( x _ { 0 } ) ) } < \infty \},$$

where x 0 is any point on H n . Using this equality Sullivan shows that

$$\delta _ { \Gamma } = \overline { \lim _ { T \rightarrow \infty } } \frac { \log N _ { T } } { R },$$

where N T is the number of orbit points y with hyperbolic distance from x 0 less than or equal than R . He further shows in [38], Corollary 10, under the additional assumption that Γ has no parabolic fixed points, that there exists constants c, C such that, as T →∞ ,

$$C e ^ { T \delta _ { \Gamma } } \leq N _ { T } \leq C e ^ { T \delta _ { \Gamma } }.$$

In particular, asymptotically, as T →∞ ,

$$N _ { T } \sim c ( T ) e ^ { T \delta _ { \Gamma } },$$

where lim T →∞ c T ( ) T = 0.

If δ Γ &gt; 1 2 ( n -1), then P. Lax and R. Phillips show that, for any geometrically finite non-elementary discrete group Γ, the function c T ( ) is a constant depending only on Γ. When Γ is of finite covolume, then δ Γ is known to be equal to n -1, and the result goes back to A. Selberg.

The assumption on δ Γ has been lifted by T. Roblin [35]. Applying (2.1), we will be able to obtain the assertion of Theorem 1.2 from the

introduction in the case when C 2 &lt; 0. To do the remaining cases where ( e, e ) is non-negative we need to replace the family of hyperbolic balls with another family of sets of growing volume.

glyph[negationslash]

Let e be a nonzero vector in R n, 1 . We assume that ( e, e ) = ± 1 if ( e, e ) = 0 and we fix an isotropic e ′ with ( e, e ′ ) = -1 if ( e, e ) = 0. The group G = SO( n, 1) 0 acts in the projective space P R ( n, 1 ) with three orbits G · [ e ]. Let G [ e ] be the stabilizer subgroup of [ e ] in G . If ( e, e ) &lt; 0, then G [ e ] is a maximal compact subgroup in K ⊂ G isomorphic to SO( n ) and the orbit G · [ e ] is the homogeneous space of left cosets G/K isomorphic to H n .

If ( e, e ) = 0, then G [ e ] is a parabolic subgroup of G . If we identify Q \{ [ e ] } with the Euclidean space E n -1 (via the projecting from [ e ]), we obtain a surjective homomorphism from G [ e ] to the group P of affine orthogonal transformations (motions) of the Euclidean space E n -1 . This homomorphism splits by a subgroup G ce of G [ e ] that stabilizes any nonzero vector ce on the line [ e ]. The kernel of the homomorphism is the group A l of hyperbolic translations along the hyperbolic line l corresponding to the subspace spanned by e and e ′ . This shows that G [ e ] = O n ( -1) · A N l · [ e ] , where O n ( -1) is the group of rotations around the line l and N [ e ] ∼ = R n -1 is the normal subgroup of translations of P , called a horospherical subgroup of G . The decomposition of G [ e ] from above is induced by the Iwasawa decomposition G = G x · A l · N [ e ] of G , where x is any point on l . The group G [ e ] acts transitively on the quadric Q and on H p .

If ( e, e ) &gt; 0, then G [ e ] ∼ = SO( n -1 1) , the homogeneous space , 0 G [ e ] \ G is known as a de Sitter space . By taking the orthogonal complement of the line spanned by e , we identify the points of the de Sitter space with oriented hyperplanes H e in H n . For example, if n = 2, the de Sitter space is the set of oriented geodesic lines in the hyperbolic plane.

Let us consider the pencil of geodesic lines P ( e ) defined by indefinite planes U ⊂ R n, 1 containing e . If ( e, e ) &lt; 0, then P ( e ) consists of geodesic lines containing the point x 0 = [ e ] ∈ H n . It is called an elliptic pencil . If ( e, e ) &gt; 0, it consists of parallel geodesic lines perpendicular to the hyperplane H e . It is called a hyperbolic pencil . Finally, if ( e, e ) = 0, it consists of geodesic lines whose closure in H n contains the point [ e ]. The pencil is called parabolic pencil with center at [ e ].

The parametric equation of the geodesic line from the pencil P ( e ) is equal to

$$\lambda ^ { \dagger \dagger \dagger \dagger \dagger \dagger \dagger \dagger } _ { \L } \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big.. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big. \Big.
 \gamma ( t ) = \begin{cases} v \sinh t + e \cosh t, ( v, e ) = 0, ( v, v ) = 1 & \text{if } ( e, e ) = - 1, \\ v \exp ( t ) + e \sinh t, ( v, e ) = 0, ( v, v ) = - 1 & \text{otherwise.} \end{cases}$$

(see [1], Chapter 4, 2.3). In the elliptic case, we have ( γ t , e ( ) ) = -cosh t , so that d γ t , ( ( ) [ e ]) = t . So, moving along the geodesic line for the distance t ≥ 0, we get the set of points in H n equidistant from [ e ]. This is a geodesic sphere

$$H _ { e } ^ { t } \coloneqq \{ [ v ] \in \mathbb { H } ^ { n } \, \colon - ( v, e ) = \cosh t \} = \{ [ v ] \in \mathbb { H } ^ { n } \, \colon d ( [ v ], [ e ] ) = t \}$$

with the center at x 0 = [ e ]. The group G [ e ] ∼ = SO( n ) acts transitively on each H t e and, projecting from [ e ] along the lines from the pencil P ( e ) identifies each geodesic sphere with the absolute as homogenous spaces with respect to G [ e ] . So, the orbits of G [ e ] in H n are Riemannian homogeneous spaces of constant positive curvature.

If ( e, e ) = 0, we have ( γ t , e ( ) ) = exp( ). t The orbits of N [ e ] are horospheres

$$H _ { e } ^ { t } \coloneqq \{ [ v ] \in \mathbb { H } ^ { n } \colon - ( v, e ) = \exp ( t ) \}.$$

In the vector model of H n , a horosphere is equal to the intersection of H n with a sphere in the Euclidean space tangent to Q at the point [ e ]. In the Klein model, they are ellipsoids tangent to the boundary of the ball. Any horosphere inherits a Riemannian metric of zero curvarure and hence isomorphic to an Euclidean space.

If ( e, e ) &gt; 0, we have ( v, e ) = ( γ (0) , e ) = 0, i.e. [ γ (0)] ∈ H e . Since ( γ t , e ( ) ) = -sinh t , after time t we move from a point on the hyperplane H e to a point [ γ t ( )] on the hypersurface

$$H _ { e } ^ { t } = \{ [ w ] \in \mathbb { H } ^ { n } \, \colon - ( w, e ) = \sinh t \}.$$

The hypersurfaces of this form are called equidistant hypersurfaces . They are orbits of the group G [ e ] and inherit a Riemannian metric of constant negative curvature. In the Klein model, they are ellipsoids whose closures in H n intersect the closure ¯ H e of the hyperplane H e along its boundary. For example, if n = 2, they are ellipses tangent to the absolute at the two points in which the closure of the geodesic line H e intersects the absolute.

In each of the three cases, let us consider the following sets B T ( e ).

$$\text{If } ( e, e ) = - 1,$$

$$B _ { T } ( e ) = \bigcup _ { t = 0 } ^ { T } H _ { e } ^ { t } = \{ [ v ] \in \mathbb { H } ^ { n } \, \colon 1 \leq - ( v, e ) \leq \cosh T \}.$$

This is just a hyperbolic closed ball with center at [ e ] and radius T . If ( e, e ) = 0,

$$B _ { T } ( e ) = \bigcup _ { t = 0 } ^ { T } H _ { e } ^ { t } = \{ [ v ] \in \mathbb { H } ^ { n } \, \colon 1 \leq | ( v, e ) | \leq \exp ( T ) \}.$$

If ( e, e ) = 1,

$$B _ { T } ( e ) = \bigcup _ { t = 0 } ^ { T } H _ { e } ^ { t } = \{ [ v ] \in \mathbb { H } ^ { n } \colon | ( v, e ) | \leq \sinh ( T ) \}.$$

Fix a point x 0 = [ v 0 ] ∈ H n and let K = G x 0 so that G/K = H n . Let H = { } 1 (resp. N [ e ] , resp. G [ e ] if ( e, e ) = -1 (resp. ( e, e ) = 0, resp. ( e, e ) = 1).

Consider the following subsets of G/H

$$\mathcal { S } _ { T } ( e, x _ { 0 } ) \coloneqq G _ { x _ { 0 } } \cdot A _ { T } \cdot H / H,$$

where A T = { a t ∈ A : 0 ≤ t ≤ log T } .

The following result from [33], Theorem 1.2 and [28], Corollary 7.14 (see a nice survey of some of these results in [32]) is crucial for our applications.

Theorem 7.1. Let Γ be a torsion free geometrically finite non-elementary discrete subgroup of SO ( n, 1) 0 . Let [ e ] ∈ P R ( n, 1 ) . If ( e, e ) ≥ 0 , assume that the orbit O Γ ([ e ]) is a discrete set and also that δ Γ &gt; 1 if ( e, e ) &gt; 0 . Then

$$\lim _ { T \to \infty } \frac { \Gamma \cdot H / H \cap \mathcal { S } _ { T } ( e, x _ { 0 } ) } { T ^ { \delta _ { \Gamma } } } = c _ { \Gamma, x _ { 0 }, [ e ] }, \\. \quad \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \ dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots\dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dts \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \sigma \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \Dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots \dots$$

where c Γ ,x 0 , [ e ] is a positive constant that depends only on Γ , x 0 , [ e ] .

Note that the results from [28] also give error terms.

Corollary 7.2. Keep the same assumptions. Fix a point x 0 ∈ H n .

Then

$$\lim _ { t \to \infty } \frac { \# \{ [ v ] \in O _ { \Gamma } ( x _ { 0 } ) \, \colon | ( e, v ) | \leq T \} } { T ^ { \delta _ { \Gamma } } } = c _ { \Gamma, x _ { 0 }, [ e ] }, \\ \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \colots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots \cdots$$

where c Γ ,x 0 , [ e ] is a positive constant that depends only on Γ , x 0 , [ e ] .

Remark 7.3 . The assumption δ Γ &gt; 1 is not needed if the group Γ ∩ N [ e ] acts on ∂H e without parabolic fixed points. Also, the assumption that Γ is torsion free made in loc.cit. is needed only for an explicit formula for the constant.

Remark 7.4 . In the case when ( e, e ) &lt; 0 (resp. ( e, e ) = 0 and n = 3), Theorem 7.1 follows from [33], Theorem 1.2 (resp. [22], Theorem 2.10) that gives an asymptotic of the number of orbit points in a ball { [ v ] ∈

H n : || v || ≤ T } , where || x || is the Euclidean norm. One has only use that if v = ( x , . . . , x 0 n ) with ∑ n i =0 x 2 i ≤ T 2 and -x 2 0 + ∑ n i =1 x 2 i = -1 (resp. = 0), then ( v, (1 , . . . , 0)) 2 = x 2 0 ≤ ( T 2 +1) 2 (resp. / ≤ T 2 ).

In the case ( e, e ) &gt; 0, and Γ and Γ ∩ G [ e ] are of finite covolume (the second condition is always satisfied if n ≥ 3 [10]), Theorem (7.1) follows from [14], [15]. In this case δ Γ = n -1.

## 8. Algebraic geometrical realization

↦

Let X be a smooth projective algebraic surface over an algebraically closed field k . The group of its automorphisms defines a natural left action γ → ( g -1 ) ∗ ( γ ) on its lattice Num( X ) of algebraic cycles modulo numerical equivalence. As usual, we denote its rank by ρ X ( ), this is the Picard number of X . By the Hodge Index Theorem, the intersection form on algebraic cycles defines a symmetric bilinear form on Num( X ) of signature (1 , n ). We denote by V X the corresponding real vector space Num( X ) R equipped with the induced symmetric bilinear form of signature (1 , n ). We assume that the image Aut( X ) ∗ of the homomorphism Aut( X ) → O(Num( X )) is infinite. It is known that this could happen only if X is an abelian surface, a K3 surface, an Enriques surface, or a rational surface. Moreover, n must be larger than or equal to 2, and in the last case, n ≥ 9). All of this is rather well-known, see, for example, [11]. Let Γ be an infinite subgroup of Aut( X ) ∗ . Since the automorphism group preserves the connected component of { v ∈ Num( X ) R : ( v, v ) &gt; 0 } containing an ample class on X , we obtain that Γ is a discrete subgroup in O( V X ) = O( ∼ n, 1). Also, since Γ preserves an integral lattice, its orbits are discrete sets. Thus, we can apply Theorem 7.1 to obtain the following.

Theorem 8.1. Assume that Γ is a geometrically finite non-elementary subgroup of Aut ( X ) ∗ . Let H,C be effective numerical divisor classes in Num ( X ) ′ such that H is ample. Then

$$( 8. 1 ) \quad \lim _ { t \to \infty } \frac { \# \{ C ^ { \prime } \in O _ { \Gamma } ( C ) \, \colon ( H, C ^ { \prime } ) / ( H, C ) \leq T \} } { T ^ { \delta _ { \Gamma } } } = c _ { \Gamma, C, H },$$

where c Γ ,H,C is a positive constant depending only on Γ [ , H , C ] [ ] but does not depend on C 2 and H 2 .

Now let us give some examples where we can apply the previous result. First, we note that we may always consider Γ up to a finite group, i.e. we either consider its subgroup of finite index, or extend it by a finite group, or do both. It affects only the constant in the asymptotic formula.

glyph[negationslash]

Let X be a complex algebraic K3 surface. It is known that the subgroup of O(Num( X )) generated by transformations g , g ∗ ∈ Aut( X ) and reflections s , r ( r, r ) = -2 is a subgroup of finite index (this is also true if the characteristic p = 2 [24]). Thus, if Num( X ) does not contain classes r with ( r, r ) = -2, then Aut( X ) is a Kleinian group of finite covolume. 9 The Hausdorff dimension δ Γ in this case is equal to n -1. So, it is more interesting to realize proper thin Kleinian groups of automorphisms.

Assume that we find a surface X and vectors α , . . . , α 1 r ∈ Num( X ) such that 2 ( α ,α i i ) α i ∈ Num( X ) ∨ , so that the reflections s α i define orthogonal transformations of Num( X ). Moreover, assume that there exist automorphisms g i ∈ Aut( X ) such that g ∗ i = s α i . Let P be the convex polyhedron with normal vectors e i = 1 ( α ,α i i ) 1 / 2 α i . Assume that it is a Coxeter polyhedron of level ≤ 2 with the reflection group Γ P generated by s α 1 , . . . , s α r . Then the limit set of Γ P is a Boyd-Maxwell sphere packing and its Hausdorff dimension δ Γ P is equal to the critical exponent (1.1). Now we can apply Theorem 8.1.

Example 8.2. Let M = A p (2) ev ∼ = ( A 1 ⊕ U )(2) be the even sublattice of the Apollonian lattice and Ap 2 be the Apollonian group. Let ω 1 = e 1 , . . . , ω 4 = e 4 be a basis in A p (2) with Gram matrix equal to G = cir(1 , -1 , -1 , -1) = ( g ij ). The lattice M ( -1 / 2) = ∼ A 1 ( -1) ⊕ U can be primitively embedded in the lattice Num( X ) = ∼ U ⊕ E 8 ( -1), where X is an Enriques surface. We identify the image with M ( -1 / 2). The group Ap 2 contains a subgroup of index 2 that consists of isometries of M ( -1 / 2) that extend to the whole lattice acting identically on the orthogonal complement. Thus the subgroup Γ = Ap ′ 2 ∩ O(Num( X ) ′ is realized by automorphisms of X . One can show that the index of Γ in Ap 2 is equal to 4. Now we can apply Theorem 8.1 with δ Γ = δ Ap 2 .

Of course, one can also realize A p (2)( -1) as the Picard lattice of a K3 surface X . Since the lattice does not represent -2, a subgroup of finite index of Ap 2 is realized as a group of automorphisms of a K3 surface. according to Vinberg's classification of 2-reflexive lattices of rank 4 [41], the 2-reflection group of the lattice A p (2) is of infinite index in the orthogonal group of the lattice. Thus Ap 2 is a thin group geometrically finite group of automorphisms of X .

The same is true for the group Ap ; 3 it contains a subgroup of finite index that can be realized as a thin geometrically finite group of automorphisms of either an Enriques or a K3 surface.

9 By Riemann-Roch theorem, such a divisor class is either effective or antieffective, and one of its irreducible component is a smooth rational curve.

Example 8.3. Consider the Coxeter group Γ( a, b, c ) of a Coxeter polyhedron P in H 2 with Gram matrix ( 1 -a -b -a 1 -c -b -c 1 ) , where a, b, c are rational numbers ≥ 1. If a = b = c = 1, the fundamental domain is an ideal triangle. If a, b, c &gt; 2, the fundamental triangle is a pair of pants as on the following picture.

Figure 7: Hyperbolic triangle with no vertices

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-a1fe4854-fff7-4a62-b2be-43a0424c96ac.jpg)

Let L be the quadratic form defined by the integral matrix A = -eG P ( ), where e is the exponent of the Coxeter matrix. Let L ∨ be the dual lattice. Suppose we can embed M = int( L ∨ ) ev in Num( X ) of some K3 surface and realize the reflection group Γ P as a group of automorphisms of X . Then we can apply Theorem 8.1 to any nef divisor class contained in the sublattice M .

For example, let us consider the case ( a, b, c ) = ( a, a, 1). The fundamental triangle has one ideal point and looks like in the following figure 8.3.

↦

Let H → 2 H 2 be the map from the upper-half plane to the unit disk given by the map z → z -i z + i . One can show that the pre-image of the sides of our triangle are the lines x = 1 , x = -1 and the upper halfcircle of radius r = 1 /a with center at the origin. Recall that a ≤ -1 so that the half-circle is between the vertical lines. Let us re-denote our group by Γ . r It was shown by C. McMullen in [25] that as r → 0, we have

$$\delta _ { \Gamma _ { r } } = \frac { r + 1 } { 2 } + O ( r ^ { 2 } ),$$

while for r → 1, we have

$$\delta _ { \Gamma _ { r } } \sim 1 - \sqrt { 1 - r }.$$

We have

$$G ( P ) ^ { - 1 } = \frac { 1 } { 4 a ^ { 2 } } \begin{pmatrix} 0 & - 2 a & - 2 a \\ - 2 a & a - 1 & - a - 1 \\ - 2 a & - a - 1 & a - 1 \end{pmatrix}$$

If a &gt; 1 is odd, the matrix 2 a G P 2 ( ) -1 defines an even integral lattice L r . The group Γ r acts on this lattice as a reflection group in the sides of the triangle. The lattice L r ( -2) is realized as the Picard lattice of a K3 surface X . Since L r ( -2) does represent -2, the surface X does not contain smooth rational curves. This implies that a subgroup of finite index of Γ r acts on X by automorphisms. So, we may apply Theorems 8.1 with the Hausdorff dimension computed by McMullen. Note that when r = a = 1, the lattice L r is isomorphic to the Apollonian lattice A p 1 and the group Γ 1 is isomorphic to the Apollonian group Ap . 1

glyph[negationslash]

Let us give another example of a realizable group Γ( a, b, c ). It is taken from [3]. Let X be a K3 surface defined over an algebraically closed field of characteristic = 2 embedded in P 2 × P 2 as a complete intersection of hypersurfaces of multi-degree (1 , 1) and (2 2). , Let p , p 1 2 : X → P 2 be the two projections. They are morphisms of degree 2 branched along a plane curve B i of degree 6. We assume that B 1 is nonsingular and B 2 has a unique double point q 0 so that the fiber p -1 2 ( q 0 ) is a smooth rational curve R that is mapped isomorphically under p 1 to a line. We assume that X is general with these properties. More precisely, we assume that Pic( X ) has a basis ( h , h 1 2 , r ), where h i = p ∗ i (line) and r is

Figure 8: Hyperbolic triangle with one ideal point

![image](https://ac.marcodsn.me/arxiv-markdown-images/1408.0374v2-figure-4e553b01-2c95-4d31-9b2e-00c868171125.jpg)

the class of R . The intersection matrix of this basis is equal to

$$\begin{pmatrix} 2 & 4 & 1 \\ 4 & 2 & 0 \\ 1 & 0 & - 2 \end{pmatrix}.$$

Let s = p ∗ 1 (( p 1 ) ∗ ( r )) -r = h 1 -r . It is a class of smooth rational curve S on X . The pre-image of the pencil of lines through q 0 is an elliptic pencil | F | on X with [ F ] = h 2 -r . The curve S is a section of the elliptic fibration defined by the linear system | F | and the curve R is its 2-section that intersects S with multiplicity 3. Consider the following three automorphisms of X . The first two Φ 1 and Φ 2 are defined by the birational deck transformations of the covers p 1 and p 2 . The third one Φ 3 is defined by the negation automorphism of the elliptic pencil with the group law defined by the choice of S as the zero section.

It is easy to compute the matrix of each Φ i in the basis ( f, s, r ) = ([ F , ] [ S , ] [ R ]) with the Gram matrix

$$\begin{pmatrix} 0 & 1 & 2 \\ 1 & - 2 & 3 \\ 2 & 3 & - 2 \end{pmatrix}.$$

We have Φ ( ) = ∗ 1 s r, Φ ( ) = ∗ 1 r s and f ′ = Φ ( ∗ 1 f ) = af + bs + cr . Since Φ 2 i is the identity and ( f, f ) = 0, we get a = -1 and b = c . Since ( f , s ′ ) = ( f, r ) = 2, we easily get b = c = 3. The matrix of Φ , and 1 similarly obtained matrices of Φ 2 and Φ 3 are as follows.

$$A _ { 1 } = \begin{pmatrix} - 1 & 0 & 0 \\ 3 & 0 & 1 \\ 3 & 1 & 0 \end{pmatrix}, \quad A _ { 2 } = \begin{pmatrix} 1 & 4 & 0 \\ 0 & - 1 & 0 \\ 0 & 1 & 1 \end{pmatrix}, \quad A _ { 3 } = \begin{pmatrix} 1 & 0 & 1 4 \\ 0 & 1 & 4 \\ 0 & 0 & - 1 \end{pmatrix}.$$

The transformations Φ ′ 1 = Φ 1 ◦ Φ 2 ◦ Φ 1 , Φ 2 , Φ 3 are the reflections with respect to the vector α i , where

$$\alpha _ { 1 } = - 4 f + 1 3 s + 1 0 r, \ \alpha _ { 2 } = 4 f - 2 s + r, \ \alpha _ { 3 } = 7 f + 2 s - r.$$

The Gram matrix of the vectors α , α , α 1 2 3 is equal to

$$G = \begin{pmatrix} - 2 2 & 1 4 3 & 2 2 0 \\ 1 4 3 & - 2 2 & 2 2 \\ 2 2 0 & 2 2 & - 2 2 \end{pmatrix} = - 2 2 \begin{pmatrix} 1 & - \frac { 1 3 } { 2 } & - 1 0 \\ - \frac { 1 3 } { 2 } & 1 & - 1 \\ - 1 0 & - 1 & 1 \end{pmatrix}.$$

So, the group generated by Φ ′ 1 , Φ 2 , Φ coincides with the triangle group 3 Γ( 13 2 , 10 1). , The fundamental triangle P has one ideal vertex. The reflection group Γ P is a subgroup of infinite index of the group Γ of automorphisms of X generated by Φ 1 , Φ 2 , Φ . 3 Baragar proves that

Γ is isomorphic to Aut( X ) (for sufficiently general X ). He finds the following bounds for δ Γ

$$. 6 5 1 5 < \delta _ { \Gamma } <. 6 5 3 8.$$

This implies that

$$\delta _ { \Gamma _ { P } } <. 6 5 3 8.$$

Example 8.4. This is again due to Baragar [4]. We consider a nonsingular hypersurface X in P 1 × P 1 × P 1 of type (2 2 2). , , It is a K3 surface whose Picard lattice contains the Apollonian lattice A p (1). If X is general, then the Picard lattice coincides with this lattice. We assume that one of the projections p ij : X → P 1 × P 1 , say p 12 , contains the whole P 1 as its fiber over some point q 0 ∈ P 1 × P 1 . All the projections are degree 2 maps. Let F , i i = 1 2 3 , , , be the general fibers of the projections p i : X → P 1 . Each F i is an elliptic curve whose image under the map p jk is a divisor of type (2 , 2). Let f , f 1 2 , f 3 , r be the classes of the curves F , F , F , R 1 2 3 . We assume that X is general with these properties so that Pic( X ) is generated by these classes. The Gram matrix of this basis is equal to

$$\begin{pmatrix} 0 & 2 & 2 & 0 \\ 2 & 0 & 2 & 0 \\ 2 & 2 & 0 & 1 \\ 0 & 0 & 1 & - 2 \end{pmatrix}.$$

It is easy to see that

$$\text{Pic} ( X ) \cong \text{U} \oplus \left ( \begin{smallmatrix} - 4 & 2 \\ - 2 & - 8 \end{smallmatrix} \right ).$$

According to Vinberg's classification of 2-reflective hyperbolic lattices of rank 4 [1], the Picard lattice is not 2-reflective. Hence the image of the group Aut( X ) in O(Pic( X )) is of infinite index.

Let Φ ij be the automorphisms of X defined by the deck transformations of the projections p ij . Let Φ ′ 4 be defined as the transformation Φ in the previous example with respect to the elliptic pencil 3 | F 3 | with section R . The transformation Φ ∗ 12 leaves the vectors f , f 1 2 , r invariant, and transforms f 3 to 2 f 1 +2 f 2 -f 3 -r . Thus Φ ∗ 12 is the reflection with respect to the vector α 1 = -2 f 1 -2 f 2 +2 f 3 + . r

The transformation Φ ∗ 13 leaves f , f 1 3 invariant and transforms r in r ′ = f 1 -r . It also transforms f 2 to some vector f ′ 2 = af 1 + bf 2 + cf 3 + dr . Computing ( f , f ′ 2 1 ) = ( f , f 2 1 ) , ( f , f ′ 2 3 ) = ( f , f 2 3 ) , ( f , r ′ 2 ) = ( f , f 2 1 -r ), we find that f ′ 2 = 2 f 1 -f 2 + 2 f 3 . Similarly, we find that Φ ∗ 23 ( f 2 ) = f , 2 Φ ∗ 23 ( f 3 ) = f , 3 Φ ∗ 23 ( r ) = f 2 -r and Φ ∗ 23 ( f 2 ) = -f 1 +2 f 2 +2 f 3 .

It follows from the definition of a group law on an elliptic curve that

$$\Phi _ { 4 } ^ { \prime } ( f _ { 3 } ) = f _ { 3 }, \ \Phi _ { 4 } ^ { \prime } ( r ) = r, \ \Phi _ { 4 } ^ { \prime } ( f _ { i } ) = - f _ { i } + 8 f _ { 3 } + 4 r, \ i = 1, 2$$

Consider the transformations

$$\Phi _ { 1 } = \Phi _ { 1 2 }, \, \Phi _ { 2 } = \Phi _ { 1 3 } \circ \Phi _ { 1 2 } \circ \Phi _ { 1 3 }, \, \Phi _ { 3 } = \Phi _ { 2 3 } \circ \Phi _ { 1 2 } \circ \Phi _ { 2 3 }, \, \Phi _ { 4 } = \Phi ^ { \prime } _ { 4 } \circ \Phi _ { 1 2 } \circ \Phi ^ { \prime } _ { 4 }.$$

′ 4 ′ 4 .

These transformations act on Pic( X ) as the reflections with respect to the vectors

$$r.$$

$$\alpha _ { 1 } \ & = \ \ - 2 f _ { 1 } - 2 f _ { 2 } + 2 f _ { 3 } + r, \\ \alpha _ { 2 } \ & = \ \ \Phi _ { 1 3 } ^ { * } ( \alpha _ { 1 } ) = - 5 f _ { 1 } + 2 f _ { 2 } - 2 f _ { 3 } - r, \\ \alpha _ { 3 } \ & = \ \ \Phi _ { 2 3 } ^ { * } ( \alpha _ { 1 } ) = 2 f _ { 1 } - 5 f _ { 2 } - 2 f _ { 3 } - r, \\ \alpha _ { 4 } \ & = \ \ \Phi _ { 4 } ( \alpha _ { 1 } ) = 2 f _ { 1 } + 2 f _ { 2 } - 3 0 f _ { 3 } - 1 5 r. \\ \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot \cdot$$

The Gram matrix of these four vectors is equal to

$$\begin{pmatrix} - 1 4 & 1 4 & 1 4 & 2 1 0 \\ 1 4 & - 1 4 & 8 4 & 1 8 2 \\ 1 4 & 8 4 & - 1 4 & 1 8 2 \\ 2 1 0 & 1 8 2 & 1 8 2 & - 1 4 \end{pmatrix} = - 1 4 \begin{pmatrix} 1 & - 1 & - 1 & - 1 5 \\ - 1 & 1 & - 6 & - 1 3 \\ - 1 & - 6 & 1 & - 1 3 \\ - 1 5 & - 1 3 & - 1 3 & 1 \end{pmatrix}$$

Let P be the Coxeter polytope defined by this matrix. The Coxeter group Γ P is generated by the reflections Φ ∗ i , i = 1 2 3 4. , , ,

Baragar proves that the automorphisms Φ ij and Φ ′ 4 generate a subgroup Γ of Aut( X ) of finite index. His computer experiments suggest that

$$1. 2 8 6 < \delta _ { \Gamma } < 1. 3 0 6.$$

Our reflection group Γ P generated by Φ 1 , . . . , Φ 4 is of infinite index in Γ. So, we obtain

$$\delta _ { \Gamma _ { P } } < 1. 3 0 6.$$

Example 8.5. We consider a general Coble rational surface [8]. The orthogonal complement of the canonical class in Num( X ) is isomorphic to U ⊕ E 8 ( -1). One can prove that the automorphism group of X is isomorphic to the automorphism group of a general Enriques surface [8] (true in any characteristic). We do the same, as in Example 8.2 to realize Γ as an automorphism group of X . Since X is rational, it gives a realization of Γ as a group of Cremona transformations of P 2 . Taking the class of a line, we obtain the asymptotic of the growth of the function deg Φ, where deg Φ is the algebraic degree of a Cremona transformation Φ from Γ.

Remark 8.6 . Taking into account Remark 6.1, we see that, for any even hyperbolic lattice of rank &gt; 4 that contains primitively the lattice M = U ⊕ A 1 ( -1) ⊕ A 1 ( -1), the orthogonal group O( L ) contains finitely generated subgroups which are not geometrically finite. Embedding M primitively in the lattice U ⊕ E 8 ( -1), we obtain that a general Enriques

surface contains finitely generated groups of automorphisms which are not geometrically finite.

## References

- [1] D. Alekseevskii, E. Vinberg, A. Solodovnikov, Geometry of spaces of constant curvature . Geometry, II, 1-138, Encyclopaedia Math. Sci., 29, Springer, Berlin, 1993.
- [2] B. Apanasov, Conformal geometry of discrete groups and manifolds . de Gruyter Expositions in Mathematics, 32. Walter de Gruyter, Berlin, 2000.
- [3] A. Baragar, Orbits of curves on certain K3 surfaces . Compositio Math. 137 (2003), 115-134.
- [4] A. Baragar, The ample cone for a K3 surface . Canad. J. Math. 63 (2011), 481-499.
- [5] Y. Benoist, P. de la Harpe, Adh´ erence de Zariski des groupes de Coxeter . Compos. Math. 140 (2004), 1357-1366.
- [6] D. Boyd, A new class of infinite sphere packings . Pacific J. Math. 50 (1974), 383-398
- [7] D. Boyd, The sequence of radii in an Apollonian packing , Math. Comp. 39 (1982) 249-254.
- [8] S. Cantat, I. Dolgachev, Rational surfaces with a large group of automorphisms . J. Amer. Math. Soc. 25 (2012), 863-905.
- [9] H. Chen, J.-Ph. Labb´ e, Lorentzian Coxeter systems and Boyd-Maxwell ball packings. arXiv:1310.8608, math.GR.
- [10] S. Dani, On invariant measures, minimal sets and a lemma of Margulis . Invent. Math. 51 (1979), 239-260.
- [11] I. Dolgachev, Reflection groups in algebraic geometry . Bull. Amer. Math. Soc. (N.S.) 45 (2008), 1-60.
- [12] I. Dolgachev, A brief introduction to Enriques surfaces , Development of moduli theory, Advanced Studies in Pure Mathematics, Math. Soc. Japan, to appear.
- [13] I. Dolgachev, B. Howard, Configuration spaces of complex and real spheres , Recent Advances in Algebraic Geometry, ed. C. Hacon, M. Mustata and M. Popa, London Math. Soc. Lect. Notes, 2014, 156-179 (to appear).
- [14] W. Duke, Z. Rudnick, P. Sarnak, P., Density of integer points on affine homogeneous varieties . Duke Math. J. 71 (1993), 143-179.
- [15] A. Eskin, C. McMullen, Mixing, counting, and equidistribution in Lie groups . Duke Math. J. 71 (1993), 181-209.
- [16] E. Fuchs, Counting problems in Apollonian packings . Bull. Amer. Math. Soc. (N.S.) 50 (2013), 229-266.
- [17] R. Graham, J. Lagarias, C. Mallows, A. Wilks, C. Yan, Apollonian circle packings: number theory . J. Number Theory 100 (2003), 1-45.
- [18] R. Graham, J. Lagarias, C. Mallows, A. Wilks, C. Yan, Apollonian circle packings: geometry and group theory. I. The Apollonian group . Discrete Comput. Geom. 34 (2005), 547-585.
- [19] R. Graham, J. Lagarias, C. Mallows, A. Wilks, C. Yan, Apollonian circle packings: geometry and group theory. II. Super-Apollonian group and integral packings . Discrete Comput. Geom. 35 (2006), 1-36.

- [20] R. Graham, J. Lagarias, C. Mallows, A. Wilks, C. Yan, Apollonian circle packings: geometry and group theory. III. Higher dimensions . Discrete Comput. Geom. 35 (2006), 37-72.
- [21] M. Kapovich, L. Potyagailo, E. Vinberg, Noncoherence of some lattices in Isom( H n ) . The Zieschang Gedenkschrift, 335-351, Geom. Topol. Monogr., 14, Geom. Topol. Publ., Coventry, 2008.
- [22] A. Kontorovich, H. Oh, Apollonian circle packings and closed horospheres on hyperbolic 3-manifolds . With an appendix by H.Oh and N. Shah. J. Amer. Math. Soc. 24 (2011), 603-648.
- [23] P. Lax, R. Phillips, The asymptotic distribution of lattice points in Euclidean and non-Euclidean spaces . J. Funct. Anal. 46 (1982), 280-350.
- [24] M. Lieblich, D. Maulik, A note on the cone conjecture for K3 surfaces in positive characteristic , math.AG arXiv:1102.3377.
- [25] C. McMullen, Hausdorff dimension and conformal dynamics. III. Computation of dimension . Amer. J. Math. 120 (1998), 691-721.
- [26] G. Maxwell, Sphere packings and hyperbolic reflection groups . J. Algebra 79 (1982), 78-97.
- [27] G. Maxwell, Wythoff's construction for Coxeter groups . J. Algebra 123 (1989), 351-377.
- [28] A. Mohammadi, H. Oh, Matrix coefficients, Counting and Primes for orbits of geometrically finite groups . arXiv:1208.4139 math.NT.
- [29] V. Nikulin, Integral quadratic forms and some of its geometric applications , Izv. Akad. Nauk SSSR, Ser. Math. 43 (1979), 111-177.
- [30] V. Nikulin, On the quotient groups of the automorphism groups of hyperbolic forms modulo subgroups generated by 2-reflections, Algebraic geometric applications , Current Problems of Mathematics, t. 18, VINITI, Moscow, 1981, pp. 3-114 [English translation:J. Soviet Math. 22 (1983), 1401-1476].
- [31] V. Nikulin, K3 surfaces with a finite group of automorphisms and a Picard group of rank three . Algebraic geometry and its applications. Trudy Mat. Inst. Steklov. 165 (1984), 119-142.
- [32] H. Oh, Harmonic analysis, Ergodic theory and Counting for thin groups, MSRI publication Vol. 61 (2014) 'Thin groups and Superstrong approximation' edited by E. Breuilliard and H. Oh,2014.
- [33] H. Oh, N. Shah, Equidistribution and counting for orbits of geometrically finite hyperbolic groups . J. Amer. Math. Soc. 26 (2013), 511-562.
- [34] J. Ratcliffe, Foundations of hyperbolic manifolds . Second edition. Graduate Texts in Mathematics, 149. Springer, New York, 2006.
- [35] T. Roblin, Ergodicit´ e et ´ equidistribution en courbure n´ egative , M´ em. Soc. Math. Fr. (N.S.) 95 (2003), vi+96 pp.
- [36] P. Sarnak, Integral Apollonian packings . Amer. Math. Monthly 118 (2011), 291-306.
- [37] D. Sullivan, Hausdorff measures old and new, and limit sets of geometrically finite Kleinian groups . Acta Math. 153 (1984), 259-277.
- [38] D. Sullivan, The density at infinity of a discrete group of hyperbolic motions . Inst. Hautes Etudes Sci. Publ. Math. ´ 50 (1979), 171-202.
- [39] P. Thomas, D. Dhar, The Hausdorff dimension of the Apollonian packing of circles . J. Phys. A 27 (1994), no. 7, 2257-2268.

42

- [40] E. Vinberg, O. Shvartsman, Discrete groups of motions of spaces of constant curvature . Geometry, II, 139-248, Encyclopaedia Math. Sci., 29, Springer, Berlin, 1993.
- [41] E. Vinberg, Classification of 2-reflective hyperbolic lattices of rank 4 . Tr. Mosk. Mat. Obs. 68 (2007), 44-76; translation in Trans. Moscow Math. Soc. 2007, 39-66.

Department of Mathematics, University of Michigan, 525 E. University Av., Ann Arbor, MI, 49109 E-mail address : idolga@umich.edu
Igor Dolgachev
Orbital counting of curves on algebraic surfaces and sphere packings
We realize the Apollonian group associated to an integral Apollonian circle
packings, and some of its generalizations, as a group of automorphisms of an
algebraic surface. Borrowing some results in the theory of orbit counting, we
study the asymptotic of the growth of degrees of elements in the orbit of a
curve on an algebraic surface with respect to a geometrically finite group of
its automorphisms.