Animated groups · reports · magic theorem

A magic theorem for clockwork animations

Five theorems about a periodic field whose spatial symmetries are allowed to cost time, their proofs, three worked examples, and a check against every animation this project has published. Nothing below assumes a differential equation, a bifurcation, or a small amplitude. Companion to the tutorial How to find a time-symmetric animation, which uses these statements but does not prove them.

What is being claimed#

Conway's Magic Theorem says that a wallpaper signature is legal exactly when its cost is 2. This report is the same shape of statement one level up. Decorate each symmetry of a wallpaper group with a fraction of a period — the time by which the animation must be shifted for that symmetry to hold — and a second closure identity appears, from the same orbifold relation, with the same proof in one line:

∑ckc/nc∈Z.\sum_c k_c/n_c \in \mathbb{Z}.

Conway's cost balances the geometry; this residue balances the choreography. It holds in all 68 entries of this project's catalogue, and — like Conway's — it enumerates: on the orbifold 333 there are 27 ways to hang offsets on the three cone classes, the identity keeps 9, relabelling collapses those to 4, and the catalogue's four p3 clockwork groups are exactly those four.

The identity does not stand alone. It is what makes two other forced structures compatible. A rotation centre whose generator costs a non-zero fraction of a period must be a zero of the first temporal harmonic — a point of the picture that does not oscillate at the fundamental frequency at all — and the phase must wind around it by an amount the offset fixes modulo the order of the centre (Theorem B). Those windings must cancel over one cell (Theorem C). B fixes their sum modulo 1 and C demands it be an integer, so B and C are simultaneously satisfiable exactly when the residue identity holds. That is Theorem D, and it is load-bearing rather than decorative.

The five statements, in one place.

Everything is a statement about one object: the first temporal harmonic of the field, one complex number per point of the plane. None of it mentions the equation that produced the animation, so none of it can be used up by changing the equation. What changes with the equation is only which of the permitted structures gets realised — and the last section of this report is about exactly that gap.

Diagram of the 333 orbifold of g227 drawn as a curved triangle with a dashed second face behind it. Each corner carries a coloured rotation marker and a label: 3 with time offset zero, 3 sub 1 with offset one third, 3 sub 2 with offset two thirds, and beside each the winding it forces: 0, minus 1, plus 1. To the right, three short equations for the Conway cost, the temporal residue and the index budget.

Shown at its own size here: drag the diagram sideways to read its labels.

The clockwork symbol, read as a decorated orbifold. Dividing the plane by g227 gives Conway's 333 turnover; the symbol r33₁3₂ then decorates each cone point with a time offset. Conway's magic theorem balances the geometry — the three costs sum to exactly 2, so the symbol tiles the plane. Applying the offset τ\tau, which is a homomorphism, to the orbifold relation balances the choreography: 0/3+1/3+2/3=10/3 + 1/3 + 2/3 = 1, an integer, and the same sum is an integer for all 68 catalogue entries (values 0, 1 and 2). The offsets then force the windings 0, −1 and +1, which sum to zero over one primitive p3 cell when each is weighted by 1/31/3. Drawn on the simulation box, which is three p3 cells, that is six pinned vortices, three of each sign.

Before any solver is started, the symbol alone decides a great deal about the picture: which points must be dark, how much phase must wind around each of them, and how the windings must balance across a cell. What follows is that statement, proved, and then tested against every animation the project has published. It is deliberately arranged the way Conway's Magic Theorem is arranged — a closure identity that the orbifold forces, and a topological budget that the identity makes satisfiable.

The data: a clockwork group#

A clockwork animation is a field u(x,t)u(x, t) on the plane, periodic in time with period TT, together with a wallpaper group GG and a map τ\tau assigning to each symmetry a fraction of the period, such that

u(gx, t+τgT)=u(x,t)u(g x,\ t + \tau_g T) = u(x, t)

for every g∈Gg \in G, every xx, every tt.

Apply the relation twice, once at gxgx and once at xx, and τhg=τh+τg\tau_{hg} = \tau_h + \tau_g drops out: τ\tau is a homomorphism G→R/ZG \to \mathbb{R}/\mathbb{Z}, the modulo coming from TT-periodicity. That single fact carries the whole theorem. Two consequences are worth naming at once. Its kernel K=ker⁡τK = \ker\tau is the subgroup of instantaneous symmetries — the ones that hold frame by frame — and the quotient G/KG/K embeds in R/Z\mathbb{R}/\mathbb{Z}. It is not automatic that it is finite: a travelling wave u(x,t)=f(k⋅x−ωt)u(x, t) = f(k \cdot x - \omega t) satisfies the defining relation for every lattice translation with τ\tau generically irrational, and then the image of τ\tau is dense. Finiteness is an assumption here, and it is the one the catalogue encodes; for a group with a rotation it is also forced, because invariance of τ\tau on the translations under an order-nn rotation clears the denominators. Granted finiteness, a finite subgroup of R/Z\mathbb{R}/\mathbb{Z} is cyclic, so G/KG/K is cyclic of some order mm.

