Follow the source behind the result.

The registered Python source is shown as text. Displaying it does not execute the code.

"""v505 -- CELEST.WP5E.BETA.01: WP5e-beta of the research contract
CELEST.SEAM.01 -- "the equivariant anomaly ledger on twistor space", the
twistor-side milestone of WP5e (the BCOV/Costello-Li twistor uplift;
preregistered kill: the anomaly inflow demands a level != 1).  WP5e-alpha
(v502) anchored the CFT side (q^{-1/3} prefactor; k = 1 forced three ways).
Here the TWISTOR side is pushed as far as exact arithmetic reaches: the
Atiyah-Bott/Lefschetz fixed-point skeleton of the one-loop box anomaly
(Costello arXiv:2111.08879: box anomaly ~ Tr_adj X^4, GS axion coupling
lambda_g/(8 pi sqrt 3); Bittleston-Sharma-Skinner arXiv:2208.12701 / CMP 403:
axion eta in Omega^{0,1}(Omega^2_cl), O(-2)-valued Poisson structure,
Poisson-CS tower (+)_s O(2s-2)) evaluated EXACTLY on C^2/Z4 with the E8
mu4-glue monodromy (v492: inner element h = (2,2,2,2,2;0^4), graded dims
(60,64,60,64)).  Exact sympy/Fraction arithmetic throughout, no floats.

[E] S1. ATIYAH-BOTT LEDGER: denominators det(1 - g_j^{-1}|_{T_0 C^2}) =
      (1-i^{-j})(1-i^j) = (2, 4, 2) for j = 1, 2, 3; sum of reciprocals
      1/2+1/4+1/2 = 5/4 = (|Z4|^2-1)/12 (the A-series Dedekind sum identity);
      the origin is the only fixed point; equivariant characters
      tr_{g_j}(Ad E8) = sum_m i^{jm} dim g_m = (248, 0, -8, 0), TWO
      independent routes (graded dims / root sum); invariant average 60 =
      the D5+A3 carrier; Frobenius sum 61568 = 4(60^2+64^2+60^2+64^2);
      tr_{g_2} = -8 = 120 - 128 (so16 adjoint-minus-spinor supertrace).
[E] S2-S3. OKUBO PER SECTOR (the honest sharpening): the INVARIANT sector
      obeys Okubo exactly -- Q^{(0)} = 36<x,x>^2 = [5/(2(248+2))](K^{(0)})^2
      with 36 = lambda~^2_{e8} (v495 re-derived through the mu4-glue
      grading); the TWISTED sectors do NOT reduce to (tr_{g_j} X^2)^2:
      exact independent-quartic content T5: (-8, 32, -8), T3: (48, -64, 48).
      In the AB-weighted sum Q^{(1)}/2 + Q^{(2)}/4 + Q^{(3)}/2 the
      D5-quartic CANCELS exactly and the A3-quartic leaves the RIGID
      residual 32 T3; no nonzero reweighting of the twisted sectors alone
      is quartic-free, and on the class blocks the only quartic-free
      weighting is the uniform Z4 average (= the invariant projector); the
      graded GS-exchange ansatz (products K^{(a)}K^{(b)}, a+b = j mod 4)
      cancels sector 0 and is exactly rank-obstructed in sectors 1, 2, 3.
      Heterotic identities drop out on the way: Tr_120 = 8 tr^4 + 3(tr^2)^2,
      Tr_128 = -8 tr^4 + 6(tr^2)^2, Tr_E8 = 9(tr^2)^2 on the D5(+)A3 Cartan.
[E] S4. INDEX BRIDGE (the centrepiece): the per-character Atiyah-Bott sum
      f(m) = (1/4) sum_j (i^{jm}-1)/det_j = (0, -3/8, -1/2, -3/8) equals
      ch2(T_m) = -(C^{-1})_mm/2 EXACTLY for m = 1, 2, 3 -- the fixed-point
      ledger (E1) and the McKay/Kronheimer intersection ledger (E2) are the
      SAME arithmetic; glue-bundle ch2 defect -78 by BOTH routes (per-class
      contributions (-24, -30, -24); -30 = -h_vee, typed [C] coincidence);
      bridge identity sum_m h_A3(m) = 5/4 = sum_j 1/det_j; sphere pairing
      <alpha,h> integer for all 240 roots, per-class Killing sums
      (320, 320, 240, 320), total 1200 = 2 h_vee (h,h).  LEVEL DIALS,
      k = 1 GEOMETRIC: lattice current count (norm-2/k vectors in the glue
      cosets) = 240 at k = 1, EXACTLY ZERO at k = 2, 3, 4 (no adjoint
      closure); embedding residual 47k^2+219k-266 = (0, 360, 814, 1362);
      h(J;k) = k integer for ALL k -- integrality alone fixes nothing
      (honest, exactly as v502); v493 lockstep periods Pi = t(i-1)(1,i,i^2)
      => ONE scale => ONE level (det(A-1) = -4: no invariant flux vector).
[E] S5. E3 REFUTATION (mandatory, prominent): weight arithmetic on PT --
      the clock-invariant modulus a0 in O(8) generates a Hamiltonian of
      weight 2+8-2 = 8, i.e. h in O(2): a0 fills the BSS GRAVITON slot
      (s = 2), NOT the axion slot; the untwisted axion eta in Omega^2_cl is
      iota_V Omega with V in T (x) O(-4), Hamiltonian weight X = -2 (the
      s = 0 tower slot O(-2)); X = -2 != +2 = Y -- "the theory brings its
      own GS axion as a0" is REFUTED at the weight level (mismatch 4 =
      -deg K_PT = |mu4|, typed [C] as a coincidence).  What DOES match:
      the three H^2(ALE) classes carry the Coxeter eigencharacters
      {i, -1, -i} = exactly the three twisted ledger sectors j = 1, 2, 3
      (bijection, no invariant class: det(A-1) = -4) -- the three sphere
      2-forms fill the three TWISTED axion slots; slot count 1 tower O(-2)
      + 3 sphere classes = 4 = |mu4| (soft count, not a construction).
[E] S6. NEGATIVE CONTROLS: (a) diag(i,i) breaks the ledger FOUR ways (AB
      denominators complex, reciprocal sum 1/4 != 5/4, omega-weight -1 not
      CY, invariant ring needs 5 generators); (b) SO(16) glue: characters
      (120, 0, +120, 0) (wrong sign AND magnitude at j = 2), odd McKay
      nodes empty, ch2 defect -30 != -78, and bulk Okubo FAILS (T5 = 16,
      T3 = -32: independent quartic Casimir, no single-axion GS -- the
      v495 sl4-style control transplanted); (c) k = 2: current count 0,
      residual 360, c-prefactor -31/48 != -1/3, while naive integrality
      dials stay integer (integrality does not kill k = 2 -- the closure
      dial does, exactly as on the v502 CFT side).
[C] Dictionaries: CS(rho) mod 1 = discriminant weight; 32 = (h,h)+(h',h');
      -30 = -h_vee; the Omega-contraction identification of the axion slot.
[O] HONEST FENCE -- what stays open: the global BCOV/Kodaira-Spencer theory
      on PT/Z4 (quantisation, non-fixed-point contributions, the twisted
      closed-string measure), the O(-2) bulk-axion construction, the
      Costello-Paquette-Sharma level-from-flux mechanism, and the pairing
      of the rigid 32 T3 twisted residual with the three sphere axions.
      VERDICT B; the preregistered kill ("inflow demands a level != 1")
      does NOT fire on the equivariant skeleton -- all computable
      twistor-side pieces are consistent with k = 1 and kill k = 2, 3, 4
      via the closure dial; what the skeleton cannot decide is the actual
      BCOV one-loop coefficient on PT/Z4.  WP5e: alpha (v502) + beta (this
      module) executed; SEAM.EQUIV.01 untouched; no marker moves anywhere.

Status: [E] exact equivariant/index/level/weight arithmetic (sympy +
Fraction); [C] the dictionaries above; [O] the global BCOV quantisation
(WP5e proper).  Python; Wolfram-mirrored (denominators/characters, index
bridge f(m) = ch2, -78 defect, current counts, slot weights), counted per
GATE.WOLFRAM.02.  Discovery provenance: experiments/tfpt-discovery/
celestial_seam_wp5e_beta_equivariant_inflow_probe.py (2026-07-21, 47/47)."""
from fractions import Fraction as F
from itertools import combinations, product

