Fractional-memory recurrent state
Implementation & benchmark of arXiv:2608.21674 — Anomalous diffusion memory factorization: Characteristic timescales and application to inverse problem
Source paper: Anomalous diffusion memory factorization: Characteristic timescales and application to inverse problem arXiv:2608.21674 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Construct an efficient recurrent or state-space layer whose impulse response follows Mittag-Leffler relaxation instead of a single exponential. A bank of stable diagonal state channels approximates the long power-law tail, allowing the layer to retain information over widely separated timescales with only \(K\) states per feature.
Formulas
Mathematical statement
The paper's temporal building block is \(e_\alpha(t)=E_\alpha(-t^\alpha)=\sum_{n=0}^{\infty}(-t^\alpha)^n/\Gamma(\alpha n+1)\), with fractional order \(\alpha\) and time \(t\). Its stated asymptotic \(e_\alpha(t)\sim t^{-\alpha}/\Gamma(1-\alpha)\) for \(0<\alpha<1\) gives a slowly decaying memory kernel. For a discrete neural sequence \(x_1,\ldots,x_T\), define the desired causal normalized filter \(y_t=\sum_{j=1}^{t}h_{t-j}x_j/\sum_{\ell=0}^{t-1}h_\ell\), where \(h_\ell=E_\alpha(-((\ell+1)/\tau)^\alpha)\), \(\tau>0\) is a characteristic timescale, and \(\ell=t-j\) is the lag. Implement it using a positive exponential mixture \(h_\ell\approx\hat h_\ell=\sum_{k=1}^{K}w_k(1-\rho_k)\rho_k^\ell\). Each component has state \(s_{t,k}=\rho_ks_{t-1,k}+(1-\rho_k)x_t\), and the output is \(y_t=\sum_kw_ks_{t,k}\); stability follows from \(0<\rho_k<1\), while logarithmically spaced \(\rho_k\) approximate the fractional power-law memory.
Implementation notes
Add the module as a drop-in replacement for the temporal mixing block in a small causal Transformer, GRU, or diagonal SSM. For an input sequence tensor \(x\in\mathbb{R}^{T\times B\times d}\), create \(K\) state tensors \(s_k\in\mathbb{R}^{B\times d}\); use \(K=8\), positive weights normalized by softmax, and decay factors \(\rho_k=\exp(-\Delta/\tau_k)\) with logarithmically spaced \(\tau_k\) from 1 to the target context length. At every timestep, compute s[k] = rho[k] * s[k] + (1-rho[k]) * x[t], then y[t] = sum(w[k] * s[k]); apply a learned input/output projection and residual connection around the filter. Optionally learn \(\alpha\) and \(\tau\) only through a differentiable initialization map that refits \(w_k,\rho_k\), but begin with fixed stable decays. Fit the bank before training by minimizing squared error between \(\hat h_\ell\) and the paper kernel \(E_\alpha(-((\ell+1)/\tau)^\alpha)\) over logarithmic lags; this fit is the empirical approximation, while the power-law asymptotic is the mathematical motivation. The first cheap experiment should compare this layer with a GRU, a single-exponential SSM, and a standard causal Transformer on synthetic delayed-copy and adding problems at sequence lengths 1k, 4k, and 16k, followed by TinyStories language modeling. Track accuracy as delay/context grows, validation perplexity at equal parameter count, and memory/runtime. Success is materially better delayed-copy accuracy and less perplexity degradation at long context with \(O(KTd)\) compute and constant recurrent state.
Verification
Stage 1 · Toy mechanism gate: Failed ✗
Stage 2 · Mechanism transferred to benchmark: Not tested
Stage 2 · Practical benchmark result: Not run
Methodology: Toy-system gate first; the benchmark stage runs only after a pass. How verification works
Stage 1 — Mechanism check agent confidence 9/10
Built a stable PyTorch fractional-memory layer using positive softmax-weighted exponential states with logarithmically spaced time constants, plus numerical verification scripts. The Mittag-Leffler power-law asymptotic matched the predicted 1/sqrt(pi) constant with relative error falling from 0.49% at lag 100 to 0.00005% at lag 1e6, and exponential stability/half-life predictions matched discrete observations within one step. However, bank approximation error improved strongly from K=1 to K=4 but then saturated rather than continuing to improve, and the fractional K=8 filter had higher random-signal MSE than the best single exponential; therefore the full claimed mechanism was not demonstrated as a win.
- Agent confidence
- 9/10
- Baseline
- Best single exponential: MSE 0.0001200 against the exact normalized fractional filter; delayed impulse at lag 4000: 1.29e-58.
- Idea
- Fractional K=8 bank: MSE 0.0006433; delayed impulse at lag 4000: 1.52e-05, approximately 1.17e53 times larger retention than the tau=32 baseline. Relative kernel-fit error by K: 1=101.99, 2=0.818, 4=0.338, 8=0.314, 16=0.313, 32=0.313.
Limitations: No trained delayed-copy model, GRU, Transformer, language-model, equal-parameter benchmark, GPU benchmark, or alpha values other than 1/2 were tested. The bank fit used fixed logarithmic time constants and a finite horizon, so improved fitting or learned decays remain unexplored.
How to run: python3 verify.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- experiment.py 5.1 KB View Raw
- fractional_memory.py 1.5 KB View Raw
- report.md 1.7 KB View
- results.json 2.5 KB View Raw
- verify.py 3.6 KB View Raw