This is the same data as a colour group whose colour permutation is cyclic: a wallpaper group with a homomorphism onto Z/m\mathbb{Z}/m, colours replaced by time offsets k/mk/m. Every entry of the catalogue is listed in exactly that form on the correspondence pages — for the worked example below, correspondence-p3.html#g227, clockwork symbol r33₁3₂, orbifold 333, mm = 3. The catalogue's own notation page defines the decorations used below — the subscripts, the tilde for a half-period mirror, the prefixes c and r: notation.html. It is also the wallpaper instance of the pair (H,K)(H, K) of spatio-temporal and spatial symmetries that Golubitsky and Stewart attach to a periodic solution of an equivariant system, with H=GH = G and H/KH/K cyclic — here cyclic by construction rather than by hypothesis.

The object every statement below is about is the first temporal harmonic of a channel,

A^j(x)=1T∫0Tu(x,t) e−2πijt/T dt,A^:=A^1.\hat A_j(x) = \frac1T \int_0^T u(x, t)\, e^{-2\pi i j t/T}\, dt, \qquad \hat A := \hat A_1.

It is one complex number per point: the amplitude and phase with which that point oscillates at the fundamental frequency. Two pieces of housekeeping, once. A rotation centre in the plane is a gyration point; its image downstairs is a cone point of the orbifold, and "cone class" below means one orbit of gyration points. And index, winding and charge are the same integer throughout — the number of turns arg⁡A^\arg \hat A makes on a small loop — called by whichever of the three names the surrounding sentence wants. On a stored movie it is one bin of a fast Fourier transform down the frame axis, which is exact for the band-limited films the repository saves. Everything that follows is a statement about A^\hat A, not about the equation that produced it — no reaction–diffusion, no near-Hopf assumption, no small amplitude.

The worked example throughout is the catalogue entry g227: p3, orbifold 333, clockwork symbol r33₁3₂, mm = 3, and — as a colour group — the entry at correspondence-p3.html#g227. The decorations are defined on the catalogue's notation page.

Theorem A — the character law#

Theorem A.

For every g∈Gg \in G and every integer jj,

A^j(gx)=e−2πijτg A^j(x).\hat A_j(g x) = e^{-2\pi i j \tau_g}\, \hat A_j(x).

That is, A^j\hat A_j transforms under the one-dimensional character χj(g)=e−2πijτg\chi_j(g) = e^{-2\pi i j\tau_g} of GG, which factors through G/K≅Z/mG/K \cong \mathbb{Z}/m.

Proof

Substituting s=t+τgTs = t + \tau_g T in the defining relation gives u(gx,s)=u(x, s−τgT)u(gx, s) = u(x,\ s - \tau_g T) for all ss. Put that into the integral defining A^j(gx)\hat A_j(gx) and shift the variable of integration by τgT\tau_g T; the integrand is TT-periodic, so the shifted interval is harmless, and the shift leaves behind the constant factor e−2πijτge^{-2\pi i j\tau_g}. ∎

The sign is a minus. It is worth saying loudly because the plus sign is the natural guess and it is wrong: on the published g227 record the minus law holds to 8.6 × 10⁻¹⁶ — on a symmetrised export, of which more below — while the plus law is off by exactly 1.7320508 = 3=∣1−e4πi/3∣\sqrt3 = |1 - e^{4\pi i/3}|, the largest error the character can produce. Flip any one of three conventions — the sense of the rotation, the sense of the winding, or the sign in the exponent of  — and the prediction flips with it. The repository's own earlier note left that sign as an unresolved ∓\mp; this fixes it.

A corollary and an honest limit. Expanding uu in harmonics, the jj-th term lives in the χj\chi_j sector: the mean field is genuinely GG-invariant and all the defects live in the sectors with jj not a multiple of mm. But A^j≡0\hat A_j \equiv 0 always satisfies the constraint, so the theorem gives no lower bound — an animation may hide its structure in a harmonic that vanishes. In the 532 record–channel pairs tested, none did: ∣A^1∣|\hat A_1| relative to the field's span stayed between 0.172 and 0.251.

Histogram on a logarithmic error axis. One population of blue bars clusters near 1e−15 and 1e−8; a separate orange population sits at and above 1, far to the right.
The sign in the character law is not a convention. Over 186 published record–channel rows with a time-shift order of at least 3, testing A^(gx)=e−2πijτ A^(x)\hat A(gx) = e^{-2\pi i j\tau}\,\hat A(x) never fails, and across all 532 rows of the sweep the worst error is 1.9e−08, while on those 186 rows the opposite sign never gets below 1.73. (Below m=3m = 3 the two signs coincide, so only those rows can tell them apart.) One honest qualification: the published exports are symmetrised by the site's own projection step, so those tiny residuals are the round-off of an identity. Re-sourced on data that was never projected, the law still holds — 7.3e−08 on an unprojected replay and 4.1e−10 on a double-precision orbit from the delayed-glue route.