import sympy as sp

from tfpt_constants import check, summary, reset, g_car, N_fam, rankE8

MU4 = N_fam + 1             # 4 = |mu4|, the seam clock order
RANK = rankE8               # 8
H_VEE = 2 * N_fam * g_car   # 30 = |Z2| N_fam g_car (v495 anchor decomposition)
HALF = F(1, 2)


def fmt(xs):
    return "(" + ", ".join(str(x) for x in xs) + ")"


# ---------------------------------------------------------------------------
# E8 roots in D5 (+) A3 glue coordinates (v128/v492/v502 construction)
# ---------------------------------------------------------------------------
def build_glue_roots():
    d5_roots, d5_v = [], []
    for i, j in combinations(range(5), 2):
        for si in (1, -1):
            for sj in (1, -1):
                v = [F(0)] * 5
                v[i], v[j] = F(si), F(sj)
                d5_roots.append(tuple(v))
    for i in range(5):
        for s in (1, -1):
            v = [F(0)] * 5
            v[i] = F(s)
            d5_v.append(tuple(v))
    d5_s, d5_c = [], []
    for signs in product((1, -1), repeat=5):
        v = tuple(HALF * s for s in signs)
        (d5_s if signs.count(-1) % 2 == 0 else d5_c).append(v)
    a3_roots = []
    for i in range(4):
        for j in range(4):
            if i != j:
                v = [F(0)] * 4
                v[i], v[j] = F(1), F(-1)
                a3_roots.append(tuple(v))

    def wclass(k):
        out = []
        for sub in combinations(range(4), k):
            v = [F(-k, 4)] * 4
            for i in sub:
                v[i] += 1
            out.append(tuple(v))
        return out

    z5, z4 = tuple([F(0)] * 5), tuple([F(0)] * 4)
    roots = {}
    for r in d5_roots:
        roots[r + z4] = 0
    for r in a3_roots:
        roots[z5 + r] = 0
    for d in d5_s:
        for w in wclass(1):
            roots[d + w] = 1
    for d in d5_v:
        for w in wclass(2):
            roots[d + w] = 2
    for d in d5_c:
        for w in wclass(3):
            roots[d + w] = 3
    return roots


# ---------------------------------------------------------------------------
# exact polynomial machinery on the D5 (+) A3 Cartan (8 parameters)
# ---------------------------------------------------------------------------
X5 = sp.symbols('x1:6')                      # D5 Cartan x1..x5
Y3 = sp.symbols('y1:4')                      # A3 Cartan (sum-zero coords)
Y4 = (*Y3, -sum(Y3))                         # y4 = -(y1+y2+y3)
ALLV = (*X5, *Y3)

S5 = sp.expand(sum(v ** 2 for v in X5))
S3 = sp.expand(sum(v ** 2 for v in Y4))
P1 = sp.expand(S5 ** 2)
P2 = sp.expand(S5 * S3)
P3 = sp.expand(S3 ** 2)
T5 = sp.expand(sum(v ** 4 for v in X5))
T3 = sp.expand(sum(v ** 4 for v in Y4))
BASIS4 = [("P1", P1), ("P2", P2), ("P3", P3), ("T5", T5), ("T3", T3)]


def lin_form(alpha):
    """<alpha, x> as a sympy linear form (A3 block in sum-zero coordinates)."""
    e = sum(sp.Rational(alpha[i].numerator, alpha[i].denominator) * X5[i]
            for i in range(5))
    e += sum(sp.Rational(alpha[5 + i].numerator, alpha[5 + i].denominator)
             * Y4[i] for i in range(4))
    return sp.expand(e)


def poly_dict(expr):
    return sp.Poly(sp.expand(expr), *ALLV, domain='QQ').as_dict()


def decompose(target, basis_exprs):
    """Exact coefficients of target in span(basis) or None if not in span."""
    dicts = [poly_dict(b) for b in basis_exprs]
    tdict = poly_dict(target)
    monos = sorted(set().union(tdict.keys(), *[d.keys() for d in dicts]))
    A = sp.Matrix([[d.get(m, 0) for d in dicts] for m in monos])
    b = sp.Matrix([tdict.get(m, 0) for m in monos])
    try:
        sol, params = A.gauss_jordan_solve(b)
    except ValueError:
        return None
    if len(params) > 0:
        sol = sol.subs({p: 0 for p in params})
    if sp.expand(target - sum(sol[i] * basis_exprs[i]
                              for i in range(len(basis_exprs)))) != 0:
        return None
    return list(sol)


def in_span(target, basis_exprs):
    return decompose(target, basis_exprs) is not None


def power_sums_by_class(roots):
    """K_m, Q_m = sum over class-m roots of <alpha,x>^2 resp. ^4 (exact)."""
    K = {m: sp.Poly(0, *ALLV, domain='QQ') for m in range(4)}
    Q = {m: sp.Poly(0, *ALLV, domain='QQ') for m in range(4)}
    for alpha, c in roots.items():
        p = sp.Poly(lin_form(alpha), *ALLV, domain='QQ')
        p2 = p ** 2
        K[c] = K[c] + p2
        Q[c] = Q[c] + p2 ** 2
    return ({m: K[m].as_expr() for m in range(4)},
            {m: Q[m].as_expr() for m in range(4)})


