Are adversarial examples inevitable?

Ali Shafahi, W. Ronny Huang, Christoph Studer, Soheil Feizi, Tom Goldstein

Introduction

A number of adversarial attacks on neural networks have been recently proposed. To counter these attacks, a number of authors have proposed a range of defenses. However, these defenses are often quickly broken by new and revised attacks. Given the lack of success at generating robust defenses, we are led to ask a fundamental question: Are adversarial attacks inevitable?

In this paper, we identify a broad class of problems for which adversarial examples cannot be avoided. We also derive fundamental limits on the susceptibility of a classifier to adversarial attacks that depend on properties of the data distribution as well as the dimensionality of the dataset.

Finally, in Section 8, we explore the causes of adversarial susceptibility in real datasets, and the effect of dimensionality. We present an example image class for which there is no fundamental link between dimensionality and robustness, and argue that the data distribution, and not dimensionality, is the primary cause of adversarial susceptibility.

Adversarial examples, first demonstrated in Szegedy et al. (2013) and Biggio et al. (2013), change the label of an image using small and often imperceptible perturbations to its pixels. A number of defenses have been proposed to harden networks against attacks, but historically, these defenses have been quickly broken. Adversarial training, one of the earliest defenses, successfully thwarted the fast gradient sign method (FGSM) (Goodfellow et al., 2014), one of the earliest and simplest attacks. However, adversarial training with FGSM examples was quickly shown to be vulnerable to more sophisticated multi-stage attacks (Kurakin et al., 2016; Tramèr et al., 2017a). More sophisticated defenses that rely on network distillation (Papernot et al., 2016b) and specialized activation functions (Zantedeschi et al., 2017) were also toppled by strong attacks (Papernot et al., 2016a; Tramèr et al., 2017b; Carlini & Wagner, 2016; 2017a).

The ongoing vulnerability of classifiers was highlighted in recent work by Athalye et al. (2018) and Athalye & Sutskever (2017) that broke an entire suite of defenses presented in ICLR 2018 including thermometer encoding (Buckman et al., 2018), detection using local intrinsic dimensionality (Ma et al., 2018), input transformations such as compression and image quilting (Guo et al., 2017), stochastic activation pruning (Dhillon et al., 2018), adding randomization at inference time (Xie et al., 2017), enhancing the confidence of image labels (Song et al., 2017), and using a generative model as a defense (Samangouei et al., 2018).

Rather than hardening classifiers to attacks, some authors have proposed sanitizing datasets to remove adversarial perturbations before classification. Approaches based on auto-encoders (Meng & Chen, 2017) and GANs (Shen et al., 2017) were broken using optimization-based attacks (Carlini & Wagner, 2017b; a).

A number of “certifiable” defense mechanisms have been developed for certain classifiers. Raghunathan et al. (2018) harden a two-layer classifier using semidefinite programming, and Sinha et al. (2018) propose a convex duality-based approach to adversarial training that works on sufficiently small adversarial perturbations with a quadratic adversarial loss. Kolter & Wong (2017) consider training a robust classifier using the convex outer adversarial polytope. All of these methods only consider robustness of the classifier on the training set, and robustness properties often fail to generalize reliably to test examples.

One place where researchers have enjoyed success is at training classifiers on low-dimensional datasets like MNIST (Madry et al., 2017; Sinha et al., 2018). The robustness achieved on more complicated datasets such as CIFAR-10 and ImageNet are nowhere near that of MNIST, which leads some researchers to speculate that adversarial defense is fundamentally harder in higher dimensions – an issue we address in Section 8.