Theorem B — what a cone point must contain#

Theorem B.

Let pp be fixed by a rotation of order n≥2n \ge 2, taken counterclockwise by +2π/n+2\pi/n, and let that rotation carry τ=k/n\tau = k/n. Then for every jj:

  1. Forced zero. If jk≢0(modn)jk \not\equiv 0 \pmod n then A^j(p)=0\hat A_j(p) = 0.
  2. Angular selection. In polar coordinates about pp, only the angular harmonics eiνθe^{i\nu\theta} with ν≡−jk(modn)\nu \equiv -jk \pmod n survive.
  3. Forced winding. If A^j\hat A_j does not vanish on a small circle about pp, its winding number there satisfies wj(p)≡−jk(modn)w_j(p) \equiv -jk \pmod n.
  4. Or no index at all. If A^j\hat A_j vanishes along a curve through pp, no winding exists there and 2 and 3 say nothing. This alternative is genuine and occurs in the data.
Proof

Write zz for the complex coordinate centred at pp; the rotation is z↦ζzz \mapsto \zeta z with ζ=e2πi/n\zeta = e^{2\pi i/n}, and Theorem A reads A^j(ζz)=e−2πijk/nA^j(z)\hat A_j(\zeta z) = e^{-2\pi i jk/n}\hat A_j(z). Expand the circle of radius rr as ∑νcν(r) eiνθ\sum_\nu c_\nu(r)\, e^{i\nu\theta}. Matching Fourier coefficients gives cν(r)(e2πiν/n−e−2πijk/n)=0c_\nu(r)\bigl(e^{2\pi i\nu/n} - e^{-2\pi i jk/n}\bigr) = 0, so cνc_\nu vanishes unless ν≡−jk\nu \equiv -jk, which is claim 2. Continuity at pp kills every cνc_\nu with ν≠0\nu \ne 0 as r→0r \to 0, and A^j(p)\hat A_j(p) is the limit of c0c_0; if ν=0\nu = 0 is excluded the value is 0, which is claim 1 (or directly: the value at a fixed point equals a non-trivial root of unity times itself). For claim 3, let Δ\Delta be the continuous increment of arg⁡A^j\arg \hat A_j across one sector of angle 2π/n2\pi/n. The relation makes the increment the same on all nn sectors and fixes it modulo 2π2\pi as −2πjk/n-2\pi jk/n, so the total winding is 2πwj=nΔ=−2πjk+2πnℓ2\pi w_j = n\Delta = -2\pi jk + 2\pi n\ell. ∎

The two halves are worth separating. Claim 1 says where the picture must be dark: any rotation centre whose generator carries a non-zero offset is a zero of the first harmonic — a point that does not oscillate at the fundamental frequency at all. Claim 3 says how the phase must run around it: the offset is a topological charge, read modulo the order of the centre. Neither statement mentions a partial differential equation.

Six panels of phase-as-hue close-ups around rotation centres, each labelled with a group, a centre type and a winding: a threefold centre with winding −1, a trio with −1, −1 and +2, a fourfold centre with −2, a sixfold centre resolved as six nodal rays with no index, the same symbol with winding +3, and a panel of crossing nodal lines.
Six cores the symbols force, side by side. The theorem fixes the winding only modulo the order of the centre, which is why g226's three 3₁ centres can read −1, −1 and +2. A 4₂ centre is driven to −2; at a 6₃ centre the tie between ±3\pm 3 need not be broken, and both outcomes occur in the published data — six nodal rays with no index, or a clean +3. In g133 the half-period mirror lines cross at the two-fold centres, which therefore lie on a nodal curve and carry no index at all, while the 4₂ cores sit off the lines and can be measured at ±2\pm 2.

Theorem C — the budget over one cell#

Theorem C.

Suppose the zeros of A^j\hat A_j on the torus are isolated. Then

  1. the indices of A^j\hat A_j over one primitive cell of the translation lattice of GG sum to 0;
  2. every orbit of the orientation-preserving subgroup carries one common index, and a cone class of order nn has ∣P+∣/n|P^{+}|/n points in that cell, where P+P^{+} is the orientation-preserving point group;
  3. therefore
    ∑cone classeswc/nc+∑free orbitswf=0.\sum_{\text{cone classes}} w_c/n_c + \sum_{\text{free orbits}} w_f = 0.
Proof

A^j\hat A_j is not in general a function on the torus: g227's translation by (1/3, 2/3)(1/3,\ 2/3) carries τ=23\tau = \tfrac23, so by Theorem A the identification multiplies A^j\hat A_j by a constant phase, and A^j\hat A_j is a section of a flat unitary line bundle. A flat bundle over a torus has first Chern number 0, which already gives the claim; but the elementary version is more instructive, and it handles the twist explicitly: cut along a fundamental parallelogram and count the boundary winding; opposite edges are identified by translations under which A^j\hat A_j is multiplied by a constant (Theorem A), and a constant factor does not change increments of argument, so the paired edges cancel exactly. For 2, an orientation-preserving symmetry multiplies A^j\hat A_j by a unimodular constant and preserves orientation, so it carries the index at xx to the index at gxgx; the orbit sizes are orbit–stabiliser. Substituting 2 into 1 and dividing by the number of orientation-preserving operations per cell gives 3. ∎