# ---------------------------------------------------------------------------
# S1 -- E1(i)+(ii): Atiyah-Bott denominators and equivariant characters
# ---------------------------------------------------------------------------
def section1(roots):
    print("  -- S1 (E1): Atiyah-Bott denominators + equivariant E8 characters")
    I = sp.I

    dens = [sp.expand((1 - I ** (-j)) * (1 - I ** j)) for j in (1, 2, 3)]
    check("S1.1 [AB DENOMINATORS] det(1 - g_j^{-1}|_{T0}) = "
          "(1 - i^{-j})(1 - i^j) = %s = (2, 4, 2) for j = 1, 2, 3 -- real "
          "and positive (the deck diag(i, i^{-1}) is in SU(2))" % fmt(dens),
          dens == [2, 4, 2])

    recip = sum(sp.Rational(1, d) for d in dens)
    check("S1.2 [DEDEKIND SUM] sum_j 1/det_j = 1/2 + 1/4 + 1/2 = %s = 5/4 = "
          "(|Z4|^2 - 1)/12 (the A-series identity sum_j 1/(4 sin^2(pi j/n)) "
          "= (n^2-1)/12 at n = 4)" % recip,
          recip == sp.Rational(5, 4)
          and sp.Rational(MU4 ** 2 - 1, 12) == sp.Rational(5, 4))

    z = sp.symbols('z')
    fixed_ok = all(sp.solve(sp.Eq(I ** j * z, z), z) == [0] for j in (1, 2, 3))
    check("S1.3 [FIXED LOCUS] i^j z = z forces z = 0 for j = 1, 2, 3 on both "
          "coordinates: the origin is the ONLY fixed point (det != 0 <=> "
          "isolated), the Lefschetz sum is a single-point ledger", fixed_ok)

    counts = [sum(1 for c in roots.values() if c == m) for m in range(4)]
    dims = [counts[0] + RANK] + counts[1:]
    norm_ok = all(sum(x * x for x in r) == 2 for r in roots)
    check("S1.4 [v492 REPLICATION] 240 roots, all norm 2, split %s = "
          "(52,64,60,64), graded dims with Cartan %s = (60,64,60,64)"
          % (counts, dims),
          len(roots) == 240 and norm_ok and counts == [52, 64, 60, 64]
          and dims == [60, 64, 60, 64])

    chi_dims = [sp.expand(sum(I ** (j * m) * dims[m] for m in range(4)))
                for j in range(4)]
    chi_roots = [sp.expand(RANK + sum(I ** (j * c) for c in roots.values()))
                 for j in range(4)]
    check("S1.5 [CHARACTERS, TWO ROUTES] tr_{g_j}(Ad E8) = sum_m i^{jm} dim "
          "g_m = %s = (248, 0, -8, 0); the root-sum route (rank + sum over "
          "240 roots) gives the same -- the glue monodromy has REAL "
          "character with a single nonzero twisted value at j = 2"
          % fmt(chi_dims),
          chi_dims == [248, 0, -8, 0] and chi_roots == chi_dims)

    inv = sp.Rational(sum(chi_dims), 4)
    check("S1.6 [INVARIANT AVERAGE] (1/4) sum_j tr_{g_j}(Ad) = (248 + 0 - 8 "
          "+ 0)/4 = %s = 60 = dim(D5 + A3 carrier + Cartan) = the v492 zero-"
          "mode algebra: the fixed-point characters know the carrier" % inv,
          inv == 60 and inv == dims[0])

    frob = sp.expand(sum(c * sp.conjugate(c) for c in chi_dims))
    check("S1.7 [FROBENIUS] sum_j |chi(g_j)|^2 = %s = 61568 = |Z4| * "
          "(60^2 + 64^2 + 60^2 + 64^2) = 4 * 15392: the multiplicity vector "
          "(60,64,60,64) is consistent as a Z4-module decomposition" % frob,
          frob == 61568
          and 4 * (60 ** 2 + 64 ** 2 + 60 ** 2 + 64 ** 2) == 61568)

    check("S1.8 [j = 2 SUPERTRACE] tr_{g_2}(Ad) = -8 = 120 - 128 = dim so16 "
          "- dim spinor: the halfway sector reads the v125 SO(16)_1 step; "
          "|tr_{g_2}| = 8 = rank = 2|mu4| (the seam winding integer)",
          chi_dims[2] == -8 and 120 - 128 == -8 and abs(-8) == RANK)
    return dims, chi_dims, dens


# ---------------------------------------------------------------------------
# S2 -- E1(iii): twisted quadratic forms
# ---------------------------------------------------------------------------
def section2(K):
    print("  -- S2 (E1): twisted quadratic forms K^{(j)}")
    I = sp.I

    coefs = {}
    ok_tab = True
    for m in range(4):
        c = decompose(K[m], [S5, S3])
        coefs[m] = c
        ok_tab &= c is not None
    tab = {m: tuple(coefs[m]) for m in range(4)}
    print("     K_m = a_m s5 + b_m s3:  %s"
          % str({m: fmt(tab[m]) for m in range(4)}))
    check("S2.1 [CLASS TABLE] (a_m, b_m) = ((16,8), (16,16), (12,20), "
          "(16,16)): class 0 = D5+A3 adjoints (2h_vee = 16, 8); class 1 = 3 "
          "= (16, s) x (4): 4*4 s5 + 16*1 s3; class 2 = (10, v) x (6): "
          "6*2 s5 + 10*2 s3; total = 60(s5+s3) = Killing 2 h_vee(E8) <x,x>",
          ok_tab and tab == {0: (16, 8), 1: (16, 16), 2: (12, 20),
                             3: (16, 16)}
          and sp.expand(sum(K.values()) - 60 * (S5 + S3)) == 0)

    Ktw = {j: sp.expand(sum(I ** (j * m) * K[m] for m in range(4)))
           for j in range(4)}
    check("S2.2 [TWISTED QUADRATICS] K^{(0)} = 60(s5+s3); K^{(1)} = K^{(3)} "
          "= 4 s5 - 12 s3 REAL (K_1 = K_3 kills the imaginary part exactly); "
          "K^{(2)} = -4(s5+s3) = -|mu4| <x,x>",
          sp.expand(Ktw[0] - 60 * (S5 + S3)) == 0
          and sp.expand(Ktw[1] - (4 * S5 - 12 * S3)) == 0
          and sp.expand(Ktw[3] - Ktw[1]) == 0
          and sp.expand(Ktw[2] + 4 * (S5 + S3)) == 0)

    so16_ad = sp.expand(K[0] + K[2])
    so16_sp = sp.expand(K[1] + K[3])
    check("S2.3 [THREE ROUTES TO -4] K_0 + K_2 = 28<x,x> (2h_vee(so16) = 28) "
          "and K_1 + K_3 = 32<x,x> (spinor-128 index 32); K^{(2)} = 28 - 32 "
          "= -4: the twisted quadratic is the so16 adjoint-minus-spinor "
          "index deficit -- matching tr_{g_2}(Ad) = -8 = 120 - 128 one "
          "level up",
          sp.expand(so16_ad - 28 * (S5 + S3)) == 0
          and sp.expand(so16_sp - 32 * (S5 + S3)) == 0
          and 28 - 32 == -4)
    return Ktw


