Unverified 2026

Noise-Induced Compact-Support Neural Field

Implementation & benchmark of arXiv:2609.01984 — Instantaneous shrinking of supports for stochastic PDEs

Usefulness6/10
Difficulty6/10
Novelty8/10

Source paper: Instantaneous shrinking of supports for stochastic PDEs arXiv:2609.01984 · analyzed Sep 3, 2026

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

Idea description

Replace ordinary dropout or soft sparsity penalties on a nonnegative spatial or token activation field with sublinear multiplicative stochastic dynamics. The field receives local diffusion or graph smoothing, while noise amplitude u^gamma vanishes at zero but is relatively strong near zero; this creates an absorbing zero state and may produce exact contiguous inactive regions. The module is suitable for feature maps, graph-node fields, token routing scores, or continuous neural operators.

Formulas

$$\partial_t u(t,x)=Lu(t,x)+\sigma(u(t,x))\xi(t,x),\qquad u(0,\cdot)=u_0,$$
$$Lu(t,x)=a(t,x)\partial_x^2u(t,x)+b(t,x)\partial_xu(t,x)+c(t,x)u(t,x),$$
$$\sigma(u)=u^\gamma\mathbf{1}_{\{u\geq0\}},\qquad 0<\gamma<1,$$
$$u_i^{k+1}=\Pi_{\geq0}\left[u_i^k+\Delta t\left(\kappa(L_Gu^k)_i-\lambda u_i^k\right)+\rho\,(u_i^k+\varepsilon)^\gamma\sqrt{\Delta t}\,\eta_i^k\right],\qquad \eta_i^k\sim\mathcal N(0,1).$$

Mathematical statement

The paper studies the nonnegative stochastic PDE partial_t u(t,x) = L u(t,x) + sigma(u(t,x)) xi(t,x), where u is a scalar field, L = a(t,x) partial_x^2 + b(t,x) partial_x + c(t,x) is a parabolic second-order operator, and xi is space-time white noise. The model coefficient is sigma(u) = u^gamma 1_{u >= 0} for 0 < gamma < 1, with sigma equal to zero at the absorbing boundary. The key boundary property is sigma(u) / u = u^(gamma-1), which diverges as u approaches zero from above, so the stochastic term is large relative to the linear signal scale near zero. The abstract states that sufficiently light-tailed initial data can yield compact support at every positive time even when the initial support is noncompact; finite initial mass suffices for 0 < gamma <= 1/2, while gamma > 1/2 requires a polynomial moment condition whose order grows as gamma approaches 1. In the neural adaptation, u_i is a nonnegative activation on a grid or graph, L_G is a graph Laplacian, eta_i^k is standard Gaussian noise, gamma controls sublinearity, and the Euler-Maruyama discretization is the implementable approximation. Exact compact-support guarantees are not automatic for a finite discretization, so support creation must be tested empirically.

Implementation notes

Integrate the module after a CNN, graph neural network, or transformer routing block whose activations are interpreted as a nonnegative field u in R_+^n. Use a 2D grid Laplacian for image feature maps, a fixed chain or token-neighbor graph for sequences, or the task graph for graph data. First create nonnegative initial values with u^0 = softplus(z), where z is the original learned activation or routing logit. Run K stochastic relaxation steps during training. Let L_G be the symmetric graph Laplacian, kappa be the diffusion strength, lambda be optional decay, rho be the noise scale, gamma be in (0,1), epsilon be a small numerical floor, and Delta t be the step size. At each step sample independent standard Gaussian eta_i^k, calculate the displayed Euler-Maruyama update, and project negative values to zero. To obtain actual computational sparsity, form a mask m_i = 1[u_i^K > tau] after relaxation and either use it in the next layer or use a straight-through estimator during training. Add a mild support cost alpha times the number of active sites, but do not use such a large alpha that the stochastic mechanism is replaced by ordinary L0 regularization. Backpropagate through sampled noise using reparameterization. Measure exact-zero fraction, active support size, connected-component count, gradient variance, and sensitivity to gamma, rho, and initial mass sum_i u_i^0. Start with a small CIFAR-10 CNN and compare against Bernoulli dropout, deterministic graph diffusion, hard-concrete gates, and L1 regularization at matched accuracy. Sweep gamma in {0.25, 0.5, 0.75}. A positive result is at least twice the exact-zero fraction of dropout at equal validation accuracy, or equal accuracy with substantially fewer active sites and lower measured inference FLOPs. A negative result is instability, accuracy collapse, or merely producing small nonzero values without a useful reduction in active computation.

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.