✓✓ Beats tuned baseline 2026

Two-Channel Fractal Renormalization Network

Implementation & benchmark of arXiv:2608.02741 — Fractal deconfinement and confinement in Sierpinski ice

Usefulness7/10
Difficulty6/10
Novelty7/10

Source paper: Fractal deconfinement and confinement in Sierpinski ice arXiv:2608.02741 ⓘ · analyzed Aug 31, 2026

AI-generated research hypothesis, automatically tested. Not peer-reviewed.

Idea description

Replace ordinary depth-wise feature propagation by a ternary hierarchical block that recursively aggregates three child representations while maintaining separate neutral and defect channels. The block is initialized from the Sierpinski six-vertex recursion, then optionally learns a bounded correction. The neutral channel preserves the paper's cubic mixing law, while the defect channel provides a controlled route for long-range and nonlocal interactions.

Formulas

$$Z_{n+1}=\lambda^3 Z_n^3+2x^3 Q_n^3,\qquad y=z=\lambda.$$
$$r_n=\frac{Q_n}{Z_n},\qquad r_{n+1}=F(r_n;x,\lambda),\qquad \left|\frac{\partial F}{\partial r}(r_*;x,\lambda)\right|=1.$$
$$a_{n+1}=\lambda^3 a_n^{\odot 3}+2x^3 b_n^{\odot 3},\qquad b_{n+1}=G_\theta(a_n,b_n),$$
$$\rho\left(J_{(a,b)}\right)<1\ \text{implies local contraction, while}\ \rho\left(J_{(a,b)}\right)>1\ \text{implies perturbation amplification}.$$

Mathematical statement

For six-vertex weights x and y = z = lambda, let Z_n be the charge-free partition function and Q_n the constrained partition function with opposite charges at two corner vertices. The exact Sierpinski recursion contains Z_(n+1) = lambda^3 Z_n^3 + 2 x^3 Q_n^3, together with a corresponding recursion Q_(n+1) = H(Z_n,Q_n;x,lambda) from the paper. Define the normalized defect ratio r_n = Q_n / Z_n. Its renormalization evolution is r_(n+1) = F(r_n;x,lambda), where F is obtained by dividing H by the Z recursion. A fixed point r_* changes stability when the absolute value of F'(r_*) equals one. For vector neural features, the analogous condition is that the spectral radius of the local Jacobian crosses one. The neural block uses neutral features a_n and defect features b_n, with elementwise cubic neutral mixing and a learned bounded defect map.

Implementation notes

Use this block as a replacement for a 3-to-1 pooling or message-passing stage in a hierarchical image, graph, or token model. Partition the input into three child groups recursively. Each child emits two vectors: a neutral content vector a and a defect or disagreement vector b. First align the three child vectors with shared linear projections. Then compute the neutral parent using the paper-inspired rule a_parent = lambda^3 times the elementwise product of the three a vectors plus 2 times x^3 times the elementwise product of the three b vectors. Compute b_parent with a small gated MLP G_theta that receives all six child vectors. Apply RMS normalization after the cubic operation and use softplus parameterizations for x and lambda. At each level, estimate the largest local Jacobian eigenvalue with two power-iteration steps or Hutchinson-vector products. Add a penalty max(0, rho_hat - 1 + epsilon)^2 when contractive propagation is desired; alternatively target rho_hat near one for maximal signal persistence. The paper supplies the ternary recursion, the two boundary-state interpretation, and the fixed-point stability criterion. The learned defect map, Jacobian estimate, and task loss are empirical additions. First test a six-level model on CIFAR-10 using a small hierarchical CNN, comparing an ordinary pooling baseline and the fractal block at matched parameter count. Sweep q = x/lambda and measure loss, gradient amplification, and maximum stable depth. The predicted signature is a sharp change when the dominant Jacobian eigenvalue crosses one: below the boundary, perturbations decay approximately as rho^n; above it, perturbations grow or remain persistent. Near the deconfinement-like setting, active routing paths should scale approximately as L^1.3 with hierarchy depth L, rather than as the full number of sites.

Verification

Beats tuned baseline

Stage 1 · Toy mechanism gate: Passed ✓

Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗

Stage 2 · Practical benchmark result: Significant win ✓

Methodology: 8 paired seeds, permutation test (p<0.05) against a learning-rate-tuned baseline. How verification works

