Packing and tiling

A width read off a staircase

Two fractions of the fourteen measured return a dip width that moves by a factor of two when the head size changes, where the others hold to three per cent. The cause is not their neighbourhood. It is that the disorder statistic changes only when the tessellation changes, so the curve a half-width is read off is a staircase, and a width narrower than the tread cannot be read at all.

Worth reading first: What a summary throws away · The six are the spirals · Packing, measured four ways.

Fourteen fractions were measured for the previous essay, twelve inside three families and two outside them. Ten returned a dip half-width that agrees between two head sizes to within three and a half per cent, once the law’s own n² scaling is applied. One agreed to twenty-one per cent, for a reason the phase before had already found.

And two disagreed by factors of 1.6 and 2.4.

They are 9/34 and 24/55, and this essay is about why — because the answer turns out to be a property of the statistic rather than of those two fractions, and it sets a resolution on every width in the thread.

The first guess, which was wrong

The obvious explanation, and the one written into this collection’s own working notes before the diagnosis was done, is the neighbourhood. This thread’s standing difficulty has been that a dip is measured against a background and a background is a stretch of neighbouring angles; the previous phase found that a single background sample 0.2° from 21/55 lands inside 13/34’s dip, reads low, and inverts the verdict.

So a fraction whose nearest rational sits unusually close would have its shoulder lifted by that neighbour’s dip, and the half-way crossing would land somewhere that has nothing to do with its own width.

That is not what is happening. 24/55’s nearest rational is 17/39, a sixth of a degree away; 9/34’s is 14/53, at a fifth of a degree. The offsets a width is read over run to five hundredths of a degree. The neighbours are three times outside the band.

The gaps close faster than the dips narrowFor each Fibonacci fraction, the distance to the nearest other rational with a denominator of 60 or less, and the half-width of its own dip at the smallest head that resolves it. The gaps fall from 0.763° at 3/8 to 0.0735° at 34/89; the dips stay between 0.0077° and 0.0155°. The dips never touch — the closest they come is a factor of 10 — so what stops the measurement is not the dips overlapping but the background between them ceasing to be flat. The clear offsets available fall from 72 to 53.-2-1.50-1-0.5000divergence angle, as a fraction of a turndegrees (logarithmic)3/85/138/2113/3421/5534/89to the nearest other rationalthe dip's own half-widthsix denominatorsgenerated from a stated rule, not drawn to look right
Fig. 1 The crowding, which is a real difficulty and is not this one. Both refused fractions have neighbours at a comparable distance to the fractions that measure cleanly.

What is actually happening

Print the profile — the second moment against the offset from the rational — and the shape is not a curve.

At 24/55 on a head of 1,279 organs, reading outwards from the exact rational, μ₂ runs

0.094 0.094 0.094 0.094 0.094 0.094 0.094 0.095 0.095 0.112 0.130 0.131 0.131 0.132 0.132 0.133 0.134 0.200 0.261

Nine identical values, a jump, two more jumps, six identical values, another jump. It is a staircase: flat treads separated by risers.

That is not a numerical accident. The second moment is a mean of squared departures from a common side count, over a finite head, and a cell’s side count is an integer. As the divergence angle is moved continuously, nothing changes until some cell in the tessellation gains or loses a neighbour, and then the statistic jumps. Between those events it is exactly constant.

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 Where the integers come from. Every cell in a head has a whole number of sides, and the statistic this thread is built on is the spread of that whole number — so it can only change when a cell’s neighbour list changes.
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. 3 The statistic itself, and what it throws away. A distribution over the integers has a second moment that moves in steps of order 1/n as one cell changes, and a head of a thousand organs has steps a thousandth of the way up its own range.

It is worth seeing the same shape at a fraction that measures cleanly, so that the staircase is not read as a property of the two failures.

At 21/55 on a head of 1,279 the profile runs 0.094 nine times, then 0.095 four times, then jumps to 0.134, 0.172, 0.173, 0.174, 0.175, and on up. Same structure: long treads, sudden risers. At 13/34 on 791 organs it runs 0.071 twelve times, then 0.072 three times, then 0.119, 0.119, 0.120, 0.120, 0.121, 0.122, then 0.192.

Every profile in this thread is a staircase. What differs between the fractions that measure and the fractions that do not is only where the half-way level happens to fall.

Why that makes a half-width unreadable, sometimes

