advanced34 min studyLesson 618 in Mathematics

Sphere packing and covering bounds

Compute lattice packing and covering densities while separating nonoverlap, coverage, and global optimality claims.

On this page

What you’ll understand

  • Translate congruent-ball packing into a separation condition on centers.
  • Compute lattice densities using the fundamental-cell volume.
  • Distinguish a shortest-vector packing radius from a worst-hole covering radius.

Packing asks how densely equal balls can be placed without overlapping interiors. Covering asks how small balls can be while reaching every point of space, allowing overlaps. The same set of centers usually needs different radii for these two tasks. A lattice makes the density calculation exact because every fundamental cell repeats, but proving that one lattice is best among all arrangements is a separate theorem. We will derive concrete planar values before making any claim about optimality.

Specify the ambient metric and center model

Use Euclidean balls of one fixed radius r in real d-space. Their interiors are disjoint exactly when every pair of distinct centers is at least 2r apart. Tangency is allowed, and boundaries have zero d-dimensional volume. For a full-rank lattice Λ = BZᵈ, the volume of a fundamental parallelepiped is the absolute determinant of B, independent of the particular lattice basis chosen. There is one center per fundamental cell.

Write κ_d for the volume of the unit Euclidean ball. The density of a lattice packing is therefore the volume κ_d rᵈ of one radius-r ball divided by the fundamental-cell volume. It cannot exceed one because packing interiors do not overlap. This volume bound is elementary and often very weak. It concerns a periodic infinite arrangement, not the fraction of a finite plotting window filled by clipped disks. For arbitrary nonperiodic arrangements, density requires an asymptotic definition; we do not replace that analysis with a screenshot.

δ(Λ,r)=κdrd∣det⁡B∣,r≤12λ1(Λ),λ1(Λ)=min⁡v∈Λ∖{0}∥v∥2\delta(\Lambda,r)=\frac{\kappa_dr^d}{|\det B|},\qquad r\le\frac12\lambda_1(\Lambda),\quad \lambda_1(\Lambda)=\min_{v\in\Lambda\setminus\{0\}}\|v\|_2
The lattice is full rank in Euclidean d-space, all balls have the same radius, and the separation condition guarantees a packing.↔ Scroll wide equations horizontally.

Covering measures the largest remaining hole

The covering radius μ(Λ) is the supremum, over all points in space, of distance to the nearest lattice point. Balls of radius R centered at the lattice cover the space exactly when R is at least μ. By periodicity it suffices to inspect one fundamental cell, although the nearest center may lie in a neighboring cell. For a rectangular planar lattice, the farthest points are rectangle centers, giving half the rectangle diagonal.

The covering density counts ball volume per cell and can exceed one because it counts overlaps with multiplicity. It must be at least one for a full covering; average multiplicity below one would leave uncovered volume. The shortest vector only detects when neighboring balls touch. It does not identify the deepest holes, so the packing radius and covering radius must not be interchanged. A stretched rectangular lattice can have a modest shortest-vector length and a very large covering radius.

μ(Λ)=sup⁡x∈Rdinf⁡v∈Λ∥x−v∥2,Θ(Λ,R)=κdRd∣det⁡B∣≥1(R≥μ(Λ))\mu(\Lambda)=\sup_{x\in\mathbb R^d}\inf_{v\in\Lambda}\|x-v\|_2,\qquad \Theta(\Lambda,R)=\frac{\kappa_dR^d}{|\det B|}\ge1\quad(R\ge\mu(\Lambda))
The covering-density lower bound follows from coverage and periodic average multiplicity. It is not the same optimization problem as maximizing packing density.↔ Scroll wide equations horizontally.

Compare constructions without overclaiming optimality

In the plane, a square lattice with nearest-center spacing two packs unit disks at density π/4. A triangular lattice with basis (2, 0) and (1, √3) has the same minimum spacing but cell area 2√3, giving density π/(2√3), approximately 0.9069. This comparison proves the triangular construction is denser than that square construction. The fact that this planar density is globally optimal is a separate classical theorem, not proved merely by comparing two examples.

