The hard disk model is one of the oldest problems in statistical mechanics: identical circles on a plane, interacting only by not overlapping. The partition function requires computing the free volume available to each disk — how much space it can explore without colliding with its neighbors. This free volume has resisted analytical treatment for decades because the excluded regions around neighboring disks create irregularly shaped voids.
Pergamenshchik, Bryk, and Trokhymchuk found the decomposition. The free volume can be expressed exactly in terms of the intersection areas of up to five exclusion circles — geometric quantities that are analytically computable. Each exclusion circle marks the region forbidden to a disk's center by one neighbor. The intersection of two circles captures pairwise crowding. Three circles capture three-body correlations. The five-circle intersection turns out to be a natural measure of local hexagonal order, because hexagonal packing is the arrangement that maximizes the simultaneous overlap of surrounding exclusion zones.
With this decomposition, the partition function factors into a product of free volumes. The formula has two exact limits: a gas-like state where free volumes are large and independent, and a liquid-like state where they're tightly coupled by neighbor correlations. Between these limits, the theory recovers the known equation of state across nearly the entire density range.
The intermediate regime reveals something new: a mixed liquid phase where local defects in hexagonal order proliferate before the system commits to crystalline ordering. The transition to the hexatic phase isn't sharp — it's preceded by a density range where hexagonal and disordered patches coexist, and the five-circle overlap area distinguishes them quantitatively.
Sixty years of intractable geometry, resolved by counting the right intersections.