# ---------------------------------------------------------------------------
# S3 -- E1(iii)+(iv): twisted quartics, Okubo per sector, AB sum
# ---------------------------------------------------------------------------
def section3(Q, Ktw):
    print("  -- S3 (E1): twisted quartic densities and the Okubo ledger")
    I = sp.I
    basis_exprs = [b for _, b in BASIS4]

    dec = {}
    ok_dec = True
    print("     class-block decomposition (P1, P2, P3, T5, T3):")
    for m in range(4):
        c = decompose(Q[m], basis_exprs)
        dec[m] = c
        ok_dec &= c is not None
        print("       Q_%d = %s" % (m, fmt(c)))
    expect = {0: [12, 0, 6, 4, 8], 1: [12, 24, 0, -8, 16],
              2: [0, 24, 30, 12, -40], 3: [12, 24, 0, -8, 16]}
    check("S3.1 [CLASS QUARTICS EXACT] all four class blocks decompose in "
          "the 5-basis {P1, P2, P3, T5, T3} with ZERO residual: Q_0 = "
          "12P1+6P3+4T5+8T3, Q_1 = Q_3 = 12P1+24P2-8T5+16T3, Q_2 = "
          "24P2+30P3+12T5-40T3 (the D5 and A3 independent quartics T5, T3 "
          "are PRESENT block by block)",
          ok_dec and all([sp.Rational(x) for x in dec[m]]
                         == [sp.Rational(x) for x in expect[m]]
                         for m in range(4)))

    tr16 = sp.expand(sum((sum(sp.Rational(s[i], 2) * X5[i] for i in range(5)))
                         ** 4
                         for s in product((1, -1), repeat=5)
                         if s.count(-1) % 2 == 0))
    tr6 = sp.expand(sum((Y4[i] + Y4[j]) ** 4
                        for i, j in combinations(range(4), 2)))
    tr4 = sp.expand(sum(w ** 4 for w in Y4))
    check("S3.2 [BLOCK IDENTITIES] spinor-16 of D5: tr_16(x^4) = 3P1 - 2T5; "
          "Lambda^2-6 of A3: tr_6(x^4) = 3P3 - 4T3; fundamental-4: tr_4(x^4) "
          "= T3 -- the classic trace identities re-derived from the explicit "
          "weight systems",
          sp.expand(tr16 - (3 * P1 - 2 * T5)) == 0
          and sp.expand(tr6 - (3 * P3 - 4 * T3)) == 0
          and sp.expand(tr4 - T3) == 0)

    Qtw = {j: sp.expand(sum(I ** (j * m) * Q[m] for m in range(4)))
           for j in range(4)}
    ok_real = (sp.expand(Qtw[1] - (Q[0] - Q[2])) == 0
               and sp.expand(Qtw[3] - Qtw[1]) == 0
               and sp.expand(Qtw[2] - (Q[0] - 2 * Q[1] + Q[2])) == 0)
    d1 = decompose(Qtw[1], basis_exprs)
    d2 = decompose(Qtw[2], basis_exprs)
    print("     twisted quartics (P1, P2, P3, T5, T3):")
    print("       Q^(1) = Q^(3) = %s" % fmt(d1))
    print("       Q^(2)         = %s" % fmt(d2))
    check("S3.3 [TWISTED QUARTICS] Q^{(1)} = Q^{(3)} = Q_0 - Q_2 = "
          "12P1 - 24P2 - 24P3 - 8T5 + 48T3 and Q^{(2)} = Q_0 - 2Q_1 + Q_2 = "
          "-12P1 - 24P2 + 36P3 + 32T5 - 64T3 (all REAL because Q_1 = Q_3)",
          ok_real
          and d1 == [12, -24, -24, -8, 48]
          and d2 == [-12, -24, 36, 32, -64])

    okubo0 = sp.expand(Qtw[0] - sp.Rational(5, 2 * (248 + 2)) * Ktw[0] ** 2)
    check("S3.4 [OKUBO, INVARIANT SECTOR] Q^{(0)} = 36 <x,x>^2 = "
          "[5/(2(248+2))] (K^{(0)})^2 = (1/100)(60<x,x>)^2 EXACTLY -- T5 and "
          "T3 vanish in the invariant sum (v495 Okubo re-derived through the "
          "mu4-glue grading; coefficient 36 = lambda~_e8^2)",
          okubo0 == 0
          and sp.expand(Qtw[0] - 36 * (S5 + S3) ** 2) == 0)

    tr120 = sp.expand(Q[0] + Q[2])
    tr128 = sp.expand(Q[1] + Q[3])
    t16v2 = sp.expand(2 * S5 + 2 * S3)
    t16v4 = sp.expand(2 * T5 + tr6)
    het_ok = (sp.expand(tr120 - (8 * t16v4 + 3 * t16v2 ** 2)) == 0
              and sp.expand(tr128 - (-8 * t16v4 + 6 * t16v2 ** 2)) == 0
              and sp.expand(tr120 + tr128 - 9 * t16v2 ** 2) == 0)
    check("S3.5 [HETEROTIC ANCHORS] on the D5(+)A3 Cartan (so16 vector = "
          "(10,1)+(1,6)): Tr_120 X^4 = 8 tr X^4 + 3 (tr X^2)^2, Tr_128 X^4 "
          "= -8 tr X^4 + 6 (tr X^2)^2, sum Tr_E8 X^4 = 9 (tr X^2)^2 -- the "
          "classic Green-Schwarz identities drop out of the glue ledger "
          "exactly", het_ok)

    ok_fail = all(not in_span(Qtw[j], [sp.expand(Ktw[j] ** 2)])
                  for j in (1, 2, 3))
    check("S3.6 [HONEST SHARPENING: PER-SECTOR OKUBO FAILS] Q^{(j)} is NOT "
          "proportional to (K^{(j)})^2 for j = 1, 2, 3: the twisted sectors "
          "carry irreducible independent-quartic content T5: (-8, 32, -8), "
          "T3: (48, -64, 48) -- the naive sector-by-sector Okubo reduction "
          "does NOT hold; only the invariant sector is a perfect square",
          ok_fail)

    Afix = sp.expand(Qtw[1] / 2 + Qtw[2] / 4 + Qtw[3] / 2)
    dA = decompose(Afix, basis_exprs)
    print("     AB fixed-point sum  A_fix = Q^(1)/2 + Q^(2)/4 + Q^(3)/2 = %s"
          % fmt(dA))
    check("S3.7 [AB LEDGER SUM] A_fix = sum_j Q^{(j)}/det_j = 9P1 - 30P2 - "
          "15P3 + 32T3 EXACTLY: the D5-quartic CANCELS in the Atiyah-Bott-"
          "weighted sum (T5: -8/2 + 32/4 - 8/2 = 0), the A3-quartic leaves "
          "the exact residual 32 T3 (32 = (h,h) + (h',h') = the inner-"
          "element norm sum -- numerical coincidence, typed [C]); the fixed-"
          "point ledger is NOT a perfect square",
          dA == [9, -30, -15, 0, 32])

    w0, w1, w2 = sp.symbols('w0 w1 w2')
    t5vec = {1: -8, 2: 32, 3: -8}
    t3vec = {1: 48, 2: -64, 3: 48}
    sols_fix = sp.solve([2 * w1 * t5vec[1] + w2 * t5vec[2],
                         2 * w1 * t3vec[1] + w2 * t3vec[2]],
                        [w1, w2], dict=True)
    only_zero = sols_fix == [{w1: 0, w2: 0}]
    t5c = [4, -8, 12]   # T5 of Q_0, Q_1(=Q_3), Q_2
    t3c = [8, 16, -40]
    solsb = sp.solve([w0 * t5c[0] + 2 * w1 * t5c[1] + w2 * t5c[2],
                      w0 * t3c[0] + 2 * w1 * t3c[1] + w2 * t3c[2]],
                     [w0, w1], dict=True)
    line_ok = (len(solsb) == 1
               and sp.simplify(solsb[0][w0] - w2) == 0
               and sp.simplify(solsb[0][w1] - w2) == 0)
    check("S3.8 [RIGIDITY OF THE RESIDUAL] no nonzero weighting of the "
          "three twisted sectors alone kills both quartics (unique solution "
          "w = 0); including the invariant block, the ONLY quartic-free "
          "weighting is the uniform one (w0 = w1 = w2 = w3) = the full Z4 "
          "orbifold average: the AB weights (1/2, 1/4, 1/2) are forced by "
          "the geometry and necessarily leave the 32 T3 residual",
          only_zero and line_ok)

    ok_span0 = in_span(Qtw[0], [sp.expand(Ktw[0] ** 2),
                                sp.expand(Ktw[1] * Ktw[3]),
                                sp.expand(Ktw[2] ** 2)])
    span1 = [sp.expand(Ktw[0] * Ktw[1]), sp.expand(Ktw[2] * Ktw[3])]
    span2 = [sp.expand(Ktw[0] * Ktw[2]), sp.expand(Ktw[1] ** 2),
             sp.expand(Ktw[3] ** 2)]
    ok_span_fail = ((not in_span(Qtw[1], span1))
                    and (not in_span(Qtw[2], span2)))
    check("S3.9 [GRADED GS-EXCHANGE TEST] the sector-0 density lies in the "
          "span of graded vertex products {K^{(a)} K^{(b)}: a+b = 0 mod 4} "
          "(cancellable by single-axion exchange, = Costello's mechanism); "
          "sectors 1 and 2 do NOT lie in their graded spans (exact rank "
          "obstruction): a K-quadratic-vertex GS mechanism alone cannot "
          "cancel the twisted fixed-point ledger -- twisted-sector closed "
          "strings must carry the rest [O]",
          ok_span0 and ok_span_fail)
    return Qtw, dec


