Jump to content

Circle packing theorem

From Wikipedia, the free encyclopedia
(Redirected from Coin graph)

A circle packing and its graph of tangencies

The circle packing theorem (also known as the Koebe–Andreev–Thurston theorem) describes the possible patterns of tangent circles among non-overlapping circles in the plane. A circle packing is a collection of circles whose union is connected and whose interiors are disjoint. The intersection graph of a circle packing, called a coin graph,[1] is the graph having a vertex for each circle, and an edge for every pair of circles that are tangent. Coin graphs are always connected, simple, and planar. The circle packing theorem states that these are the only requirements for a graph to be a coin graph: every finite connected simple planar graph has a circle packing in the plane whose intersection graph is isomorphic to .

A stronger form of the circle packing theorem applies to any polyhedral graph and its dual graph, and proves the existence of a primal–dual packing, circle packings for both graphs that cross at right angles. Circle packings and their tangencies, and the circle packing theorem, have been extended to arbitrary Riemannian surfaces including the sphere, the hyperbolic plane, and to surfaces of bounded genus. More generally, intersection graphs of interior-disjoint geometric objects are called tangency graphs[2] or contact graphs.[3] As a special case, the coin graphs of packings of unit circles are called penny graphs.[4]

Circle packings have applications in conformal mapping, the construction of polyhedra, planar separator theorems, graph drawing, and the theory of random walks. The study of the tangencies of circle packings, for which the circle packing theorem is central, should be distinguished from the study of circle packings within fixed shapes but without specified tangencies, typically studied with respect to their density (the fraction of area covered by circles).

Theorem statement and proofs

[edit]

A maximal planar graph is a finite simple planar graph to which no more edges can be added while preserving planarity. It may be embedded in the plane with different choices of the outer face, but all such embeddings share the same set of faces (including the outer face) which must all be triangles. The circle packing theorem guarantees the existence of a circle packing with finitely many circles whose intersection graph is isomorphic to . As the following theorem states more formally, every maximal planar graph can have at most one packing.[5]

Koebe–Andreev–Thurston circle packing theoremIf is a finite connected planar graph, then a circle packing in the plane whose tangency graph is isomorphic to exists.[5] If has a fixed planar embedding, the circle packing can be chosen so that the cyclic order of tangencies around each circle corresponds to the order of edges around each vertex in the embedding.[6] If is maximal planar, this packing is unique, up to reflections in lines and Möbius transformations.[5]

The geometric transformations in the uniqueness part of the theorem, reflections and Möbius transformation, preserve circles, and therefore the tangency graph of any system of circles. Uniqueness means that, for any two packings of the same maximal planar graph, there exists a transformation of these types that takes one packing into the other. For every finite planar graph , it is possible to add vertices to to construct a larger maximal planar graph in which is an induced subgraph. Constructing a circle packing from this larger graph, and then removing the circles for the added vertices, produces a circle packing for . This allows the existence of circle packings for arbitrary planar graphs to be derived from the special case for maximal planar graphs.[7]

Existence

[edit]

Paul Koebe's original proof of the existence of a circle packing for any planar graph is based on his conformal uniformization theorem saying that a finitely connected planar domain is conformally equivalent to a circle domain.[8] A circle domain is a domain, a connected open subset of the plane, whose boundary components are circles or points. It is finitely connected when it has finitely many boundary components, and hence finite first Betti number (the number of generators of its fundamental group). The circle packing is a limiting case of Koebe's result, for domains complementary to the union of disks of the desired circle packing. Koebe conjectured that the finite connectivity assumption was unnecessary in his theorem but was unable to prove it. He & Schramm (1993) extend Koebe's theorem to countably connected domains, and to certain packings of countably many circles.[9]

Thurston's proof is based on the Brouwer fixed point theorem.[10] Another proof uses a discrete variant of Perron's method of constructing solutions to the Dirichlet problem.[11] Yves Colin de Verdière proved the existence of the circle packing as a minimizer of a convex function on a certain configuration space,[12] and Chow & Luo (2003) find another proof using a combinatorial analogue of Ricci flow. As MathSciNet reviewer Igor Rivin writes, this is also the gradient flow of a function that "may well be the one introduced by Colin de Verdière".[13]

Another proof starts with a circle packing with the correct number of disks but an incorrect pattern of tangencies, for which existence is easy to prove, such as an Apollonian packing. It then corrects the tangencies, one at a time, by removing one edge from the corresponding maximal planar graph and replacing it by the other diagonal of the resulting quadrilateral, at the same time adjusting the packing to match this change, in effect walking in a flip graph of triangulations, until reaching the desired graph of tangencies.[14]

Uniqueness

[edit]

