Packing and tiling

The second moment is the measurement

The mean number of sides in a cellular tissue is six, and Euler's formula leaves it no choice — so it takes the same value on a golden-angle head, a whorled head and a set of random points. The mean squared departure from six varies by a factor of eighty across the same three, and almost nobody reports it.

Worth reading first: Why the average cell has six sides · Packing, measured four ways · Lewis's law wants disorder.

An earlier essay established why the average cell in a planar tiling has six sides: Euler’s formula for a planar graph, counted two ways, forces it as the tiling grows, and no preference of the cells enters. It is a theorem about graphs, not a fact about tissue, and the honeycomb efficiency argument is a different claim about a different thing.

That essay stopped there, which was the right place to stop for the question it was asking. This one asks the next question: if the mean is forced, what should a measurement of a tissue report instead?

Six arrangements, one mean

Build six sets of nine hundred points. Two irrational divergence angles that produce spiral lattices, two rational angles a hair away from them, one rational angle that produces a whorled pattern, and one set with no rule in it at all. Tessellate each, cut the rim off, and count sides.

The mean number of sides per cell:

arrangement mean sides
whorled, 144° 5.99
golden, 137.508° 5.99
Lucas, 99.502° 6.04
rational, 137.5° 5.99
137.0° 5.99
Poisson 5.97

They agree to about one per cent. A set of points with no structure whatever averages the same number of sides as a lattice built to five decimal places, because both are planar tilings and a planar tiling cannot average anything else.

The statistic everybody reports is the one that cannot varySix arrangements of 900 points, from a whorled lattice to a set with no rule in it. The mean number of sides per cell is 5.97–6.04 on all six, because Euler's formula forces it. The mean squared departure from six runs from 0.023 to 1.83 — a factor of 79 — and the most hexagonal tissue in the set is the whorled one, at a rational angle.mean sides per cell(forced to six)mean squared departure from six(not forced)whorled, 144°5.9860.023golden, 137.508°5.9900.253rational, 137.5°5.9900.255Lucas, 99.502°6.0360.255137.0°5.9900.291Poisson5.9691.830six433–637 interior cells each, inside 86% of the radiussame cells, same cut, two statistics
Fig. 1 The two statistics side by side on the same six arrangements. The left column is the mean, which is flat because a theorem holds it there. The right is the mean squared departure from six, on the same cells with the same rim cut.

The same six, one statistic that varies

Now measure μ2\mu_2, the mean squared departure from six — the standard disorder measure for a cellular tissue, and the quantity Aboav–Weaire is a statement about:

arrangement μ2\mu_2 hexagons
whorled, 144° 0.023 99%
golden, 137.508° 0.253 75%
Lucas, 99.502° 0.255 75%
rational, 137.5° 0.255 75%
137.0° 0.291 71%
Poisson 1.830 27%

A factor of eighty between the ends, on arrangements whose means agree to a per cent.

That is the essay in one comparison. The statistic a tissue paper reports is the one the topology has already decided; the statistic that distinguishes a lattice from a random scatter of points is one step further into the same distribution, and it costs nothing extra to compute because it is the same side counts.

Why the mean cannot help, said precisely

It is worth being exact about what “forced” means, because it is not that the mean happens to be six.

The derivation is the earlier essay’s and is not repeated here. What matters for this one is its shape: counting edges two ways in a planar graph where three cells meet at each vertex gives a mean of 612/F6 - 12/F over FF cells, which goes to six as the tiling grows and is six to a per cent by a few hundred cells.

Nothing in that derivation mentions where the points are. Move them anywhere — cluster them, spread them, put them on a lattice — and as long as three cells meet at each vertex the arithmetic is unchanged. The mean is a property of the plane and the number of cells.

What the arrangement controls is how the side counts are distributed about that mean, and μ2\mu_2 is the first number that sees it.

How many sides the cells actually haveThe mean is 5.908, which Euler's formula forces. The spread around it is not noise — a tiling of only hexagons cannot close up on a finite patch.5 sides13420%6 sides45869%7 sides7311%665 bounded cellsmean 5.908 sidessix is forced, not chosen
Fig. 2 The distributions the two statistics are computed from. Every one of these has the same mean; what differs is how much of the weight sits away from six, which is the quantity the mean is constructed to be blind to.