# ---------------------------------------------------------------------------
# S4 -- E2: McKay/Kronheimer index bridge and the level dial
# ---------------------------------------------------------------------------
def d5_coset_norms(maxnorm):
    out = {k: {} for k in ('0', 'v', 's', 'c')}
    for v in product(range(-2, 3), repeat=5):
        n = F(sum(x * x for x in v))
        if n > maxnorm:
            continue
        key = '0' if sum(v) % 2 == 0 else 'v'
        out[key][n] = out[key].get(n, 0) + 1
    rng_half = [F(k, 2) for k in (-3, -1, 1, 3)]
    for v in product(rng_half, repeat=5):
        n = sum(x * x for x in v)
        if n > maxnorm:
            continue
        key = 's' if (sum(v) - F(5, 2)) % 2 == 0 else 'c'
        out[key][n] = out[key].get(n, 0) + 1
    return out


def a3_coset_norms(maxnorm):
    out = {k: {} for k in range(4)}
    for k in range(4):
        shift = F(-k, 4)
        rng = [m + shift for m in range(-3, 4)]
        for v in product(rng, repeat=4):
            if sum(v) != 0:
                continue
            n = sum(x * x for x in v)
            if n > maxnorm:
                continue
            out[k][n] = out[k].get(n, 0) + 1
    return out


