Path 2 — Selection · Lesson 8 of 11

How MDL Selects the Polynomial

A five-step funnel eliminates 10290 candidates. The orbit interpolation then pins the last constraint. Exactly one multilinear GF(7) rule survives.

The Bridge from Path 1

In L06 (GTE Orbits) you saw that the orbit {73, 42, 275} encodes the three charged lepton generations, and that the structure MDL forces is Z₇×Z₃. But how does the actual polynomial formula emerge from this?

The answer is a two-stage process. MDL acts twice in sequence:

Selection 1
Selects the Lepton Seed

MDL picks the arithmetic triple (1, 73, 823) as the unique starting point consistent with all four UGP invariants.

Selection 2
Selects the polynomial

The orbit's parity shadow, together with the vacuum transparency condition, uniquely interpolates the polynomial p(L,C,R) = C+R−CR−LCR.

This lesson focuses on Selection 2: how the orbit data pins the polynomial and why no other rule is consistent.

The Five-Step Elimination Funnel

MDL doesn't search through candidates one by one. It applies five structural constraints in sequence, each collapsing the space by an enormous factor:

All k=7, r=1 CA rules 7343 ≈ 10290
↓ Step 1: Binary restriction must be Rule 110
Rules whose binary restriction is Rule 110 7335 remain
↓ Step 2: Binary field is insufficient — 3 generations need k≥4
Rules over a non-binary prime field primes k ≥ 3
↓ Step 3: Derivability criterion — Z₃ must cost zero extra bits
Rules with derivable color symmetry primes k with 3 | k−1
↓ Step 4: GF(7) is the unique minimal prime
GF(7) rules consistent with orbit data one prime, 78 = 5,764,801 candidates
↓ Step 5: Direct-Interpolation Lift pins the polynomial
Unique multilinear GF(7) rule 1 rule
p(L,C,R) = C + R − CR − LCR (mod 7)
This lesson covers Step 5

Steps 1–4 establish that only GF(7) rules need be considered. L09 explains those steps in detail. Here we focus on Step 5: how the orbit data (together with vacuum transparency) uniquely interpolates the polynomial among all 78 = 5,764,801 remaining candidates.

The Direct-Interpolation Lift

After Steps 1–4, we know the rule must be a multilinear polynomial over GF(7). Such a polynomial is determined by its values on the 8 binary input triples (L,C,R) ∈ {0,1}³ — the 8 "binary corner" points. That gives exactly 8 degrees of freedom, each a value in {0,1,2,3,4,5,6}: so 78 = 5,764,801 candidates.

The orbit's parity shadow pins 7 of those 8 corners. The 8th is pinned by the vacuum transparency condition: p(0,0,0) = 0. Together, 8 constraints determine 8 unknowns. Exactly one polynomial survives.

What is the "parity shadow"?

The GTE orbit is a sequence of triples with values in GF(7). When you project an orbit state onto {0,1} by the map v → v mod 2, you get the parity shadow — a binary state. The parity shadow of the orbit {(1,73,823) → (9,42,1023) → (5,275,65535)} generates specific binary neighborhoods as the orbit evolves.

Reading off the (L,C,R) neighborhoods visited and their outputs from the orbit evolution gives us 7 input-output pairs for the polynomial. Each pair is a constraint: p(L,C,R) = output.

The 10 constraints, step by step

The full orbit produces 10 constraints (some redundant). After removing redundancies, exactly 7 are independent. Use step buttons to walk through them:

Orbit state g=1: (1, 73, 823)

The parity projection is (1 mod 2, 73 mod 2, 823 mod 2) = (1, 1, 1).

Parity neighborhood: (L, C, R) = (1, 1, 1)

The orbit map sends this state to g=2: (9, 42, 1023).
Parity of output b-value: 42 mod 2 = 0.
→ Constraint: p(1, 1, 1) = 0

The Rule 110 truth table also gives p(1,1,1) = 0 — this is consistent with Step 1.

Orbit state g=2: (9, 42, 1023)

Parity projection: (9 mod 2, 42 mod 2, 1023 mod 2) = (1, 0, 1).

Parity neighborhood: (L, C, R) = (1, 0, 1)