This paper uses well-known results from high-dimensional geometry, specifically isoperimetric inequalities, to provide bounds on the robustness of classifiers. Several other authors have investigated adversarial susceptibility through the lens of geometry. Wang et al. (2018) analyze the robustness of nearest neighbor classifiers, and provide a more robust NN classifier. Fawzi et al. (2018) study adversarial susceptibility of datasets under the assumption that they are produced by a generative model that maps random Gaussian vectors onto images. Gilmer et al. (2018) do a detailed case study, including empirical and theoretical results, of classifiers for a synthetic dataset that lies on two concentric spheres. Simon-Gabriel et al. (2018) show that the Lipschitz constant of untrained networks with random weights gets large in high dimensions. Shortly after the original appearance of our work, Mahloujifar et al. (2018) presented a study of adversarial susceptibility that included both evasion and poisoning attacks. Our work is distinct in that it studies adversarial robustness for arbitrary data distributions, and also that it rigorously looks at the effect of dimensionality on robustness limits.

2 Notation

Problem setup

We also consider a “classifier” function C:Ω→{1,2,…,m}\mathcal{C}:\Omega\to\{1,2,\ldots,m\} that partitions Ω\Omega into disjoint measurable subsets, one for each class label. The classifier we consider is discrete valued – it provides a label for each data point but not a confidence level.

With this setup, we can give a formal definition of an adversarial example.

Consider a point x∈Ωx\in\Omega drawn from class c,c, a scalar ϵ>0,\epsilon>0, and a metric d.d. We say that xx admits an ϵ\epsilon-adversarial example in the metric dd if there exists a point x^∈Ω\hat{x}\in\Omega with C(x^)≠c,\mathcal{C}(\hat{x})\neq c, and d(x,x^)≤ϵ.d(x,\hat{x})\leq\epsilon.

In plain words, a point has an ϵ\epsilon-adversarial example if we can sneak it into a different class by moving it at most ϵ\epsilon units in the distance dd.

We also consider sparse adversarial examples in which only a small subset of pixels are manipulated. This corresponds to the metric d0,d_{0}, in which case the constraint ∥x−x^∥0≤ϵ\|x-\hat{x}\|_{0}\leq\epsilon means that an adversarial example was crafted by changing at most ϵ\epsilon pixels, and leaving the others alone.

The simple case: adversarial examples on the unit sphere

We begin by looking at the case of classifiers for data on the sphere. While this data model may be less relevant than the other models studied below, it provides a straightforward case where results can be proven using simple, geometric lemmas. The more realistic case of images with pixels in $$ will be studied in Section 4.

The idea is to show that, provided a class of data points takes up enough space, nearly every point in the class lies close to the class boundary. To show this, we begin with a simple definition.

The ϵ\epsilon-expansion of a subset A⊂Ω\mathcal{A}\subset\Omega with respect to distance metric d,d, denoted A(ϵ,d),\mathcal{A}(\epsilon,d), contains all points that are at most ϵ\epsilon units away from A\mathcal{A}. To be precise

We sometimes simply write A(ϵ)\mathcal{A}(\epsilon) when the distance metric is clear from context.

Our result provides bounds on the probability of adversarial examples that are independent of the shape of the class boundary. This independence is a simple consequence of an isoperimetric inequality. The classical isoperimetric inequality states that, of all closed surfaces that enclose a unit volume, the sphere has the smallest surface area. This simple fact is intuitive but famously difficult to prove. For a historical review of the isoperimetric inequality and its variants, see Osserman et al. (1978). We will use a special variant of the isoperimetric inequality first proved by Lévy & Pellegrino (1951) and simplified by Talagrand (1995).

The classical isoperimetric inequality is a simple geometric statement, and frequently appears without absolute bounds on the size of the ϵ\epsilon-expansion of a half-sphere, or with bounds that involve unspecified constants (Vershynin, 2017). A tight bound derived by Milman & Schechtman (1986) is given below. The asymptotic blow-up of the ϵ\epsilon-expansion of a half sphere predicted by this bound is shown in Figure 3.

The geodesic ϵ\epsilon-expansion of a half sphere has normalized measure at least

Lemmas 1 and 2 together can be taken to mean that, if a set is not too small, then in high dimensions almost all points on the sphere are reachable within a short ϵ\epsilon jump from that set. These lemmas have immediate implications for adversarial examples, which are formed by mapping one class into another using small perturbations. Despite its complex appearance, the result below is a consequence of the (relatively simple) isoperimetric inequality.

