Research
LipNeXt: Scaling up Lipschitz-based Certified Robustness to Billion-parameter Models
Overview Research area: Certified adversarial robustness, specifically Lipschitz-based deterministic certification for image classifiers, combined with scalable architecture design. Technical level: A

- arXiv
- 2601.18513
- Published
- 2026-01-26
- Authors
- Kai Hu, Haoqi Hu, Matt Fredrikson
AI summary
Overview
Research area: Certified adversarial robustness, specifically Lipschitz-based deterministic certification for image classifiers, combined with scalable architecture design.
Technical level: Advanced. The paper assumes familiarity with Lipschitz constants, orthogonal/Stiefel manifolds, Riemannian gradient descent, and certified robust accuracy metrics.
Scope: The paper introduces a 1-Lipschitz architecture (LipNeXt) that removes both the weight-constraint re-parameterization and the convolutional layers used by prior Lipschitz-certified networks, and shows it scales to 1–2 billion parameters on ImageNet.
What This Paper Is About
Lipschitz-based certification gives deterministic robustness guarantees (unlike probabilistic randomized smoothing), but it has been stuck with small VGG-style models and ≤32M parameters, and it degrades badly on ImageNet. The bottleneck is that building tight 1-Lipschitz networks requires orthogonal weights, and every existing way of learning them (matrix exponential, Cayley transform, LOT-Orth, Cholesky-Orth, AOL/CPL/SLL layers, power-iteration regularization) either adds large computational overhead or fails to scale. The goal is to make Lipschitz certification scale like modern vision backbones do, without giving up determinism or efficiency.
Key Contributions
- A constraint-free manifold optimization scheme for learning orthogonal parameters at scale, with a norm-adaptive Taylor approximation of the matrix exponential (
FastExp) and two stabilization mechanisms (periodic SVD polar retraction and a manifold-adapted Lookahead optimizer). Per-update overhead is reduced to at most five matrix multiplications. - A convolution-free Spatial Shift Module, motivated by a theorem stating that a depthwise circular-padded convolution is norm-preserving (tight 1-Lipschitz isometric) if and only if its kernel contains exactly one non-zero element with value ±1.
- The LipNeXt architecture — the first constraint-free and convolution-free 1-Lipschitz certified model — built from orthogonal projections, spatial shifts, a 1-Lipschitz β-Abs nonlinearity, and L2 spatial pooling.
- State-of-the-art certified robust accuracy (CRA) and clean accuracy on CIFAR-10, CIFAR-100, Tiny-ImageNet, and ImageNet, including the first demonstration that Lipschitz certification scales to billion-parameter models.
Main Findings
- CIFAR-10 (no synthetic data): LipNeXt L32W1024 (64M params) reaches 81.5% clean accuracy and CRA of 71.2 / 59.2 / 45.9 at ε = 36/255, 72/255, 108/255. LipNeXt L32W2048 (256M) reaches 85.0% clean and 73.2 / 58.8 / 43.3. The only metric where a baseline wins is CRA at 108/255 on CIFAR-10, where AOL Large (136M) reports 49.0 versus 45.9 for L32W1024 — and the paper notes AOL pays for this with a substantial drop in clean accuracy and in CRA at other radii.
- CIFAR-100 (no synthetic data): LipNeXt L32W1024 (64M) reaches 53.3% clean and 41.3 / 30.5 / 21.8; L32W2048 (256M) reaches 57.4% clean and 44.1 / 31.9 / 22.2, compared with 54.3% clean for BRONet (68M) and 53.0% for LiResNet (83M).
- Tiny-ImageNet: LipNeXt L32W1024 (64M) reaches 42.5% clean and 32.0 / 21.8 / 15.2; L32W2048 (256M) reaches 45.5% clean and 35.0 / 25.9 / 18.0, versus 41.2% clean for BRONet (75M) and 32.1% clean for SLL X-Large (1.1B).
- Synthetic data helps at larger size: With diffusion-generated data, LipNeXt L32W2048 (256M) reaches 88.2% clean and 79.2 / 68.0 / 54.9 on CIFAR-10 and 62.1% clean and 51.2 / 38.5 / 27.5 on CIFAR-100; L32W2896 (512M) reaches 92.7% clean and 81.7 / 68.6 / 55.8 on CIFAR-10 and 63.6% clean and 55.2 / 39.2 / 28.3 on CIFAR-100.
- Robust overfitting: Scaling LipNeXt without extra data improves clean accuracy but decreases CRA, which the authors attribute to the small dataset size and address with diffusion-generated data.
- ImageNet scale-up: LipNeXt L32W4096 (1B parameters) reaches 40.2% clean / 21.1% CRA at ε = 1 and 55.9% clean / 40.3% CRA at ε = 36/255. L32W5792 (2B) reaches 41.0% clean / 22.4% CRA at ε = 1 and 57.0% clean / 41.2% CRA at ε = 36/255. The paper reports gains over prior Lipschitz models of up to +8% CRA at ε = 1 and 3% at ε = 36/255.
- Training stability and low precision: LipNeXt trains in bfloat16 using only orthogonal matrix multiplications and spatial shifts. LiResNet requires float32 on ImageNet because bfloat16 causes occasional numerical overflow, and BRONet's FFT-based complex convolutions are limited to complex32, incurring float64-like memory and compute overhead.
- Throughput: On a single 8×H100 machine, LipNeXt L32W4096 (1B) trains at 8.9 minutes per epoch and L32W5792 (2B) at 17.8 minutes per epoch (L32W5792 was trained on two such machines and the speed is normalized to equivalent single-node throughput). LiResNet (51M) reports 5.3 minutes per epoch and BRONet (86M) reports 10.5 minutes per epoch on the same hardware description.
- Scaling behavior (ImageNet, 400 randomly sampled classes, ε = 1): Fixing depth at 32, widening from 1024 to 4096 raises clean accuracy from 40.5% to 51.7% and CRA from 22.9% to 30.0%. Fixing width at 2048, deepening from 8 to 128 layers raises clean accuracy from 30.7% to 47.5% and CRA from 22.4% to 26.8% (CRA peaks at depth 64 with 26.9%). Under a fixed 1B parameter budget, the (32, 4096) configuration is best at 51.7% clean / 30.0% CRA.
- Padding finding: Empirically, zero-padding shifts outperform circular shifts in the convolution-free architecture, but the authors hypothesize this is because zero-padding injects position information. Adding explicit positional encoding while keeping circular padding gives superior performance, consistent with the norm-preservation guarantee — a hypothesis they state their experiments confirm.
- Datasets used are described only by name and class-sampling scheme; the number of training or test images is not reported in the provided content.
Methodology in Plain English
The core problem is that a network's certified radius depends on how tightly you can bound its Lipschitz constant, and that bound is only tight when weight matrices are orthogonal. Rather than parameterizing weights so they happen to be orthogonal (the "constrained" approach used by prior work), the authors let the weights live directly on the orthogonal manifold and constrain the updates instead. They compute a projected (Riemannian) gradient, then move along the manifold using a matrix exponential of a skew-symmetric update.
The expensive part is the matrix exponential. Their observation is that in large-model training, learning rates are small (on the order of 10⁻³), so the update matrix has a small Frobenius norm and the exponential can be replaced by a short Taylor series. The truncation is chosen by norm: second order below 0.05, third order below 0.25, fourth order below 1, and the exact exponential above 1. Truncation error is cleaned up by an SVD-based polar retraction once per epoch, and by a Lookahead optimizer adapted to the manifold that interpolates the skew-symmetric updates in the tangent space (at a 0.5 mixing factor) rather than interpolating the orthogonal matrices themselves, which would break orthogonality.
For spatial mixing, they avoid convolution entirely. Their theorem shows that under circular padding the only depthwise convolution that preserves norms exactly is a shifted identity with a single ±1 kernel entry, so they simply implement that shift: split each token's channels into partitions and circularly shift some partitions forward and others backward by one position (five partitions for 2D data, with shift ratios of 1/8 or 1/16). This mixes neighboring positions at essentially zero cost. A LipNeXt block wraps this by adding a learnable positional embedding, projecting channel features with an orthogonal matrix R, applying the shift, projecting back with Rᵀ, applying a second orthogonal matrix M and a 1-Lipschitz β-Abs activation (β-Abs generalizes MinMax, which is recovered at β = 0.5). The overall network follows LiResNet's macro-structure but swaps in LipNeXt blocks and replaces the Neck with L2 spatial pooling so the whole network stays 1-Lipschitz. Training uses the EMMA loss and the LiResNet++ recipe, without diffusion-generated data by default.
Why This Matters
Impact on research: This work argues that deterministic Lipschitz certification is not inherently unscalable — the previous ceiling came from re-parameterization overhead and second-best architecture choices, not from theory. It gives the field a viable path to billion-parameter provably robust models and shows they can be trained on low-precision hardware, which is the first time Lipschitz methods have plausibly tracked the same scaling curve as mainstream vision models. It also contrasts explicitly with randomized smoothing, which the appendix criticizes for being probabilistic (false positives, typically below 0.1% but still nonzero) and computationally expensive from noise averaging.
Real-world applications (as motivating contexts named in the paper):
- Autonomous driving, where certified behavior under input perturbation is safety-critical.
- Medical image processing, where a misclassified or manipulable input can have severe clinical consequences.
- Malware classification, where adversaries can deliberately perturb inputs.
- Any deployment that requires an absolute, mathematical robustness guarantee rather than an empirical defense that may fail against future attack strategies.
Industry relevance: The constraint-free optimizer avoids power-iteration penalties and the numerical fragility of prior methods, and the architecture is stable in bfloat16 while still achieving throughput on par with much smaller prior models. That means LipNeXt can keep benefiting from hardware advances (low-precision accelerators) that LiResNet and BRONet currently cannot fully exploit.
Future Directions
- Training LipNeXt on large-scale image-text pair datasets is explicitly stated as future work, which would test whether the scaling trend continues beyond the 2B parameters reported here.
- The CRA loss at larger model sizes without synthetic data (robust overfitting) is addressed only by adding diffusion-generated images; better remedies for small-data certified scaling remain open.
- Whether the 1B-parameter-depth-width optimum (depth 32) holds at larger budgets or on full ImageNet-1K rather than the 400-class subset used for the scaling study.
- Transferring the approach to the safety-critical domains used as motivation (autonomous driving, medical imaging, malware classification), which the paper does not evaluate.
Target Audience
Researchers and engineers working on adversarial robustness certification, provable deep learning, and Lipschitz-constrained or norm-preserving network design. It is most useful to readers already comfortable with manifold optimization and certified accuracy metrics, and to practitioners evaluating whether deterministic certification can realistically be applied to billion-parameter image models.
Authors’ abstract
Lipschitz-based certification offers efficient, deterministic robustness guarantees but has struggled to scale in model size, training efficiency, and ImageNet performance. We introduce \emph{LipNeXt}, the first \emph{constraint-free} and \emph{convolution-free} 1-Lipschitz architecture for certified robustness. LipNeXt is built using two techniques: (1) a manifold optimization procedure that updates parameters directly on the orthogonal manifold and (2) a \emph{Spatial Shift Module} to model spatial pattern without convolutions. The full network uses orthogonal projections, spatial shifts, a simple 1-Lipschitz $β$-Abs nonlinearity, and $L_2$ spatial pooling to maintain tight Lipschitz control while enabling expressive feature mixing. Across CIFAR-10/100 and Tiny-ImageNet, LipNeXt achieves state-of-the-art clean and certified robust accuracy (CRA), and on ImageNet it scales to 1-2B large models, improving CRA over prior Lipschitz models (e.g., up to $+8\%$ at $\varepsilon{=}1$) while retaining efficient, stable low-precision training. These results demonstrate that Lipschitz-based certification can benefit from modern scaling trends without sacrificing determinism or efficiency.