Derangetropy Rank Warp
Implementation & benchmark of arXiv:2607.24705 — Derangetropy Operators
Source paper: Derangetropy Operators arXiv:2607.24705 ⓘ · analyzed Aug 31, 2026
AI-generated research hypothesis, automatically tested. Not peer-reviewed.
Idea description
Insert a distribution-free rank warp before selected MLP or attention projections. For each scalar activation, replace its empirical rank u by the cumulative interval map induced by the Type-III derangetropy kernel, optionally followed by Gaussian or affine output calibration. The transform is invariant to strictly increasing reparameterizations of the feature and contracts the marginal toward central ranks, potentially reducing sensitivity to heavy tails and outliers.
Formulas
Mathematical statement
For an absolutely continuous scalar random variable X with density f and CDF F(x)=Pr[X<=x], the paper's Type-III operator is rho_III[f](x)=2 sin^2(pi F(x)) f(x). The rank-only kernel is w(u)=2 sin^2(pi u) for u in [0,1]. It is a valid density reweighting because the integral from 0 to 1 of w(u) du equals 1. If Y has density rho_III[f], then its CDF is the interval map G(u)=integral from 0 to u of 2 sin^2(pi v) dv = u - sin(2 pi u)/(2 pi), so the transformed rank is G(F(x)). The construction depends only on rank: for every strictly increasing coordinate change h, F_(h(X))(h(x))=F_X(x). In a neural layer, use r=softCDF(a) and r'=G(r); either output r' directly or calibrate it with y=mu+sigma Phi^(-1)(clip(r',epsilon,1-epsilon)), where a is an activation, Phi^(-1) is the standard-normal quantile, and mu and sigma are learned or running affine parameters.
Implementation notes
Integrate the method as a per-channel scalar normalization layer after an MLP activation or before a transformer projection; initially test it in a small MLP rather than inside attention. For a batch tensor A in R^(B times d), maintain one running CDF approximation per feature using K=128 fixed quantile knots q[j,k], where k ranges from 0 to K-1. Estimate each activation rank by piecewise-linear interpolation between adjacent running quantiles. Update the quantile estimates with an exponential moving average during training, but detach the knot values in the forward pass to avoid unstable gradients through sorting. Compute r_prime = r - sin(2 pi r)/(2 pi), clip it to [1e-4,1-1e-4], and either use r_prime directly or map it through the inverse standard-normal CDF followed by learned scale and bias. A differentiable batch-only variant can use r_i = (1/B) sum_j sigmoid((a_i-a_j)/tau), with tau annealed from 0.2 to 0.02, but this costs O(B squared times d), so the running-quantile version should be the main implementation. The kernel and antiderivative are computed exactly; CDFs and quantiles are estimated empirically. First test on CIFAR-10 with a three-layer width-512 MLP and on a small transformer language model, comparing against LayerNorm, RMSNorm, BatchNorm, and identity at equal parameter count. Record validation accuracy or perplexity, activation kurtosis, gradient-norm variance, and robustness to monotone rescaling and heavy-tailed perturbations. Success means lower gradient variance and equal-or-better clean validation performance, with improved accuracy under feature rescaling or outlier corruption and no more than 10 percent training-time overhead.
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.