Lorentzian coefficient regularization
Implementation & benchmark of arXiv:2608.12266 — Normalized skew Schur polynomials are Lorentzian
Source paper: Normalized skew Schur polynomials are Lorentzian arXiv:2608.12266 ⓘ · analyzed Sep 1, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Represent a nonnegative neural output as a homogeneous polynomial with coefficients indexed by count vectors, and penalize violations of the Lorentzian Hessian signature after factorial normalization. Add an M-convex support penalty so mass can move between coordinates through valid exchange operations rather than forming disconnected or brittle coefficient patterns.
Formulas
Mathematical statement
The paper defines the normalization operator for f(x)=\sum_{\alpha\in\mathbb{Z}_{\geq0}^n}c_\alpha x^\alpha by \mathcal{N}(f)=\sum_\alpha(c_\alpha/\alpha!)x^\alpha, where \alpha=(\alpha_1,\ldots,\alpha_n), x^\alpha=\prod_i x_i^{\alpha_i}, and \alpha!=\prod_i\alpha_i!. Its support is \operatorname{supp}(f)=\{\alpha:c_\alpha
eq0\}. A support J is M-convex when for any \alpha,\beta\in J with \alpha_i>\beta_i, there exists j with \alpha_j<\beta_j such that both exchanged vectors \alpha-\mathbf e_i+\mathbf e_j and \beta-\mathbf e_j+\mathbf e_i remain in J; \mathbf e_i is the i-th standard basis vector. The Lorentzian characterization can be operationalized for a homogeneous degree-d polynomial g=\mathcal N(f): for every multi-index \gamma with |\gamma|=d-2, the Hessian H_\gamma=
abla^2(\partial^\gamma g) has at most one positive eigenvalue, while \operatorname{supp}(g) is M-convex. Here \partial^\gamma=\prod_i\partial_{x_i}^{\gamma_i}, and |\gamma|=\sum_i\gamma_i. Penalizing the second-largest positive eigenvalue transfers this concrete signature condition into a differentiable neural regularizer.
Implementation notes
Integrate this at the output head of a small mixture-count or routing model. Let the network emit logits u_alpha for a fixed finite set A of multi-indices alpha with |alpha|=d; define nonnegative coefficients c_alpha=softplus(u_alpha)+epsilon and normalized polynomial coefficients g_alpha=c_alpha/alpha!. For each minibatch, sample K multi-indices gamma with |gamma|=d-2. Construct the degree-two polynomial q_gamma(x)=partial^gamma g(x) analytically from the coefficient table, then form its n-by-n Hessian H_gamma at x=1, with H_gamma[i,j]=partial_{x_i}partial_{x_j}q_gamma(1). Compute the two largest eigenvalues using a symmetric eigensolver for n<=32, or estimate lambda_2 with two-step Lanczos for larger n, and add max(0,lambda_2)^2 to the task loss. For the support term, use a soft support mask m_alpha=sigmoid((c_alpha-tau)/temperature). Sample pairs alpha,beta and coordinates i,j satisfying alpha_i>beta_i and alpha_j<beta_j; penalize missing exchange points with m_alpha*m_beta*(1-m_{alpha-e_i+e_j}) and the analogous term for beta-e_j+e_i. The normalization and Hessian are computed exactly from the coefficient table; only gamma sampling, eigenvalue estimation, and thresholding are approximate. Start with a 4-expert, degree-4 router on CIFAR-10 or a tiny language-model MoE on WikiText-2, comparing ordinary softmax routing against the regularized router at equal parameters and FLOPs. Success means lower validation loss, fewer expert-collapse events, and smoother per-expert load trajectories without reducing routing sparsity; also test robustness under noisy or shifted inputs.
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.