def section4(roots, dims, chi_dims, dens):
    print("  -- S4 (E2): McKay/Kronheimer index bridge and the level dial")
    I = sp.I

    C = sp.Matrix([[2, -1, 0], [-1, 2, -1], [0, -1, 2]])
    Cinv = C.inv()
    ch2T = [-Cinv[m, m] / 2 for m in range(3)]
    hA3 = [sp.Rational(j * (4 - j), 8) for j in (1, 2, 3)]
    check("S4.1 [INTERSECTION / c1 ARITHMETIC] intersection form = -C(A3), "
          "det C = %s = 4 = |Z4|; tautological bundles deg(T_m|S_l) = "
          "delta_ml give int c1(T_m)^2 = -(C^{-1})_mm and ch2(T_m) = %s = "
          "(-3/8, -1/2, -3/8) = MINUS the A3 discriminant weights j(4-j)/8 "
          "= the v502 sector Casimir energies: the geometric face of the "
          "discriminant form" % (C.det(), fmt(ch2T)),
          C.det() == 4
          and ch2T == [sp.Rational(-3, 8), sp.Rational(-1, 2),
                       sp.Rational(-3, 8)]
          and [-c for c in ch2T] == hA3)

    adj = sp.zeros(4, 4)
    for m in range(4):
        adj[m, (m + 1) % 4] += 1
        adj[m, (m - 1) % 4] += 1
    affine = 2 * sp.eye(4) - adj
    ker = affine.nullspace()
    check("S4.2 [McKAY] Q (x) chi_m = chi_{m+1} + chi_{m-1}: adjacency = "
          "the affine A3 4-cycle; 2*1 - Adj = affine Cartan with det %s = 0 "
          "and kernel (1,1,1,1) = tautological ranks / regular-rep "
          "multiplicities (Kronheimer quiver data)" % affine.det(),
          affine.det() == 0 and len(ker) == 1
          and list(ker[0] / ker[0][0]) == [1, 1, 1, 1])

    fvals = [sp.nsimplify(sp.expand(
        sp.Rational(1, 4) * sum((I ** (j * m) - 1) / dens[j - 1]
                                for j in (1, 2, 3)))) for m in range(4)]
    print("     per-character AB sums f(m) = (1/4) sum_j (i^{jm}-1)/det_j "
          "= %s" % fmt(fvals))
    check("S4.3 [THE INDEX BRIDGE, per character] f(m) = %s = (0, -3/8, "
          "-1/2, -3/8) = ch2(T_m) EXACTLY for m = 1, 2, 3: the Atiyah-Bott "
          "fixed-point sum per Z4 character equals the intersection-form "
          "charge of the McKay tautological bundle -- the E1 ledger and the "
          "E2 ledger are the SAME arithmetic" % fmt(fvals),
          fvals == [0] + ch2T)

    defect_ab = sp.nsimplify(sp.expand(
        sp.Rational(1, 4) * sum((chi_dims[j] - 248) / dens[j - 1]
                                for j in (1, 2, 3))))
    defect_geo = sum(dims[m + 1] * ch2T[m] for m in range(3))
    contrib = [dims[m + 1] * ch2T[m] for m in range(3)]
    print("     glue-bundle ch2 defect: AB route %s, intersection route %s, "
          "per class %s" % (defect_ab, defect_geo, fmt(contrib)))
    check("S4.4 [GLUE BUNDLE DEFECT] equivariant ch2 defect of Ad = (+)_m "
          "g_m (x) T_m: AB route (1/4) sum_j (chi(g_j) - 248)/det_j = %s "
          "and intersection route sum_m dim g_m ch2(T_m) = 64(-3/8) + "
          "60(-1/2) + 64(-3/8) = %s: both EXACTLY -78, an INTEGER (the E8 "
          "glue extension carries integral charge); per-class contributions "
          "(-24, -30, -24) -- the middle one -30 = -h_vee(E8) (numerical "
          "coincidence, typed [C])" % (defect_ab, defect_geo),
          defect_ab == -78 and defect_geo == -78
          and contrib == [-24, -30, -24])

    bridge = sum(hA3)
    check("S4.5 [BRIDGE IDENTITY] sum_m h_A3(m) = 3/8 + 1/2 + 3/8 = %s = "
          "5/4 = sum_j 1/det_j (S1.2): the total discriminant weight IS the "
          "Dedekind sum of the AB denominators -- E1(i) and E2 talk to each "
          "other" % bridge, bridge == sp.Rational(5, 4))

    h9 = (F(2),) * 5 + (F(0),) * 4
    ip = lambda a, b: sum(p * q for p, q in zip(a, b))
    pair_int = all(ip(r, h9).denominator == 1 for r in roots)
    per_class = {m: sum(F(ip(r, h9)) ** 2
                        for r, c in roots.items() if c == m)
                 for m in range(4)}
    tot = sum(per_class.values())
    check("S4.6 [SPHERE PAIRING] <alpha,h> INTEGER for all 240 roots (the "
          "inner grading is integral -- the root line bundles have integer "
          "degrees on the spheres); per-class Killing sums K_m(h) = %s = "
          "(320, 320, 240, 320), total %s = 2 h_vee (h,h) = 60 * 20 = 1200 "
          "(the v495 Killing identity evaluated at x = h)"
          % (fmt(per_class.values()), tot),
          pair_int and tot == 1200
          and per_class == {0: 320, 1: 320, 2: 240, 3: 320})

    ks = sp.symbols('k', positive=True)
    hs = sp.together((ks ** 2 * sp.Rational(5, 4) + 2 * ks * 5)
                     / (2 * (ks + 8)))
    ha1 = sp.together((ks ** 2 * sp.Rational(3, 4)
                       + 2 * ks * sp.Rational(3, 2)) / (2 * (ks + 4)))
    ok_forms = (sp.simplify(hs - 5 * ks / 8) == 0
                and sp.simplify(ha1 - 3 * ks / 8) == 0)
    hJ = {k: F(5 * k, 8) + F(3 * k, 8) for k in range(1, 5)}
    check("S4.7 [LEVEL DIAL i, honest] affine closed forms h(k om_s; D5) = "
          "5k/8, h(k om_1; A3) = 3k/8, glue diagonal h(J; k) = k: INTEGER "
          "for ALL k -- integrality alone fixes nothing (v502 honest "
          "sharpening replicated); the current condition h(J) = 1 forces "
          "k = 1", ok_forms and all(hJ[k] == k for k in range(1, 5))
          and [k for k in range(1, 5) if hJ[k] == 1] == [1])

    d5 = d5_coset_norms(4)
    a3 = a3_coset_norms(4)
    pairing = {0: ('0', 0), 1: ('s', 1), 2: ('v', 2), 3: ('c', 3)}
    counts_k = {}
    for k in (1, 2, 3, 4):
        target = F(2, k)
        tot_k = 0
        for j, (dk, ak) in pairing.items():
            for n5, c5 in d5[dk].items():
                for n3, c3 in a3[ak].items():
                    if n5 + n3 == target:
                        tot_k += c5 * c3
        counts_k[k] = tot_k
    check("S4.8 [LEVEL DIAL ii, lattice] current-slot vectors (norm 2/k in "
          "the glue cosets, i.e. norm 2 after the level-k rescaling): "
          "counts %s for k = 1..4: 240 root currents at k = 1 (+ 8 Cartan = "
          "248 = adjoint closure), ZERO at k = 2, 3, 4 -- the seam lattice "
          "realises the (E8)_1 extension at k = 1 ONLY (geometric "
          "counterpart of the v502 current condition)" % str(counts_k),
          counts_k == {1: 240, 2: 0, 3: 0, 4: 0})

    resid = {k: 47 * k * k + 219 * k - 266 for k in (1, 2, 3, 4)}
    check("S4.9 [LEVEL DIAL iii] conformal-embedding residual 47k^2 + 219k "
          "- 266 = %s for k = 1..4: zero exactly at k = 1 (v502 S5.4 "
          "replication)" % str(resid),
          resid == {1: 0, 2: 360, 3: 814, 4: 1362})

    t = sp.symbols('t', positive=True)
    roots4 = [t, I * t, -t, -I * t]
    Pi = [sp.simplify(roots4[(j + 1) % 4] - roots4[j]) for j in range(3)]
    dft = sp.Matrix([1, I, I ** 2])
    ok_pi = sp.simplify(sp.Matrix(Pi) - t * (I - 1) * dft) == sp.zeros(3, 1)
    mods = [sp.simplify(sp.Abs(p)) for p in Pi]
    A = sp.Matrix([[0, 0, -1], [1, 0, -1], [0, 1, -1]])
    check("S4.10 [ONE SCALE => ONE LEVEL] v493 replication: periods Pi = "
          "t(i-1)(1, i, i^2), all three sphere volumes = sqrt(2) t "
          "(lockstep); det(A - 1) = %s != 0: the Coxeter clock admits NO "
          "invariant flux vector on H2 -- a level assignment cannot vary "
          "over the three spheres, k is a single scalar; with S4.7-S4.9 the "
          "geometric side reproduces k = 1" % (A - sp.eye(3)).det(),
          ok_pi and all(sp.simplify(m - sp.sqrt(2) * t) == 0 for m in mods)
          and (A - sp.eye(3)).det() == -4)


