{
 "artifacts": null,
 "category": "architecture",
 "description": "Replace an Erlang delay or exponential smoothing cascade in a recurrent or state-space layer by a positive rational kernel of the form \\(\\kappa(u)=C e^{-a u}p(u)^2\\). Choose the degree-\\(m\\) polynomial by deleting the adjacent pair of Laguerre zeros with smallest relative gap from \\(L_{m+2}\\), then rescale the resulting density to unit mean. This preserves a nonnegative impulse response while reducing temporal jitter relative to Erlang filters.",
 "formulas_latex": [
  "$$\\kappa_N^{\\mathrm{Er}}(u)=\\frac{N^{N}}{(N-1)!}u^{N-1}e^{-Nu},\\qquad K_N^{\\mathrm{Er}}(s)=\\left(\\frac{N}{N+s}\\right)^N,\\qquad \\operatorname{Var}(U)=\\frac{1}{N}.$$",
  "$$\\kappa(u)=C e^{-a u}p(u)^2,\\qquad \\int_0^\\infty \\kappa(u)\\,du=1,\\qquad \\int_0^\\infty u\\kappa(u)\\,du=1.$$",
  "$$r_1\u003c\\cdots\u003cr_{m+2},\\qquad g_j=\\frac{r_{j+1}-r_j}{r_j},\\qquad j^*=\\arg\\min_{1\\le j\\le m+1}g_j,\\qquad q(x)=\\prod_{i\\notin\\{j^*,j^*+1\\}}(x-r_i).$$",
  "$$\\operatorname{Var}_{\\mathrm{opt}}(m)=\\min_{1\\le j\\le m+1}\\frac{r_{j+1}-r_j}{r_j},\\qquad \\operatorname{Var}(U)=\\int_0^\\infty (u-1)^2\\kappa(u)\\,du.$$"
 ],
 "id": 2809,
 "implementation": "Integrate this at the impulse-response level of a causal SSM or recurrent smoothing layer, rather than in the token-mixing projection. For sequence step size \\(\\Delta t\\), choose polynomial degree \\(m\\), compute the roots of the generalized Laguerre polynomial \\(L_{m+2}\\) using a symmetric tridiagonal Jacobi-matrix eigensolver, and calculate all relative gaps \\(g_j=(r_{j+1}-r_j)/r_j\\). Delete the pair attaining the minimum gap and construct \\(q(x)=\\prod_{i\\ne j^*,j^*+1}(x-r_i)\\). Numerically find \\(a\\) by bisection so that the normalized function \\(e^{-au}q(u)^2\\) has mean one; obtain \\(C\\) from the unit-integral constraint. Sample the kernel at \\(u_k=k\\Delta t\\), set \\(w_k=e^{-a u_k}q(u_k)^2\\), truncate at a fixed horizon \\(K\\), and renormalize with \\(w_k\\leftarrow w_k/\\sum_{k=0}^K w_k\\). The forward operation is `y[t] = sum_k w[k] * x[t-k]`; for an SSM implementation, fit a stable finite-dimensional realization to this impulse response or retain it as a depthwise convolution. The paper provides the positivity guarantee and Laguerre root-selection rule; quadrature, discretization, truncation, and realization are estimated numerically. Begin with degrees 4, 8, and 16 on Copying Memory, Adding Problem, and sequential MNIST. Compare against equal-state Erlang cascades and learned unconstrained FIR filters at matched receptive field, parameter count, and FLOPs. Record impulse-response variance, validation loss, gradient norm, and recurrent rollout stability. The key falsifiable signal is lower validation loss or longer effective memory at equal state budget, together with variance decreasing faster with degree than the Erlang baseline and no negative or oscillatory impulse weights.",
 "math_summary": "The paper's benchmark is the mean-one Erlang kernel \\(\\kappa_N^{\\mathrm{Er}}(u)=\\frac{N^N}{(N-1)!}u^{N-1}e^{-Nu}\\), with Laplace transform \\(K_N^{\\mathrm{Er}}(s)=\\left(\\frac{N}{N+s}\\right)^N\\) and variance \\(1/N\\). The proposed class is a nonnegative kernel \\(\\kappa(u)=C e^{-a u}p(u)^2\\), where \\(u\\ge 0\\), \\(p\\) is a real polynomial of degree \\(m\\), \\(a\u003e0\\), and \\(C\\) normalizes the integral to one. Positivity follows from the square. After imposing unit mean, the paper states that the minimum variance equals the smallest relative gap between adjacent zeros of the Laguerre polynomial \\(L_{m+2}\\), and every optimizer is obtained by deleting the corresponding pair of zeros. Let \\(r_1\u003c\\cdots\u003cr_{m+2}\\) be the positive zeros of \\(L_{m+2}\\), define \\(g_j=(r_{j+1}-r_j)/r_j\\), and choose \\(j^*=\\arg\\min_j g_j\\). Form \\(q(x)=\\prod_{i\\notin\\{j^*,j^*+1\\}}(x-r_i)\\), then determine the exponential scale \\(a\\) and normalization \\(C\\) so that the integral and mean equal one. The resulting variance is \\(\\operatorname{Var}(U)=\\int_0^\\infty (u-1)^2\\kappa(u)\\,du\\).",
 "math_tags": [
  "approximation-theory",
  "spectral-theory",
  "linear-algebra",
  "probability"
 ],
 "ml_areas": [
  "ssm",
  "rnn",
  "training-dynamics"
 ],
 "paper": {
  "arxiv_id": "2608.29143",
  "arxiv_url": "https://arxiv.org/abs/2608.29143",
  "summary_what_math_gives_to_ml": "The paper gives a constructive way to build nonnegative rational temporal kernels that approximate a deterministic delay with substantially lower variance than Erlang/gamma cascades at the same polynomial order. The transferable asset is the polynomial-square form, which guarantees positivity while remaining implementable as a finite-dimensional repeated-pole state-space filter, together with the Laguerre-zero characterization that selects the least-variable kernel. This suggests replacing Erlang delay lines or positive exponential smoothing cascades in state-space and recurrent modules with a spectrally designed polynomial-square kernel. The first test should compare equal-state causal temporal filters on long-range sequence prediction, measuring delay fidelity, validation loss, and stability under recurrent rollout.",
  "title": "Least Variability in a Polynomial-Square Class of Rational Kernels",
  "year": "2026"
 },
 "ratings": {
  "difficulty": 6,
  "novelty": 7,
  "usefulness": 6
 },
 "solves": [
  "accuracy",
  "stability",
  "scalability"
 ],
 "title": "Laguerre-Optimal Positive Delay Filter",
 "url": "https://synthcore.org/idea/2809/laguerre-optimal-positive-delay-filter",
 "verification": {
  "peer_reviewed": false,
  "status": "unverified",
  "status_label": "Unverified",
  "verdict_source": "deterministic test code (paired-seed permutation statistics)"
 }
}