Let VcV_{c} denote the magnitude of the supremum of ρc\rho_{c} relative to the uniform density. This can be written Vc:=sn−1⋅sup⁡xρc(x).V_{c}:=s_{n-1}\cdot\sup_{x}\rho_{c}(x).

Let fc=μ1{x∣C(x)=c}f_{c}=\mu_{1}\{x|\mathcal{C}(x)=c\} be the fraction of the sphere labeled as cc by classifier C\mathcal{C}.

Choose some class cc with fc≤12f_{c}\leq\frac{1}{2}. Sample a random data point xx from ρc.\rho_{c}. Then with probability at least

xx is misclassified by C,\mathcal{C}, or

xx admits an ϵ\epsilon-adversarial example in the geodesic distance.

Choose a class cc with fc≤12.f_{c}\leq\frac{1}{2}. Let R={x∣C(x)=c}\mathcal{R}=\{x|\mathcal{C}(x)=c\} denote the region of the sphere labeled as class cc by C\mathcal{C}, and let R‾\overline{\mathcal{R}} be its complement. R‾(ϵ)\overline{\mathcal{R}}(\epsilon) is the ϵ\epsilon-expansion of R‾\overline{\mathcal{R}} in the geodesic metric. Because R‾\overline{\mathcal{R}} covers at least half the sphere, the isoperimetric inequality (Lemma 1) tells us that the epsilon expansion is at least as great as the epsilon expansion of a half sphere. We thus have

Now, consider the set Sc\mathcal{S}_{c} of “safe” points from class cc that are correctly classified and do not admit adversarial perturbations. A point is correctly classified only if it lies inside R,\mathcal{R}, and therefore outside of R‾\overline{\mathcal{R}}. To be safe from adversarial perturbations, a point cannot lie within ϵ\epsilon distance from the class boundary, and so it cannot lie within R‾(ϵ)\overline{R}(\epsilon). It is clear that the set Sc\mathcal{S}_{c} of safe points is exactly the complement of R‾(ϵ).\overline{R}(\epsilon). This set has normalized measure

The probability of a random point lying in Sc\mathcal{S}_{c} is bounded above by the normalized supremum of ρc\rho_{c} times the normalized measure μ1[Sc].\mu_{1}[\mathcal{S}_{c}]. This product is given by

We then subtract this probability from 1 to obtain the probability of a point lying outside the safe region, and arrive at equation 1. ∎

It is easily observed that, for any two points xx and yy on a sphere,

where d∞(x,y)d_{\infty}(x,y), d2(x,y)d_{2}(x,y), and dg(x,y)d_{g}(x,y) denote the l∞l_{\infty}, Euclidean, and geodesic distance, respectively. From this, we see that Theorem 1 is actually fairly conservative; any ϵ\epsilon-adversarial example in the geodesic metric would also be adversarial in the other two metrics, and the bound in Theorem 1 holds regardless of which of the three metrics we choose (although different values of ϵ\epsilon will be appropriate depending on the norm).

What about the unit cube?

The above result about the sphere is simple and easy to prove using classical results. However, real world images do not lie on the sphere. In a more typical situation, images will be scaled so that their pixels lie in $$, and data lies inside a high-dimensional hypercube (but, unlike the sphere, data is not confined to its surface). The proof of Theorem 1 makes extensive use of properties that are exclusive to the sphere, and is not applicable to this more realistic setting. Are there still problem classes on the cube where adversarial examples are inevitable?

This question is complicated by the fact that geometric isoperimetric inequalities do not exist for the cube, as the shapes that achieve minimal ϵ\epsilon-expansion (if they exist) depend on the volume they enclose and the choice of ϵ\epsilon (Ros, 2001). Fortunately, researchers have been able to derive “algebraic” isoperimetric inequalities that provide lower bounds on the size of the ϵ\epsilon-expansion of sets without identifying the shape that achieves this minimum (Talagrand, 1996; Milman & Schechtman, 1986). The result below about the unit cube is analogous to Proposition 2.8 in Ledoux (2001), except with tighter constants. For completeness, a proof (which utilizes methods from Ledoux) is provided in Appendix A.