The half-width is defined as a level crossing: find where the profile has climbed half way from the dip’s floor to its background, and report the offset.

On a staircase, a level crossing has two cases.

The level falls inside a riser. The crossing is where the riser is, the interpolation between the two sampled offsets either side is a good estimate of it, and the answer is reproducible.

The level falls on a tread. There is no crossing there at all. The interpolation returns the offset of the next riser’s foot, which is a property of where the tessellation happens to change rather than of the dip’s shape — and when the head size changes, the tessellation changes at different angles, so the treads move and the reported width moves with them.

At 24/55 the half-way level is 0.125. The profile goes 0.112 at an offset of 0.0025° and 0.130 at 0.003°, so the level falls inside a riser and is reported at 0.00285°. On the larger head the level is 0.090 and the profile goes 0.084 at 0.0037° and 0.097 at 0.0045°, giving 0.00405°. Both are honest readings of where a riser is; they are not readings of the same quantity.

For comparison, 21/55 on the same two heads gives 0.00554° and 0.00234°, which after the n² scaling agree to one per cent.

Four fractions of 34, one widthThe half-width of the second-moment dip at four divergence angles with the same denominator, on heads of 791 organs — 23 in each of 34 rows. 13/34 is the Fibonacci convergent, the fraction every earlier measurement of this law was made at and the one whose neighbourhood is emptiest; 9/34, 15/34, 11/34 are not. The four agree within a factor of 1.14, well inside the factor of two the law is stated to, so the width is a function of the denominator and not of how well the fraction approximates its neighbours. What is left over is ordered by the note beside each bar: the more crowded the neighbourhood, the wider the dip comes out, which is the direction a neighbour's own shoulder would push it.half-width of the dip, in degrees, at 791 organsnote: how far the nearest other rational sits13/341.10e-221/55 at 0.193° · convergent9/34unmeasurable156% apart at two head sizes15/341.04e-226/59 at 0.179°11/349.68e-312/37 at 0.286°q = 34 · 791 organsgenerated from a stated rule, not drawn to look right
Fig. 4 The family the refusal appears in. Three bars and one refusal, and the refusal is stated on the figure rather than left as a gap — a measurement that quietly omits its failures is reporting a subset chosen by its own defect.

There is a tempting misreading to head off. It is not that the two refused fractions have narrower dips than the others — 24/55’s floor and background are ordinary, and its width, whatever it is, is somewhere in the range its family occupies. It is that the half-way level lands badly for them. A fraction with a perfectly typical dip can be unmeasurable by this method and a fraction with an unusual one can measure cleanly; which happens is decided by the arithmetic of where the tessellation flips, and that has nothing to do with the dip.

Which means the refusal carries no information about the fraction. It is not a finding about 9/34 or 24/55. It is a finding about the method, and the two fractions are the two that happened to expose it.

The resolution this puts on the whole thread

The useful form of the finding is not “two fractions failed”. It is that the method has a resolution, it can be stated, and every width in the thread should be read against it.

The tread spacing near the half-way level is what sets it. On these heads the risers near the crossing are a thousandth to a few thousandths of a degree apart, which is the same order as the widths being measured — 5.5 × 10⁻³ at 21/55, 1.1 × 10⁻² at 13/34. So the thread has been operating within a factor of a few of its own resolution throughout, and it did not know.

Two consequences follow, and one of them is reassuring.

The measured widths are good to about a tread. Which is to say a few tenths of themselves for the narrow ones and a few per cent for the wide ones — and that is consistent with the agreement actually observed between head sizes, which is what makes the account believable rather than merely available.

And the law is safe, because the law is a scaling over a factor of eleven in denominator and a factor of two and a half in head size. A resolution of a few per cent to a few tenths does not touch a factor-of-twelve trend, which is why the previous phase’s result stands and only its precision is at issue.

Divided by the denominator, the coefficient is one numbern²·w divided by q, measured at the two larger of each denominator's three head sizes and averaged, for six Fibonacci fractions. Undivided it runs 1188, 1285, 2999, 7005, 9125, 14297 — a spread of 12.0. Divided by the denominator it runs 149, 99, 143, 206, 166, 161 — a spread of 2.1. So the law is w ≈ 150·q/n² to within a factor of two, and the previous phase's constant of 1,400 to 3,300 was three measurements of the small-denominator end of it.22.102.202.30divergence angle, as a fraction of a turnn²·w divided by the denominator (logarithmic)149991432061661613/85/138/2113/3421/5534/89spread 2.1 across the sixsix denominators · n²·w/qgenerated from a stated rule, not drawn to look right
Fig. 5 The law with the resolution in mind. The trend spans a factor of twelve; the resolution costs a few tenths at worst. What cannot be extracted from these measurements is a coefficient to better than the factor of two it is already quoted to.

