Source-linked AI summary

Measuring Model Robustness via Fisher Information: Spectral Bounds, Theoretical Guarantees, and Practical Algorithms

Chong Zhang, Xiang Li, Jia Wang, Qiufeng Wang, Xiaobo Jin

arXiv:2606.04767v1cs.LGcs.CV

TL;DR

Robustness evaluation for deep neural networks remains costly, attack-dependent, and difficult to interpret, especially when connecting geometric sensitivity with uncertainty. The paper proposes an FIM spectral-norm framework with theoretical architectural bounds and scalable estimation algorithms. Experiments across datasets and architectures show correlation with adversarial vulnerability, while the authors position the metric as a complementary diagnostic tool rather than a replacement for benchmark attacks.

  • Problem

    Existing robustness evaluations are often attack-dependent, computationally expensive, or insufficiently interpretable and scalable, lacking a unified connection between geometric sensitivity and statistical uncertainty.

  • Method

    The paper uses the maximum eigenvalue of the input-space FIM, or its inverse, as a robustness metric, establishes its connection to input-Jacobian variance, derives architectural spectral bounds, and develops scalable estimation algorithms.

  • Results

    Experiments across CIFAR-10, ImageNet, and medical images show correlation between the metric and adversarial vulnerability across datasets and architectures.

  • Takeaways & Limitations

    The framework provides an interpretable diagnostic tool for robustness-aware model selection and complements attack-based evaluations with architectural sensitivity analysis.

  • Takeaways & Limitations

    The theoretical architectural ranking may differ from dataset-specific rankings because of weights, training dynamics, data-dependent curvature, and spectral-estimation approximation errors.

Abstract

from arXiv · show

The robustness of deep neural networks is crucial for safety-critical deployments, yet existing evaluation methods are often attack-dependent and lack interpretability. We propose a principled, attack-agnostic robustness metric based on the spectral norm of the Fisher Information Matrix (FIM), which quantifies the worst-case sensitivity of the model's output distribution to input perturbations. Theoretically, we establish that the FIM equals the variance of the input Jacobian and derive closed-form spectral bounds for common architectures, including VGG, ResNet, DenseNet, and Transformer, providing the first theoretical robustness ranking. To enable scalable evaluation, we develop efficient algorithms, including power iteration and Hutchinson-based estimation, that support both white-box and black-box settings. Extensive experiments across multiple datasets, including CIFAR, ImageNet, and medical images, and across multiple architectures show a strong correlation between our metric and adversarial vulnerability. Our framework serves as an interpretable diagnostic tool that complements attack-based evaluations, offering insights into architectural sensitivity and guiding the design of more robust models. Code is available at: https://github.com/franz-chang/SRP/.

1. Introduction

The paper addresses limitations of attack-dependent and difficult-to-interpret robustness evaluation by proposing an attack-agnostic Fisher Information Matrix metric. It connects probabilistic sensitivity to architecture-level theory and scalable evaluation, with experiments supporting correlation with adversarial vulnerability.

  • Motivation: Attack-dependent metrics are computationally expensive and sensitive to attack choices, while theoretical bounds often lack probabilistic interpretation and scale poorly to modern architectures.The paper identifies a missing framework connecting geometric sensitivity with statistical uncertainty.
  • Core proposal: The authors establish that the FIM equals the variance of the input Jacobian, linking information-theoretic geometry with practical robustness evaluation.This connection supplies the theoretical basis for interpreting FIM sensitivity in the input space.
  • Core proposal: The paper introduces the first FIM spectral norm-based robustness metric, using the largest eigenvalue or its reciprocal to measure worst-case output-distribution sensitivity.The metric is described as attack-agnostic, interpretable, and probabilistically grounded.
  • Theoretical analysis: Closed-form spectral bounds for ReLU, convolution, and attention support theoretical robustness rankings of VGG, ResNet, DenseNet, and Transformer architectures.The ranking is intended to expose sensitivity differences arising from architectural design choices.
  • Practical algorithms: Power iteration and Hutchinson-based methods enable scalable FIM spectral-norm estimation in white-box and black-box settings.These algorithms are designed to extend evaluation to large-scale models.
  • Empirical validation: Experiments across CIFAR-10, ImageNet, and medical images verify correlation between the metric and adversarial vulnerability and demonstrate usefulness for robustness-aware model selection.The framework is positioned as a diagnostic complement to attack-based benchmarks.