Consider a measurable subset of the cube A⊂n,\mathcal{A}\subset^{n}, and a p-norm distance metric dp(x,y)=∥x−y∥pd_{p}(x,y)=\|x-y\|_{p} for p>0.p>0. Let Φ(z)=(2π)−12∫−∞ze−t2/2dt,\Phi(z)=(2\pi)^{-\frac{1}{2}}\int_{-\infty}^{z}e^{-t^{2}/2}dt, and let α\alpha be the scalar that satisfies Φ(α)=vol⁡[A].\Phi(\alpha)=\operatorname{vol}[\mathcal{A}]. Then

where p∗=min⁡(p,2).p^{*}=\min(p,2). In particular, if vol⁡(A)≥1/2,\operatorname{vol}(\mathcal{A})\geq 1/2, then we simply have

Using this result, we can show that most data samples in a cube admit adversarial examples, provided the data distribution is not excessively concentrated.

Consider a classification problem with mm classes, each distributed over the unit hypercube n^{n} with density functions {ρc}c=1m\{\rho_{c}\}_{c=1}^{m}. Choose a classifier function C:n→{1,2,…,m}{\mathcal{C}:^{n}\to\{1,2,\ldots,m\}} that partitions the hypercube into disjoint measurable subsets. Define the following scalar constants:

Let UcU_{c} denote the supremum of ρc\rho_{c}.

Let fcf_{c} be the fraction of hypercube partitioned into class cc by C\mathcal{C}.

xx is misclassified by C,\mathcal{C}, or

xx has an adversarial example x^,\hat{x}, with ∥x−x^∥p≤ϵ\|x-\hat{x}\|_{p}\leq\epsilon.

for p≥2,p\geq 2, where Φ^(z)=12π∫z∞e−t2/2dt≥12πze−z2/2\hat{\Phi}(z)=\frac{1}{\sqrt{2\pi}}\int_{z}^{\infty}e^{-t^{2}/2}dt\geq\frac{1}{\sqrt{2\pi}z}e^{-z^{2}/2} (for z>0z>0), and α=Φ−1(1−fc)\alpha=\Phi^{-1}(1-f_{c}). For this bound to be meaningful with ϵ<1\epsilon<1, we need fcf_{c} to be relatively small, and ϵ\epsilon to be roughly fcf_{c} or smaller. This is realistic for some problems; ImageNet has 1000 classes, and so fc<10−3f_{c}<10^{-3} for at least one class.

What about sparse adversarial examples?

If a point xx has an ϵ\epsilon-adversarial example in this norm, then it can be perturbed into a different class by modifying at most ϵ\epsilon pixels (in this case ϵ\epsilon is taken to be a positive integer).

Theorem 2 is fairly tight for p=1p=1 or 22. However, the bound becomes quite loose for small p,p, and in particular it fails completely for the important case of p=0.p=0. For this reason, we present a different bound that is considerably tighter for small pp (although slightly looser for large pp).

The case p=0p=0 was studied by Milman & Schechtman (1986) (Section 6.2) and McDiarmid (1989), and later by Talagrand (1995; 1996). The proof of the following theorem (appendix B) follows the method used in Section 5 of Talagrand (1996), with modifications made to extend the proof to arbitrary pp.

Consider a measurable subset of the cube A⊂n,\mathcal{A}\subset^{n}, and a p-norm distance metric d(x,y)=∥x−y∥pd(x,y)=\|x-y\|_{p} for any p≥0.p\geq 0. We have

Using this result, we can prove a statement analogous to Theorem 2, but for sparse adversarial examples. We present only the case of p=0,p=0, but the generalization to the case of other small pp using Lemma 4 is straightforward.

Consider the problem setup of Theorem 2. Choose some class cc with fc≤12f_{c}\leq\frac{1}{2}, and sample a random data point xx from the class distribution ρc.\rho_{c}. Then with probability at least

xx is misclassified by C,\mathcal{C}, or

