Unverified 2026

Curved-Critical Residual Lattice

Implementation & benchmark of arXiv:2608.02903 — Two dimensional inhomogeneous classical systems at criticality

Usefulness6/10
Difficulty6/10
Novelty8/10

Source paper: Two dimensional inhomogeneous classical systems at criticality arXiv:2608.02903 ⓘ · analyzed Aug 31, 2026

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

Idea description

Construct a 2D recurrent or residual neural lattice with slowly varying local couplings, while parameterizing those couplings so that an anisotropy invariant remains constant across all spatial and depth locations. The network obtains controlled local propagation velocities rather than arbitrary inhomogeneous amplification, enabling depth-dependent receptive fields while preserving near-critical signal propagation.

Formulas

$$\Delta_{jk}=\frac{a_{jk}^{2}+b_{jk}^{2}-c_{jk}^{2}}{2a_{jk}b_{jk}}=\Delta,\qquad a_{jk}=\sin(\lambda_j-\mu_k+\gamma),\quad b_{jk}=\sin(\lambda_j-\mu_k),\quad c_{jk}=\sin\gamma.$$
$$A_{jk}=\frac{a_{jk}}{s_{jk}},\quad B_{jk}=\frac{b_{jk}}{s_{jk}},\quad C_{jk}=\frac{c_{jk}}{s_{jk}},\quad s_{jk}=|a_{jk}|+|b_{jk}|+|c_{jk}|.$$
$$x_{jk}^{\mathrm{new}}=\phi\!\left(A_{jk}x_{j-1,k}+B_{jk}x_{j+1,k}+C_{jk}x_{j,k-1}+u_{jk}\right),\qquad ds^{2}=dk^{2}+v_{jk}^{2}dj^{2},\quad v_{jk}=\left|\frac{A_{jk}}{B_{jk}}\right|.$$
$$\mathcal{L}_{\mathrm{total}}=\mathcal{L}_{\mathrm{task}}+\alpha\sum_{j,k}\left[(\lambda_{j+1}-\lambda_j)^2+(\mu_{k+1}-\mu_k)^2\right]+\beta\sum_{j,k}(\Delta_{jk}-\bar\Delta)^2.$$

Mathematical statement

The six-vertex construction keeps the local anisotropy invariant \(\Delta_{jk}=(a_{jk}^{2}+b_{jk}^{2}-c_{jk}^{2})/(2a_{jk}b_{jk})\) constant, with \(a_{jk}=\sin(\lambda_j-\mu_k+\gamma)\), \(b_{jk}=\sin(\lambda_j-\mu_k)\), and \(c_{jk}=\sin\gamma\). Here \(j\) and \(k\) index the two lattice directions, \(\lambda_j\) and \(\mu_k\) are slowly varying spectral parameters, and \(\gamma\) fixes the interaction and therefore \(\Delta\). Transfer this by defining normalized neural couplings \(A_{jk}=a_{jk}/s_{jk}\), \(B_{jk}=b_{jk}/s_{jk}\), and \(C_{jk}=c_{jk}/s_{jk}\), where \(s_{jk}=|a_{jk}|+|b_{jk}|+|c_{jk}|\). The recurrent state update is \(x_{jk}^{\mathrm{new}}=\phi(A_{jk}x_{j-1,k}+B_{jk}x_{j+1,k}+C_{jk}x_{j,k-1}+u_{jk})\). The invariant constrains local anisotropy, while the slow-variation condition \(|\lambda_{j+1}-\lambda_j|\ll 1\) and \(|\mu_{k+1}-\mu_k|\ll 1\) suppresses abrupt signal scattering. An engineering approximation to the induced local metric is \(ds^{2}=dk^{2}+v_{jk}^{2}dj^{2}\), with propagation-speed ratio \(v_{jk}=|A_{jk}/B_{jk}|\).

Implementation notes

1. Exact integration point: implement the model as a recurrent 2D lattice, a deep-equilibrium layer, or a stack of residual blocks indexed by horizontal position \(j\) and depth/time \(k\). Replace ordinary convolution or recurrent mixing with the four-neighbor update above. Use trainable vectors \(\lambda_j\) and \(\mu_k\), constrain their increments with the smoothness penalty, and fix \(\gamma\) initially. Restrict parameters to a range where the sine weights have consistent signs, or retain absolute values in the normalization. 2. Pseudocode: compute \(q_{jk}=\lambda_j-\mu_k\); set \(a=\sin(q+\gamma)\), \(b=\sin q\), and \(c=\sin\gamma\); normalize by \(s=|a|+|b|+|c|\); update every lattice state with the weighted neighbor sum; add the smoothness and invariant penalties to the task loss. 3. Computed versus estimated: \(a,b,c,\Delta_{jk}\), and \(v_{jk}\) are computed exactly at every step. The continuum metric is an empirical diagnostic; estimate it by injecting impulses at selected lattice sites and fitting the displacement of the largest response over depth. Estimate the largest Jacobian singular value with ten-step power iteration. 4. First cheap experiment: use a 32-layer, 128-token sequence model represented as a \(32\times128\) lattice on permuted-MNIST, comparing a standard residual RNN, an unconstrained inhomogeneous lattice, and this invariant parameterization at equal parameter count. Sweep coupling amplitude and record gradient norms, impulse-response width, validation accuracy, and estimated Jacobian singular values. The predicted signature is a broad near-critical regime where the per-layer Jacobian singular value remains near one and impulse correlations decay much more slowly than exponentially. Unconstrained inhomogeneity should show exploding or frozen propagation earlier. The transition to frozen behavior should occur when the measured local response gain falls below approximately one, while the constrained model should preserve its transition boundary within 20 percent as variation increases until the parameter increments become order one.

Verification

This idea has not been verified yet.

Verification happens in two stages: Stage 1 — a mechanism check on a toy system confirms the claimed mathematical phenomenon reproduces; Stage 2 — a benchmark implements the idea on a real (small) neural network task and compares it against a tuned baseline over 8 paired seeds with a permutation test.

Artifacts

Artifacts unavailable.