The lift-counting in claim 2 is the step that makes the budget an orbifold statement rather than a torus statement: a cone point of order nn is one point downstairs but ∣P+∣/n|P^{+}|/n points upstairs, so it enters the balance with weight 1/n1/n. The multiplicity law was checked against all 68 catalogue entries and is exact in every one.

Two panels over several lattice cells. The left shows the phase of the first temporal harmonic as hue with the nine threefold centres marked by circles and up and down triangles, each labelled with a winding of 0, −1 or +1. The right shows the magnitude, dark at the marked screw centres.
What the symbol forces, measured on the published orbit. The phase of the first temporal harmonic winds by 0 at the three plain threefold centres, by −1 at the three 3₁ centres and by +1 at the three 3₂ centres — exactly −k mod 3-k \bmod 3 at a centre whose generator carries the offset k/3k/3. The magnitude vanishes at the six screw centres and peaks at the three τ=0\tau = 0 centres. The charges sum to zero over one cell, as the index budget demands.

Theorem D — the temporal magic theorem#

Conway's Magic Theorem says a wallpaper signature is legal exactly when its cost is 2 — equivalently, when the orbifold Euler characteristic vanishes. The time offsets obey their own closure identity, with the same source.

Theorem D.

For an orientation-preserving group with cone classes of orders n1,…,nsn_1, \ldots, n_s and offsets kc/nck_c/n_c,

∑ckc/nc∈Z.\sum_c k_c/n_c \in \mathbb{Z}.
Proof

Take Conway's orbifold presentation. For the sphere-with-cone-points cases 2222, 333, 442 and 632 the generators are the primitive counterclockwise rotations, one per cone class, with γc nc=1\gamma_c^{\,n_c} = 1 and γ1γ2⋯γs=1\gamma_1\gamma_2\cdots\gamma_s = 1. Apply the homomorphism τ\tau to that last relation: each generator contributes kc/nck_c/n_c, and the product relation says the sum is 00 in R/Z\mathbb{R}/\mathbb{Z}. For the torus case there are no cone points and the sum is empty. The target is abelian, so the conjugacy ambiguity in the presentation is harmless and the argument survives verbatim on higher-genus orbifolds. ∎

Theorems B and C are compatible if and only if D holds. B forces ∑wc/nc≡−∑kc/nc\sum w_c/n_c \equiv -\sum k_c/n_c modulo 1; C demands that the same sum be an integer; the two can hold together exactly when ∑kc/nc\sum k_c/n_c is an integer, which is D. That is the shape of Conway's theorem repeated one level up: Conway's cost balances the geometry; this residue balances the choreography.

What the residue means for a group with mirrors. The proof uses the orientation-preserving presentation, so for the forty-odd catalogue entries that have mirrors or glides the sum below is taken over the cone classes of the group's rotation subgroup, which is orientation-preserving and to which the theorem applies verbatim. That is what the residue column of the table records, and it is why a ∗-entry and its rotation subgroup share a residue.

Checked that way on all 68 catalogue entries, the residue is an integer in every one — the value is 0 for 27 entries, 1 for 32 and 2 for 9. But a consistency check is the least of it. Conway's theorem is famous because it enumerates, and so does this one. Take 333: there are 333^3 = 27 ways to hang offsets on its three cone classes, Theorem D keeps the 9 with k1+k2+k3≡0(mod3)k_1 + k_2 + k_3 \equiv 0 \pmod 3, and relabelling the classes collapses those to four multisets — {0,0,0}\{0,0,0\}, {1,1,1}\{1,1,1\}, {2,2,2}\{2,2,2\} and {0,1,2}\{0,1,2\}. The catalogue's p3 entries are 333, 3₁3₁3₁, 3₂3₂3₂ and r33₁3₂: exactly those four and no others.

What D does not do is guarantee existence. 4₁4₃2 has residue 14+34+0\tfrac14 + \tfrac34 + 0 = 1, an integer, and no group on a primitive lattice realises it; it exists only on a centred space–time lattice, as c4₁4₃2 — catalogue entry g99. The residue decides what may be decorated; the lattice decides what can be built.

Theorem E — mirrors#