Why the departure from six rather than the variance

There are two second moments available and they are not quite the same, so the choice is worth a line.

The variance is the mean squared departure from the sample mean. μ2\mu_2 is the mean squared departure from six. They differ by the square of the distance between the sample mean and six — which the theorem says is small, and which the table above confirms is at most 0.04, so the squared difference is under two thousandths.

On the numbers here the two are interchangeable to three decimal places, and the table would read the same either way. The reason to prefer μ2\mu_2 is not numerical: six is the value the topology forces, so a departure from it is the physical quantity, and a departure from the sample mean is a departure from an estimate of a quantity already known exactly. Estimating something a theorem supplies, and then measuring deviations from that estimate, throws away a degree of freedom for nothing.

It also matters at the edges. On a small or heavily rim-cut sample the mean can sit appreciably away from six, and in that case the difference between the two moments is a measure of how badly the tessellation was extracted rather than of anything about the tissue. μ2\mu_2 carries that information; the variance hides it.

How many cells it takes

The statistic is worth a sample-size calculation, in the currency this collection prices everything else in.

μ2\mu_2 is a mean of squared departures over nn cells, so its standard error is the standard deviation of those squared departures over n\sqrt{n}. On a spiral lattice about a quarter of the cells are non-hexagonal and contribute 1 each, so the squared departures are mostly zeros and ones and their standard deviation is about 0.45 — giving a standard error near 0.45/n0.45/\sqrt{n}.

To separate a spiral lattice’s 0.25 from a Poisson set’s 1.83 needs essentially no cells at all: a dozen would do it. To separate the golden angle’s 0.253 from the 137.0° head’s 0.291 — a difference of 0.038 — needs about 1,100 cells, which is a head of fifteen hundred primordia measured to its rim.

That asymmetry is the useful output. The statistic settles “is this tissue ordered” on a fragment and cannot settle “which ordered angle is this” on a whole head. Since the second question is answered instantly by counting the spirals, the division of labour is clear enough: count for the angle, measure μ2\mu_2 for the disorder, and do not expect either to do the other’s work.

The most hexagonal tissue is the least interesting one

The row that should stop a reader is the first. At 144° — two fifths of a turn exactly — the tessellation is 99% hexagons and μ2\mu_2 is 0.023, an order of magnitude tidier than the golden-angle head.

That angle produces a whorled pattern: five straight radial rows, with everything lined up. It is the arrangement the whole subject treats as the degenerate case, the one where the parastichy numbers share a factor and no divergence angle is recoverable from the counts.

So if “hexagonality” were a measure of how well a tissue is arranged, the best-arranged tissue in this set would be the one with the least structure in it.

That is the same inversion the site’s foundation phase found for area-evenness, where the criterion was won by rational angles whose sliver cells are near-identical while the gaps between the rays are enormous. Two criteria, two phases apart, the same answer: a regular-looking local statistic is optimised by a pattern that is globally degenerate.

One packing criterion across the angles, with the others' winners markedThe three criteria pick 137.5°, 138.0° and 135.0°. The golden angle is near the top of all three and the exact winner of none at this size.00.2500.5000.7501120130140150divergence angle (°)closest pair, as a fraction of the mean spacing (higher is better)137.508°400 points per anglethree criteria, three winners
Fig. 3 The earlier version of the same lesson. Four packing criteria on the same heads, three winners, and the criterion that rewards local uniformity rewarded the arrangement with the largest gaps.

And it does not single out the golden angle

The other row worth pausing on is the pair in the middle. The golden angle gives 0.253. A rational angle eight thousandths of a degree away gives 0.255.

No distinction. The head at 137.0°, half a degree away, gives 0.291 — a real difference, and not an optimum: it is simply a different arrangement with a slightly broader side distribution.

This is the seventh criterion this collection has measured against the claim that the golden angle is special at a fixed head size, and like the other six it does not single it out. The finding is now well enough established to state as a general result rather than as a series of surprises: no local geometric criterion measured on a finite head distinguishes the golden angle from its rational neighbours.