2. Related Work

Prior robustness research relies on attacks, sensitivity bounds, or information-theoretic quantities that each leave important gaps in cost, scalability, probabilistic interpretation, or input-space geometry. This paper positions FIM spectral analysis as a bridge between these perspectives and adapts randomized algorithms for scalable computation.

  • Adversarial attack-based metrics: Adversarial attacks such as PGD and C&W expose vulnerabilities but are computationally expensive and attack-dependent, limiting generalization to unknown threat models or real-world noise.These methods create perturbations intended to induce misclassification.
  • Sensitivity-based metrics: Lipschitz and Jacobian-norm methods bound model sensitivity but lack probabilistic interpretation and are difficult to scale to complex architectures such as Transformers.The paper attributes these difficulties to fuzzy boundaries or exponential computational complexity.
  • Information-theoretic perspective: Earlier information-theoretic robustness metrics did not connect KL divergence and mutual information to input-space geometry.The paper addresses this gap by linking KL divergence to Fisher information.
  • Fisher information foundations: FIM has traditionally supported statistical estimation, natural-gradient optimization, and uncertainty quantification, primarily through parameter-space rather than input-space analysis.The paper distinguishes its robustness focus from these established uses.
  • FIM for adversarial robustness: Prior FIM studies examined adversarial detection, robust training, or manifold geometry without establishing FIM eigenvalues as direct measures of inherent model robustness.The paper identifies the largest FIM eigenvalue as encoding worst-case sensitivity.
  • Spectral methods: Spectral normalization studies mainly target generative models, whereas this work analyzes discriminative CNN and Transformer architectures through FIM spectral properties.The comparison highlights a difference in both model scope and analytical perspective.
  • Randomized algorithms: The authors adapt Hutchinson estimation and power iteration to FIM structure, providing efficient λmax(F) estimation with provable convergence for modern architectures.These randomized methods address the computational scale of large matrix operations.

3. Methodology

The method links output-distribution geometry and predictive uncertainty through the Fisher Information Matrix, using its spectral norm to quantify worst-case input sensitivity. It derives architectural bounds and scalable estimators, then compares theoretical trends with empirical robustness while noting dataset- and model-dependent deviations.

  • Problem Formulation: The framework models input perturbations through KL divergence between output distributions, approximated by a Mahalanobis distance governed by the Fisher Information Matrix.The perturbation direction is approximately the principal eigenvector of the FIM.
  • Problem Formulation: The FIM equals the variance of the input Jacobian and weights gradient directions by the model’s predictive probabilities.This connects geometric sensitivity with statistical confidence rather than treating the Jacobian as purely deterministic.
  • Robustness Metric: The largest FIM eigenvalue, or spectral norm, bounds worst-case perturbation impact and defines an attack-agnostic robustness measure.The principal eigenvector identifies the perturbation direction that maximizes the local KL divergence.
  • Theoretical Analysis: Theoretical bounds reduce deep-network analysis to component Jacobian norms, covering ReLU, pooling, normalization, softmax, convolution, attention, and common architectures.ReLU and max pooling have spectral norm 1, while average pooling can smooth perturbations and normalization layers can amplify or suppress them.
  • Empirical Evaluation: The theoretical architectural trend is largely consistent with empirical rankings, with DenseNet121 least robust and ViT-B-16 among the most robust.Actual rankings can differ because of weights, training dynamics, data-dependent curvature, and spectral-norm estimation errors.
  • Practical Algorithms: Power iteration and Hutchinson-based procedures estimate the FIM spectral norm efficiently by exploiting its special structure.Indirect estimation through P can reduce space complexity and approximation error, while convergence probability improves as the number of samples increases.

4. Experiments