Theorem E.

  1. Only two offsets. A reflection squares to the identity, so 2τ≡02\tau \equiv 0 and τ∈{0,12}\tau \in \{0, \tfrac12\}.
  2. A half-period mirror is a nodal line. If τ=12\tau = \tfrac12 on a reflection in a line ℓ\ell, then on ℓ\ell Theorem A reads A^j=e−πijA^j\hat A_j = e^{-\pi i j}\hat A_j, so A^j≡0\hat A_j \equiv 0 along the whole line for every odd jj. The zero set is one-dimensional, not a point.
  3. A neutral mirror carries no charge. If τ=0\tau = 0, a reflection reverses orientation while fixing A^j\hat A_j, so the index at a fixed point is its own negative and vanishes; off the line, vortices come in mirror pairs ±w\pm w.
  4. A dihedral half-offset centre cannot have an isolated zero. If a centre of even order nn with k=n/2k = n/2 also lies on a mirror, B demands a non-zero index and E.3 demands zero; the only escape is B's fourth case — the index is undefined and the zero set is a curve of nn rays.
  5. Glides are free. A glide reflection has no fixed point, so it imposes no pointwise condition.
Proof

All five are one-line substitutions into Theorem A, together with the elementary fact that an orientation-reversing map negates a winding number. ∎

Clause 2 is the familiar antisymmetry statement in a new setting: when m=2m = 2 a clockwork group is a two-colour group, a half-period mirror is a primed mirror, and the statement is that the order parameter vanishes on a primed mirror. What is new here is that the same argument runs for every mm and every harmonic.

What this does and does not explain about the gaps in the atlas. Of the 51 catalogue entries with a non-zero offset, 35 have a Ginzburg–Landau record and 29 a Brusselator record. Every one of the 16 entries with no Ginzburg–Landau record, and every one of the 22 with no Brusselator record, is forced by the theorem to carry either a nodal line or a core of charge ≥2\ge 2 — 16/16 and 22/22, no counter-examples (Fisher exact, pp = 7.8 × 10⁻⁴ and 1.0 × 10⁻⁵). Three qualifications, all of which matter. First, that composite predicate was chosen after seeing which one reached 16/16 and 22/22: the scripts compute eight predicates against two models, sixteen tests in all, and the quoted pp values are not corrected for that choice. Second, a predicate with no theorem content at all — "the group has reflections" — reaches pp = 6.3 × 10⁻⁴ on the Brusselator gap, so the composite is not uniquely good; the sharpest single predictor is the forced nodal line alone, pp = 7.0 × 10⁻⁸, and that one the repository's own notes had already reported for Ginzburg–Landau. Third, pushing the other way, the search budget was uniform: the Brusselator batch gave exactly 4 jobs to each of the 51 entries, and all 88 jobs of the 22 entries without a record failed, which removes the most obvious confounder. The honest summary is a necessary condition satisfied with no exceptions, plus a plausible mechanism — a one-dimensional zero set is fragile under a perturbation that does not respect it, and a multiply charged core tends to split. It is not a demonstration that no orbit exists for the missing entries. The decisive experiment, running the solver on those entries and asking whether they fail for the predicted reason, has not been done.

Three worked examples#

g227 — r33₁3₂, the published spiral lattice

p3, orbifold 333, mm = 3, residue 0+13+230 + \tfrac13 + \tfrac23 = 1, forced windings (0, −1, +1), and the budget 0−13+13=00 - \tfrac13 + \tfrac13 = 0 closing on its own — so no free vortices are needed. Per p3 cell that is two charged centres plus one neutral one where ∣A^∣|\hat A| is instead maximal; on the simulation box, which is three p3 cells, six pinned vortices, three of charge −1 and three of charge +1, plus three neutral centres, total charge 0.

3₁ · τ = ⅓ · w = −1 3₂ · τ = ⅔ · w = +1 3 · τ = 0 · w = 0 one p3 cell — the box is three of these
The forced structure of g227, schematically, on one primitive p3 cell — the simulation box drawn everywhere else on this page is three of these, with nine centres rather than three. The three classes of threefold centre of a 333 orbifold sit at the corners of the cell and at the two triangle centroids. Which centroid class is labelled 3₁ and which 3₂ is a choice of basis; what the symbol forces is that the three classes carry windings 0, −1, +1, that the two charged ones are exact zeros of A^\hat A, and that the total over the cell is 0. τ=0\tau = 0, no charge, ∣A^∣|\hat A| maximal τ=13\tau = \tfrac13, w=−1w = -1 τ=23\tau = \tfrac23, w=+1w = +1

All of that is legible in the published animation. The pattern turns; the cores do not move, because Theorem B pins them to the group and not to the frame. Every third of a period the picture returns to itself rotated by 120° about a 3₁ centre — the plain centres are the ones that cost no time, and a turn about one of those is a symmetry of each frame on its own.

The pair of orbifolds is worth reading off, because it says what a simulator has to work with. In colour-group form g227 is 333/333: the kernel is again a p3, on a lattice three times coarser, so every frame separately still has threefold symmetry. Its neighbour g226 is 333/◦ — every rotation carries τ=13\tau = \tfrac13, the kernel is the translation lattice alone, and no single frame has any rotational symmetry at all. That is what the catalogue's mm column measures, and it is why the entries with a trivial kernel are the ones where a statement about the kernel-symmetric subspace has nothing to say.

g226 — 3₁3₁3₁, where the group offers a choice