Thurston observes that this uniqueness is a consequence of the Mostow rigidity theorem. To see this, let be represented by a circle packing. Then the plane in which the circles are packed may be viewed as the boundary of a halfspace model for three-dimensional hyperbolic space; with this view, each circle is the boundary of a plane within the hyperbolic space. One can define a set of disjoint planes in this way from the circles of the packing, and a second set of disjoint planes defined by the circles that circumscribe each triangular gap between three of the circles in the packing. These two sets of planes meet at right angles, and form the generators of a reflection group whose fundamental domain can be viewed as a hyperbolic manifold. By Mostow rigidity, the hyperbolic structure of this domain is uniquely determined, up to isometry of the hyperbolic space; these isometries, when viewed in terms of their actions on the Euclidean plane on the boundary of the half-plane model, translate to Möbius transformations.[15]

Schramm (1991) generalizes the uniqueness of circle packings to certain packings of infinitely many circles on a sphere or open disk. His uniqueness theorem applies to circle packings in which the connected components of the space outside all circles are either circular triangles or individual singular points. The graph of tangencies of such a packing and the triangles determined by the non-singular components form a triangulation of the sphere or disk, punctured at the singular points. Schramm proves that when any two packings have at most a countably infinite set of singularities and have combinatorially equivalent triangulations, they are equivalent under Möbius transformations. This uniqueness argument does not, however, characterize the triangulations for which such a packing exists.[16]

Generalizations

[edit]
Infinite periodic locally finite circle packing
A Doyle spiral, not locally finite at the spiral center

A version of the circle packing applies to some infinite graphs. In particular, an infinite planar triangulation of an open disk has a locally finite packing in either the Euclidean plane or the hyperbolic plane (but not both). Here, locally finite means that every point of the plane has a neighborhood that intersects only finitely many circles. In the Euclidean case, the packing is unique up to similarity; in the hyperbolic case, it is unique up to isometry.[17]

The circle packing theorem generalizes to graphs that are not planar. If is a graph that can be embedded on a surface (more precisely, an orientable 2-manifold), then there is a Riemannian metric of constant curvature on and a circle packing on whose contact graph is isomorphic to . If is closed (compact and without boundary) and is a triangulation of , then and the packing are unique up to conformal equivalence. If is the sphere, then its geometry is spherical geometry and the equivalence is up to Möbius transformations. If it is a torus, the geometry is locally Euclidean (as a flat torus), and the equivalence is up to scaling by a constant and isometries. If has genus at least 2, then the geometry is locally hyperbolic and the equivalence is up to isometries.[18]

Another generalization of the circle packing theorem involves replacing the condition of tangency with a specified intersection angle between circles corresponding to neighboring vertices. One version is as follows. Suppose that is a finite 3-connected planar graph (that is, a polyhedral graph), then there is a pair of primal–dual circle packings: one packing whose intersection graph is isomorphic to , and another whose intersection graph is isomorphic to the planar dual of , such that for every vertex in and face adjacent to it, the circle in the first packing corresponding to the vertex intersects orthogonally with the circle in the second packing corresponding to the face.[1] For instance, applying this result to the graph of the tetrahedron gives, for any four mutually tangent circles, a second set of four mutually tangent circles each of which is orthogonal to three of the first four.[19] The circle packing for a maximal planar graph is automatically part of a pair of primal–dual circle packings, because each triple of tangent circles in the primal packing has a dual circle that passes at right angles through its three points of tangency. Like the circle packing for a maximal planar graphs, primal–dual circle packings are unique up to reflections and Möbius transformations.[20] A further generalization, replacing intersection angle with inversive distance, allows the specification of packings in which some circles are required to be disjoint from each other rather than crossing or being tangent.[21]

Yet another variety of generalizations allow shapes that are not circles.[3] Suppose that each vertex of a finite planar graph corresponds to a smooth convex set . Then there is a packing in the plane whose tangencies correspond one-for-one with the edges of where each set is obtained from by translating and scaling. This result, from Oded Schramm's 1990 Ph.D. thesis,[22] is a special case of Schramm's "monster packing theorem". The "monster" of the theorem is a spanning tree of the given graph, with a function for each vertex specifying how the shape it should be packed with depends on the relative positions of tangencies around all the shapes, satisfying certain consistency conditions. To prove that monsters can always be packed, Schramm followed Thurston's proof of the circle packing theorem, using the Brouwer fixed point theorem. He wrote that "One can just see the terrible monster swinging its arms in sheer rage, the tentacles causing a frightful hiss, as they rub against each other."[23]

Other properties

[edit]

Relative radii

[edit]
Construction for rings of circles achieving the most extreme ratio of inner radius to surrounding radius in the ring lemma