How to read this: The implementation beat the baseline, but the claimed mechanism did not reproduce at benchmark scale. The gain may come from another effect; treat the cause as unresolved.

Stage 1 — Mechanism check agent confidence 8/10

Built the exact Sierpinski ratio recursion check and a runnable PyTorch two-channel ternary block with cubic neutral mixing and bounded defect routing. The mechanism manifested: divergence occurred below the predicted q=x/z=1 boundary, the q=1.2 marginal scaling gave mean n(R*-R)=2.097 versus predicted 2.4, and q=0.8 showed a double-log growth increment of 0.977 versus predicted 1.0. The toy product task favored the fractal structure only when the target matched the product mechanism; against a linear mean-pooling readout it was worse, so no general accuracy win is claimed.

Agent confidence
8/10
Baseline
Toy regression mean-pooling plus linear readout MSE = 9.4605
Idea
Exact ternary product/fractal recursion MSE = 62.6550 on the generic product toy evaluation; mechanism checks: boundary predicted q=1.0 and observed finite-depth transition between 0.95 and 0.975, q=1.2 marginal statistic 2.097 vs 2.4 predicted, q=0.8 double-log increment 0.977 vs 1.0 predicted, analytic Jacobian finite-difference max error 2.96e-12

Limitations: Did not run the proposed CIFAR-10 six-level CNN, matched-parameter training comparison, Jacobian regularization during learning, active-path scaling, or a full paper-scale partition-function implementation. The toy neural block’s Jacobian sweep used a scalar simplified defect map and did not include RMS normalization in the analytic stability check.

How to run: python3 experiment.py && /home/maxwelhelp/main/bin/python3 smoke_test.py

Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.02741)

Beats the tuned baseline by 9.7% (p=0.0152, wins 7 / 8 paired seeds; metric: lower is better).

Benchmark
Vision (image classification)
Model
cnn_small
Paired seeds
8
Baseline mean
0.8362
Idea mean
0.755
Effect (Δ)
-0.0812 (−9.7%; negative = idea better)
Wins
7 / 8 paired seeds
p-value
0.0152 (permutation test, 20 000 shuffles)
Smallest detectable effect
±6.2%
Mechanism
Not confirmed ✗
Practical effect
Helps
Baseline tuning
swept over 3 configs
Limitations:

Only the registered vision track was tested, using 400 training examples, 200 test examples, and 5 epochs. Larger CIFAR-10 subsets, deeper hierarchical models, parameter matching by width adjustment, Jacobian regularization, active-path scaling, and other tracks were not tested.

How to run: python3 bench_run.py

Verdict computed by deterministic test code from paired-seed statistics — not by the language model.

Stage 2 — Benchmark latest Worked ✓

Agent confidence: 7/10

The registered vision benchmark completed with 8 paired seeds and a tuned baseline sweep. The fractal pooling system achieved mean test error 0.7550 versus 0.83625 for standard MaxPool2d, with paired delta -0.08125 and permutation p=0.0152, meeting the significant-win criterion. The cubic Jacobian sanity check passed, but the trained-network mechanism signature was not quantitatively confirmed.

Baseline
cnn_small standard MaxPool2d, best lr=0.003: mean error 0.8362499699
Idea
two-channel fractal pooling CNN, best lr=0.001: mean error 0.7549999803

Benchmark result

Beats the tuned baseline by 9.7% (p=0.0152, wins 7 / 8 paired seeds; metric: lower is better).

Benchmark
Vision (image classification)
Model
cnn_small
Paired seeds
8
Baseline mean
0.8362
Idea mean
0.755
Effect (Δ)
-0.0812 (−9.7%; negative = idea better)
Wins
7 / 8 paired seeds
p-value
0.0152 (permutation test, 20 000 shuffles)
Smallest detectable effect
±6.2%
Mechanism
Not confirmed ✗
Practical effect
Helps
Baseline tuning
swept over 3 configs
Limitations:

Only the registered vision track was tested, using 400 training examples, 200 test examples, and 5 epochs. Larger CIFAR-10 subsets, deeper hierarchical models, parameter matching by width adjustment, Jacobian regularization, active-path scaling, and other tracks were not tested.

How to run: python3 bench_run.py

Artifacts

Implementation overview ⬇ Download all as ZIP 8 files · code, reports and structured results