Delay-Aware Frequency-Preserving Recurrent Coupling
Implementation & benchmark of arXiv:2607.18913 — How network perturbations distort agreement trajectories in LTI multi-agent systems
Source paper: How network perturbations distort agreement trajectories in LTI multi-agent systems arXiv:2607.18913 ⓘ · analyzed Aug 30, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
For coupled recurrent or state-space modules that represent oscillatory or periodic signals, explicitly account for communication or attention delay in the characteristic equation. Tune the coupling gain or add a phase-lead compensator so that the desired latent frequency remains a closed-loop mode instead of being shifted by small delays.
Formulas
Mathematical statement
Let \(x_i(t)\in\mathbb{R}^d\) be the state of module \(i\), with local dynamics \(A\), coupling matrices \(B,K\), graph Laplacian \(L\), and transmission delay \(\tau\). The delayed coupled system is \(\dot x_i(t)=Ax_i(t)-BK\sum_jL_{ij}x_j(t-\tau)\). After diagonalizing \(L\), each graph mode with eigenvalue \(\lambda_j\) has characteristic equation \(\det[sI-A+\lambda_jBK e^{-s\tau}]=0\). The consensus mode \(\lambda_1=0\) is delay-independent, but nonzero modes acquire a frequency-dependent phase factor. At a target oscillation \(s=\mathrm{i}\omega_0\), delay changes the coupling phase by \(-\omega_0\tau\), so periodic synchronization is generally lost. A phase-lead compensator \(C(s)\) should satisfy \(C(\mathrm{i}\omega_0)\approx e^{\mathrm{i}\omega_0\tau}\).
Implementation notes
Integrate this into a recurrent or state-space architecture whose hidden state is updated in several coupled branches, for example four diagonal SSM blocks or four small GRU-like modules connected by a graph. Use a continuous-time-inspired Euler update x_i[t+1] = x_i[t] + h*(A*x_i[t] - B*K*sum_j L[i,j]*x_j[t-d]), where \(h\) is the integration step and \(d\) is an integer delay. First choose \(A\) with a lightly damped complex pair at target frequency \(\omega_0\), and train only a readout on synthetic mixtures of sinusoids so the effect of coupling is isolated. Implement two variants: ordinary delayed coupling and phase-compensated coupling using a two-tap predictor \(\widehat{x}_j[t+d]=x_j[t]+d(x_j[t]-x_j[t-1])\), or a learnable FIR filter constrained to have frequency response near \(e^{\mathrm{i}\omega_0\tau}\). The paper contributes the modal delay characteristic equation and the distinction between delay-invariant zero-frequency consensus and delay-sensitive periodic modes. Estimate \(A,B,K\) from module weights and estimate dominant roots by scanning complex \(s=\alpha+\mathrm{i}\omega\) for minima of the characteristic determinant; measure the output spectrum with an FFT. Sweep \(\tau\), coupling gain \(\gamma\), and target frequency. The falsifiable prediction is a loss of phase locking or a dominant-frequency shift once \(\omega_0\tau\) becomes appreciable, with instability beginning when a characteristic root reaches \(\operatorname{Re}s=0\). The compensated model should preserve peak frequency within 5 percent up to delays for which the uncompensated model exhibits a measurable shift, while constant inputs should remain nearly delay-invariant.
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 9/10
Built and ran a delayed modal recurrent-coupling MVP with ordinary delayed feedback and a two-tap phase-lead predictor. The three mechanism predictions were confirmed: delay phase slope was -1.000 versus -1.000 predicted, the lambda=0 mode was exactly delay-invariant while a nonzero mode changed by relative RMS 0.825, and the small-gain frequency-shift slope differed from perturbation theory by 3.9%. In the four-branch ring toy experiment, compensation reduced mean nonzero-mode frequency error from 0.03425 rad/time to 0.000563 rad/time, a 98.4% reduction; this is a real signal for the stated mechanism, though not evidence of a general ML-training win.
- Agent confidence
- 9/10
- Baseline
- Delayed coupling: mean nonzero-mode frequency error 0.0342477 rad/time; lambda=2 error 0.0251704; lambda=4 error 0.0524022.
- Idea
- Two-tap predictor compensation: mean nonzero-mode frequency error 0.0005633 rad/time; lambda=2 error 0.0007321; lambda=4 error 0.0002257; relative reduction 98.36%. Phase slope -1.000 predicted/observed; frequency-shift slope 0.31074 observed vs 0.29909 predicted, a 3.90% relative error.
Limitations: Only a deterministic complex single-mode Euler DDE and a four-branch ring modal equivalent were tested. No learned GRU/SSM, synthetic readout training, multi-frequency mixture, broad delay/gain stability-boundary sweep, wall-clock comparison, or GPU experiment was included; FFT output-spectrum evaluation was replaced by a continuous phase-slope estimator.
How to run: python3 delay_coupling_experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2607.18913)
No significant difference vs. the tuned baseline (−-0.0%, p=1, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0021
- Idea mean
- 0.0021
- Effect (Δ)
- 0 (−-0.0%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 1 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.0%
- Mechanism
- Confirmed ✓
- Practical effect
- No effect
- Baseline tuning
- swept over 9 configs
Only the registered built-in dynamics track was tested. The experiment used a one-step delay, short 12-epoch training, 400 training and 200 test examples, and did not test longer horizons, multiple delays, broader oscillatory frequency tasks, larger recurrent architectures, or wall-clock/FLOP effects. The behavioral signature confirms delayed cross-branch alignment in trained hidden states but not preservation of a known target frequency on the benchmark task.
How to run: python3 delay_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
The delay-aware compensated coupled GRU was evaluated on the registered dynamics/control track with the required paired-seed protocol. The mechanism signature was confirmed behaviorally, but the best idea setting selected zero coupling and exactly matched the tuned baseline; nonzero compensated coupling was worse. Therefore the bench verdict is no measurable effect and worked=false.
- Baseline
- best_cfg={lr:0.006,coupling:0.0}; full mean test MSE=0.002108921129547525, std=0.001215784166403164, per_seed=[0.0009661344229243696,0.0021135066635906696,0.0017933575436472893,0.0040876418352127075,0.0011596226831898093,0.0011802533408626914,0.004168401472270489,0.0014024510746821761]
- Idea
- best_cfg={lr:0.006,coupling:0.0}; mean test MSE=0.002108921129547525, std=0.001215784166403164, per_seed=[0.0009661344229243696,0.0021135066635906696,0.0017933575436472893,0.0040876418352127075,0.0011596226831898093,0.0011802533408626914,0.004168401472270489,0.0014024510746821761]; compensated lr=0.001,coupling=0.08 mean=0.046488794789183885; compensated lr=0.006,coupling=0.16 mean=0.0030070416687522084
Benchmark result
No significant difference vs. the tuned baseline (−-0.0%, p=1, wins 0 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- rnn_small
- Paired seeds
- 8
- Baseline mean
- 0.0021
- Idea mean
- 0.0021
- Effect (Δ)
- 0 (−-0.0%; negative = idea better)
- Wins
- 0 / 8 paired seeds
- p-value
- 1 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±0.0%
- Mechanism
- Confirmed ✓
- Practical effect
- No effect
- Baseline tuning
- swept over 9 configs
Only the registered built-in dynamics track was tested. The experiment used a one-step delay, short 12-epoch training, 400 training and 200 test examples, and did not test longer horizons, multiple delays, broader oscillatory frequency tasks, larger recurrent architectures, or wall-clock/FLOP effects. The behavioral signature confirms delayed cross-branch alignment in trained hidden states but not preservation of a known target frequency on the benchmark task.
How to run: python3 delay_bench.py
Artifacts
- bench_report.json 5.1 KB View Raw
- delay_bench.py 5.3 KB View Raw
- delay_coupling_experiment.py 4.8 KB View Raw
- report.md 1.8 KB View
- report_bench_2026-08-31T123143.md 7.6 KB View
- results.json 2.6 KB View Raw