✗ Mechanism failed 2026

Excitation-Controlled Recurrent Learning

Implementation & benchmark of arXiv:2608.09722 — Decoding Gene Regulatory Networks from Single-Cell RNA Velocity

Usefulness7/10

Difficulty6/10

Novelty8/10

Dynamics Sample efficiency Accuracy Stability Rnn Ssm World model Curriculum Training dynamics Dynamical systems Control theory Linear algebra Optimization

Source paper: Decoding Gene Regulatory Networks from Single-Cell RNA Velocity arXiv:2608.09722 ⓘ · analyzed Sep 1, 2026

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

Idea description

Expose a recurrent model to deliberately designed input pulses or latent-state perturbations instead of training only on passive trajectories. Choose perturbations that maximize the smallest eigenvalue of the accumulated feature Gramian, making otherwise indistinguishable recurrent couplings recoverable and reducing uncertainty in long-horizon predictions.

Formulas

$$\dot{x}(t)=F_{\theta}(x(t))+Bv(t),\qquad \frac{\partial\Psi(t,\tau)}{\partial t}=J(t)\Psi(t,\tau),\quad J(t)=\frac{\partial F_{\theta}}{\partial x}(x(t)).$$

$$W_T=\int_{0}^{T}\Psi(T,\tau)BB^{\top}\Psi(T,\tau)^{\top}\,d\tau.$$

$$v^{\star}=\arg\max_{v\in\mathcal{V}}\lambda_{\min}(G(v)+\epsilon I),\qquad G(v)=\frac{1}{N}\Phi(v)^{\top}\Phi(v).$$

Mathematical statement

For hidden dynamics $\dot{x}=F_{\theta}(x)+Bv(t)$, let $J(t)=\partial F_{\theta}/\partial x$ along a nominal trajectory and let $\Psi(t,\tau)$ satisfy $\partial_t\Psi(t,\tau)=J(t)\Psi(t,\tau)$, $\Psi(\tau,\tau)=I$. The finite-horizon controllability or excitation Gramian is $W_T=\int_0^T\Psi(T,\tau)BB^{\top}\Psi(T,\tau)^{\top}d\tau$. Perturbations are informative when $\lambda_{\min}(W_T)$ is bounded away from zero. For neural identification, use an empirical feature Gram matrix $G(v)=\Phi(v)^{\top}\Phi(v)/N$, where $\Phi(v)$ contains trajectory features or parameter sensitivities under control sequence $v$, and select interventions by maximizing $\lambda_{\min}(G+\epsilon I)$ or $\log\det(G+\epsilon I)$.

Implementation notes

Integrate this as an active-data-collection controller or curriculum module for an RNN, neural state-space model, or world model. During an initial passive rollout, estimate the hidden-state Jacobian using automatic differentiation and identify weakly excited directions. Inject bounded pulses, colored noise, or channel-specific interventions into observable inputs or selected latent coordinates. For each candidate control sequence $v^{(m)}$, roll out the current model, compute sensitivities $\partial x_j/\partial\theta_S$ or candidate regulator features, and construct $G(v)=\Phi(v)^\top\Phi(v)/N$. Select the candidate maximizing $\lambda_{\min}(G+\epsilon I)$ subject to $\|v\|_2\leq c$, collect that trajectory, and train with the integral sparse loss. The paper contributes the identifiability and Gramian criterion; the implementation estimates $J$, $\Phi$, and $G$ with autodiff and minibatch rollouts. First test a 20-dimensional sparse tanh RNN with only five observed coordinates, comparing passive Gaussian inputs with greedy pulse design at equal trajectory budgets. The prediction is that interventions increase $\lambda_{\min}(G)$ by at least 2x and reduce parameter error; when $\lambda_{\min}(G)$ is near zero, edge recovery and long-horizon prediction should fail sharply, while designed perturbations move the model above this boundary.

Verification

Mechanism failed

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 reproducible linear recurrent identification MVP with empirical feature Gramians, greedy pulse selection, parameter recovery, and fixed-seed sweeps. The mechanism checks passed: Gramian minimum eigenvalue scaled as amplitude squared, vanished at an unexcited channel, and approximately followed b squared scaling. However, under equal input energy the designed pulse set had lower lambda_min than passive excitation (0.05425 vs 0.07168), so the central promised Gramian improvement was not observed; active data nevertheless reduced parameter error.

Agent confidence: 9/10

Baseline: Equal-energy passive trajectories: lambda_min(G)=0.071684; parameter error 1.638 at noise 0.02 and 1.637 at noise 0.1.
Idea: Greedy equal-energy bounded pulses: lambda_min(G)=0.054250, ratio 0.757; parameter error 0.630 at noise 0.02 and 0.634 at noise 0.1. Amplitude sweep observed normalized [0.0625000000, 0.2500000000, 1, 4.0000000000] versus predicted [0.0625, 0.25, 1, 4].

Limitations: This is a 2D linear recurrent toy, not the proposed 20D sparse tanh RNN with five observed coordinates. No online neural autodiff Jacobian estimation, latent interventions, sparse loss, long-horizon prediction, FLOP accounting, or robustness to model misspecification was tested.

How to run: python3 experiment.py

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

Artifacts

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

🐍 experiment.py 2.5 KB View Raw
📄 report.md 1.5 KB View
📊 results.json 1.5 KB View Raw