A result known as the ring lemma controls the sizes of adjacent circles in a Euclidean circle packing. According to the ring lemma, if one circle in a packing is surrounded by a ring of others, then the ratios of the inner circle's radius to the radius of each surrounding circle is at most exponential in . More precisely, this ratio is at most , where denotes the th Fibonacci number.[24]

Monotonicity of angles

[edit]

In the triangle connecting the centers of three mutually tangent circles, the angle formed at the center of one of the circles is monotone decreasing in its radius and monotone increasing in the two other radii.[25] Beardon & Stephenson (1991) described related monotonicity properties of circle packings in the hyperbolic plane that they describe as being a discrete analogue to the Schwarz–Pick theorem in complex differential geometry. They consider any finite circle packing in the hyperbolic plane whose coin graph has the structure of a triangulated disk, and they compare this packing to another packing representing the same graph, obtained by treating the line at infinity of the hyperbolic plane as another circle, tangent to all the circles on the outer boundary of the disk (which become horocycles of infinite radius in the second packing). They prove that in the second packing, each circle's radius is non-decreasing, as is the hyperbolic distance between each pair of centers of tangent circles. If any of the radii or distances is equal, then they all are, and the two packings are congruent.[26]

Ordering by radii

[edit]

A Koebe ordering is the sequence of vertices of a planar graph, ordered by their radii in a circle packing (larger to smaller). If there are no ties (which can be achieved by performing a Möbius transformation) then in the resulting ordering each vertex has at most five earlier neighbors, matching the worst-case degeneracy of planar graphs. This property has been generalized as the -admissability of the ordering, the maximum number of disjoint paths of at most steps that all start at the same vertex , continue through vertices later than in the ordering, and end at a vertex earlier than . Any Koebe ordering has -admissability , best possible for planar graphs to within a constant factor.[27]

Triangulated packings of surfaces

[edit]
A periodic triangulated packing of the plane with three radii within 0.651 of each other, conjectured by Connelly & Zhang (2025) to have the closest ratio to 1 of any non-hexagonal packing

Instead of starting with a topological embedding of a graph, and searching for a surface with uniform geometry on which that graph can be represented by a circle packing, several authors have studied a reversed version of this question: which uniform surfaces admit triangulated packings, circle packings for which each gap between circles is a circular horn triangle surrounded by three tangent circles? By the uniqueness part of the circle packing theorem, each packing is uniquely determined (up to symmetries of the surface) by its graph of tangencies, but for this problem, the graph of tangencies is determined from the geometry, rather than vice versa.[28]

Among periodic infinite triangulated packings of the Euclidean plane, there is only one triangulated packing with a single radius, the hexagonal packing. There are nine distinct pairs of radii (normalized so the maximum radius is one), and 164 triples of radii, that form triangulated packings. In general the number of -tuples of radii that can form a triangulated packing is finite, but with an unknown growth rate as a function of .[29] The ratio of the smallest to largest radius is at most 0.701 for any triangulated packing other than the hexagonal packing, for which this ratio is 1. There is a packing with radius ratio approximately 0.651, with three distinct radii; it has been conjectured by Connelly & Zhang (2025) to have the largest possible radius ratio of any non-hexagonal triangulated packing.[30]

A circular rectangle with Brooks invariant 3 + 1/(2 + 1/4) = 31/9

Brooks (1986) describes a method of continuously associating a real number as an invariant of a circular rectangle, the region bounded by a cycle of four tangent circles (which may be given the symmetries of a rectangle by a Möbius transformation), in such a way that when the invariant is a rational number, the rectangle has a triangulated packing. His method first packs a chain of circles tangent to the top and bottom sides of , with the leftmost circle in the chain tangent to the left side of ; let be the number of circles that can be packed in this way without overlapping the right side of . This chain leaves a smaller circular rectangle between its rightmost circle and the right side of the starting rectangle; pack this rectangle recursively in the same way with another chain of circles, tangent to its left and right sides, and let be the number of circles in this second chain. This packing process may continue forever, giving an infinite sequence of positive integers , the numbers of circles in each chain, or it may terminate after a finite number of stages with a triangulated packing. Regardless, Brooks defines his invariant to be the simple continued fraction This is a rational number if and only if this packing process terminates with a finite triangulated packing.[28][31]

Not every uniform surface admits a triangulated packing. For instance, in the case of flat tori (the quotient spaces of the Euclidean plane by a lattice, obtained by gluing opposite edges of a parallelogram), the geometry can vary continuously but, by the uniqueness of toroidal circle packings, only countably many distinct tori (up to scaling) have triangulated packings, and in particular they must have shapes that can be described by algebraic numbers.[32] Nevertheless, Brooks (1986) showed that, for each possible genus of compact surfaces, the surfaces that have triangulated packings form a dense set in the Teichmuller space of all surfaces with that genus. This means that, for any given uniform compact surface, there exists another uniform compact surface, with almost the same geometry, on which there exists a triangulated packing. To prove this, Brooks first considers packings in which all gaps between triangles have either three or four sides, which can be obtained by placing circles to break up gaps with more sides into pieces with fewer sides. He then deforms each circular rectangle obtained in this way into a nearby rectangle whose invariant is rational to obtain a triangulated packing within each of the deformed rectangles. This produces a triangulation of the whole surface which, by the circle packing theorem, can be represented by a circle packing on a different uniform surface, slightly deformed from the given surface.[28]