Output maps to g=3: (5, 275, 65535).
Parity of 275: 275 mod 2 = 1.
→ Constraint: p(1, 0, 1) = 1

Rule 110 also gives p(1,0,1) = 1 — all orbit constraints are consistent with the binary floor.

Orbit state g=3: (5, 275, 65535)

Parity projection: (5 mod 2, 275 mod 2, 65535 mod 2) = (1, 1, 1).

Parity neighborhood: (L, C, R) = (1, 1, 1)

This repeats the constraint from g=1: p(1,1,1) = 0.
→ Redundant.

The full orbit analysis (including all neighborhoods visited across the ring during evolution) produces 10 input-output pairs total. Removing redundancies leaves 7 independent constraints.

The 7 independent constraints from the orbit

(L,C,R)
p output
Rule 110?
Source
(1,1,1)
0
✓
Orbit g=1
(1,1,0)
1
✓
Orbit ring
(1,0,1)
1
✓
Orbit g=2
(1,0,0)
0
✓
Orbit ring
(0,1,1)
1
✓
Orbit ring
(0,1,0)
1
✓
Orbit ring
(0,0,1)
1
✓
Orbit ring

The missing corner: (0,0,0). This is the 8th input triple, not yet determined by the orbit.

Pinning the 8th: Vacuum Transparency

The vacuum transparency condition says: when all three input cells are in the vacuum state (0,0,0), the output must also be the vacuum (0). Equivalently: the vacuum is a fixed point of the update map.

This is not an arbitrary choice. MDL demands it: any rule where p(0,0,0) ≠ 0 would spontaneously generate excitations from empty space — an uncontrollable source of description cost. A rule that creates something from nothing costs strictly more bits to specify than one that preserves the vacuum.

The 8th constraint

p(0,0,0) = 0

Vacuum transparency is a derived requirement, not an assumption. MDL forces it on any physically realizable substrate.

All 8 constraints together

(L,C,R)
p output
Rule 110?
Source
(1,1,1)
0
✓
Orbit
(1,1,0)
1
✓
Orbit
(1,0,1)
1
✓
Orbit
(1,0,0)
0
✓
Orbit
(0,1,1)
1
✓
Orbit
(0,1,0)
1
✓
Orbit
(0,0,1)
1
✓
Orbit
(0,0,0)
0
✓
Vacuum

Green = pinned by vacuum transparency condition.

The unique survivor

8 constraints, 8 unknowns. Over GF(7), there is exactly one multilinear polynomial that satisfies all 8 simultaneously. The Lagrange interpolation formula over GF(7) computes it:

p(L,C,R) = C + R − CR − LCR (mod 7)

Every other multilinear GF(7) rule fails at least one of the 8 constraints. This is a theorem, machine-certified in Lean 4:

Candidates going in
78
5,764,801
→
Survivors
1
p = C+R−CR−LCR

Verification: p Satisfies All 8 Constraints

The formula p(L,C,R) = C + R − CR − LCR can be checked against all 8 binary inputs by hand. Here are three spot-checks:

p(1,1,1)

= 1 + 1 − 1·1 − 1·1·1
= 2 − 1 − 1
= 0 ✓

p(0,1,1)

= 1 + 1 − 1·1 − 0·1·1
= 2 − 1 − 0
= 1 ✓

p(0,0,0)

= 0 + 0 − 0 − 0
= 0 ✓

The full table is given in L03 (The Polynomial). Every row matches. This agreement is not a design coincidence — it is a consequence of the interpolation theorem.

Why Every Other Rule Fails

The proof of uniqueness is by exhaustive interpolation over GF(7). Any multilinear degree-≤3 polynomial over GF(7) in three variables is determined by 8 values: one for each binary corner (L,C,R) ∈ {0,1}³. The 8 constraints above pin each of those values exactly.

Any other choice of values for any corner would mean:

Either way, any competing rule is MDL-penalized relative to the unique survivor. The Lean 4 proof formalizes this by machine-verifying that the Lagrange interpolant is unique over GF(7) and equals C+R−CR−LCR:

Machine certification

 — zero sorry, zero custom axioms. The orbit data uniquely determines the polynomial.

Key Takeaways

See Also

Lean 4 proofs (ugp-lean)