✗ Failed on benchmark 2026

Hessian-Spectrum Transition Monitor and Beta Controller

Implementation & benchmark of arXiv:2608.06597 — Cascading Through the Hierarchy: Regularizer-Induced Feature Detection as Phase Transitions in Deep Linear Neural Networks

Usefulness7/10
Difficulty5/10
Novelty6/10

Source paper: Cascading Through the Hierarchy: Regularizer-Induced Feature Detection as Phase Transitions in Deep Linear Neural Networks arXiv:2608.06597 ⓘ · analyzed Aug 31, 2026

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

Idea description

Use the paper's explicit Hessian dependence on learned singular values to detect when a feature mode approaches a curvature transition, then adapt weight decay or learning rate before the mode destabilizes. This turns regularization from a static hyperparameter into feedback control based on mode-wise curvature and feature amplitude.

Formulas

lambda_l = 2[n sigma_l^(2n-2) + (2-n) beta]
lambda_l = 0 if and only if sigma_l,0 = ((n-2) beta / n)^(1/(2n-2)), for n > 2
eta_t < 2 / lambda_max,t for local stability of gradient descent on a positive-curvature quadratic mode
beta_(t+1) = clip(beta_t exp(kappa (lambda_target - lambda_min,t) / lambda_target), beta_min, beta_max)

Mathematical statement

For mode l in an n-layer deep linear network with end-to-end singular value sigma_l, the paper gives the aligned-subspace Hessian eigenvalue lambda_l = 2[n sigma_l^(2n-2) + (2-n) beta]. Here n is depth, beta is weight decay, and sigma_l is the end-to-end singular value. For n greater than 2, curvature changes sign at sigma_l,0 = ((n-2) beta / n)^(1/(2n-2)). Negative lambda_l indicates local negative curvature in that mode, while large positive lambda_l imposes a local gradient-descent step-size ceiling approximately eta < 2/lambda_max. The proposed controller keeps active modes in a target curvature band by adjusting beta and optionally the learning rate.

Implementation notes

1. Integration point: add a low-frequency controller around the optimizer for a deep linear baseline, then apply it to factorized low-rank layers such as W = U V or to individual transformer MLP and projection blocks. Every K steps compute end-to-end singular values sigma_l of each factorized block, or estimate the largest values using randomized Lanczos or power iteration. 2. Pseudocode: estimate lambda_l = 2[n sigma_l^(2n-2) + (2-n) beta_t]; set lambda_max to the largest positive estimate and lambda_min to the smallest estimate; replace the optimizer learning rate by min(eta_0, 1.5 / max(lambda_max, epsilon)); update beta_t using the displayed multiplicative controller and clip it to the allowed range. If several modes have negative lambda_l, temporarily increase beta or reduce the learning rate, and decrease beta only after negative curvature disappears. 3. Taken from the paper: the mode-wise Hessian formula and zero-curvature boundary. Estimated empirically: singular values and actual Hessian curvature, which can be checked with Hessian-vector products. 4. First cheap experiment: train depth-4 linear networks on synthetic diagonal regression and compare SGD with fixed beta and eta against the controller. Sweep eta and plot divergence probability against the predicted boundary eta_c approximately equal to 2 divided by lambda_max. Separately sweep beta and verify that each measured Hessian eigenvalue crosses zero near sigma_l,0, within 20 percent. The controller should reduce oscillation or divergence when the uncontrolled run crosses eta_c.

Verification

Failed on benchmark

Stage 1 · Toy mechanism gate: Passed ✓

Stage 2 · Mechanism transferred to benchmark: Not reproduced ✗

Stage 2 · Practical benchmark result: Significantly worse ✗

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 reproducible depth-4 aligned deep-linear monitor/controller in experiment.py. The curvature transition matched the predicted beta boundary to machine precision, and the gradient-descent stability boundary matched eta_c=2/Lambda across three Lambda values. The controller reduced predicted maximum curvature but worsened this toy training loss, so the mechanism manifested but no optimization win was demonstrated.

