△ Mechanism confirmed, baseline not beaten 2026

Quotient Spectral Positional Encoding

Implementation & benchmark of arXiv:2607.08987 — Group Invariant Spectral Embedding

Usefulness8/10
Difficulty6/10
Novelty6/10

Source paper: Group Invariant Spectral Embedding arXiv:2607.08987 ⓘ · analyzed Aug 30, 2026

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

Idea description

Construct a graph and its spectral positional features using affinities between inputs after optimally aligning one input over the known symmetry group. Feed these quotient-space eigenvectors to a transformer or graph neural network, so symmetry-equivalent examples receive the same structural coordinates without storing augmented copies.

Formulas

$$[x]=\{g\cdot x:g\in G\},\qquad \mathcal{N}=\mathcal{M}/G$$
$$W_{ij}^{\mathrm{quot}}=\exp\left(-\frac{\min_{g\in G}d(x_i,g\cdot x_j)^2}{\varepsilon}\right),\qquad D_{ii}=\sum_{j=1}^{n}W_{ij}^{\mathrm{quot}}$$
$$L_{\mathrm{sym}}=I-D^{-1/2}W^{\mathrm{quot}}D^{-1/2}$$
$$\text{standard error}=O(\varepsilon)+O_P\left(n^{-1/2}\varepsilon^{-1/2-\dim(\mathcal{M})/4}\right)$$

Mathematical statement

The paper considers a Riemannian manifold \(\mathcal{M}\) with a compact Lie group \(G\) acting smoothly, freely, and isometrically, and forms the quotient \(\mathcal{N}=\mathcal{M}/G\). For observations \(x_i\), define the orbit \([x_i]=\{g\cdot x_i:g\in G\}\). A minimization-based invariant affinity can be implemented as \(W_{ij}=\exp[-\min_{g\in G}d(x_i,g\cdot x_j)^2/\varepsilon]\), where \(d\) is the input metric and \(\varepsilon>0\) is the bandwidth. With degree \(D_{ii}=\sum_jW_{ij}\), the normalized graph Laplacian is \(L=I-D^{-1/2}WD^{-1/2}\); its low-eigenvalue eigenvectors are quotient-space coordinates. The paper states that the ordinary graph-Laplacian convergence error has the form \(O(\varepsilon)+O_P(n^{-1/2}\varepsilon^{-1/2-\dim(\mathcal{M})/4})\), while invariant kernels converge to differential operators on \(\mathcal{M}/G\), replacing the ambient manifold dimension by the lower effective quotient dimension under the theorem's regularity assumptions. The transferable property is orbit identification before spectral estimation, which removes variation along group directions rather than asking a downstream network to learn it from augmentation.

Implementation notes

1. Integration point: use this as a frozen or periodically refreshed positional-encoding module before a transformer encoder, or as structural node features in a graph neural network. Inputs are a batch or memory bank of vectors \(x_1,\ldots,x_n\), together with a known differentiable group action \(g\cdot x\), such as 2D rotations of images, 3D rotations of point clouds, or permutations represented by a finite group. 2. Pseudocode: choose bandwidth \(\varepsilon\); for each pair \((i,j)\), sample group elements \(g_1,\ldots,g_K\), compute \(r_{ij}=\min_k\|x_i-g_k\cdot x_j\|_2^2\), set \(W_{ij}=\exp(-r_{ij}/\varepsilon)\), set \(D_{ii}=\sum_jW_{ij}\), form \(L=I-D^{-1/2}WD^{-1/2}\), and compute the smallest \(m\) nontrivial eigenvectors using Lanczos. Concatenate these \(m\)-dimensional coordinates to token or node embeddings. 3. The paper supplies the quotient-manifold construction and convergence motivation; estimate the group minimization numerically, choose \(K\), \(\varepsilon\), and \(m\) empirically, and refresh the graph every few epochs if embeddings change. For continuous rotations, use a coarse angular grid followed by local refinement; for finite groups enumerate all elements. 4. First experiment: train a small ViT on rotated MNIST or CIFAR-10 with random SO(2) rotations. Compare ordinary Laplacian positional encodings, explicit 8-fold augmentation, and quotient spectral encodings at equal stored examples and FLOPs. Success is higher accuracy on unseen rotation angles, lower train-test gap under rotation, or matching augmentation accuracy with substantially less data and memory.

Verification

Mechanism confirmed, baseline not beaten

Stage 1 · Toy mechanism gate: Passed ✓

Stage 2 · Mechanism transferred to benchmark: Confirmed ✓

Stage 2 · Practical benchmark result: No measurable effect

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

Stage 1 — Mechanism check agent confidence 8/10

