Wearie-Phelan Bubbles

Updated 19 Mar 2021

The Problem

How would bubbles pack together, to give the least possible amount of surface film between them? This question has not been answered yet.

The question is actually a little more precise than that, as the obvious solution is one giant spherical bubble, but that is clearly not what we have in mind. What about bubbles of a given size? Well, why not allow bubbles of different sizes if that turns out to be more efficient? Really, what we mean is bubbles at a given density, where the bubble density is the number of bubbles per unit volume.

In the nineteenth century Lord Kelvin considered the simplest case of identically sized and shaped bubbles. He discovered a packing of 14-sided bubbles or tetrakaidecahedra, based on the truncated octahedron but with slightly curved faces. The faces and edges have to curve a little so that they meet at the correct angles to minimise the surface energy, as described by Plateau's laws. These angles are always the same for all bubbles; the bubbles always meet three to an edge and four to a vertex. The angle of a bubble at a vertex is known as the tetrahdral angle. The faces are minimal surfaces stretching between the edges.

the truncated octahedron
Truncated octahedron
A group of Kelvin bubbles
A group of Kelvin bubbles

Kelvin conjectured that this was the most efficient solution, but he could neither prove nor disprove his conjecture.

The Bubbles

A hundred years went by before Denis Weaire and Robert Phelan disproved the Kelvin conjecture by finding a better solution. They asked the slightly less strict question, how would equal sized, but not necessarily identically shaped, bubbles pack? The best answer they could find was a mixture of 12- and 14-sided bubbles, also having slightly curved surfaces.

In the Weaire-Phelan packing, two irregular pentagonal dodecahedra (12-sided) and six tetrakaidecahedra (of a rather different shape from Kelvin's, which also has 14 faces) form a translation unit with a lattice periodicity which is simple cubic. The illustration shows the dodecahedra as wire frames and the tetrakaidecahedra as solid. The dodecahedra do not touch each other, but are entirely surrounded by tetrakaidecahedra.

the w-p polyhedra
   Dodecahedron          Tetrakaidecahedron
A translation unit
Translation unit

The dodecahedron has pentagonal faces with sides of unequal length, so it is "pentagonal" to distinguish it from the rhombic dodecahedron and "irregular" to distinguish it from the regular (Platonic) dodecahedron.

Another way to visualise the packing is to note that the hexagonal faces of the tetrakaidecahedra are truly flat and meet to stack them in long rods. A set of rods lies parallel to each of the three orthogonal axes of space. These three sets of rods interlace, touching on their eight smaller pentagonal faces but leaving dodecahedral voids between the larger pentagons.

More about bubble packings can be found in [1].

Since then no better solution has been found, though at least one has come very close. Actual Weaire-Phelan bubbles also resisted experimental demonstration for more than a decade, but a foam was eventually made. Proof of the optimal solution remains as elusive as ever, we still do not know if this is it.

The Polyhedra

We can flatten out all the bubble faces to form conventional polyhedra. These faces lie orthogonal to the lines joining adjacent cell centroids, and to retain equal cell volumes some vertices must move very slightly. The geometry of the resulting spacefilling is in all other respects identical to Weaire and Phelan's: the differences are, in any case, barely noticeable.

The tetrakaidecahedron may now be seen an example of the Goldberg polyhedron, which helps to distinguish it from the one associated with Kelvin, the truncated octahedron. The dodecahedron has the same symmetry seen in crystals of iron pyrites and so is sometimes called a pyritohedron, although its angles are different.

Vertex coordinates and printable nets for card models of the two polyhedra are given below. The coordinate scale gives a lattice period, i.e. the size of the repeating pattern, of 20.


    x         y         z
 3.14980   0         6.29961
-3.14980   0         6.29961
 4.19974   4.19974   4.19974
 0         6.29961   3.14980
-4.19974   4.19974   4.19974
-4.19974  -4.19974   4.19974
 0        -6.29961   3.14980
 4.19974  -4.19974   4.19974
 6.29961   3.14980   0
-6.29961   3.14980   0
-6.29961  -3.14980   0
 6.29961  -3.14980   0
 4.19974   4.19974  -4.19974
 0         6.29961  -3.14980
-4.19974   4.19974  -4.19974
-4.19974  -4.19974  -4.19974
 0        -6.29961  -3.14980
 4.19974  -4.19974  -4.19974
 3.14980   0        -6.29961
-3.14980   0        -6.29961

(Goldberg polyhedron)

    x         y         z
 3.14980   3.70039   5
-3.14980   3.70039   5
-5         0         5
-3.14980  -3.70039   5
 3.14980  -3.70039   5
 5         0         5
 4.19974   5.80026   0.80026
-4.19974   5.80026   0.80026
-6.85020   0         1.29961
-4.19974  -5.80026   0.80026
 4.19974  -5.80026   0.80026
 6.85020   0         1.29961
 5.80026   4.19974  -0.80026
 0         6.85020  -1.29961
-5.80026   4.19974  -0.80026
-5.80026  -4.19974  -0.80026
 0        -6.85020  -1.29961
 5.80026  -4.19974  -0.80026
 3.70039   3.14980  -5
 0         5        -5
-3.70039   3.14980  -5
-3.70039  -3.14980  -5
 0        -5        -5
 3.70039  -3.14980  -5

The Spacefilling

The polyhedral spacefilling is related to what crystallographers call a Frank-Kasper phase. These phases are characterised by the arrangement of atoms (here represented by the cell centroids) at the vertices of planar arrays of triangles and hexagons.

If you place a small sphere at each atom in the array and expand them all together until they meet to form a boundary, these boundaries will be flat faces. If you keep expanding the spheres until no intervening space is left, the flat faces all join up and the spheres become polyhedra known as Voronoi cells. The Weaire-Phelan polyhedra are the Voronoi cells of an A15 Frank-Kasper phase.

spheres meeting
Equally expanding spheres meet in a plane surface (shaded).

A natural example may be found in certain gas hydrates, in which the gas molecules lie at the centres of the cells while the oxygen atoms of water form the vertices of the lattice and the hydrogen atoms lie along the edges. It is described by Williams as the 12 Å cubic gas hydrate, where the angstrom Å is a unit of length equal to 10−10 m.[2]

Since it has flat faces, the polyheral spacefilling is better suited to making card models than the actual bubbles are. Downloadable nets are provided above.

The polyhedra are also easier to design for 3D printing than the bubbles! Tip: it can be worth making them hollow with a fairly thin wall, so as to use less material. This saves cost and time, especially if you are printing off large numbers, and also makes them lighter.


  1. Weaire, D.; "Froths, Foams and Heady Geometry", New Scientist, 21 May 1994.
  2. Robert Williams; The Geometrical Foundation of Natural Structure, 2nd Edn, Dover, 1979, pp.186-9. (1st Edn. published as Natural Structure, Eudaemon, 1972. Further editions have been published under other titles.)