Again p3 and 333, but now all three cone classes carry offset 13\tfrac13; the residue is 13+13+13\tfrac13 + \tfrac13 + \tfrac13 = 1, so the entry is legal. Theorem B allows −1-1 at each class, but then ∑w/n=−1\sum w/n = -1, and Theorem C needs the total to be 00. There are exactly two cheap ways out, and the group does not choose between them:

Both branches occur in the published data. The Brusselator records and three of the four Gray–Scott rows take the doubled core; both Ginzburg–Landau records of g226 take the minimal windings with three free +1 vortices. The mirror-image statement holds for g225, whose Gray–Scott row also takes the free-vortex branch. It is tempting — and it was, at first, tempting here — to say that the doubled core is arithmetically forced. It is not: it is merely the cheaper of two admissible structures, and which one appears is a question for the dynamics, not for the symbol. The same caution applies to g244, where nothing is forced at all — its minimal structure 6₂[−2] 3₂[+1] 3₂[+1] already balances, and the one Gray–Scott record that departs from it simply chooses −2 at the threefold classes and pays with free charge.

g133 — 4₂∗~2, a mirror example

Orbifold 4*2, mm = 2. Its boundary is a single circle with one corner, so there is one class of mirror line — the two perpendicular families one sees in the plane are conjugate by the fourfold rotation — and it carries τ=12\tau = \tfrac12, which is the single tilde in the catalogue's symbol 4₂∗~2. By Theorem E.2 the first harmonic vanishes identically along that whole family of lines: the picture is a grid of dark lines, not a lattice of isolated cores. The geometry is worth getting the right way round, because it is easy to invert. The lines cross at the twofold centres, which therefore lie on a nodal curve and have no index at all — Theorem B's fourth case, and exactly what the measurements record. The fourfold centres, which carry kk = 2 = n/2n/2, sit off the lines, at distance 0.354 in units of the cell edge, which is why they can carry the measurable charges ±2 that the budget needs. At such a chiral centre +n/2+n/2 and −n/2-n/2 are the same residue, so the group cannot fix the sign: this is the one place where the theorem yields a dichotomy rather than a value, and both branches — an isolated core of charge ±n/2\pm n/2, or a crossing of nn nodal rays with no index — are realised in the catalogue.

What was tested, and what the test can carry — 532 record–channel rows, and why the big number is weaker evidence than it looks

The experiment scripts applied every prediction above to 532 (record, channel) pairs: all 116 published equation orbits — 65 Ginzburg–Landau records and 51 Brusselator records — plus 150 rows of the atlas produced by a third equation, Gray–Scott, a stratified sample covering all 68 catalogue entries, both channels, in 203 seconds of local processor time and with zero evaluations of any right-hand side. Across those rows, 2990 rotation centres were examined, 2032 of them had a well-defined index, and the congruence w≡−k(modn)w \equiv -k \pmod n held at every single one: 0 failures. The generalisation to higher harmonics was checked separately on 13 groups at j=1…4j = 1\ldots4 — 52 pairs, 244 measurable centres, again 0 failures.

That number is weaker evidence than it looks, and the reason is instructive. Every movie the repository exports is passed through a space–time projection before it is encoded, which averages over the group's action so that the defining relation holds exactly. On the g227 record the nine operations agree bitwise. So residuals at the 10⁻¹⁵ level in that sweep are the round-off of an identity — a check on the implementation and on the sign conventions, not a measurement of the dynamics. The independent evidence comes from fields no projector ever touched:

Theorems A and B measured on unprojected data; two channels where two figures are given, and residuals relative to the field's own scale. On the double-precision field the plus sign is off by 1.7320508, exactly as on the exports. The table scrolls sideways on a narrow screen.

FieldCharacter law, minus sign ∣A^∣|\hat A| at a forced-zero centreWindings
Unprojected fourth-order Runge–Kutta replay of the g227 orbit, 336 steps from the published first frame 7.3 × 10⁻⁸ / 5.7 × 10⁻⁸ 3.4 × 10⁻⁸ / 2.3 × 10⁻⁸ 0, −1, +1 — unchanged
A double-precision g227 orbit from a different experiment, never projected (its own symmetry residual is 1.1 × 10⁻⁹) 4.1 × 10⁻¹⁰ / 2.6 × 10⁻¹⁰ 2.6 × 10⁻¹¹ three at 0, three at +1, three at −1; total 0

In other words: the theorem holds on unprojected trajectories to the accuracy the orbit itself has, which is 10⁻⁸ for float32-stored fields and 10⁻¹⁰ for double precision. That is the number to quote. What the big sweep genuinely adds is different and more interesting: it records which of the structures the theorem permits each equation actually selects — the sign at a half-offset centre, minimal windings against a doubled core, free vortices present or absent.

Two predictions in the sweep failed, and both failures were deliberate over-reaches rather than theorem failures.