The claim that does survive is arithmetic rather than geometric — resistance to rational approximation, which is a statement about limits and cannot be tested by building one head — and this measurement does not touch it either way.

The spiral counts four different divergence angles produceFibonacci counts come from one angle. The Lucas angle — which the same dynamical model reaches on a different branch — gives 47 and 76, and neither number is a Fibonacci number.golden 137.51°55 and 89FibonacciLucas 99.50°47 and 76Lucas151.14°50 and 81neither77.96°37 and 60neithercounted from the pointsone sequence per branch
Fig. 4 The same negative result from the counting side. Sweep the divergence angle and measure a property of the resulting head; the golden angle is not at an extremum of it, and its rational neighbours are indistinguishable from it.

Why the rim cut is not cosmetic

One methodological point, because the site has been burned by it here before.

The measurements above use only cells inside 86% of the head’s radius. That cut was introduced in an earlier essay for a reason worth repeating: sixty boundary cells once inverted the sign of Lewis’s law on this site, giving a fitted slope of −0.63 where the true value is small and positive. The population responsible was a few dozen enormous five-sided cells one ring inside the convex hull.

Side counts have the same defect in a milder form. A cell near the boundary has fewer neighbours because the neighbours are missing, not because the arrangement put fewer there — so a rim cell’s side count is a measurement of the edge of the sample.

Including the rim would raise every μ2\mu_2 in the table, and it would raise them unequally, because the arrangements differ in how much of the disc the outermost ring occupies. The comparison would then be partly a comparison of rim effects.

Voronoi cells of a head at 137.51°171 bounded cells, averaging 5.87 sides. Cells with six sides are shaded; the others are what a tiling has to contain to close up.170 bounded cellsmean 5.87 sides
Fig. 5 The tessellation the numbers come from, with the rim visible. The cells outside the cut are bounded, have areas, and are not measurements of the arrangement.

What a tissue paper should report

Three numbers, and the first two are free once the cells are drawn.

The mean side count — as a check that the tessellation is well formed, not as a result. If it is not six, something is wrong with the tiling or the rim cut, and that is exactly what it is good for.

μ2\mu_2 — the disorder, which is the quantity that varies and the quantity every relation in this field is implicitly about.

The distribution itself — because μ2\mu_2 is a second moment and two very different distributions can share one. A golden-angle head is 13% pentagons, 75% hexagons, 12% heptagons; a Poisson set runs from three sides to eleven. Both are summarised by a single number and only one of them is well summarised by it.

That is a modest recommendation and it is not novel — the disorder measure is standard in the cellular-tissue literature and has been for decades. What is worth saying is the reason: not that μ2\mu_2 is a better statistic in general, but that the mean is determined by a theorem and therefore cannot be evidence about anything, which is a stronger objection than “it is uninformative”.

Cell area against side count, at 0% disorderThe dashed line is Lewis's law, (n−2)/4. The fitted slope here is 0.009 against his 0.25, and the side count accounts for 16% of the variation in area.00.500155.5066.507sides of the cellarea, as a multiple of the mean cell areaLewis607 interior cellsslope 0.009 against 0.25
Fig. 6 Where the distribution matters rather than the mean, in the site’s earlier work. Lewis’s law is a relation between side count and cell area, and how much of the variance it explains is a question about the joint distribution that no pair of means can answer.

The head half a degree away is the most disordered lattice

One row has been passed over and it deserves its own paragraph, because it is the only lattice in the table that is worse than the golden angle rather than indistinguishable from it.

The head at 137.0° gives μ2=0.291\mu_2 = 0.291 against the golden angle’s 0.253 — about 15% more disorder — and 71% hexagons against 75%. The difference is small and it is larger than the difference between the golden angle and its own rational neighbour, so it is the only comparison in the table that is doing anything.

The reason is visible in the contact network rather than in the side counts. At 137.0° the pattern’s families are 8, 13, 21, 29, 50, 71, 92, 113 — a sequence that adds 21 over and over rather than adding the two previous — and the essay two fields away sets out why. A lattice whose family orders grow slowly has many families of comparable strength at a given radius, so its cells are drawn from a broader population of neighbours, and a broader population of neighbours is a broader population of side counts.

