Onsager–Casimir Response Regularizer
Implementation & benchmark of arXiv:2608.00344 — Memory with Onsager-Casimir symmetry: Rotating particle in a viscoelastic fluid
Source paper: Memory with Onsager-Casimir symmetry: Rotating particle in a viscoelastic fluid arXiv:2608.00344 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Train a sequence model so that measured perturbation responses and spontaneous hidden-state correlations satisfy the paper's off-diagonal fluctuation–response identity. This discourages arbitrary non-reciprocal dynamics while preserving a controlled antisymmetric response that can encode directional temporal dependencies.
Formulas
Mathematical statement
Let $\mathbb{C}(t,t')$ be the two-time covariance matrix of a vector hidden representation and let $\boldsymbol{\chi}(t)$ be its causal linear response to an impulse. For $t>t'$, the paper relates the covariance derivative to the difference between delayed and current responses; for $t<t'$, it uses the transpose of the reversed-time response. In a neural model, replace $k_BT$ by an effective noise scale $\tau>0$. Estimate $\chi_k$ with a Jacobian-vector product or finite difference. For a signed orientation variable $s\in\{-1,+1\}$, enforce Onsager–Casimir reversal on the antisymmetric response $\chi_k^A=(\chi_k-\chi_k^T)/2$ by requiring $\chi_k^A(s)\approx-\chi_k^A(-s)$.
Implementation notes
Apply the regularizer to hidden states of an RNN, SSM, transformer with a recurrent memory cache, or world model. During each training batch, retain hidden vectors $h_0,\ldots,h_T$ and compute empirical lagged covariances $\mathbb{C}_{k,j}=B^{-1}\sum_b(h_k^{(b)}-\bar h_k)(h_j^{(b)}-\bar h_j)^T$, using covariance shrinkage when the batch is small. Estimate the response $\chi_q$ by injecting a small impulse $\delta u$ into the input or control channel and calculating $[h_q(u+\delta u)-h_q(u)]/\delta u$; use Jacobian-vector products to avoid constructing a full Jacobian. Add $\lambda_{\mathrm{FDR}}\mathcal{L}_{\mathrm{FDR}}+\lambda_{\mathrm{OC}}\mathcal{L}_{\mathrm{OC}}$ to the task loss. The paper supplies the correlation-response structure; $\tau$ is an effective neural temperature estimated from injected input-noise variance or tuned on validation data. The first cheap experiment should use a two-dimensional synthetic rotating linear system with known transition $A=-\gamma I+\omega J$, followed by delayed-MNIST with a compact GRU. Compare an unregularized baseline with regularized models at equal compute. The prediction is that normalized FDR residual decreases monotonically for small-to-moderate $\lambda_{\mathrm{FDR}}$, while antisymmetric lag correlation changes sign under $s\to-s$ and approaches zero as $\omega\to0$. Excessive regularization should cause a measurable U-shaped validation-loss curve.
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 reproducible 2D rotating linear-system experiment with OC reversal, antisymmetric-response scaling, Euler stability, FDR residual, and a small equal-data transition-fit comparison. OC reversal was exact, antisymmetric response vanished at omega=0 and matched predicted small-omega slopes within 0.06%, and the stability boundary was observed at omega=0.65 versus predicted 0.6614. However, the supplied FDR formula failed for this discrete noisy system: predicted tau=0.49 was not optimal (grid optimum=0.25), with normalized residuals near 1, so the overall mechanism criterion is not met.
- Agent confidence
- 9/10
- Baseline
- Transition-fit test MSE 2.287e-6; OC residual 1.437e-3
- Idea
- At lambda=1.0, test MSE 1.480e-6 and OC residual 7.229e-4; OC residual decreased monotonically across tested lambdas, but this does not offset the failed FDR check
Limitations: Only the analytic 2D linear toy system and tiny linear fitting comparison were tested; no GRU, delayed-MNIST, Jacobian-vector-product implementation, covariance shrinkage, continuous-time integrator, or neural sequence benchmark was evaluated. The FDR expression was tested in the stated discrete-time Euler setup and appears incompatible without a more careful discrete-time response/covariance derivation.
How to run: python3 oc_experiment.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- oc_experiment.py 5.8 KB View Raw
- report.md 1.5 KB View
- results.json 4.6 KB View Raw