The two checks that were not consequences of the theorem.
Stronger guessPassFailWhat it means
the index is always the minimal representative of −k mod n-k \bmod n 39816 All in g225, g226, g244. For g225 and g226 the minimal windings do not sum to zero, so the budget must be closed by a doubled core or by a free orbit — both are legal and both occur in the published records. For g244 nothing is forced at all: one record simply chose a non-minimal winding and paid for it.
cone indices plus free vortices sum to zero when counted by elementary mesh loops 1664 A mesh artefact. Beside a core of charge ≥2\ge 2 the phase jump between adjacent nodes reaches 0.96π–1.00π and the principal-value loop map aliases; on a mesh refined twofold and fourfold the free charge is exactly zero. The band-limited contour index used for every congruence test is stable on six radii.
All 68 catalogue entries: clockwork symbol, kernel, residue, forced minimum

These are the catalogue's own symbols, not a notation invented here: every one has Conway cost exactly 2, so it still names a wallpaper group, and the decoration is what turns that group into a clock. 3₁ is a threefold gyration point whose counterclockwise generator carries τ=13\tau = \tfrac13; 2₁ a half-turn at half a period; ~ a mirror arc that holds only after half a period; ×½ a glide that also waits half a period; and the prefixes c and r a centring of the space–time lattice — a half-cell shift at half a period, a third-cell shift at a third of a period. They are defined in full on the catalogue's notation page. The Kernel column is the orbifold of the subgroup with τ=0\tau = 0 — the symmetry a single frame still has, and the domain a fundamental-domain simulator would be left with. The residue column is ∑kc/nc\sum k_c/n_c over the cone classes of the group's rotation subgroup, an integer in every row — that is Theorem D. The last two columns give the cheapest defect structure compatible with Theorems B and C, one winding per cone class, and its cost in cores per primitive cell of the rotation subgroup's translation lattice; a cheapest structure is not a forecast of what an equation will produce, only a lower bound on what it must pay. Nine columns: the table scrolls sideways on a narrow screen.

EntryFamilyOrbifoldmm Kernel Clockwork symbol ∑k/n\sum k/n Cheapest windingscores / cell
g1p1◦1◦◦0—0
g5p2222212222222200, 0, 0, 00
g6p222222◦2₁2₁2₁2₁2−1, −1, +1, +14
g7p2222222222c222₁2₁10, 0, −1, +12
g8cm∗×1∗×∗×0—0
g9cm∗×2◦~∗×½0—0
g10pm∗∗1∗∗∗∗0—0
g11pg××1××××0—0
g54pmm∗22221∗2222∗222200, 0, 0, 00
g55pmm∗222222222∗~2~2~2~200, 0, 0, 00
g56pmg22∗122∗22∗00, 0, 0, 00
g57pmg22∗2222222~∗00, 0, 0, 00
g58pgg22×122×22×00, 0, 0, 00
g59pgg22×2222222×½00, 0, 0, 00
g60pmm∗22222∗∗∗2₁~2₁2₁~2₁2−1, −1, +1, +14
g61pmg22∗2××2₁2₁~∗2−1, −1, +1, +14
g62pmg22∗2∗∗2₁2₁∗2−1, −1, +1, +14
g63pgg22×2××2₁2₁×2−1, −1, +1, +14
g64pmm∗22222∗2222c∗222₁2₁ᵃ10, 0, −1, +12
g65pmm∗2222222∗c∗222₁2₁ᵇ10, 0, −1, +12
g66pmg22∗222∗c22₁∗ᵃ10, 0, −1, +12
g67pmg22∗222×c22₁∗ᵇ10, 0, −1, +12
g68cmm2∗2212∗222∗2200, 0, 0, 00
g69cmm2∗22222222∗~2~200, 0, 0, 00
g70cmm2∗222∗×2₁∗2₁~2₁2−1, −1, +1, +14
g71cmm2∗222∗2222c2₁∗22ᵃ10, 0, −1, +12
g72cmm2∗22222×c2₁∗22ᵇ10, 0, −1, +12
g73cmm2∗22222∗c2∗2₁2₁10, 0, −1, +12
g74pmm∗222222∗22c∗22₁22₁10, 0, −1, +12
g75pgg22×42222c22₁×¼10, 0, −1, +12
g94p4442144244200, 0, 00
g95p4442222224₂4₂21−2, +2, 04
g96p44424◦4₁4₁2₁1−1, −1, +14
g97p44424◦4₃4₃2₁2+1, +1, −14
g98p44422442c44₂2₁10, −2, +14
g99p444242222c4₁4₃21−1, +1, 02
g128p4m∗4421∗442∗44200, 0, 00
g129p4m∗4422∗2222∗4₂~4₂21−2, +2, 04
g130p4m∗44222∗22∗~4₂4₂~21−2, +2, 04
g131p4m∗4422442∗~4~4~200, 0, 00
g132p4g4∗214∗24∗200, 0, 00
g133p4g4∗2222×4₂∗~21−2, +2, 04
g134p4g4∗222∗224₂∗21−2, +2, 04
g135p4g4∗224424∗~200, 0, 00
g136p4m∗4422∗442c∗44₂2₁ᵃ10, −2, +14
g137p4g4∗24∗2222c4₃∗21−1, +1, 02
g138p4m∗44224∗2c∗44₂2₁ᵇ10, −2, +14
g139p4g4∗2422×c4₁∗21−1, +1, 02
g224p3333133333300, 0, 00
g225p33333◦3₂3₂3₂2+1, +1, −24
g226p33333◦3₁3₁3₁1−1, −1, +24
g227p33333333r33₁3₂10, −1, +12
g230p3m1∗3331∗333∗33300, 0, 00
g231p3m1∗3332333∗~3~3~300, 0, 00
g232p31m3∗313∗33∗300, 0, 00
g233p31m3∗323333∗~300, 0, 00
g234p31m3∗33∗333r3₂∗310, −1, +12
g235p31m3∗36333r3₂∗~310, −1, +12
g243p6632163263200, 0, 00
g244p6632322226₂3₂21−2, +1, 04
g245p6632322226₄3₁21+2, −1, 04
g246p663223336₃32₁1−3, 0, +16
g247p66326◦6₅3₂2₁2+1, +1, −16
g248p66326◦6₁3₁2₁1−1, −1, +16
g268p6m∗6321∗632∗63200, 0, 00
g269p6m∗63223∗3∗6₃~3~2₁1−3, 0, +16
g270p6m∗6322632∗~6~3~200, 0, 00
g271p6m∗6322∗333∗~6₃32₁1−3, 0, +16