Applications

[edit]

The circle packing theorem is a useful tool to problems including conformal maps, polyhedral realization of graphs, the planar separator theorem, graph drawing, random walks on planar graphs, and the visualization of the human brain.

Conformal mapping

[edit]
Circle packings can be used to approximate conformal mappings between specified domains. Each circle on the left corresponds to a circle on the right.

A conformal map between two open sets in the plane or in a higher-dimensional space is a continuous function from one set to the other that preserves the angles between any two curves. The Riemann mapping theorem, formulated by Bernhard Riemann in 1851, states that, for any two open topological disks in the plane, there is a conformal map from one disk to the other. Conformal mappings have applications in mesh generation, map projection, and other areas. However, it is not always easy to construct a conformal mapping between two given domains in an explicit way.[33]

A construction of William Thurston uses circle packings to approximate conformal mappings.[34] More precisely, this construction maps an arbitrary open disk to the interior of a circle; the mapping from one topological disk to another disk could then be found by composing the map from to with the inverse of the map from to .[33] Thurston's idea was to pack circles of some small radius in a hexagonal tessellation of the plane, within region , leaving a narrow region near the boundary of , of width , where no more circles of this radius can fit. He then constructs a maximal planar graph from the intersection graph of the circles, together with one additional vertex adjacent to all the circles on the boundary of the packing. By the circle packing theorem, this planar graph can be represented by a circle packing within circle in which all the edges (including the ones incident to the boundary vertex) are represented by tangencies of circles. The circles from the packing of correspond one-for-one with the circles within , with the circle that forms the boundary of corresponding to the boundary of .[33]

This correspondence of circles can be used to construct a continuous function from to the interior of in which each circle and each gap between three circles is mapped from one packing to the other by a Möbius transformation. As goes to infinity, the function determined using Thurston's method from hexagonal packings of radius- circles converges uniformly on compact subsets of to a conformal map from to .[35][33] More generally, this holds for sequences of packings that are not necessarily hexagonal packings, as long as the sequence of maximum radii of the packings converges to zero.[36] An alternative method hexagonally packs circles of small radius into , and constructs a combinatorially equivalent packing of circles into , tangent to its boundary; it converges in the same way.[37] Practical applications of this method have been hindered by the difficulty of computing circle packings and by its relatively slow convergence rate.[38] However, it has some advantages when applied to non-simply-connected domains and in selecting initial approximations for numerical techniques that compute Schwarz–Christoffel mappings, a different technique for conformal mapping of polygonal domains.[33]

This construction of approximate conformal maps can also be described in terms of supplying a uniform geometry to a topological surface, and it can be extended from the surface itself to other objects embedded within the surface. An example comes from dessins d'enfant, a certain type of graph embedding used in algebraic geometry to study Riemann surfaces and the absolute Galois group. These automatically come with a triangulation, within which the dessin is a subset of the edges, but these triangulations may form multigraphs rather than simple graphs. By applying the circle packing theorem to subdivisions of these triangulations, an approximation of the dessin on a surface of uniform geometry may be constructed.[33][39]

Quasiconformal mapping

[edit]

When a conformal mapping is not possible, it may still may be possible to obtain a quasiconformal mapping. Conformal mappings are smooth mappings from one Riemannian surface to another that, in the limit of small neighborhoods around any point, map circles to circles. Instead, quasiconformal mappings may, in the same limit, map circles to ellipses of bounded eccentricity, with smaller deformation as the bound on eccentricity approaches one. As one application of this idea, Hurdal & Stephenson (2004) apply circle packing to the problem of constructing a flattened map of the functional regions of the human brain. To do this, they triangulate a three-dimensional model of the brain surface, find a circle packing representing the triangulation, and construct a piecewise-linear mapping from the triangles of the triangulated brain to the triangles of circle centers in the packing. They argue that this produces a quasiconformal mapping and therefore that it preserves the approximate shape of brain structures.[40]