Experiments evaluate whether the FIM-based metrics capture relative robustness trends across architectures, datasets, attack settings, and access regimes. Results support correlation with attack-based measures while highlighting diagnostic value and practical trade-offs.

  • Scope boundaries: The study does not aim to produce absolute rankings under the strongest attacks or reproduce a full RobustBench leaderboard, partly because large-scale estimation is costly.
  • Experimental scope: The evaluation spans CNNs and Transformers, datasets from MNIST to ImageNet, and medical images to reduce dependence on one architecture or domain.
  • Experimental objective: The experiments target reliable relative robustness trends and interpretable insights complementary to attack-based evaluation.
  • Metric validation: The FIM spectral norm tracks CW-based reductions closely, while PGD changes less because training and testing use different attacks.
  • Metric stability: FIM spectral norm, L(x), and CLEVER remain stable as dataset size changes, whereas CW and PGD fluctuate more because their accuracy functions are nondifferentiable.
  • Model ranking: On CIFAR10, FIM spectral norm matches the rankings of L(x) and CLEVER, while Rspec matches CW, supporting their use as complementary ranking metrics.
  • Access regimes: Black-box and white-box evaluations yield consistent dataset comparisons, although black-box estimated eigenvalues are an order of magnitude lower.
  • Practical trade-offs: The metric supports gradient-inaccessible auditing and model selection, but costs more computation than Lipschitz or CLEVER estimates because it provides richer curvature and uncertainty information.

5. Conclusion

The paper presents Fisher information as a unified framework for quantifying neural-network robustness. It combines theoretical spectral analysis with scalable white-box and black-box computation, and experiments validate the framework across datasets and architectures.

  • The framework uses the maximum FIM eigenvalue, or its inverse, as an interpretable robustness metric grounded in information theory.
  • The analysis derives spectral bounds for components and compares theoretical robustness across VGG, ResNet, DenseNet, and Transformer architectures.
  • Three spectral-norm algorithms, including Hutchinson estimation with finite differences, support scalable evaluation across white-box and black-box settings.
  • Experiments across datasets of varying sizes and types validate the theoretical results and support the metric’s attack-independent applicability.

6. Impact Statement

The framework is positioned as a tool for safer, more interpretable robustness evaluation, architecture comparison, and model selection. Its scalable algorithms are intended to reduce evaluation barriers, while the impact claims remain connected to attack-independent diagnostics.

  • Safety-critical applications: The framework targets safety-critical applications by quantifying robustness without relying on adversarial attacks.The stated application areas include autonomous systems, healthcare, and finance.
  • Transparency and interpretability: Its Fisher-information, Jacobian-variance, and robustness connections provide interpretable insights into model behavior and vulnerabilities.The paper links this interpretability to explainable AI and regulated deployment contexts.
  • Model selection and benchmarking: The spectral norm ∥F(x)∥2 or its reciprocal supports architecture comparison and model selection with less reliance on empirical adversarial testing.The paper specifically frames comparisons such as VGG versus Transformer as a model-selection use case.
  • Efficiency of robustness evaluation: Power iteration and Hutchinson approximation enable scalable evaluation of large models and can reduce the computational barrier relative to attack-based evaluation.The stated benefit is improved accessibility for resource-constrained researchers and practitioners.

A.4. Summary on the relationship between the three metrics

The paper relates the FIM spectral norm to Lipschitz, CLEVER, and randomized-smoothing measures through local sensitivity and probability-distribution geometry. These relationships connect the proposed metric to both gradient-based robustness estimates and certification radius bounds.

  • Metric relationships: All three metric families are related to gradient norms used to measure local model sensitivity and stability.The comparison covers the FIM spectral norm, the model Lipschitz constant, and the CLEVER score.
  • Metric relationships: The FIM spectral norm measures model-output sensitivity, whereas CLEVER estimates the Lipschitz constant of the cross-entropy loss.The former concerns the model f(x), while the latter concerns the loss function.
  • Randomized smoothing: The FIM spectral norm limits an upper bound on the expected certification radius produced by randomized smoothing.The relationship is stated for the certification radius r under the smoothing framework.
  • Derivation: The theoretical connection follows by expanding log p(y|x′) to second order and substituting the expansion into the KL divergence.A third-order remainder is then bounded to characterize the approximation error.
  • Derivation: The approximation bounds can be expressed for both l2 and l∞ perturbations, with dataset-level analysis replacing perturbation norms by an upper bound θ.The paper discusses both norms explicitly and applies the bound across an entire dataset.