xx can be adversarially perturbed by modifying at most ϵ\epsilon pixels, while still remaining in the unit hypercube.

What if we just show that adversarial examples exist?

Tighter bounds can be obtained if we only guarantee that adversarial examples exist for some data points in a class, without bounding the probability of this event.

Let supp⁡(ρc)\operatorname{supp}(\rho_{c}) denote the support of ρc.\rho_{c}. Then there is a point xx with ρc(x)>0\rho_{c}(x)>0 that admits an ϵ\epsilon-adversarial example if

The bound for the case p=0p=0 is valid only if ϵ≥nlog⁡2/2.\epsilon\geq\sqrt{n\log 2/2}.

Discussion: Can we escape fundamental bounds?

There are a number of ways to escape the guarantees of adversarial examples made by Theorems 1-4. One potential escape is for the class density functions to take on extremely large values (i.e., exponentially large UcU_{c}); the dependence of UcU_{c} on nn is addressed separately in Section 8.

In practice, image datasets might lie on low-dimensional manifolds within the cube, and the support of these distributions could have measure zero, making the density function infinite (i.e., Uc=∞U_{c}=\infty). The arguments above are still relevant (at least in theory) in this case; we can expand the data manifold by adding a uniform random noise to each image pixel of magnitude at most ϵ1.\epsilon_{1}. The expanded dataset has positive volume. Then, adversarial examples of this expanded dataset can be crafted with perturbations of size ϵ2\epsilon_{2}. This method of expanding the manifold before crafting adversarial examples is often used in practice. Tramèr et al. (2017a) proposed adding a small perturbation to step off the image manifold before crafting adversarial examples. This strategy is also used during adversarial training (Madry et al., 2017).

Adding a “don’t know” class

The analysis above assumes the classifier assigns a label to every point in the cube. If a classifier has the ability to say “I don’t know,” rather than assign a label to every input, then the region of the cube that is assigned class labels might be very small, and adversarial examples could be escaped even if the other assumptions of Theorem 4 are satisfied. In this case, it would still be easy for the adversary to degrade classifier performance by perturbing images into the “don’t know” class.

Feature squeezing

If decreasing the dimensionality of data does not lead to substantially increased values for UcU_{c} (we see in Section 8 that this is a reasonable assumption) or loss in accuracy (a stronger assumption), measuring data in lower dimensions could increase robustness. This can be done via an auto-encoder (Meng & Chen, 2017; Shen et al., 2017), JPEG encoding (Das et al., 2018), or quantization (Xu et al., 2017).

Computational hardness

It may be computationally hard to craft adversarial examples because of local flatness of the classification function, obscurity of the classifier function, or other computational difficulties. Computational hardness could prevent adversarial attacks in practice, even if adversarial examples still exist.

Experiments & Effect of Dimensionality

In this section, we discuss the relationship between dimensionality and adversarial robustness, and explore how the predictions made by the theorems above are reflected in experiments.

It is commonly thought that high-dimensional classifiers are more susceptible to adversarial examples than low-dimensional classifiers. This perception is partially motivated by the observation that classifiers on high-resolution image distributions like ImageNet are more easily fooled than low resolution classifiers on MNIST (Tramèr et al., 2017a). Indeed, Theorem 2 predicts that high-dimensional classifiers should be much easier to fool than low-dimensional classifiers, assuming the datasets they classify have comparable probability density limits Uc.U_{c}. However, this is not a reasonable assumption; we will see below that high dimensional distributions may be more concentrated than their low-dimensional counterparts.

We study the effects of dimensionality with a thought experiment involving a “big MNIST” image distribution. Given an integer expansion factor b,b, we can make a big MNIST distribution, denoted bb-MNIST, by replacing each pixel in an MNIST image with a b×bb\times b array of identical pixels. This expands an original 28×2828\times 28 image into a 28b×28b28b\times 28b image. Figure 4a shows that, without adversarial training, a classifier on big MNIST is far more susceptible to attacks than a classifier trained on the original MNIST Only the fully-connected layer is modified to handle the difference in dimensionality between datasets..