Brooks (1986) uses this idea in a different way, to study Kleinian groups, discrete subgroups of the symmetries of three-dimensional hyperbolic space. A Kleinian group is geometrically finite if it has a fundamental domain that is the intersection of finitely many half-spaces. Brooks shows that any geometrically finite group can be transformed by a quasiconformal mapping of arbitrarily small distortion into a subgroup of a Kleinian group with a fundamental domain of finite volume. If, moreover, the given Kleinian group has no cusps (either points at infinity where two bounding planes of the fundamental domain meet tangentially, or isolated points at infinity) it can be transformed by a quasiconformal mapping of arbitrarily small distortion into a subgroup of a Kleinian group with a compact fundamental domain, without points at infinity. The proof involves finding triangulated packings of surfaces whose geometry approximates the Riemann surfaces at the ends of the fundamental domain (the connected components of its points at infinity). The circles of these packings generate another Kleinian group, within which the circles coming from the boundaries of the fundamental domain generate a subgroup, close to the originally given group.[28] Williams (2019) uses similar ideas to construct quasiconformal maps from a specified domain to the unit disk, having a specified distortion function.[31]

Random walks

[edit]

Several applications of the circle packing theorem use it to study random walks on planar graphs, based on the intuition that the approximation to conformal mapping of a circle packing can be used to relate these walks to Brownian motion, a geometric random process that is invariant under conformal maps.[41] One result in this direction is that unbiased limit graphs of bounded-degree planar rooted graphs are almost surely recurrent, meaning that random walks on these graph limits almost surely return to the root infinitely often. The proof involves finding a circle packing that represents the limit graph, and using the ring lemma to bound the sizes of circles a small number of steps from the start of the walk, allowing the walk to be compared to a random walk on an integer grid.[42] Jonnason & Schramm (2000) used circle packings to prove that the expected cover time of a random walk on every -vertex planar graph (the average number of steps of the walk before all vertices are visited) is .[43]

Instead of using methods from Brownian motion to understand random walks on graphs, circle packings can also be used to find graphs whose random walks approximate Brownian motion. For the same sequences of refined circle packings used to approximate conformal mappings of a domain, a random walk on the associated coin graph will reach a boundary vertex with a probability that approaches, in the limit as the circle packings become more fine, the probability that a Brownian motion will reach a nearby point on the boundary of the domain.[2]

Polyhedral realization and planar separators

[edit]
A polyhedron and its midsphere. The red circles on the sphere (its horizons as viewed from the polyhedron vertices), and the blue domes where the midsphere protrudes through the polyhedron, form orthogonal spherical circle packings of the graph of the polyhedron and its dual graph.

The primal–dual packing of any polyhedral graph can be lifted from the plane to a sphere by stereographic projection. It can then be used to construct a convex polyhedron that has the given graph as its vertices and edges and that has a midsphere, a sphere tangent to all of the edges of the polyhedron. Each vertex of the polyhedron lies outside the sphere, at the apex of a cone that is tangent to the sphere along the circle corresponding to that vertex. This circle forms the horizon on the sphere, as viewed from the vertex. Each face of the polyhedron lies on the plane through its corresponding circle. Each edge of the polyhedron passes from one cone apex to another through the point of tangency of two circles, where both cones are tangent to the sphere. Conversely, if a polyhedron has a midsphere, then the circles formed by the intersections of the sphere with the polyhedron faces and the circles formed by the horizons on the sphere as viewed from each polyhedron vertex form a dual packing of this type.[44]

Möbius transformations of the sphere lead to geometrically different polyhedron realizations. One particular realization, called the canonical polyhedron, is obtained by using a unit sphere and choosing a Möbius transformation that causes the center of the sphere to coincide with the centroid of the points of tangency of the polyhedron edges. All combinatorially equivalent polyhedra will produce in this way the same canonical polyhedron, up to rotations of the sphere.[45] An alternative proof of the planar separator theorem, originally proved in a different way by Lipton and Tarjan,[46] uses a similar idea of applying a Möbius transformation to a spherical circle packing. The construction obtains a circle packing on a sphere, representing the given graph. Next, a Möbius transformation of the sphere is used to transform a centerpoint of the circle centers, a point in space such that every plane through the centerpoint divides the centers into two subsets that each have size , to coincide with the center of the sphere. A uniformly random plane through the center of the sphere intersects of the circles of any packing of circles, in expectation. For the circle packing with the centerpoint at the sphere center, the vertices corresponding to these intersected circles form a separator that, when removed from the graph, leaves connected subgraphs of at most vertices.[47] Circle packing has also been used to explain the success of spectral methods for constructing separators for graphs of bounded genus and bounded degree.[48]

Graph drawing

[edit]

Circle packing has been used as a frequent tool in graph drawing, the study of methods for visualizing graphs. Fáry's theorem, that every graph that can be drawn without crossings in the plane using curved edges can also be drawn without crossings using straight line segment edges, follows as a simple corollary of the circle packing theorem: by placing vertices at the centers of the circles and drawing straight edges between them, a straight-line planar embedding is obtained.[49]