Why the check was two head sizes and not two anything else

The disagreement that produced this essay came from measuring each fraction at two head sizes and comparing through the law’s own scaling. That choice is worth defending, because it is the only one of the available checks that would have worked.

Comparing two fractions would not: the whole question of the previous essay is whether two fractions should agree, so a disagreement is ambiguous between a defect and a result.

Comparing two background methods would not: the previous phase already did that and fixed the background, and a staircase in the profile survives any choice of background level.

Comparing two sampling grids on the offsets would not, for the reason given below: the risers are already resolved.

Two head sizes work because the law says exactly what should happen — n²w is constant — and because changing the head changes where the tessellation flips without changing the dip. It is a comparison with a prediction in it, which is what makes a disagreement diagnostic rather than merely a spread.

One number is worth recording for whoever does the repair, because it is the thing that makes the area version attractive. The dip’s depth — the distance from its floor to its background — is a difference of two numbers each of which is stable to a per cent or so, and it is reproducible at every fraction measured here, including the two whose widths are not. 24/55’s floor is 0.094 at 1,279 organs and 0.060 at 1,980, and its background 0.218 and 0.179; the ratio of the two falls in the same proportion as every other fraction’s.

So the depth is measurable where the width is not, and an area is a depth times a width. A quantity built from the profile as a whole rather than from one level crossing on it would inherit the depth’s stability, and the two refusals in this essay would become measurements.

What would fix it

Three repairs, none of them run in this phase, and it is worth being explicit about why not.

Read the width from an area rather than from a crossing. The integral of the profile’s departure from its background is a smooth function of the head size even when the profile is a staircase, because integration does not care where the risers are. This is the right repair and it changes the definition of the quantity, so it cannot be done without re-measuring every width in the previous phase — which would mean this phase reporting a different law from the one it set out to test.

Average over head sizes. Cheap and partial: it reduces the tread’s influence without removing it, and it costs a factor of three in computation on the thread’s most expensive machinery.

Or sample the profile densely enough to find the risers. This does not help. The problem is not that the risers are missed — the sampling already resolves them — but that the level being sought does not always fall on one.

The first is the one a later phase should do, and doing it would settle whether the residual ordering by crowding in the previous essay is real or is a second face of this same effect.

The neighbourhood of 21/55, and where its background was taken fromμ₂ across nine tenths of a degree either side of 21/55, on a head of 825 organs — 15 in each of its 55 rows. The deep notch at the centre is the dip. The shaded columns are the offsets this site now takes a background from: those at least 0.03° from every rational with a denominator of 60 or less, of which there are 54 here. The marked sample at 0.2° is the one the previous phase used, and it lands 0.007° from 13/34 — inside that rational's own dip. It reads 0.115 against a floor of 0.148, so measured that way the dip is an inversion. The clear offsets give 0.339.00.2000.400-0.400-0.300-0.200-0.10000.1000.2000.3000.400degrees from 21/55, at 137.4545°μ₂, the second moment of the side-count distribution13/34the background from the clear offsets: 0.339the old single sample: 0.11521/55 · head of 825generated from a stated rule, not drawn to look right
Fig. 6 The background the level is measured against, which has its own uncertainty. A crossing between two numbers that each carry a resolution is worse than either, and this is the other half of the width’s error budget.

What a staircase says about the dip itself

Set the measurement aside for a paragraph, because the staircase is not only an obstacle. It is a statement about what a rational divergence does to a head, and it is a sharper one than the dip.

At an exactly rational angle the organs fall into q radial rows, and the tessellation is as regular as that arrangement allows. Move the angle by a millionth of a degree and nothing happens at all — not a small change, no change: the same cells have the same neighbours and the statistic is identical to the last digit. The arrangement’s combinatorics are unchanged over a finite interval of angles.

That interval is the first tread, and its width is the dip’s true floor: the range of angles over which a head of n organs is combinatorially the rational head. Everything the dip’s width is trying to capture is a coarse-grained version of that.