However, each curve in Figure 4a only shows the attack susceptibility of one particular classifier. In contrast, Theorems 1-4 describe the fundamental limits of susceptibility for all classifiers. These limits are an inherent property of the data distribution. The theorem below shows that these fundamental limits do not depend in a non-trivial way on the dimensionality of the images in big MNIST, and so the relationship between dimensionality and susceptibility in Figure 4a results from the weakness of the training process.

Likewise, if all bb-MNIST classifiers have bϵb\epsilon-adversarial examples with probability pp for some b≥1b\geq 1, then all classifiers on the original MNIST distribution have ϵ\epsilon-adversarial examples with probability pp.

We get a better picture of the fundamental limits of MNIST by considering classifiers that are hardened by adversarial trainingAdversarial examples for MNIST/CIFAR-10 were produced as in Madry et al. (2017) using 100-step/20-step PGD. (Figure 4b). These curves display several properties of fundamental limits predicted by our theorems. As predicted by Theorem 5, the 112×112112\times 112 classifer curve is twice as wide as the 56×5656\times 56 curve, which in turn is twice as wide as the 28×2828\times 28 curve. In addition, we see the kind of “phase transition” behavior predicted by Theorem 2, in which the classifier suddenly changes from being highly robust to being highly susceptible as ϵ\epsilon passes a critical threshold. For these reasons, it is reasonable to suspect that the adversarially trained classifiers in Figure 4b are operating near the fundamental limit predicted by Theorem 2.

Theorem 5 shows that increased dimensionality does not increase adversarial susceptibility in a fundamental way. But then why are high-dimensional classifiers so easy to fool? To answer this question, we look at the concentration bound UcU_{c} for object classes. The smallest possible value of UcU_{c} is 1, which only occurs when images are “spread out” with uniform, uncorrelated pixels. In contrast, adjacent pixels in MNIST (and especially big MNIST) are very highly correlated, and images are concentrated near simple, low-dimensional manifolds, resulting in highly concentrated image classes with large UcU_{c}. Theory predicts that such highly concentrated datasets can be relatively safe from adversarial examples.

We can reduce UcU_{c} and dramatically increase susceptibility by choosing a more “spread out” dataset, like CIFAR-10, in which adjacent pixels are less strongly correlated and images appear to concentrate near complex, higher-dimensional manifolds. We observe the effect of decreasing UcU_{c} by plotting the susceptibility of a 56×5656\times 56 MNIST classifier against a classifier for CIFAR-10 (Figure 4, right). The former problem lives in 31363136 dimensions, while the latter lives in 3072,3072, and both have 10 classes. Despite the structural similarities between these problems, the decreased concentration of CIFAR-10 results in vastly more susceptibility to attacks, regardless of whether adversarial training is used. The theory above suggests that this increased susceptibility is caused at least in part by a shift in the fundamental limits for CIFAR-10, rather than the weakness of the particular classifiers we chose.

Informally, the concentration limit UcU_{c} can be interpreted as a measure of image complexity. Image classes with smaller UcU_{c} are likely concentrated near high-dimensional complex manifolds, have more intra-class variation, and thus more apparent complexity. An informal interpretation of Theorem 2 is that “high complexity” image classes are fundamentally more susceptible to adversarial examples, and Figure 4 suggests that complexity (rather than dimensionality) is largely responsible for differences we observe in the effectiveness of adversarial training for different datasets.

So…are adversarial examples inevitable?

The question of whether adversarial examples are inevitable is an ill-posed one. Clearly, any classification problem has a fundamental limit on robustness to adversarial attacks that cannot be escaped by any classifier. However, we have seen that these limits depend not only on fundamental properties of the dataset, but also on the strength of the adversary and the metric used to measure perturbations. This paper provides a characterization of these limits and how they depend on properties of the data distribution. Unfortunately, it is impossible to know the exact properties of real-world image distributions or the resulting fundamental limits of adversarial training for specific datasets. However, the analysis and experiments in this paper suggest that, especially for complex image classes in high-dimensional spaces, these limits may be far worse than our intuition tells us.