Built a quotient spectral positional-encoding prototype using a torus toy manifold with a nuisance S1 rotation group, normalized graph Laplacian eigenvectors, exact or grid-sampled orbit alignment, and a downstream classifier. The mechanism manifested: quotient affinity was exactly invariant (maximum change 0), angular grid error scaled as K^-1.997 versus the predicted K^-2, and operator RMSE showed a finite-bandwidth optimum near epsilon=3.0 with quotient RMSE 0.0728 versus ordinary 0.0794. The sample-size sweep gave slope -0.149 rather than the ideal -0.5, so asymptotic rate improvement was not established; downstream accuracy nevertheless improved from 0.945 ordinary to 0.967 exact quotient and 0.978 with K=8.

Agent confidence
8/10
Baseline
Ordinary spectral features: operator RMSE at eps=[0.04,0.12,0.35,0.85,1.5,3.0,6.0] = [6.994,1.879,0.623,0.276,0.189,0.0794,0.312]; classifier accuracy 0.9451.
Idea
Quotient spectral features: operator RMSE = [1.386,0.614,0.290,0.194,0.167,0.0728,0.312]; exact orbit invariance max affinity-distance change 0; K-grid median distance errors [0.1518,0.0383,0.00960,0.00240,0.000599,0.000151] with slope -1.997; classifier accuracy 0.9670 exact and 0.9780 at K=8.

Limitations: This is a small synthetic torus experiment, not rotated MNIST/CIFAR or a transformer/GNN. The quotient distance uses an analytic S1 alignment for the exact case and does not benchmark runtime, memory, refresh schedules, learned embeddings, or real image noise. The sample-size scaling sweep did not verify the theorem's ideal n^-1/2 stochastic exponent, and the bandwidth/operator test is a finite-sample proxy rather than a proof of manifold convergence.

How to run: python3 quotient_spectral_experiment.py

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

Ahead of the baseline (-0.0212 metric points) but not statistically significant (p=0.7959, wins 5 / 8 paired seeds).

Benchmark
Rotating_shapes
Model
custom_cnn
Paired seeds
8
Baseline mean
0
Idea mean
0
Effect (Δ)
-0.0212 (n/a; negative = idea better)
Wins
5 / 8 paired seeds
p-value
0.7959 (permutation test, 20 000 shuffles)
Smallest detectable effect
n/a
Mechanism
Confirmed ✓
Practical effect
No effect
Baseline tuning
swept over 0 configs
Limitations:

The built-in CIFAR track was not used because it has no known rotation group, so a local custom rotating-shapes track was used. The implementation used a custom CNN rather than bench.make_model and manually ran the learning-rate grid rather than bench.sweep_baseline. Exact continuous SO(2) alignment, larger real-image data, runtime/memory scaling, graph refreshes, and transformer/GNN integration were not tested.

How to run: python3 quotient_bench.py

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

Stage 2 — Benchmark latest Failed ✗

Agent confidence: 9/10

Implemented a local rotation-symmetric image track and matched CNN systems differing in ordinary versus quotient spectral positional features. The idea mean error was 0.16500 versus baseline 0.18625, but the paired delta was -0.02125 with permutation p=0.7959, so there was no significant win; the mechanism check showed small rotation-induced feature change but is not fully compliant as an independently measured trained-model prediction signature.

Baseline
8-seed test errors [0.120, 0.045, 0.420, 0.095, 0.285, 0.160, 0.095, 0.270], mean 0.18625; best lr=0.003
Idea
8-seed test errors [0.060, 0.145, 0.140, 0.185, 0.050, 0.115, 0.455, 0.170], mean 0.16500; best lr=0.003
Paired delta
-0.02125 idea-minus-baseline; p=0.7959; 5/8 wins

Benchmark result

Ahead of the baseline (-0.0212 metric points) but not statistically significant (p=0.7959, wins 5 / 8 paired seeds).

Benchmark
Rotating_shapes
Model
custom_cnn
Paired seeds
8
Baseline mean
0
Idea mean
0
Effect (Δ)
-0.0212 (n/a; negative = idea better)
Wins
5 / 8 paired seeds
p-value
0.7959 (permutation test, 20 000 shuffles)
Smallest detectable effect
n/a
Mechanism
Confirmed ✓
Practical effect
No effect
Baseline tuning
swept over 0 configs
Limitations:

The built-in CIFAR track was not used because it has no known rotation group, so a local custom rotating-shapes track was used. The implementation used a custom CNN rather than bench.make_model and manually ran the learning-rate grid rather than bench.sweep_baseline. Exact continuous SO(2) alignment, larger real-image data, runtime/memory scaling, graph refreshes, and transformer/GNN integration were not tested.

How to run: python3 quotient_bench.py

Artifacts

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