C.2. Experimental estimation

The experimental-estimation section derives the FIM from model-output probabilities and input gradients, then evaluates approximation error in sampled model and dataset settings. The reported approximation error and proportionality coefficient are both very small.

  • Experimental setup: The experiments sample 500 CIFAR10 examples from four classic models trained with CW adversarial training, and evaluate ResNet18 across three datasets.The CIFAR10 perturbations satisfy ∥δ∥∞≤θ.
  • Experimental findings: The approximation error and proportionality coefficient are both reported as very small, so the approximation error can be ignored in practice.This conclusion is based on the results in the sampled CIFAR10 models and the ResNet18 multi-dataset evaluation.
  • FIM estimation: For softmax classifiers, the FIM is expressed through the Jacobian of the softmax input and the variance of the Jacobian matrix.The Jacobian is taken with respect to the input, while variance denotes the variance of the matrix random variable.
  • FIM estimation: The stochasticity used in FIM estimation comes from the predictive distribution p(y|x), which probabilistically reweights deterministic Jacobian directions.This produces a covariance-like summary of gradient directions.
  • FIM estimation: The classification FIM is computed from class probabilities and their input gradients, yielding an estimated Fisher information matrix.The construction uses the model’s output probabilities and the gradient of each output with respect to x.

E. General Analysis of Model Robustness

The paper develops spectral bounds for the FIM through output-probability covariance and network Jacobians, then specializes the analysis to common layers and architectures. The resulting bounds expose how operators, depth, and residual structure affect sensitivity.

  • General bounds: The FIM spectral norm is bounded using the output-probability covariance matrix and the spectral norm of the model-output Jacobian.The bound is framed for classifiers of the form h(x)=softmax(f(x)).
  • General bounds: When the principal eigenvector of the probability covariance aligns with the principal left singular vector of the Jacobian, the bound becomes an exact expression.The resulting expression is 2 max_k p_k(1−p_k)∥Jf(x)∥2.
  • Limitations of the bound: The general spectral bound depends on alignment that is usually unlikely, limiting when the exact equality condition can be expected.The paper states that the principal eigenvector of the probability covariance is usually unlikely to be the principal left singular vector of the Jacobian.
  • Convolution layer: For convolution operators, the Jacobian spectral norm is approximated by the spectral norm of the convolution kernel.The analysis covers circular, zero-padded, and valid convolution settings.
  • Pooling and normalization: Pooling and normalization layers receive separate Jacobian analyses based on their window structure and inference-stage parameters.Average pooling yields eigenvalue-based conclusions, while batch normalization can risk gradient explosion or vanishing depending on γ relative to σ.
  • Architecture analysis: VGG’s Jacobian bound can grow or decay exponentially with depth because it lacks residual connections, while ResNet grows more modestly.The VGG discussion attributes the behavior to the product of layer spectral norms.
  • Transformer: The transformer analysis combines previously derived component bounds for classification without a decoding layer.The paper explicitly notes that the decoding layer is not used for transformer classification.

H.2. Alignment of randomly sampled vectors with the principal eigenvector of a matrix

The section analyzes how a random unit vector aligns with a fixed principal eigenvector, showing that typical alignment decreases rapidly as dimension increases. Concentration and variance arguments yield an O(1/√n) scale for the inner product.

  • For uT v, symmetry gives zero expectation because opposite directions produce opposite inner products.
  • The variance of uT v is derived from the uniform distribution of u on the unit sphere.
  • |uT v| is usually no more than O(1/√n), with 95% probability within [−2σ,+2σ].
  • Applying Levy’s lemma to the 1-Lipschitz function f(u)=uT v provides the concentration argument behind the alignment bound.
  • A random unit vector’s alignment with a fixed principal eigenvector decays exponentially as the dimension increases.

