Volume-Threshold Contracting State Layer
Implementation & benchmark of arXiv:2607.27089 — Geometric Properties of Higher Dimensional Solenoidal Attractors
Source paper: Geometric Properties of Higher Dimensional Solenoidal Attractors arXiv:2607.27089 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Construct a recurrent or state-space layer as a skew product: an expanding bounded feature coordinate drives a linearly contracting hidden state. Constrain the hidden transition matrix A to have spectral radius below one, and monitor the predicted transition ell times the absolute determinant of A equals one: below it, hidden trajectories should occupy a thin or fractal set, while above it they should have substantially higher-dimensional state coverage without losing local contraction.
Formulas
Mathematical statement
The paper studies T(x,y)=(ell x,Ay+phi(x)), where x belongs to the circle, y belongs to R^d, ell is at least 2, A is an invertible d by d matrix, rho(A) is its spectral radius, and phi is a smooth forcing map. Fiber stability is controlled by rho(A)<1. The geometric phase boundary is ell times the absolute determinant of A equal to one: below this value, generic attractors have dimension below the ambient value and dimension described by an affinity-pressure root; above it, the SRB measure is generically absolutely continuous. In a neural layer, x_t is a bounded scalar or low-dimensional driver, h_t is the hidden state, and phi_theta(x_t,u_t) is a learned forcing map.
Implementation notes
Integrate this as an experimental recurrent or state-space block between an embedding layer and the output head. Represent A as Q diag(tanh(a_1),...,tanh(a_d)) Q^{-1}, or use a normal or orthogonal parameterization initially so that its spectral radius and determinant are cheap to compute; initialize eigenvalue magnitudes below one. Add a bounded driver x_t=sigmoid(w_x^T u_t+b_x), optionally followed by x_{t+1}=ell x_t modulo one, and compute h_{t+1}=A h_t+phi_theta(x_t,u_t). Pseudocode: estimate rho(A) by power iteration; compute g=log(ell)+log(abs(det(A))); if rho(A)>=1, project eigenvalue magnitudes below 1-epsilon; run the state update; record hidden-state singular values, pairwise distances, and covariance rank. The theorem-level ingredients taken from the paper are the fiber-contraction requirement rho(A)<1 and the volume exponent g; neural training must estimate effective dimension empirically. The first cheap experiment should use a width-8 recurrent layer on sequential MNIST or permuted sequential MNIST, comparing a standard RNN against layers with ell times abs(det(A)) equal to 0.5, 0.9, 1.1, and 2.0, while keeping rho(A)=0.95. The quantitative prediction is a change near log(ell)+log(abs(det(A)))=0: below the boundary, covariance participation ratio and estimated correlation dimension should remain submaximal and often collapse with sequence length; above it, state-space occupancy and covariance rank should increase, while gradients remain bounded because rho(A)<1. Test whether the empirical transition lies within 20 percent of ell times abs(det(A))=1.
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 skew-product toy experiment and an 8-unit contracting recurrent layer. The exact contraction prediction held: observed log decay -0.22314355 matched log(0.8)=-0.22314355, and the 30-step ratio matched 0.8^30. However, the predicted geometric transition at ell|det(A)|=1 was not observed: empirical box/correlation dimensions showed no clear discontinuity across r=1/sqrt(2), so the required mechanism was not confirmed. In the learning task, the standard RNN also outperformed the contracting layer (accuracy 0.8833 vs 0.8600; loss 0.2537 vs 0.3203).
- Agent confidence
- 9/10
- Baseline
- Standard tanh RNN, 130 parameters: test accuracy 0.8833, test loss 0.2537, runtime 14.96 s
- Idea
- Contracting state layer with r=0.95, bounded driver, 255 parameters: test accuracy 0.8600, test loss 0.3203, runtime 25.26 s; toy sweep box dimensions 1.823-1.935 below/above threshold and correlation dimensions 2.421-2.339, with no transition
Limitations: Only a low-dimensional isotropic A=rI toy system and a synthetic sequence-sum classification task were tested; no sequential MNIST, non-normal matrices, learned A parameterization, affinity-pressure computation, covariance-rank metric, or large-scale/generalization study was performed. Empirical dimension estimation is finite-sample and does not establish the theorem's genericity assumptions.
How to run: python3 volume_threshold_experiment.py && /home/maxwelhelp/main/bin/python3 mini_train.py
Verdict computed by deterministic test code from paired-seed statistics — not by the language model.
Artifacts
- mini_results.json 0.4 KB View Raw
- mini_train.py 3.1 KB View Raw
- report.md 1.8 KB View
- toy_results.json 2.8 KB View Raw
- volume_threshold_experiment.py 3.8 KB View Raw