Malitz & Papakostas (1994) apply circle packing to angular resolution, a measure of the quality of a graph drawing defined by the sharpest angle any two edges form at a shared vertex. They show that planar graphs of bounded maximum degree have drawings of bounded angular resolution, using circle packing. To do this, they derive a version of the ring lemma and use it to prove that the angular resolution of the drawing in which each vertex is placed at the center of its circle is at least an inverse-exponential function of the degree.[49] Keszegh, Pach & Pálvölgyi (2013) use a more complicated strategy for placing vertices, near the circle centers on subdivisions of integer lattices, in order to find drawings where the number of distinct edge slopes (the slope number) is at most an exponential function of the degree.[50]

Primal–dual circle packings in the plane can be used to obtain simultaneous orthogonal drawings, straight-line planar drawings of any polyhedral graph and its dual graph, with one dual vertex at infinity, so that each primal and dual pair of edges cross at right angles. Such a drawing can be obtained from a primal–dual circle packing for which one dual circle surrounds the others, with all vertices (except the outer dual vertex) placed at the center of the corresponding circle. The use of circle packings to construct these drawings resolves a problem posed by W. T. Tutte in 1963.[51] A pointed drawing of a given planar graph, in which the edges are drawn as smooth curves that, at each vertex, leave it in the same direction as each other, may be obtained by placing each vertex on the circle of a circle packing (rather than at its center) and drawing each edge as a biarc formed by two circular arcs, perpendicular to the circles of the packing, through the point of tangency of two circles.[52] Circle packings have also been used to obtain strongly monotone drawings of polyhedral graphs. These are straight-line drawings in which each two vertices are connected by a polygonal path that is monotone with respect to the line segment between the two vertices (its orthogonal projection onto this line segment is one-to-one).[53] Another construction involving circle packing produces a Lombardi drawing for any subcubic planar graph, a drawing resembling a two-dimensional soap bubble foam in which the edges are drawn as circular arcs that meet at equal angles at each vertex.[54]

Algorithmic aspects

[edit]

When a circle packing exists, it exists with algebraic numbers for its center and radii, solving a system of quadratic equations that represent the requirements that each two centers of tangent circles be at a distance from each other equal to the sum of their radii.[32] Therefore, in principle, it is possible to find circle packings using methods from symbolic computation for solving systems of polynomial equations, such as the Gröbner basis method. However, this approach can require the solution of high-degree polynomial equations, with unsolvable Galois groups. Instead, researchers have focused on more efficient numerical methods for finding accurate approximations to the centers and radii of circle packings.[55]

One of the earliest and most versatile of these approximation results, announced by Bojan Mohar in 1993 with details in two subsequent papers, can be used to find primal–dual circle packings on surfaces of arbitrary genus, in time polynomial in the number of circles and in the number of bits of precision used to specify the output to the desired accuracy. Mohar's algorithm is an iterative method applied to the system of radii of a primal–dual circle packing, starting from an initial inaccurate system of radii. It considers the orthodiagonal quadrilaterals formed by each pair of crossing primal and dual edges. From a system of radii, one can calculate the angles of these quadrilaterals (if the circles could be packed with these radii): for a system of radii to form a packing, these angles should sum to around each vertex. For systems of radii that do not give the correct sums, Mohar partitions the radii into a subset of radii that are too large (the angles around that vertex sum to less than and a complementary subset of radii that are too small, and uses the bisection method to find two parameters and such that multiplying the large radii by and the small radii by reduces the mean squared error between the angle sums and . After sufficiently many repetitions of this process lead to an accurate system of radii, Mohar uses these radii to place the circle centers, starting from one primal–dual quadrilateral and then using the already-placed centers of each quadrilateral and the radii of its circles to determine the locations of the remaining centers. Although Mohar states that the algorithm runs in polynomial time, he does not state its polynomial time bound explicitly.[56] Dong, Lee & Quanrud (2020) suggest that the computational complexity of Mohar's packing algorithm is at least quintic in the number of circles to be packed.[57]

Collins & Stephenson (2003) provide a simpler iterative method that modifies only one circle radius at a time, again computing the system of circle radii before using the radii to place the circle centers. The version of the circle packing problem that they solve produces a circle packing in the Euclidean plane for a maximal planar graph, with a specified outer triangle and with specified radii for the three circles of this triangle. They begin with a set of tentative radii that do not correspond to a valid packing, and then repeatedly perform the following steps:

  1. Choose an internal vertex of the input graph.
  2. Calculate the total angle that its neighboring circles would cover around the circle for , if the neighbors were placed tangent to each other and to the central circle using their tentative radii.
  3. Determine a representative radius for the neighboring circles, such that circles of radius would give the same covering angle as the neighbors of give.
  4. Set the new radius for to be the value for which circles of radius would give a covering angle of exactly .