I. Overview and Analysis of Algorithms

The paper develops indirect spectral-norm estimators for the FIM that exploit its low-rank structure, reducing computational demands while supporting practical robustness evaluation across models and datasets. Power iteration and Hutchinson estimation are presented alongside comparison metrics and experimental settings.

  • Algorithmic framework: Indirect estimation replaces F(x)=QΛQT with B(x)=Λ^1/2QTQΛ^1/2, exploiting matrix identities and lower-dimensional computation.The approach is used in both power iteration and Hutchinson-based estimation.
  • Algorithmic framework: The proposed algorithms generally reduce time and space complexity, making spectral robustness metrics more feasible for large-scale applications.
  • Power iteration: Power iteration initializes a gradient-based vector and repeatedly applies Q(Λ(QTbt)) until the estimated eigenvalue converges.
  • Hutchinson estimation: Hutchinson estimation repeatedly samples Gaussian vectors and applies the low-rank representation to estimate the principal eigenvalue.
  • Evaluation setup: Experiments evaluate VGG16, ResNet18, DenseNet121, and ViT-B/16 on CIFAR10, MNIST, and Tiny-ImageNet after resizing inputs to 224×224.
  • Evaluation setup: The evaluation compares spectral robustness with Lipschitz, CLEVER, PGD, and C&W metrics using clean test inputs and perturbation samples.

J.3. Analysis of Spectral Robustness Measure

Experiments verify the FIM–Jacobian variance relationship, examine architecture and metric rankings, and evaluate estimator accuracy, runtime, randomness, and cross-dataset behavior. Results generally support the spectral robustness measure while identifying dataset-specific effects.

  • Theoretical verification: As sample size increases, the estimated FIM and Jacobian variance converge, consistent with the law of large numbers.The comparison uses relative Frobenius error between F(x) and the estimated Jacobian variance.
  • Architecture analysis: For models sharing components, robustness Rspec decreases as the number of layers increases, illustrated by the ResNet18–ResNet34 comparison.The reported relationship is approximately O(1/L).
  • Architecture analysis: ViT has the highest reported robustness metric among the architectures discussed in the corresponding comparison.
  • Metric comparison: Rspec and the Lipschitz-based robustness measure show linear correlation and consistent evaluation in the reported Pearson analysis.
  • Algorithm evaluation: When output dimensionality is small, direct eigenvalue decomposition is usually faster; medium and larger output scales favor iterative estimators.The runtime comparison varies with model output scale.
  • Algorithm evaluation: Hutchinson approximation error decreases as the number of samples increases, while Gaussian sampling performs better than Rademacher sampling.
  • Algorithm evaluation: Using Λ^1/2QTQΛ^1/2 makes Hutchinson error depend on the number of categories rather than parameter dimension, unlike QΛQT.
  • Robustness evaluation: Across repeated experiments and datasets, metric variance is generally small except for DenseNet121 on MNIST, and ResNet18 ranks below ViT-B/16 on medical data and CIFAR100.Excluding MNIST outliers, Rspec agrees with L(x), CLEVER, and CW in the reported comparisons.

K. Limitations and Future Work

The paper identifies limitations involving data dependence, computational cost for extremely large models, task scope, and incomplete validation on RobustBench models. It proposes normalization, efficiency, task-extension, and benchmark-validation directions.

  • Because F(x) depends on the input distribution, cross-dataset robustness comparisons reflect both model and data properties.This dependence challenges universal ranking and motivates dataset-agnostic bounds or normalization.
  • Estimating F(x) remains expensive for extremely large models such as billion-parameter LLMs despite the proposed scalable algorithms.Future directions include more efficient gradient estimators and distributed Hutchinson sampling.
  • The current analysis focuses on classification with softmax outputs, so other task types require reformulating F(x) for continuous or structured outputs.
  • Top RobustBench models have not yet been evaluated with the proposed metric, and applying it to them is identified as future work.
Loading 2606.04767v1…