Acknowledgements

T. Goldstein and his students were generously supported by DARPA Lifelong Learning Machines (FA8650-18-2-7833), the Office of Naval Research (N00014-17-1-2078), the DARPA Young Faculty Award program (D18AP00055), and the Sloan Foundation. The work of C. Studer was supported in part by Xilinx Inc., and by the US NSF under grants ECCS-1408006, CCF-1535897, CCF-1652065, and CNS-1717559.

References

Appendix A Proof of lemma 3

We now prove Lemma 3. To do this, we begin with a classical isoperimetric inequality for random Gaussian variables. Unlike the case of a cube, tight geometric isoperimetric inequalities exist in this case. We then prove results about the cube by creating a mapping between uniform random variables on the cube and random Gaussian vectors.

which is the cumulative density of a Gaussian curve.

The following Lemma was first proved in Sudakov & Tsirelson (1974), and an elementary proof was given in Bobkov et al. (1997).

Using this result we can now give a proof of Lemma 3.

Since ∂∂ziΦ(z)≤12π,\frac{\partial}{\partial z_{i}}\Phi(z)\leq\frac{1}{\sqrt{2\pi}}, we also have

where we have used the identity ∥u∥p≤n1/min⁡(p,2)−1/2∥u∥2.\|u\|_{p}\leq n^{1/\min(p,2)-1/2}\|u\|_{2}.

Now, consider any set A\mathcal{A} in the cube, and let B=Φ−1(A).\mathcal{B}=\Phi^{-1}(\mathcal{A}). From equation 10, we see that

where α=Φ−1(σ[A]).\alpha=\Phi^{-1}(\sigma[\mathcal{A}]).

To obtain the simplified formula in the theorem, we use the identity

which is valid for x>0,x>0, and can be found in Abramowitz & Stegun (1965).

Appendix B Proof of Lemma 4

Consider a probability space Ω\Omega with measure μ.\mu. For g:Ω→,g:\Omega\to, we have

Our proof of Lemma 3 follows the three-step process of Talagrand illustrated in Talagrand (1995). We begin by proving the bound

where f(x,A)=min⁡y∈A∑i∣xi−yi∣p=dpp(x,A)f(x,\mathcal{A})=\min_{y\in\mathcal{A}}\sum_{i}|x_{i}-y_{i}|^{p}=d_{p}^{p}(x,\mathcal{A}) is a measure of distance from A\mathcal{A} to xx, and α,t\alpha,t are arbitrary positive constants. Once this bound is established, a Markov bound can be used to obtain the final result. Finally, constants are tuned in order to optimize the tightness of the bound.

We start by proving the bound in equation 12 using induction on the dimension. The base case for the induction is n=1,n=1, and we have

We now prove the result for nn dimensions using the inductive hypothesis. We can upper bound the integral by integrating over “slices” along one dimension. Let A⊂n.\mathcal{A}\subset^{n}. Define

Clearly, the distance from (ω,z)(\omega,z) to A\mathcal{A} is at most the distance from zz to Aω,\mathcal{A}_{\omega}, and so

We also have that the distance from xx to A\mathcal{A} is at most one unit greater than the distance from xx to B\mathcal{B}. This gives us

Now, we apply lemma 6 to equation 13 with g(ω)=α[Aω]/α[B]g(\omega)=\alpha[\mathcal{A}_{\omega}]/\alpha[\mathcal{B}] to arrive at equation 12.

The second step of the proof is to produce a Markov inequality from equation 12. For the bound in equation 12 to hold, we need

The third step is to optimize the bound by choosing constants. We minimize the right hand side by choosing t=4αϵpn(α+1)t=\frac{4\alpha\epsilon^{p}}{n(\alpha+1)} to get

Now, we can simply choose α=1\alpha=1 to get the simple bound

or we can choose the optimal value of α=2ϵ2pnlog⁡(1/σ)−1,\alpha=\sqrt{\frac{2\epsilon^{2p}}{n\log(1/\sigma)}}-1, which optimizes the bound in the case ϵ2p≥n2log⁡(1/σ(A)).\epsilon^{2p}\geq\frac{n}{2}\log(1/\sigma(\mathcal{A})). We arrive at

