Positive-envelope stability for complex state updates
Implementation & benchmark of arXiv:2608.27939 — Entropy and domination for quasi-Hitchin representations
Source paper: Entropy and domination for quasi-Hitchin representations arXiv:2608.27939 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
For a complex-valued recurrent or state-space layer, construct a positive envelope by replacing each factor matrix with its entrywise modulus. The envelope provably upper-bounds every entry of the complex product and therefore gives a cheap conservative estimate of worst-case amplification, while a learned phase-cancellation term can exploit complex interference without allowing unstable growth.
Formulas
Mathematical statement
The paper represents a matrix entry as a sum of weighted paths in a planar network. If $B$ is the complex matrix and $B_0$ is obtained by replacing every network parameter by its modulus, then each entry satisfies $|(B)_{ij}|\leq (B_0)_{ij}$ because $B_{ij}=\sum_{p:i\to j}w_p$ and $(B_0)_{ij}=\sum_{p:i\to j}|w_p|$. For a product $A=A_L\cdots A_1$, define the positive envelope $P=|A_L|\cdots |A_1|$, where $|A_k|$ is entrywise modulus. Repeated application of the triangle inequality gives $|A|\leq P$ entrywise. Consequently, for any vector $x$, $|Ax|\leq P|x|$ componentwise, and any induced monotone norm obeys $\|A\|\leq\|P\|$. Here $A_k$ are complex transition factors, $P$ is a nonnegative matrix, $x$ is an input or hidden state, and $\rho(P)$ is the spectral radius of the envelope. Strict inequality can occur when different complex path weights have incompatible phases, so the envelope can distinguish safe cancellation from unsafe magnitude growth.
Implementation notes
Integrate this into a complex-valued SSM or GRU-style recurrent layer whose transition over one chunk is factored as $A=A_L\cdots A_1$, preferably using sparse triangular, diagonal, or low-rank factors so the envelope is cheap. At every optimizer step, compute the usual complex transition $A$ and, in parallel, compute $P=|A_L|\cdots|A_1|$ using differentiable modulus operations. Use the following pseudocode: P=I; A=I; for k in 1..L: A=A_k@A; P=abs(A_k)@P; y=A@x; v=random_positive_vector(); for j in 1..3: v=P@v; v=v/(norm(v,2)+eps); rhoP=(v@(P@v))/(v@v+eps); loss=task_loss(y,target)+lambda_stab*relu(log(rhoP)-log(rho_target))**2 + beta*(norm(P,1)-norm(A,1))/(norm(P,1)+eps); backprop(loss). The paper-derived quantities are the modulus product $P$, the componentwise inequality $|A|\leq P$, and the spectral-radius upper-bound heuristic; the target radius, penalty weights, and number of power iterations are empirical choices. Begin with a 2-layer complex diagonal-plus-low-rank SSM on sequential MNIST or the adding problem, comparing unconstrained complex training, ordinary spectral normalization, and the envelope penalty at equal hidden size. Measure exploding-gradient frequency, validation loss, long-horizon rollout error, and wall-clock overhead. Success is fewer unstable runs and lower long-horizon error with less than 10% training overhead.
Verification
Stage 1 · Toy mechanism gate: Passed ✓
Stage 2 · Mechanism transferred to benchmark: Confirmed ✓
Stage 2 · Practical benchmark result: Inconclusive
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 a complex two-factor recurrent transition with a differentiable positive envelope, numerical verification, and a matched CUDA training benchmark. The mechanism manifested: across 5,000 random products the maximum entrywise-bound violation was 0, phase-cancellation ratios matched |cos(theta/2)| to floating-point precision, and scalar growth thresholds matched predictions. The envelope model also showed slightly lower long-rollout RMSE with controlled envelope radius, though this is only a small single-seed signal.
- Agent confidence
- 9/10
- Baseline
- train loss 9.113e-09; long-rollout RMSE 1.019e-04; unstable steps 0; rho(A)=0.481; rho(P)=1.366; 14.04 s
- Idea
- train loss 8.538e-09; long-rollout RMSE 9.614e-05; unstable steps 0; rho(A)=0.602; rho(P)=0.636; 14.91 s; 6.2% overhead
Limitations: Only one small synthetic complex recurrence and one seed were used. Sequential MNIST, the adding problem, ordinary spectral normalization, multi-seed statistics, stress tests, and a larger SSM were not evaluated. The benchmark uses a row-sum envelope penalty as a proxy rather than the full power-iteration spectral-radius loss.
How to run: python3 experiment.py
Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.27939)
No significant difference vs. the tuned baseline (+6.7%, p=0.7833, wins 6 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- factor_rnn
- Paired seeds
- 8
- Baseline mean
- 0.0211
- Idea mean
- 0.0225
- Effect (Δ)
- 0.0014 (+6.7%; negative = idea better)
- Wins
- 6 / 8 paired seeds
- p-value
- 0.7833 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±47.2%
- Mechanism
- Confirmed ✓
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
The fixed built-in rnn_small interface is real-valued, so the experiment used a real two-factor recurrent analogue of the complex transition rather than a fully complex GRU/SSM. The run tested only the built-in 400-sample dynamics task, 18 epochs, three learning rates, and three envelope penalty settings; ordinary spectral normalization, longer horizons, and a larger complex SSM were not tested.
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
The positive-envelope intervention was implemented and evaluated on the structurally matched built-in dynamics track using end-to-end trained, architecture-matched recurrent systems. The trained-model envelope mechanism was quantitatively confirmed, but the best idea configuration was slightly worse on independent test MSE and not statistically significant; therefore the idea did not win this bench round.
- Baseline
- Best baseline lr=0.006, 18 epochs; test MSE mean 0.02109241997823119 over 8 paired seeds.
- Delta
- idea - baseline = +0.001420103246346116 MSE; p=0.78325; no significant win.
- Idea
- Best envelope setting lr=0.006, beta=0.02, target=0.92, 18 epochs; test MSE mean 0.022512523224577308 over 8 paired seeds.
- Mechanism signature
- confirmed=true; maximum observed componentwise ratio |Ax|/(P|x|)=0.7172038555; maximum violation=-0.35021981597.
Benchmark result
No significant difference vs. the tuned baseline (+6.7%, p=0.7833, wins 6 / 8 paired seeds).
- Benchmark
- Dynamics forecasting
- Model
- factor_rnn
- Paired seeds
- 8
- Baseline mean
- 0.0211
- Idea mean
- 0.0225
- Effect (Δ)
- 0.0014 (+6.7%; negative = idea better)
- Wins
- 6 / 8 paired seeds
- p-value
- 0.7833 (permutation test, 20 000 shuffles)
- Smallest detectable effect
- ±47.2%
- Mechanism
- Confirmed ✓
- Practical effect
- Inconclusive
- Baseline tuning
- swept over 3 configs
The fixed built-in rnn_small interface is real-valued, so the experiment used a real two-factor recurrent analogue of the complex transition rather than a fully complex GRU/SSM. The run tested only the built-in 400-sample dynamics task, 18 epochs, three learning rates, and three envelope penalty settings; ordinary spectral normalization, longer horizons, and a larger complex SSM were not tested.
How to run: python3 bench_experiment.py
Artifacts
- bench_experiment.py 6.3 KB View Raw
- bench_report.json 5.5 KB View Raw
- experiment.py 7.0 KB View Raw
- report.md 1.4 KB View
- report_bench_2026-09-02T005838.md 6.0 KB View
- results.json 3.1 KB View Raw