Variance-budgeted stochastic momentum
Implementation & benchmark of arXiv:2608.04245 — Strong error analysis for the stochastic momentum optimizer
Source paper: Strong error analysis for the stochastic momentum optimizer arXiv:2608.04245 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Replace fixed momentum with an online controller that selects the momentum coefficient from an upper bound on the next-step momentum second moment. The controller lowers momentum when minibatch noise dominates and permits higher momentum when the gradient estimate is stable.
Formulas
Mathematical statement
The paper assumes the one-point convexity inequality <θ-ϑ,g(θ)> ≤ -c max{||θ-ϑ||²,||g(θ)||²}, where θ is the current parameter vector, ϑ is a target point, g(θ)=E[G(θ,U)] is the population stochastic-gradient field, and c>0 is the one-point convexity constant. Momentum is m_n=αm_{n-1}+(1-α)Ḡ_n, where α∈[0,1), Ḡ_n=J^{-1}∑_{j=1}^J G(Θ_{n-1},U_{n,j}), and J is the minibatch size. Lemma 2.3 gives E[||m_n||²] ≤ αE[||m_{n-1}||²] + ((1-α)(J-1+α)/(Jc²))E[||Θ_{n-1}-ϑ||²] + (1-α)²𝔠²/J, assuming E[||G(θ,U)||²]≤𝔠². The proof also gives E[||ε_n||²] = J^{-1}(E[||G(Θ_{n-1},U)||²]-E[||g(Θ_{n-1})||²]) ≤ 𝔠²/J, where ε_n=Ḡ_n-g(Θ_{n-1}) is minibatch noise. In a neural network, ϑ and c are unavailable, so estimate local displacement, gradient variance, and gradient alignment empirically, then use the recursion as a clipped heuristic budget rather than as a literal theorem guarantee.
Implementation notes
(1) Integrate the method into the momentum-buffer update of SGD, or replace AdamW's fixed first-moment coefficient β1 while leaving its second-moment estimator unchanged. Maintain a momentum buffer m_t, an exponential estimate of its squared norm, an estimate of minibatch gradient variance, and a short parameter displacement. (2) Pseudocode: compute two independent microbatch gradients g1 and g2; set g=(g1+g2)/2 and sigma_hat²=||g1-g2||²/2. Maintain M_hat=EMA(||m_t||²), displacement d=theta_t-theta_{t-K}, and D_hat=||d||²/(K²+epsilon). Estimate local alignment with c_hat=clip(abs(<d,g_t-g_{t-K}>)/(||d||²+epsilon),c_min,c_max). Set a target momentum-energy budget B_t=q(||g||²+epsilon). For each candidate alpha in a grid from alpha_min to alpha_max, compute B_hat(alpha) using the displayed adapted inequality and choose the largest candidate satisfying B_hat(alpha)≤B_t. Update m_{t+1}=alpha_t m_t+(1-alpha_t)g, then theta_{t+1}=theta_t-eta_t m_{t+1}. Smooth sigma_hat, c_hat, and D_hat with EMAs and use conservative clipping because the theorem's target point and convexity constant are not known. (3) The paper supplies the recursion and the 1/J noise scaling; all estimates of c, distance to the target, and the noise bound are empirical proxies. (4) First experiment: train ResNet-18 on CIFAR-10 and a small Transformer on WikiText-2, comparing fixed-momentum SGD and AdamW against the controller at identical batch size, learning-rate schedule, and update count. Track loss at equal FLOPs, gradient-noise estimates, momentum norm, divergence or overflow frequency, and final validation accuracy or perplexity. The expected signal is fewer loss spikes and fewer unstable runs during noisy early training, with equal or faster loss reduction and no late-training accuracy loss.
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.