# ---------------------------------------------------------------------------
# S5 -- E3: the axion slot (weight/degree arithmetic)
# ---------------------------------------------------------------------------
def section5():
    print("  -- S5 (E3): the axion slot -- exact weight arithmetic")
    I = sp.I

    w_vol = 2 + 1 + 1
    check("S5.1 [VOLUME FORM] Omega = <lambda dlambda> dmu^1 dmu^2 has "
          "homogeneity 2 + 1 + 1 = %d = 4 = -deg K_PT (K_PT = O(-4)): the "
          "canonical O(4)-twist of twistor space" % w_vol,
          w_vol == 4)

    w_omega_a3 = 4 + 4 + 2 - 8
    w_omega_eh = 2 + 2 + 2 - 4
    check("S5.2 [FIBRE SYMPLECTIC FORM] Poincare residue omega = "
          "Res(dX dY dZ / F): A3 family (X, Y in O(4), Z in O(2), F = XY - "
          "Z^4 - a0 in O(8)): weight 4+4+2-8 = %d; EH: 2+2+2-4 = %d; both "
          "= 2 = the O(2)-valued fibre symplectic form of the nonlinear "
          "graviton -- the graviton condition is k-independent"
          % (w_omega_a3, w_omega_eh),
          w_omega_a3 == 2 and w_omega_eh == 2)

    wts = {name: 8 - 2 * m for name, m in
           (('a0', 0), ('a1', 1), ('a2', 2))}
    check("S5.3 [MODULI WEIGHTS] a_m Z^m demands a_m in O(8 - 2m): (a0, a1, "
          "a2) in (O(8), O(6), O(4)); the clock-invariant modulus a0 sits "
          "at the TOP weight O(8) = O(2|Z4|): the Z2 -> Z4 doubling O(4) -> "
          "O(8) (EH c^2 is a quartic section, BHS) with ratio 8/4 = 2 = "
          "|Z4|/|Z2| (v493 S3.8 replication)",
          wts == {'a0': 8, 'a1': 6, 'a2': 4}
          and F(8, 4) == F(MU4, 2))

    Zs, a0s, a1s, a2s, a3s = sp.symbols('Z a0 a1 a2 a3')
    P = Zs ** 4 + a3s * Zs ** 3 + a2s * Zs ** 2 + a1s * Zs + a0s
    Pnew = sp.expand(P.subs(Zs, Zs / I))
    wts_clock = {s: sp.simplify(Pnew.coeff(Zs, k) / s)
                 for s, k in ((a3s, 3), (a2s, 2), (a1s, 1), (a0s, 0))}
    check("S5.4 [CLOCK WEIGHTS] the mu4 clock acts on (a3, a2, a1, a0) with "
          "weights (i, -1, -i, 1) = the regular representation; a0 is the "
          "UNIQUE clock-invariant modulus (v493 S2.11 replication)",
          wts_clock[a3s] == I and wts_clock[a2s] == -1
          and wts_clock[a1s] == -I and wts_clock[a0s] == 1)

    check("S5.5 [a0 = GRAVITON SLOT] the Poisson bracket maps O(a) x O(b) "
          "-> O(a+b-2) (O(-2)-valued Poisson structure, BSS), so the "
          "Hamiltonian generating the a0-deformation ({h, .} with {h, F} "
          "matching delta F = delta a0 in O(8)) has weight 2 + 8 - 2 = 8, "
          "i.e. h in O(2): a0 is a zero mode of the BSS GRAVITON field h in "
          "Omega^{0,1}(O(2)) = the s = 2 slot of the Poisson-CS tower "
          "(+)_s O(2s-2)",
          2 + 8 - 2 == 8 and 2 * 2 - 2 == 2)

    w_axion_vec = 0 - 4
    w_axion_ham = w_axion_vec + 2
    check("S5.6 [AXION SLOT] the BSS/Costello axion eta lives in "
          "Omega^{0,1}(Omega^2_cl) UNTWISTED; via Omega-contraction an "
          "untwisted 2-form is iota_V Omega with V in T (x) O(-4), i.e. "
          "Hamiltonian weight -4 + 2 = %d: the axion is the s = 0 slot "
          "O(2s-2) = O(-2) of the tower (BSS: the tower spectrum 'includes "
          "a scalar'; H^1(PT, O(-2)) = massless scalar = the 4th-order "
          "axion)" % w_axion_ham,
          w_axion_ham == -2 and 2 * 0 - 2 == -2)

    X_w, Y_w = -2, 2
    check("S5.7 [E3 VERDICT: X != Y] axion Hamiltonian weight X = %d vs a0 "
          "Hamiltonian weight Y = %d: the identification 'a0 = GS axion' "
          "FAILS at the weight level; the mismatch Y - X = 4 = -deg K_PT = "
          "|mu4| is exactly one canonical twist (numerical coincidence, "
          "typed [C]); a0 fills the GRAVITON slot instead -- 'the theory "
          "brings its own axion' is NOT established by a0" % (X_w, Y_w),
          X_w == -2 and Y_w == 2 and Y_w - X_w == 4 == MU4)

    A = sp.Matrix([[0, 0, -1], [1, 0, -1], [0, 1, -1]])
    xs = sp.symbols('x')
    eigs = set(sp.roots(A.charpoly(xs).as_expr(), xs).keys())
    sector_chars = {I ** j for j in (1, 2, 3)}
    check("S5.8 [TWISTED SLOTS MATCH] the three H2(ALE) classes carry the "
          "Coxeter eigencharacters %s = {i, -1, -i} = exactly the "
          "characters i^j of the three twisted ledger sectors j = 1, 2, 3 "
          "(bijection); det(A - 1) = -4 != 0: NO invariant class -- the "
          "three sphere 2-forms fill the three TWISTED axion slots, the "
          "bulk axion must come from the O(-2) tower field itself "
          "(construction [O])" % sorted(eigs, key=str),
          eigs == {-1, I, -I} and eigs == sector_chars
          and (A - sp.eye(3)).det() == -4)

    check("S5.9 [SLOT COUNT, soft] the equivariant ledger calls for 1 bulk "
          "axion (cancelling the Okubo square 36<x,x>^2, coupling "
          "lambda~_e8 = 6 per v495) + 3 twisted-sector partners (S3.9 "
          "obstruction); available: 1 tower slot O(-2) + 3 sphere classes = "
          "1 + 3 = 4 = |mu4| (count matches, N_fam = 3) -- a COUNT "
          "statement, not a construction",
          1 + 3 == 4 and N_fam == 3)


