Fourth-Mass Regularization for Signed Projections
Implementation & benchmark of arXiv:2608.17802 — Fourth-Moment Geometry of Rademacher Sums
Source paper: Fourth-Moment Geometry of Rademacher Sums arXiv:2608.17802 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Add a differentiable fourth-order-mass penalty to coefficient vectors used by randomized signed projections, sign-noise layers, or stochastic quantizers. The penalty controls the effective number of active coefficients and therefore the distribution of the injected random fluctuation: diffuse coefficients generate nearly Gaussian perturbations, whereas concentrated coefficients generate larger non-Gaussian deviations.
Formulas
Mathematical statement
For independent Rademacher signs \(\varepsilon_i\in\{-1,+1\}\) and coefficients \(a=(a_1,\ldots,a_n)\in\mathbb{R}^n\), define the signed sum \(S=\sum_{i=1}^n a_i\varepsilon_i\). Under the paper's normalization \(\sum_i a_i^2=1\), the classical Gaussian comparison is \(\mathbb{E}|S|^p\leq \mu_p\) for \(p>2\), where \(\mu_p=\mathbb{E}|G|^p\) and \(G\sim\mathcal{N}(0,1)\). The exact fourth-moment identity is \(\mathbb{E}S^4=3-2\sum_i a_i^4\). Let \(m_4(a)=\sum_i a_i^4\); then \(m_4\in[1/n,1]\), and the effective dimension is \(d_{\mathrm{eff}}=1/m_4\). The fourth moment decreases from \(1\) for a one-hot coefficient vector to \(3-2/n\) for a flat vector. Use \(m_4\) as a differentiable proxy for the coefficient-sensitive higher-moment and tail behavior studied in the paper. Its gradient is \(\partial m_4/\partial a_i=4a_i^3\).
Implementation notes
Integrate this at a stochastic signed linear projection or quantization module, not at an arbitrary dense layer. For an input tensor \(x\in\mathbb{R}^{B\times n}\), use a coefficient vector \(a\in\mathbb{R}^n\) to form \(z_b=\sum_i a_i x_{b,i}\varepsilon_{b,i}\), with fresh independent signs \(\varepsilon_{b,i}\) per batch element or per training step. Parameterize \(a_i=u_i/\sqrt{\sum_j u_j^2+10^{-8}}\) so the required \(\ell_2\) normalization is automatic. Choose a target effective dimension \(d_0\), set \(\tau=1/d_0\), and add \(\lambda(\sum_i a_i^4-\tau)^2\) to the task loss. Pseudocode: compute a = u / sqrt(sum(u*u)+eps); compute m4 = sum(a**4); sample eps_sign = 2*Bernoulli(0.5)-1; compute z = sum(a*x*eps_sign, dim=-1); optimize task_loss + lambda*(m4-tau)**2. The paper supplies the exact fourth-moment dependence; the target \(d_0\), penalty weight \(\lambda\), and whether signs are per-example or shared are empirical choices. Estimate realized behavior by measuring batch excess kurtosis and comparing it with \(3-2m_4\) after normalizing each projection to unit variance. First test on CIFAR-10 with a small MLP or ResNet-18 containing one signed projection, comparing no penalty, an \(\ell_2\)-only baseline, and the fourth-mass penalty at matched parameter count and FLOPs. Sweep \(d_0\in\{4,8,16,n\}\). Success is lower activation-gradient variance, fewer unstable runs, and equal or better validation accuracy at the same injected-noise level; a secondary test should compare quantization error or perturbation robustness against a dense Gaussian-noise baseline.
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.