Agent confidence
9/10
Baseline
Depth-4 diagonal SGD, 160 steps, 3 seeds: final losses 0.06019, 0.06018, 0.06017; mean 0.06018; no divergence; beta=0.018; maximum predicted lambda about 45.10.
Idea
Feedback beta/learning-rate controller, 160 steps, 3 seeds: final losses 0.31410, 0.31422, 0.31434; mean 0.31422; no divergence; beta saturated at 0.1; maximum predicted lambda about 31.76. Curvature transition: predicted beta*=0.60801334 and observed beta*=0.60801334, relative error 1.8e-16. Stability: predicted eta_c=2/Lambda and observed transition between 0.99 eta_c and 1.01 eta_c for Lambda=1,3,10.

Limitations: The stability sweep is an exact scalar quadratic recurrence, and the Hessian check uses a stationary aligned scalar mode rather than a full matrix-network Hessian-vector-product test. The controller was tested on only one small diagonal regression setup, with hand-chosen controller gains and beta bounds; no randomized Lanczos estimates, transformer blocks, FLOP-normalized comparisons, or broad hyperparameter sweep were tested.

How to run: python3 experiment.py

Stage 2 — Benchmark vs. tuned baseline (arXiv:2608.06597)

Significantly WORSE than the tuned baseline (+5.2%, p=0.0081, wins 0 / 8 paired seeds).

Benchmark
Dynamics forecasting
Model
rnn_small
Paired seeds
8
Baseline mean
0.002
Idea mean
0.0021
Effect (Δ)
0.0001 (+5.2%; negative = idea better)
Wins
0 / 8 paired seeds
p-value
0.0081 (permutation test, 20 000 shuffles)
Smallest detectable effect
±2.8%
Mechanism
Not confirmed ✗
Practical effect
Hurts
Baseline tuning
swept over 9 configs
Limitations:

Only the registered dynamics track was tested. No vision, tabular, sequence, custom track, speed/FLOP analysis, randomized Lanczos estimate, or full Hessian spectrum was evaluated. The GRU is not the paper's exact deep-linear aligned factorization, so its monitored curvature is an empirical spectral proxy rather than a direct application of the theorem.

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: 10/10

Implemented and ran the Hessian-spectrum monitor/controller on the structurally matched registered dynamics track with rnn_small, eight paired seeds, equal-budget baseline tuning, and shared learning-rate search space. The tuned baseline achieved mean test MSE 0.00198587 versus 0.00208846 for the controller; paired delta was +0.00010259 with permutation p=0.0081, so the idea was significantly worse. The trained-model mechanism signature was not confirmed: predicted curvature proxy and HVP Rayleigh curvature had correlation 0.2723.

Baseline
dynamics/rnn_small; sweep lr={0.001,0.003,0.006} and weight_decay={0,0.0001,0.0005}; best cfg {lr:0.006,weight_decay:0.0005}; 8-seed mean test MSE 0.001985870083444752.
Idea
Feedback controller; best cfg {lr:0.006,weight_decay:0.0005}; 8-seed mean test MSE 0.0020884609839413315; paired delta +0.0001025909004965797; permutation p=0.0081; verdict idea worse (significant).

Benchmark result

Significantly WORSE than the tuned baseline (+5.2%, p=0.0081, wins 0 / 8 paired seeds).

Benchmark
Dynamics forecasting
Model
rnn_small
Paired seeds
8
Baseline mean
0.002
Idea mean
0.0021
Effect (Δ)
0.0001 (+5.2%; negative = idea better)
Wins
0 / 8 paired seeds
p-value
0.0081 (permutation test, 20 000 shuffles)
Smallest detectable effect
±2.8%
Mechanism
Not confirmed ✗
Practical effect
Hurts
Baseline tuning
swept over 9 configs
Limitations:

Only the registered dynamics track was tested. No vision, tabular, sequence, custom track, speed/FLOP analysis, randomized Lanczos estimate, or full Hessian spectrum was evaluated. The GRU is not the paper's exact deep-linear aligned factorization, so its monitored curvature is an empirical spectral proxy rather than a direct application of the theorem.

How to run: python3 bench_experiment.py

Artifacts

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