Each of these steps may be performed with simple trigonometric calculations. Collins and Stephenson prove that, applied only to a single vertex, this method rapidly converges to a radius that gives covering angle . Further, they argue that the system also globally converges to a consistent system of radii, and they provide computational experiments suggesting that the rapid local convergence of this method extends to rapid global convergence of the entire system of radii. Once the system has converged, the circles may be placed one at a time, similarly to the algorithm of Mohar, at each step using the positions and radii of two neighboring circles to determine the center of each successive circle.[58]

Another algorithm by Orick, Stephenson & Collins (2017) applies to the same planar packings as Collins & Stephenson (2003), and approximates both the system of radii and circle centers simultaneously and iteratively. It alternates between steps that use a weighted Tutte embedding to place the centers, and that use the resulting circle centers to determine improved radii. Although they provide no theoretical guarantees on the performance of this algorithm, Orick et al. claim that their experiments show linear convergence to a consistent packing in practice, which when true would lead to a polynomial time bound.[59]

Dong, Lee & Quanrud (2020) again address the case of finding planar packings for maximal planar graphs, but with an explicit time bound. Although focused only on a primal circle packing, they work with the system of radii of a primal–dual packing, for maximal planar graphs only, with a dual circle for each triangular face of the graph. Like Mohar, they use orthodiagonal quadrilaterals determined by crossing pairs of primal and dual edges, with angles calculated from a given system of radii, which for a valid system of radii should sum to at each vertex. Following Colin de Verdière (1991), they set up a convex optimization problem in which the variables to be optimized are the logarithms of the radii of the primal–dual packing, and the convex function to be minimized combines terms that are linear in these variables and in the antiderivatives of the arctangents of ratios of radii. They show that the minimizer of this function provides a system of radii that meet the conditions of a circle packing in that the angles of the quadrilaterals around each vertex sum to . By finding a spanning tree of the graph of orthogonal circles in the packing with the property that each spanning tree edge connects circles whose radii are not too far from each other, they prove that this function is strongly convex near its minimum, allowing the application of powerful methods from convex programming. Based on these techniques, they give a time bound for finding a packing of where the is a form of big O notation that suppresses terms that are logarithmic in its argument, is the number of circles to be packed, is the ratio between the largest and smallest radius in the resulting packing (which may be at worst exponential in ), and is a parameter that controls how accurate the numerical approximation to the packing should be. All computed radii are within a factor of the radius of an exact circle packing, and all computed centers are at a distance from the centers of an exact circle packing that is at most times the radius of the smallest circle.[57]

Yu et al. (2011) develop a numerical distributed algorithm based on a discrete version of Ricci flow to construct approximate primal–dual circle packings on a sphere. Their motivation is to use the midsphere polyhedra constructed from these packings for load-balanced position-based routing, following a scheme of Papadimitriou & Ratajczak (2005).[60] According to this scheme, on a polyhedron with a midsphere, messages can be forwarded from vertex to a neighboring vertex of this polyhedron in a way that always increases the inner product of the three-dimensional coordinates of the current location and destination. The existence of a midsphere is used to prove that each message eventually reaches its destination.[61]

History

[edit]
An Apollonian gasket, mentioned by Leibniz in 1706
1911 image of a Doyle spiral from an article on plant growth

Ken Stephenson writes that "there is indeed a long tradition behind our story", and attributes the familiar hexagonal packing of congruent circles to folklore.[62] Other early works on circle packings with specific patterns of tangency include:

The circle packing theorem itself was first stated and proved by Paul Koebe in 1936.[67] In the 1970s, William Thurston rediscovered the circle packing theorem,[62] and noted that it followed from the work of E. M. Andreev.[68] He mentioned these results in a 1978 talk at the International Congress of Mathematicians; in the same year, German mathematician Gerd Wegner also conjectured that the circle packing theorem holds.[69] At the Bieberbach conference in 1985, Thurston proposed a scheme for using the circle packing theorem to obtain a homeomorphism of a simply connected proper subset of the plane onto the interior of the unit disk, and he conjectured that this scheme converges to the Riemann mapping as the radii of the circles tend to zero.[34] Thurston's conjecture was proven in 1987 by Burton Rodin and Dennis Sullivan.[35] Thurston's work in this area brought the topic of circle packings to the attention of mathematicians and led to the formation of a research community centered on this idea.[62]

See also

[edit]

Notes