This latter bound is stronger than we need to prove Lemma 3, but it will come in handy later to prove Theorem 4.

Appendix C Proof of Theorems 2 and 3

We combine the proofs of these results since their proofs are nearly identical. The proofs closely follow the argument of Theorem 1.

The set R‾(ϵ;h)‾\overline{\overline{\mathcal{R}}(\epsilon;h)} contains all points that are correctly classified and safe from adversarial perturbations. This region has volume at most δ\delta, and the probability of a sample from the class distribution ρc\rho_{c} lying in this region is at most Ucδ.U_{c}\delta. We then subtract this from 1 to obtain the mass of the class distribution lying in the “unsafe” region R‾c.\overline{\mathcal{R}}_{c}.

Appendix D Proof of Theorem 4

Let A\mathcal{A} denote the support of pc,p_{c}, and suppose that this support has measure vol⁡[A]=η.\operatorname{vol}[\mathcal{A}]=\eta. We want to show that, for large enough ϵ,\epsilon, the expansion A(ϵ,dp)\mathcal{A}(\epsilon,d_{p}) is larger than half the cube. Since class cc occupies less than half the cube, this would imply that A(ϵ,dp)\mathcal{A}(\epsilon,d_{p}) overlaps with other classes, and so there must be data points in A\mathcal{A} with ϵ\epsilon-adversarial examples.

We start with the case p>0,p>0, where we bound A(ϵ,dp)\mathcal{A}(\epsilon,d_{p}) using equation 2 of Lemma 3. To do this, we need to approximate Φ−1(η).\Phi^{-1}(\eta). This can be done using the inequality

which holds for α<0.\alpha<0. Rearranging, we obtain

Now, if α=Φ−1(η),\alpha=\Phi^{-1}(\eta), then Φ(α)=η,\Phi(\alpha)=\eta, and equation 19 gives us α≥−log⁡14η2.\alpha\geq-\sqrt{\log\frac{1}{4\eta^{2}}}. Plugging this into equation 2 of Lemma 3, we get

The quantity on the left will be greater than 12,\frac{1}{2}, thus guaranteeing adversarial examples, if

This can be re-arranged to obtain the desired result.

In the case p=0,p=0, we need to use equation 17 from the proof of Lemma 3 in Appendix B, which we restate here

This bound is valid, and produces a non-vacuous guarantee of adversarial examples, if

Appendix E Proof of Theorem 5

Assume that any MNIST classifier can be fooled by perturbations of size at most ϵ\epsilon with probability at least pp. To begin, we put a bound on the susceptibility of any bb-MNIST classifier (for b≥1b\geq 1) under this assumption. We can classify MNIST images by upsampling them to resolution 28b×28b28b\times 28b and feeding them into a high-resolution “back-end” classifier. After upsampling, an MNIST image with perturbation of norm ϵ\epsilon becomes a 28b×28b28b\times 28b image with perturbation of norm bϵ.b\epsilon. The classifier we have constructed takes low-resolution images as inputs, and so by assumption it is fooled with probability at least pp. However, the low-resolution classifier is fooled only when the high-resolution “back-end” classifier is fooled, and so the high-resolution classifier is fooled with probability at least pp at well. Note that we can build this two-scale classifier using any high-resolution classifier as a back-end, and so this bound holds uniformly over all high-resolution classifiers.

Likewise, suppose we classify bb-MNIST images (for integer b≥1b\geq 1) by downsampling them to the original 28×2828\times 28 resolution (by averaging pixel blocks) and feeding them into a “back-end” low-resolution classifier. After downsampling, a 28b×28b28b\times 28b image with perturbation of norm bϵb\epsilon becomes a 28×2828\times 28 image with perturbation of norm at most ϵ.\epsilon. Whenever the high-resolution classifier is fooled, it is only because the back-end classifier is fooled by a perturbation of size at most ϵ,\epsilon, and this happens with probability at least pp.