# ---------------------------------------------------------------------------
# S6 -- E4: negative controls
# ---------------------------------------------------------------------------
def section6(Q):
    print("  -- S6 (E4): negative controls")
    I = sp.I

    dens_f = [sp.expand((1 - I ** (-j)) ** 2) for j in (1, 2, 3)]
    recip_f = sp.simplify(sum(1 / d for d in dens_f))
    w_omega_f = sp.simplify(I * I)
    n_gens = 5
    check("S6.1 [NC-a: diag(i,i)] AB denominators (1 - i^{-j})^2 = %s = "
          "(2i, 4, -2i): COMPLEX (not in SU(2)); reciprocal sum = %s = 1/4 "
          "!= 5/4 (Dedekind identity broken); omega-weight i*i = %s != 1 "
          "(not Calabi-Yau: no O(2) symplectic slot, no ALE resolution); "
          "invariant ring needs %d generators (all degree-4 monomials; cone "
          "over the rational normal quartic, embedding dim 5): NOT a "
          "hypersurface -- the ledger breaks four independent ways"
          % (fmt(dens_f), recip_f, w_omega_f, n_gens),
          dens_f == [2 * I, 4, -2 * I] and recip_f == sp.Rational(1, 4)
          and w_omega_f == -1 and n_gens == 5)

    so16_dims = [45 + 15, 0, 60, 0]
    chi_so16 = [sp.expand(sum(I ** (j * m) * so16_dims[m] for m in range(4)))
                for j in range(4)]
    C = sp.Matrix([[2, -1, 0], [-1, 2, -1], [0, -1, 2]])
    Cinv = C.inv()
    defect_so16 = sum(so16_dims[m + 1] * (-Cinv[m, m] / 2) for m in range(3))
    check("S6.2 [NC-b: SO(16) GLUE, classes] so16 = so10 (+) so6 (+) (10,6) "
          "under the same h: Z4-decomposition 120 = 60 + 0 + 60 + 0 "
          "(classes {0, 2} ONLY -- a mere Z2 monodromy); characters %s = "
          "(120, 0, 120, 0): at j = 2 the value +120 != -8 (wrong sign AND "
          "magnitude); odd McKay nodes EMPTY (no Z4 quiver pairing, the "
          "64/64 slots are gone); ch2 defect %s = -30 != -78: a DIFFERENT "
          "ledger, quantified" % (fmt(chi_so16), defect_so16),
          so16_dims[0] + so16_dims[2] == 120
          and chi_so16 == [120, 0, 120, 0] and chi_so16[2] != -8
          and defect_so16 == -30)

    Qso = sp.expand(Q[0] + Q[2])
    Kso = sp.expand(28 * (S5 + S3))
    ok_okubo_fail = not in_span(Qso, [sp.expand(Kso ** 2)])
    dso = decompose(Qso, [b for _, b in BASIS4])
    check("S6.3 [NC-b: SO(16) BULK OKUBO FAILS] Tr_adj so16 X^4 = %s in "
          "(P1, P2, P3, T5, T3): T5 = 16, T3 = -32 NONZERO -- so16 has an "
          "independent quartic Casimir, its box anomaly is NOT a perfect "
          "square, no single-axion GS exists (so16 is not on Costello's "
          "list without matter; the v495 sl4-style control transplanted); "
          "D8-vector variant: h_v = 1/2 non-integer => no current closure "
          "(v502 S6.1)" % fmt(dso),
          ok_okubo_fail and dso == [12, 24, 36, 16, -32])

    cvals = {k: F(248 * k, k + H_VEE) for k in (1, 2)}
    n_k = {k: 12 * k * k for k in (1, 2)}
    check("S6.4 [NC-c: k = 2] the closure/current dial kills k = 2 three "
          "ways (S4.8 count 0, S4.9 residual 360, c(E8,2) = %s => prefactor "
          "-c/24 = %s != -1/3), while naive flux-integrality-type dials "
          "(e.g. 12k^2 = %s) stay INTEGER -- integrality alone does NOT "
          "kill k = 2; the kill is the closure dial, exactly as on the "
          "v502 CFT side (honest)" % (cvals[2], -cvals[2] / 24, n_k[2]),
          cvals[2] == F(31, 2) and -cvals[2] / 24 == F(-31, 48)
          and n_k == {1: 12, 2: 48})


# ---------------------------------------------------------------------------
# S7 -- E5: honest fences + kill-test status
# ---------------------------------------------------------------------------
def section7():
    print("  -- S7 (E5): honest fences and kill-test status")
    check("S7.1 [WHAT IS EXACT] (a) the Atiyah-Bott skeleton of the twistor "
          "inflow: denominators (2,4,2), characters (248,0,-8,0), twisted "
          "quadratics/quartics, invariant Okubo 36<x,x>^2, heterotic "
          "identities; (b) the index bridge f(m) = ch2(T_m) and the -78 "
          "glue defect (two routes); (c) the level dials (current count "
          "240/0/0/0, residual (0,360,814,1362), one-scale-one-level); "
          "(d) the weight arithmetic of the axion slot (X = -2 vs Y = +2)",
          True)
    check("S7.2 [HONEST SHARPENINGS] (i) per-sector Okubo FAILS (twisted "
          "quartics T5: (-8,32,-8), T3: (48,-64,48)); (ii) the AB sum "
          "cancels the D5 quartic but leaves 32 T3 exactly -- the fixed-"
          "point ledger alone is not GS-square, and no admissible "
          "reweighting fixes it; (iii) a0 is the GRAVITON modulus, not the "
          "axion (weight mismatch 4); (iv) h(J;k) = k integer for all k -- "
          "the k = 1 selection lives in the closure/current dial, "
          "geometrically exactly as on the CFT side", True)
    check("S7.3 [KILL-TEST STATUS] preregistered kill 'anomaly inflow "
          "demands a level != 1': the computable twistor-side pieces "
          "(fixed-point characters, index bridge, integral -78 defect, "
          "lattice current count, embedding residuals, one-scale argument) "
          "are ALL consistent with k = 1 and kill k = 2, 3, 4 via the "
          "closure dial: the kill branch does NOT fire on the equivariant "
          "skeleton; what the skeleton CANNOT decide: the actual BCOV "
          "one-loop coefficient on PT/Z4 -- the 32 T3 residual must be "
          "matched by twisted closed-string exchange, and if it is NOT the "
          "inflow would fail for reasons ORTHOGONAL to the level", True)
    check("S7.4 [O-FENCE] open: the global BCOV/Kodaira-Spencer theory on "
          "PT/Z4 (quantisation, non-fixed-point contributions, the twisted "
          "closed-string measure), the O(-2) bulk-axion construction, the "
          "CPS level-from-flux mechanism, and the pairing of the 32 T3 "
          "residual with the three sphere axions; SEAM.EQUIV.01 untouched; "
          "NO marker moves anywhere", True)


# ---------------------------------------------------------------------------
def run():
    reset()
    print("v505  CELEST.WP5E.BETA.01: the equivariant anomaly ledger on "
          "twistor space (WP5e-beta of CELEST.SEAM.01 -- the twistor-side "
          "milestone)")
    roots = build_glue_roots()
    dims, chi_dims, dens = section1(roots)
    K, Q = power_sums_by_class(roots)
    Ktw = section2(K)
    section3(Q, Ktw)
    section4(roots, dims, chi_dims, dens)
    section5()
    section6(Q)
    section7()

    return summary("v505 CELEST.WP5E.BETA.01: the Atiyah-Bott/Lefschetz "
                   "skeleton of the WP5e twistor inflow is exact -- "
                   "denominators (2,4,2) with Dedekind sum 5/4 = "
                   "(|Z4|^2-1)/12, equivariant characters (248,0,-8,0) by "
                   "two routes, invariant average 60 = carrier, Frobenius "
                   "61568; the invariant sector obeys Okubo exactly "
                   "(36<x,x>^2, 36 = lambda~^2_e8), the twisted sectors do "
                   "NOT, and the AB sum cancels the D5 quartic leaving the "
                   "rigid residual 32 T3; the index bridge f(m) = ch2(T_m) "
                   "= -(C^{-1})_mm/2 holds exactly with glue defect -78 by "
                   "both routes; the level dials say k = 1 geometrically "
                   "(current count 240/0/0/0, residual 0/360/814/1362, one "
                   "scale => one level); E3 REFUTED: a0 fills the GRAVITON "
                   "slot O(2), not the axion slot O(-2) (X = -2 != +2 = Y), "
                   "while the three H2(ALE) classes fill the three twisted "
                   "axion slots (Coxeter characters = sector characters); "
                   "controls diag(i,i)/SO(16)/k = 2 break the ledger in "
                   "quantified ways.  Verdict B; the kill branch does not "
                   "fire on the equivariant skeleton; the global BCOV "
                   "quantisation stays [O]; SEAM.EQUIV.01 untouched")


if __name__ == "__main__":
    raise SystemExit(1 if run() else 0)

b0ad20db2bdeda29c98604ca76d3ebe9d74acb60399d66b643b1f2c9f9295703

The public registry includes 116 document views, 6,693 section records and 985 physical LaTeX sources. Active publication documents, research material and archives retain their recorded roles. This is source navigation, not a new semantic review or proof run.