△ Mechanism confirmed, baseline not beaten 2026

Cross-Channel Scattering Front End

Implementation & benchmark of arXiv:2607.08855 — Spatial Neighboring Scattering Transform: A Cross-Channel Amplitude Coupling Measure for EEG Connectivity

Usefulness6/10
Difficulty4/10
Novelty7/10

Source paper: Spatial Neighboring Scattering Transform: A Cross-Channel Amplitude Coupling Measure for EEG Connectivity arXiv:2607.08855 ⓘ · analyzed Aug 30, 2026

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

Idea description

Add a differentiable SNST layer before an EEG classifier or sequence model. For every local channel pair and wavelet band, compute the magnitude of the complex cross-channel analytic response, then average it over a controllable temporal window and concatenate it with ordinary channelwise features. This gives the model an explicit, phase-robust amplitude-coupling representation that is especially useful when labeled training data are scarce.

Formulas

$$S_{1}^{\mathrm{WST}}(t,\lambda)=\big(|x*\psi_{\lambda}|*\phi_{J}\big)(t)$$
$$U_{1}(t;m,n,\lambda)=\left|(x_m*\psi_{\lambda})(t)\cdot(x_n*\psi_{\lambda})^{*}(t)\right|$$
$$S_{1}(t;m,n,\lambda)=\big(U_1(\cdot;m,n,\lambda)*\phi_J\big)(t)$$
$$\widetilde{U}_{1}(t;m,n,\lambda)=\sqrt{\operatorname{Re}(z_m^\lambda(t)z_n^\lambda(t)^*)^2+\operatorname{Im}(z_m^\lambda(t)z_n^\lambda(t)^*)^2+\epsilon},\quad z_m^\lambda=x_m*\psi_\lambda$$

Mathematical statement

The paper's first-order single-channel scattering coefficient is $S_{1}^{\mathrm{WST}}(t,\lambda)=\big(|x*\psi_{\lambda}|*\phi_J\big)(t)$, where $x$ is a scalar signal, $\psi_{\lambda}$ is a band-pass wavelet centered at scale or frequency $\lambda$, $\phi_J$ is a low-pass averaging window with invariance scale $T=2^J$, and $*$ denotes convolution. The multichannel SNST replaces the pointwise modulus with $U_1(t;m,n,\lambda)=\left|(x_m*\psi_{\lambda})(t)\cdot(x_n*\psi_{\lambda})^*(t)\right|$, where $x_m,x_n$ are two distinct channels and $^*$ is complex conjugation. The pooled feature is $S_1(t;m,n,\lambda)=\big(U_1(\cdot;m,n,\lambda)*\phi_J\big)(t)$, using the same temporal averaging idea as scattering. Since $|ab^*|=|a||b|$, this feature measures joint band-limited amplitude envelopes while discarding the phase of the cross-product. The neural adaptation uses learnable or fixed complex wavelets $\psi_\lambda$, local channel pairs $(m,n)$ from a sensor graph, and a differentiable magnitude $\sqrt{\operatorname{Re}(q)^2+\operatorname{Im}(q)^2+\epsilon}$.

Implementation notes

(1) Integration point: place this layer immediately after input normalization and before a temporal CNN, graph neural network, or Transformer. The input is a batch $X\in\mathbb{R}^{B\times C\times T}$, with $C$ sensor channels. Choose 4-8 wavelet bands and a sensor adjacency list containing spatially neighboring pairs; do not form all $C^2$ pairs initially. Use complex Morlet filters implemented as two real 1-D convolutions, or use fixed filters for the first experiment. For each band $\lambda$, compute complex responses $Z[:,m,\lambda,:]=(X_m*\psi_\lambda)$ and pairwise products only for $(m,n)$ in the adjacency list. Concatenate the pooled SNST features with the original input or with channelwise scattering features. (2) Pseudocode: Z = complex_conv1d(X, psi); for lambda in bands: for (m,n) in edges: q = Z[:,m,lambda,:] * conj(Z[:,n,lambda,:]); u = sqrt(real(q)2 + imag(q)2 + eps); s = depthwise_avg_pool(u, kernel_size=2**J); features.append(s); H = concat(features, baseline_features); loss = cross_entropy(classifier(H), y). The displayed $U_1$ formula is the computation of u, and the displayed $S_1$ formula is the averaging operation. (3) Mathematical quantities are the complex wavelet responses, conjugate products, magnitude, and low-pass pooling. Estimate no statistics from the paper; select $J$ so that $2^J$ is between 50 and 250 ms of data, and optionally learn a positive scalar multiplier per band after initialization. Normalize each pair feature using training-set mean and variance, with the normalization statistics frozen for validation and test. To limit redundancy, share a linear projection across all edge features or retain only the top-$K$ spatial edges by validation performance. (4) First experiment: use BCI Competition IV-2a, 4-class motor-imagery EEG, with a small 2-layer temporal CNN or EEGNet baseline. Compare raw-input EEGNet, channelwise scattering plus classifier, and channelwise scattering plus SNST edge features under identical parameter budgets and five subject-wise seeds. The primary success signal is higher held-out accuracy or macro-F1 at equal labels and FLOPs; a useful secondary signal is a larger accuracy gain in the first 10-20% of the training data and lower seed-to-seed variance. Also ablate conjugation, edge locality, and temporal pooling to verify that the gain comes from cross-channel amplitude coupling rather than simply adding filter-bank capacity.

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 7/10