Reading of the whole column: the residue is 0 in 27 entries, 1 in 32 and 2 in 9. Of the 51 entries with a non-zero offset, none requires a free vortex orbit — every clockwork group in the catalogue can be balanced by putting all of its topological charge on cone points, which is a statement about the arithmetic and not a promise that any particular equation will produce it — and 15 require at least one core of charge ≥2\ge 2, because their minimal windings do not sum to zero. Nothing forbids free vortices; the second reciprocal shell of g227 produces 24 cores in the box: the same six charged cone points, with the same charges, plus 18 free ones in two orbits of nine, whose charges must therefore sum to zero among themselves. They do — nine at +1 and nine at −1.

One period of the published orbit with the forced cores marked: the pattern turns, the cores do not move. Every vortex is pinned to a rotation centre of the group, and the marker at each centre is the winding the symbol predicted before the equation was solved.

What it buys, before any solver starts#

The statements above are not only a check on an answer; they are enough to draw one. Project a plane wave onto the character of Theorem A and most waves cancel identically — for g227, three of the six shortest reciprocal vectors survive and three read exactly zero, and the whole of the next shell cancels. Symmetrising one survivor gives a field with three Fourier components: three plane waves at 120°, which is the minimal p3 clockwork field.

Three panels of the seed field over two lattice cells: the real part in the ember palette showing bright hexagonal blobs, the modulus in greyscale showing dark points on a honeycomb, and the phase in hue with black and white diamonds marking vortices of opposite sign.
The seed, built from the clockwork symbol alone and before a single step of the reaction–diffusion equation. Symmetrising the three surviving waves gives a star with exactly three Fourier components — three plane waves at 120°, ∣k∣2|k|^2 = 0.032815 — that satisfies the character law to 7.0e−16. Its modulus already has exact zeros at the six 3₁ and 3₂ screw centres of the box with windings ∓1\mp 1, and the three τ=0\tau = 0 centres sit at maxima. The picture is right before any computation of the dynamics.

That picture already has the structure Theorems B and C forced: exact zeros at the six screw centres with windings ∓1, maxima at the three τ=0\tau = 0 centres, total charge 0, and the character law to 7.0e−16. Lifting it through the oscillatory eigenvector of a 2 × 2 matrix and handing it to Newton's method on a twisted shooting equation reaches the published orbit in 5 steps. The construction, the solver and two live searches are in the companion tutorial, How to find a time-symmetric animation.

A structure the theorem permits is not a structure an equation will produce. Of the 51 catalogue entries with a non-zero offset, 29 have a Brusselator record; every one of the 22 without one is forced to carry a nodal line or a core of charge at least 2. That is an exceptionless necessary condition with a plausible mechanism — a one-dimensional zero set is fragile, a multiply charged core tends to split — and it is not a proof that no orbit exists. The qualifications are stated in full under Theorem E.

What is proved, what is measured, what is open#

Proved, from the defining relation alone, with no equation and no small-amplitude assumption: Theorems A, B, C, D, E. Each was independently re-derived during review. Measured, on unprojected data: the character law and the forced zeros and windings of g227, to 7.3 × 10⁻⁸ and 4.1 × 10⁻¹⁰ on two independent fields; and, on the symmetrised exports, the full catalogue sweep as a check of conventions and implementation. Open: whether a dynamically realised orbit can hide its structure in a vanishing harmonic — permitted by the theorem, never observed here; which branch of the half-offset tie a given equation selects; whether a forced structure is attainable at all, which is the stability question the coverage gap gestures at without settling; and harmonics beyond the fourth.

Where this sits#