Higher dimensions introduce increasingly complex lattices and nonlattice periodic arrangements. The determinant controls volume per center, but a short vector can violate separation even when the determinant looks favorable. Covering has its own competing arrangements and bounds. In an engineering application, first decide whether the task is collision avoidance, guaranteed spatial service, or another metric entirely. Changing the norm, allowing unequal radii, or constraining centers to a finite container defines a different problem and changes the relevant bounds.

One square lattice, two radii

For the planar lattice with basis (2, 0) and (0, 2), find the maximum packing radius and the minimum covering radius, with their respective densities.

  1. The shortest nonzero lattice vectors have length two, so the largest nonoverlapping radius is one. The fundamental square has area four, giving packing density π/4.

  2. Every point in a fundamental square is within distance √2 of a corner, and the square center attains that distance to its nearest lattice centers. Hence the covering radius is √2.

  3. At that radius the disk area is 2π, so covering density is 2π/4 = π/2. Its value exceeds one because the covering disks overlap; this does not contradict the packing-density bound.

The result

The packing radius is 1 with density π/4; the covering radius is √2 with density π/2.

Explore the model

Interactive lab

Interactive lab · matrix

A fundamental cell exposes two different distances

Pause and Reset, then set diagonal scale a = 2 and shear b = 0. Predict the final rectangular cell is 2 by 1, with shortest lattice vector length one. Its packing radius is 0.5 and its covering radius is √5/2, obtained from the cell center. Play to inspect the final shape, remembering that readouts describe the selected matrix. The lab draws a cell, not balls or a search for optimal lattices.

A fundamental cell exposes two different distancesThe selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Current readings: Area scale · det A: 1.4; Singular values · σmax, σmin: 1.597, 0.876; Condition number · κ₂(A): 1.823. All parameters are adjustable with the labeled controls below.Matrix A with first row 1.4, 0.6 and second row zero, one; display interpolates between identity and A while final determinant remains 1.4.−4−2024−2−1012x · equal geometric scaleyA = [ 1.4 0.6 ; 0 1 ]Display blend: 0% · final values below
Interpolated geometryOriginal unit squareSelected final transform
Area scale · det A
1.4
Singular values · σmax, σmin
1.597, 0.876
Condition number · κ₂(A)
1.823
t = 0 s

The selected matrix is A = [[a,b],[0,1]]. Playback interpolates geometry between the identity and A; readouts always describe A. Its determinant is a and its singular values come from AᵀA. Positive a makes every supported matrix invertible. This two-dimensional family is not a general eigensolver or a learned network. Parameter changes pause playback and recompute the model at the current elapsed time. Reset restores both parameters and returns time to zero. Playback starts paused. Adjust either parameter, then choose Play or Step.

Look forNearest-neighbor separation and the farthest uncovered location are different geometric measurements.

Check your understanding

What does a lattice covering density greater than one mean?

Try it yourself

For the rectangular lattice with basis (4, 0) and (0, 2), find its packing radius, covering radius, and both densities. Compare with the same lattice scaled by one-half.

Take a moment to make a prediction or work it out on paper. Open a hint when you need one.

Give me a hint

Use the shorter side for the packing radius, half the diagonal for the covering radius, and the cell area for each density.

Compare with the solution

The cell area is eight and its shortest nonzero vector has length two, so packing radius is one and packing density is π/8. The farthest cell-center distance is √(2² + 1²) = √5, giving covering density 5π/8. Scaling the entire lattice by one-half gives cell area two, packing radius one-half, and covering radius √5/2. Both numerator ball areas and denominator cell area scale by one-quarter, so the two densities remain π/8 and 5π/8.

Take it with you

  • Packing imposes minimum center separation; covering bounds every point’s distance to a center.
  • Lattice density divides one ball’s volume by a fundamental-cell volume.
  • Comparing constructions does not by itself establish global optimality.

Read further

Original sources and deeper explanations to follow at your own pace.

Optional progress, saved in this browser when available. Revisit any lesson at any time.

Lesson not marked complete.