Cross-Channel Scattering Front End
Implementation & benchmark of arXiv:2607.08855 — Spatial Neighboring Scattering Transform: A Cross-Channel Amplitude Coupling Measure for EEG Connectivity
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
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
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
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
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
- bench_experiment.py 6.1 KB View Raw
- bench_report.json 4.2 KB View Raw
- experiment.py 6.1 KB View Raw
- report.md 1.4 KB View
- report_bench_2026-09-04T201907.md 5.0 KB View