Built a differentiable fixed-Morlet Cross-Channel Scattering front end with local edge pairing, complex conjugate products, differentiable magnitude, and temporal averaging. The numerical checks confirmed |z_m z_n*|=|z_m||z_n|, phase invariance, and exact agreement with an independent implementation to floating-point precision. On the scarce-label synthetic amplitude-coupling task, SNST achieved 0.5406 +/- 0.0181 accuracy versus 0.5104 +/- 0.0053 for the raw CNN across three seeds, a promising but modest signal; this is not evidence of a real-EEG win.

Agent confidence
7/10
Baseline
Raw temporal CNN: held-out accuracy [0.5156, 0.5125, 0.5031], mean 0.5104, std 0.0053
Idea
SNST front end plus classifier: held-out accuracy [0.5156, 0.5484, 0.5578], mean 0.5406, std 0.0181; +0.0302 absolute mean accuracy

Limitations: Tested only on a small synthetic four-channel amplitude-coupling task, not BCI Competition IV-2a or real EEG. Only three seeds were used; no channelwise-scattering baseline, ablations, learned wavelets, strict FLOP/parameter matching, or subject-wise evaluation was performed. CUDA convolution failed on the available runtime and the script correctly fell back to CPU.

How to run: python3 experiment.py

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

No significant difference vs. the tuned baseline (+0.9%, p=0.6465, wins 3 / 8 paired seeds).

Benchmark
Multichannel_amplitude_coupling
Model
matched_temporal_cnn
Paired seeds
8
Baseline mean
0.4831
Idea mean
0.4875
Effect (Δ)
0.0044 (+0.9%; negative = idea better)
Wins
3 / 8 paired seeds
p-value
0.6465 (permutation test, 20 000 shuffles)
Smallest detectable effect
±5.0%
Mechanism
Confirmed ✓
Practical effect
No effect
Baseline tuning
swept over 6 configs
Limitations:

This is a deterministic synthetic multichannel sequence track rather than real EEG or BCI Competition IV-2a. The SNST front end is fixed rather than learned, only three local edges and two bands were tested, and parameter/FLOP equality was not formally counted. The mechanism signature verifies trained-model implementation consistency, not that coupling improves prediction; no channelwise-scattering ablation or real-subject evaluation was run.

How to run: python3 bench_experiment.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

Built a local multichannel amplitude-coupling sequence track and matched end-to-end temporal CNN systems, with fixed Morlet cross-channel scattering as the only intervention. The math sanity check passed to approximately 1.2e-7 error and the implementation signature matched exactly, but SNST test error was 0.48750 versus 0.48312 for the tuned raw baseline (paired delta +0.00438, p=0.64655), so there is no observed transfer win.

Baseline
8-seed test error mean 0.483125, std 0.031319; per-seed [0.4850, 0.4850, 0.4950, 0.4850, 0.5000, 0.5400, 0.4400, 0.4350]; best sweep lr=0.006, weight_decay=0.001, epochs=18
Idea
8-seed test error mean 0.487500, std 0.018708; per-seed [0.4750, 0.5000, 0.4800, 0.5000, 0.5000, 0.5000, 0.5000, 0.4450]; best idea setting lr=0.006, weight_decay=0.001, epochs=18
Paired comparison
idea minus baseline = +0.004375 mean; 3/8 idea wins; permutation p=0.64655; verdict no measurable effect

Benchmark result

No significant difference vs. the tuned baseline (+0.9%, p=0.6465, wins 3 / 8 paired seeds).

Benchmark
Multichannel_amplitude_coupling
Model
matched_temporal_cnn
Paired seeds
8
Baseline mean
0.4831
Idea mean
0.4875
Effect (Δ)
0.0044 (+0.9%; negative = idea better)
Wins
3 / 8 paired seeds
p-value
0.6465 (permutation test, 20 000 shuffles)
Smallest detectable effect
±5.0%
Mechanism
Confirmed ✓
Practical effect
No effect
Baseline tuning
swept over 6 configs
Limitations:

This is a deterministic synthetic multichannel sequence track rather than real EEG or BCI Competition IV-2a. The SNST front end is fixed rather than learned, only three local edges and two bands were tested, and parameter/FLOP equality was not formally counted. The mechanism signature verifies trained-model implementation consistency, not that coupling improves prediction; no channelwise-scattering ablation or real-subject evaluation was run.

How to run: python3 bench_experiment.py

Artifacts

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