It also says what a real plant would have to do to sit in a dip. Not to have a rational divergence — no measurement could establish that and no plant could hold it — but to have one within the first tread, which for a head of a thousand organs at denominator 55 is a few thousandths of a degree. The dip is not a region of near-rationality; it is a region of exactly the same tessellation, and its edges are the angles at which the first cell flips.

That reading was available from the previous phase’s data and nobody looked. It is the more interesting half of what this essay found while diagnosing a failure.

Where this fits in the collection

This is the fourth defect on this site whose symptom was a plausible number.

A spiral counter returned the two smallest index offsets instead of the two shortest surface hops, and was believed until a recovery step refused. Ticks vanished from an axis with a descending domain, and two figures carried no gridlines for three phases while every gate passed. A single background sample landed inside a neighbour’s dip and inverted a verdict. And a width has been read off a staircase at a resolution comparable to itself.

All four share a shape: the output was well-formed and about something else. No assertion failed, because nothing had been asserted about the thing that was wrong.

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. 7 The thread’s founding observation, which turns out to apply to its own machinery: a summary statistic discards the structure that would have shown it was not measuring what it was taken to. A mean side count hides the distribution; a level crossing hides the staircase.

The habit that catches them is not more assertions of the same kind. It is measuring the same quantity a second way and requiring the two to agree — a recovery against a count, a second window against a first, two head sizes against each other. The refusal in this essay exists because somebody made the width measure itself twice, and the whole of the finding is the disagreement.

One more thing is worth recording, because it is the sort of detail that gets lost between phases and then costs a day. The two refusals were not found by a gate. There is no check on this site that asks whether a width is measurable; what exists is a check that asks whether the reported widths agree, and it exists because the previous essay needed a stability criterion for a different reason — to decide which fractions could enter a family comparison.

So the refusal is a by-product of a threshold introduced for bookkeeping. That is worth noticing rather than celebrating: a defect whose discovery depended on somebody needing an unrelated criterion is a defect that could as easily have gone another phase unnoticed, and the collection’s other three of the same shape were all found the same way.

The dip at 5/13 — 138.4615° — at three head sizesWalking the divergence angle off an exact rational, at 300, 600, 1200 points. The floor falls as the head grows — 0.059 at 300, 0.029 at 600, 0.014 at 1200, halving for each doubling — and the dip narrows faster: half-widths of 0.0162°, 0.0048°, 0.0010°, a factor of four for each doubling rather than two. Half-width times the square of the head size is 1456, 1726, 1402, which is what makes the width a property of the sample rather than of the angle.00.2000.4000.600-3-2.50-2-1.50-1-0.700offset from the rational, log₁₀ of degreesμ₂300 points600 points1200 points5/13 of a turn · 12 offsets · vertical rules are the half-widthsgenerated from a stated rule, not drawn to look right
Fig. 8 One fraction at three head sizes, which is the check that produced this essay — at 5/13, where the older background method still works. Three points that scale as they should are evidence; two points that do not are a refusal.
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. 9 The two laws about a head’s cells this thread separates. Both are statements about integers per cell, and both therefore have a staircase under them; only one of them has had a width read off it.
The disorder of a head against its divergence angle, 300 pointsμ₂ is the mean squared departure of a Voronoi cell's side count from six, measured on 300 points inside 86% of the radius. Swept across 1.60° it is a staircase: flat over stretches of a few hundredths of a degree, with sharp steps between them and narrow deep dips wherever a rational falls. At 137.144° — which is 360 × 8/21 — it is 0.078; At 137.648° — which is 360 × 13/34 — it is 0.197; At 138.460° — which is 360 × 5/13 — it is 0.060. The ticks along the top are the fractions p/q, placed from arithmetic rather than from the curve.00.2000.4000.600137138138139divergence angle, degreesμ₂, the mean squared departure of a cell's side count from six5/138/2113/34401 angles · 0.0040° apart · 300 points eachmarks are the fractions, placed from arithmetic
Fig. 10 And the object under all of it. Every dip in this curve is measured by the method this essay puts a resolution on, and the curve itself is the staircase at a coarser scale.

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.

ArtefactCensusConvergentsDiscretisationFalsifiabilityHonest limitsHow many sidesIdentifiabilityMeasurementMeasurement errorRational angleSamplingSummary statisticToleranceVoronoi cells