[edit]
  1. ^ a b Brightwell & Scheinerman (1993).
  2. ^ a b Dubejko (1997).
  3. ^ a b de Fraysseix, Ossona de Mendez & Rosenstiehl (1994).
  4. ^ Pisanski & Randić (2000).
  5. ^ a b c Thurston (2002), p. 330, Corollary 13.6.2.
  6. ^ Nachmias (2020).
  7. ^ Thurston (2002, p. 331) and Nachmias (2020) justify this claim by adding one vertex per face, adjacent to all vertices in the face. This is not a correct construction: It fails when a vertex has more than one incidence with a face, as happens in plane graphs that are not 2-vertex-connected. de Fraysseix, Ossona de Mendez & Rosenstiehl (1994, p. 241, Theorem 3.2) claim that two added vertices per face suffice, but do not provide details. A more complicated augmentation strategy that does not require 2-connectivity, adding a binary tree within each face and then triangulating the new faces created by this tree, is given by Alam et al. (2014).
  8. ^ Koebe (1936).
  9. ^ He & Schramm (1993).
  10. ^ Rohde (2011), p. 1628.
  11. ^ Beardon & Stephenson (1991); Carter & Rodin (1992).
  12. ^ Colin de Verdière (1991).
  13. ^ Chow & Luo (2003).
  14. ^ Connelly & Gortler (2021).
  15. ^ Thurston (2002), pp. 331–332.
  16. ^ Schramm (1991).
  17. ^ He & Schramm (1995).
  18. ^ Beardon & Stephenson (1990).
  19. ^ Coxeter (2006).
  20. ^ Schramm (1992); Brightwell & Scheinerman (1993); Sachs (1994)
  21. ^ Bowers & Stephenson (2004), pp. 78–82, Section 8.2 Inversive distance packings.
  22. ^ Schramm (1990).
  23. ^ Rohde (2011), pp. 1627–1628.
  24. ^ Rodin & Sullivan (1987, p. 352); Hansen (1988); Aharonov (1997)
  25. ^ Stephenson (2005), pp. 63–64, Lemma 6.2 (Monotonicity in Triangles).
  26. ^ Beardon & Stephenson (1991).
  27. ^ Nederlof, Pilipczuk & Węgrzycki (2023).
  28. ^ a b c d Brooks (1986); Kapovich (2001)
  29. ^ Fernique (2025).
  30. ^ Connelly & Zhang (2025).
  31. ^ a b Williams (2019).
  32. ^ a b Louder, Mishchenko & Souto (2014).
  33. ^ a b c d e f Stephenson (1999).
  34. ^ a b Thurston (1985).
  35. ^ a b Rodin & Sullivan (1987).
  36. ^ He & Schramm (1996).
  37. ^ Carter & Rodin (1992).
  38. ^ Stephenson (1999): "Circle packing certainly cannot compete with the classical numerical methods for speed and accuracy."
  39. ^ Bowers & Stephenson (2004).
  40. ^ Hurdal & Stephenson (2004).
  41. ^ Nachmias (2020), p. 7.
  42. ^ Benjamini & Schramm (2001).
  43. ^ Jonnason & Schramm (2000).
  44. ^ Schramm (1992); Sachs (1994). Schramm states that the existence of an equivalent polyhedron with a midsphere was claimed by Koebe (1936), but that Koebe only proved this result for polyhedra with triangular faces. Schramm credits the full result to William Thurston, but the relevant portion of Thurston's lecture notes (Thurston 2002, pp. 331–332) again only states the result explicitly for triangulated polyhedra.
  45. ^ Ziegler (1995).
  46. ^ Lipton & Tarjan (1979).
  47. ^ Miller et al. (1997).
  48. ^ Kelner (2006).
  49. ^ a b Malitz & Papakostas (1994).
  50. ^ Keszegh, Pach & Pálvölgyi (2013).
  51. ^ Tutte (1963); Brightwell & Scheinerman (1993); Felsner & Rote (2019)
  52. ^ Aichholzer et al. (2012).
  53. ^ Felsner et al. (2016).
  54. ^ Eppstein (2014).
  55. ^ Bannister et al. (2015).
  56. ^ Mohar (1993); Mohar (1997a); Mohar (1997b)
  57. ^ a b Dong, Lee & Quanrud (2020).
  58. ^ Collins & Stephenson (2003).
  59. ^ Orick, Stephenson & Collins (2017).
  60. ^ Yu et al. (2011).
  61. ^ Papadimitriou & Ratajczak (2005).
  62. ^ a b c Stephenson (2003).
  63. ^ Bos (2010).
  64. ^ Michiwaki (2008).
  65. ^ Leibniz (1706).
  66. ^ Emch (1910).
  67. ^ Koebe (1936).
  68. ^ Andreev (1970a, 1970b)
  69. ^ Sachs (1994).

References

[edit]
[edit]
  • CirclePack (free software for constructing circle packings from graphs, by Kenneth Stephenson, Univ. of Tennessee)
  • CirclePackings, open source software for constructing circle packings from graphs, by Benjamin and Brice Loustau