So μ2\mu_2 is not measuring “distance from the golden angle”. It is measuring something about the continued fraction of the divergence angle, and 137.0° scores worse than 137.5° for the same reason it produces a whorled-looking pattern: a large partial quotient. Whether that generalises — whether μ2\mu_2 tracks the continued fraction across many angles — is a sweep this collection has not run, and it is the obvious next measurement on this statistic.

What this leaves open

The measurement here is of arrangements the site builds. Whether a real tissue’s μ2\mu_2 distinguishes anything is a question about specimens, and the honest answer is that this collection has none.

What can be said is what a specimen would have to show. A real capitulum, if it is a golden-angle lattice, should give μ2\mu_2 near 0.25 and about three quarters hexagons. Substantially more disorder than that is either a genuinely irregular pattern or a tessellation taken from a poor image, and the two are separable by remeasurement. Substantially less is a whorled arrangement, and the counts would show it.

That is a check with a number attached, which is more than the mean has ever been able to offer.

Two laws, two tilings, and they disagree about which tiling is tissueLewis's law wants disorder: its slope is 0.231 on the random set and 0.009 on the golden head. Aboav's relation wants order: a = 1.18 on the head, 0.59 on the random set.Lewis slope — golden head0.009the law says 0.25Lewis slope — random set0.231the law says 0.25Aboav a — golden head1.177the law says 1.2Aboav a — random set0.593the law says 1.2filled where the tiling obeys the law it is being judged by900 points in each tilingLewis explains 32% of the area spread at best
Fig. 7 The other two relations this field cites about the same cells. Both are statements about the distribution rather than about its mean, and both were measured here to want opposite kinds of tissue.

The general shape, which is not about tissue

Strip the subject away and the argument has a form worth carrying elsewhere.

A quantity is reported everywhere in a field. It is easy to compute, it is stable across specimens, and its stability is read as evidence that the systems being measured have something in common. Then someone works out that a theorem fixes it, and the stability stops being evidence of anything: it was never a measurement of the specimens at all.

The diagnostic is cheap and it is the question this essay is really recommending: what would have to be true of the sample for this number to come out differently? If the answer is “nothing that could happen in this system”, the number is a check on the apparatus rather than a result, and it should be reported as one.

Applied here it takes a minute. Move the points anywhere; as long as three cells meet at each vertex, the mean is six. There is no arrangement of primordia that averages five. So the mean cannot distinguish arrangements, and anything a paper concludes from its constancy across species is a conclusion about planar graphs.

This collection has now found the same shape three times in different clothes — a statistic invariant to shuffling, a statistic fixed by a theorem, and a parameter that stopped mattering above a threshold. In each case the reported quantity was real, correctly computed, and incapable of varying with the thing it was being used to describe.

How the largest gap behaves as the head fillsThe rational angle's gap grows by a factor of 3.8 over this range; the golden angle's stays within 1.30. This is the claim about 137.5° that survives measurement.02.5057.50105001e+31.5e+32e+3number of primordialargest empty gap, in units of the mean spacing137.508° (golden)137.3°135° = 3/8 of a turngaps measured from the triangulationan asymptotic claim, not a contest
Fig. 8 The largest empty gap as a head grows, which is a third statistic on the same cells and varies where the mean cannot.
How nearly each angle is a simple fraction of a turnA dip means a rational approximation and a lattice with visible rows. The golden angle's floor is 0.008; the rational angle reaches exactly zero.00.2000.400204060denominator qhow nearly q × angle is a whole turn (0 = exactly)golden 137.508°137.3°135° = 3/8distance to the nearest whole turnno dips means no rows
Fig. 9 The claim about the golden angle that survives measurement, and it is arithmetic rather than geometric — which is why a seventh geometric criterion failing to single it out changes nothing.

Shares its objects with

Essays that name at least two of the same things, and that neither author linked.

Named objects

A flat tag is an object no other essay names yet.

Aboav–WeaireCell areaDisorderEuler's formulaGolden angleLatticeMeasurementNull modelPackingRational angleRim effectSummary statisticVoronoi cellsWhorled