arXiv is now an independent nonprofit! Learn more
License: arXiv.org perpetual non-exclusive license
arXiv:1809.02104v3 [cs.LG] 03 Feb 2020

Are adversarial examples inevitable?

Ali Shafahi    Ronny Huang    Christoph Studer    Soheil Feizi & Tom Goldstein
Abstract

A wide range of defenses have been proposed to harden neural networks against adversarial attacks. However, a pattern has emerged in which the majority of adversarial defenses are quickly broken by new attacks. Given the lack of success at generating robust defenses, we are led to ask a fundamental question: Are adversarial attacks inevitable? This paper analyzes adversarial examples from a theoretical perspective, and identifies fundamental bounds on the susceptibility of a classifier to adversarial attacks. We show that, for certain classes of problems, adversarial examples are inescapable. Using experiments, we explore the implications of theoretical guarantees for real-world problems and discuss how factors such as dimensionality and image complexity limit a classifier’s robustness against adversarial examples.

1 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.

Original

Refer to caption

egyptian cat (28%)

2\ell_{2}-norm=10 Refer to caption traffic light (97%)

\ell_{\infty}-norm=0.05 Refer to caption traffic light (96%)

0\ell_{0}-norm=5000 (sparse) Refer to caption traffic light (80%)

Figure 1: Adversarial examples with different norm constraints formed via the projected gradient method (Madry et al. 2017) on Resnet50, along with the distance between the base image and the adversarial example, and the top class label.

Adversarial examples occur when a small perturbation to an image changes its class label. There are different ways of measuring what it means for a perturbation to be “small”; as such, our analysis considers a range of different norms. While the \ell_{\infty}-norm is commonly used, adversarial examples can be crafted in any p\ell_{p}-norm (see Figure 1). We will see that the choice of norm can have a dramatic effect on the strength of theoretical guarantees for the existence of adversarial examples. Our analysis also extends to the 0\ell_{0}-norm, which yields “sparse” adversarial examples that only perturb a small subset of image pixels (Figure 2).

As a simple example result, consider a classification problem with nn-dimensional images with pixels scaled between 0 and 1 (in this case images live inside the unit hypercube). If the image classes each occupy a fraction of the cube greater than 12exp(πϵ2)\frac{1}{2}\exp(-\pi\epsilon^{2}), then images exist that are susceptible to adversarial perturbations of 2\ell_{2}-norm at most ϵ\epsilon. Note that ϵ=10\epsilon=10 was used in Figure 1, and larger values are typical for larger images.

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.

original Refer to caption

\ell_{\infty}-norm Refer to caption

0\ell_{0}-norm (sparse) Refer to caption

sparse perturbation Refer to caption

Figure 2: Sparse adversarial examples perturb a small subset of pixels and can hide adversarial “fuzz” inside high-frequency image regions. The original image (left) is classified as an “ox.” Under \ell_{\infty}-norm perturbations, it is classified as “traffic light”, but the perturbations visibly distort smooth regions of the image (the sky). These effects are hidden in the grass using 0\ell_{0}-norm (sparse) perturbations limited to a small subset of pixels.

1.1 Background: A brief history of adversarial examples

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; Carlini & Wagner 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; Carlini & Wagner 2017a).

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.

1.2 Notation

We use [0,1]n[0,1]^{n} to denote the unit hypercube in nn dimensions, and vol(𝒜)\vol(\mathcal{A}) to denote the volume (i.e., n-dimensional Lebesgue measure) of a subset 𝒜[0,1]n\mathcal{A}\subset[0,1]^{n}. We use 𝕊n1={xn|x2=1}\mathbb{S}^{n-1}=\{x\in\mathbb{R}^{n}|\,\|x\|_{2}=1\} to denote the unit sphere embedded in n,\mathbb{R}^{n}, and sn1s_{n-1} to denote its surface area. The size of a subset 𝒜𝕊n1\mathcal{A}\in\mathbb{S}^{n-1} can be quantified by its (n1n-1 dimensional) measure μ[𝒜],\mu[\mathcal{A}], which is simply the surface area the set covers. Because the surface area of the unit sphere varies greatly with nn, it is much easier in practice to work with the normalized measure, which we denote μ1[𝒜]=μ[𝒜]/sn1\mu_{1}[\mathcal{A}]=\mu[\mathcal{A}]/s_{n-1}. This normalized measure has the property that μ1[𝕊n1]=1,\mu_{1}[\mathbb{S}^{n-1}]=1, and so we can interpret μ1[𝒜]\mu_{1}[\mathcal{A}] as the probability of a uniform random point from the sphere lying in 𝒜.\mathcal{A}. When working with points on a sphere, we often use geodesic distance, which is always somewhat larger than (but comparable to) the Euclidean distance. In the cube, we measure distance between points using p\ell_{p}-norms, which are denoted

zp\displaystyle\|z\|_{p} =(i|zi|p)1/p if p>0, and z0=card{zi|zi0}.\displaystyle=\left(\sum_{i}|z_{i}|^{p}\right)^{1/p}\text{ if }p>0\text{, and }\|z\|_{0}=\card\{z_{i}\,|\,z_{i}\neq 0\}.

Note that p\|\cdot\|_{p} is not truly a norm for p<1,p<1, but rather a semi-norm. Such metrics are still commonly used, particularly the “0\ell_{0}-norm” which counts the number of non-zero entries in a vector.

2 Problem setup

We consider the problem of classifying data points that lie in a space Ω\Omega (either a sphere or a hypercube) into mm different object classes. The mm object classes are defined by probability density functions {ρc}c=1m,\{\rho_{c}\}_{c=1}^{m}, where ρc:Ω\rho_{c}:\Omega\to\mathbb{R}. A “random” point from class cc is a random variable with density ρc.\rho_{c}. We assume ρc\rho_{c} to be bounded (i.e., we don’t allow delta functions or other generalized functions), and denote its upper bound by Uc=supxρc(x).U_{c}=\sup_{x}\rho_{c}(x).

We also consider a “classifier” function 𝒞:Ω{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.

Definition 1.

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 𝒞(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 consider adversarial examples with respect to different p\ell_{p}-norm metrics. These metrics are written dp(x,x^)=xx^p.d_{p}(x,\hat{x})=\|x-\hat{x}\|_{p}. A common choice is p=,p=\infty, which limits the absolute change that can be made to any one pixel. However, 2\ell_{2}-norm and 1\ell_{1}-norm adversarial examples are also used, as it is frequently easier to create adversarial examples in these less restrictive metrics.

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 xx^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.

3 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 [0,1][0,1] 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.

Definition 2.

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

𝒜(ϵ,d)={xΩ|d(x,y)ϵ for some y𝒜}.\mathcal{A}(\epsilon,d)=\{x\in\Omega\,|\,d(x,y)\leq\epsilon\text{ for some }y\in\mathcal{A}\}.

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

Refer to caption
Refer to caption
Figure 3: The ϵ\epsilon-expansion of a half sphere nearly covers the whole sphere for small ϵ\epsilon and large nn. Visualizations show the fraction of the sphere captured within ϵ\epsilon units of a half sphere in different dimensions. Results from a near-exact experimental method are compared to the theoretical lower bound in Lemma 2.

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.

Lemma 1 (Isoperimetric inequality).

Consider a subset of the sphere 𝒜𝕊n1n\mathcal{A}\subset\mathbb{S}^{n-1}\subset\mathbb{R}^{n} with normalized measure μ1(𝒜)1/2\mu_{1}(\mathcal{A})\geq 1/2. When using the geodesic metric, the ϵ\epsilon-expansion 𝒜(ϵ)\mathcal{A}(\epsilon) is at least as large as the ϵ\epsilon-expansion of a half sphere.

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.

Lemma 2 (ϵ\epsilon-expansion of half sphere).

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

1(π8)12exp(n12ϵ2).1-\left(\frac{\pi}{8}\right)^{\frac{1}{2}}\exp\left(-\frac{n-1}{2}\epsilon^{2}\right).

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.

Theorem 1 (Existence of Adversarial Examples).

Consider a classification problem with mm object classes, each distributed over the unit sphere 𝕊n1n\mathbb{S}^{n-1}\subset\mathbb{R}^{n} with density functions {ρc}c=1m\{\rho_{c}\}_{c=1}^{m}. Choose a classifier function 𝒞:𝕊n1{1,2,,m}\mathcal{C}:\mathbb{S}^{n-1}\to\{1,2,\ldots,m\} that partitions the sphere into disjoint measurable subsets. Define the following scalar constants:

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

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

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

1Vc(π8)12exp(n12ϵ2)\displaystyle 1-V_{c}\left(\frac{\pi}{8}\right)^{\frac{1}{2}}\exp\left(-\frac{n-1}{2}\epsilon^{2}\right) (1)

one of the following conditions holds:

  1. 1.

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

  2. 2.

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

Proof.

Choose a class cc with fc12.f_{c}\leq\frac{1}{2}. Let ={x|𝒞(x)=c}\mathcal{R}=\{x|\mathcal{C}(x)=c\} denote the region of the sphere labeled as class cc by 𝒞\mathcal{C}, and let ¯\overline{\mathcal{R}} be its complement. ¯(ϵ)\overline{\mathcal{R}}(\epsilon) is the ϵ\epsilon-expansion of ¯\overline{\mathcal{R}} in the geodesic metric. Because ¯\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

μ1[R¯(ϵ)]1(π8)12exp(n12ϵ2).\mu_{1}[\overline{R}(\epsilon)]\geq 1-\left(\frac{\pi}{8}\right)^{\frac{1}{2}}\exp\left(-\frac{n-1}{2}\epsilon^{2}\right).

Now, consider the set 𝒮c\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 ,\mathcal{R}, and therefore outside of ¯\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 𝒮c\mathcal{S}_{c} of safe points is exactly the complement of R¯(ϵ).\overline{R}(\epsilon). This set has normalized measure

μ1[𝒮c](π8)12exp(n12ϵ2).\mu_{1}[\mathcal{S}_{c}]\leq\left(\frac{\pi}{8}\right)^{\frac{1}{2}}\exp\left(-\frac{n-1}{2}\epsilon^{2}\right).

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

Vc(π8)12exp(n12ϵ2).V_{c}\left(\frac{\pi}{8}\right)^{\frac{1}{2}}\exp\left(-\frac{n-1}{2}\epsilon^{2}\right).

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

In the above result, we measure the size of adversarial perturbations using the geodesic distance. Most studies of adversarial examples measure the size of perturbation in either the 2\ell_{2} (Euclidean) norm or the \ell_{\infty} (max) norm, and so it is natural to wonder whether Theorem 1 depends strongly on the distance metric. Fortunately (or, rather unfortunately) it does not.

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

d(x,y)d2(x,y)dg(x,y),d_{\infty}(x,y)\leq d_{2}(x,y)\leq d_{g}(x,y),

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 ll_{\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).

4 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 [0,1][0,1], 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.

Lemma 3 (Isoperimetric inequality on a cube).

Consider a measurable subset of the cube 𝒜[0,1]n,\mathcal{A}\subset[0,1]^{n}, and a p-norm distance metric dp(x,y)=xypd_{p}(x,y)=\|x-y\|_{p} for p>0.p>0. Let Φ(z)=(2π)12zet2/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[𝒜].\Phi(\alpha)=\vol[\mathcal{A}]. Then

vol[𝒜(ϵ,dp)]Φ(α+2πnn1/pϵ)\displaystyle\vol[\mathcal{A}(\epsilon,d_{p})]\geq\Phi\left(\alpha+\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon\right) (2)

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

vol[𝒜(ϵ,dp)]1exp(2πn12/pϵ2)2πϵn1/21/p.\displaystyle\vol[\mathcal{A}(\epsilon,d_{p})]\geq 1-\frac{\exp(-2\pi n^{1-2/p^{*}}\epsilon^{2})}{2\pi\epsilon n^{1/2-1/p^{*}}}. (3)

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

Theorem 2 (Adversarial examples on the cube).

Consider a classification problem with mm classes, each distributed over the unit hypercube [0,1]n[0,1]^{n} with density functions {ρc}c=1m\{\rho_{c}\}_{c=1}^{m}. Choose a classifier function 𝒞:[0,1]n{1,2,,m}{\mathcal{C}:[0,1]^{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 𝒞\mathcal{C}.

Choose some class cc with fc12,f_{c}\leq\frac{1}{2}, and select an p\ell_{p}-norm with p>0p>0. Define p=min(p,2).p^{*}=\min(p,2). Sample a random data point xx from the class distribution ρc.\rho_{c}. Then with probability at least

1Ucexp(2πn12/pϵ2)2πϵn1/21/p.\displaystyle 1-U_{c}\frac{\exp(-2\pi n^{1-2/p^{*}}\epsilon^{2})}{2\pi\epsilon n^{1/2-1/p^{*}}}. (4)

one of the following conditions holds:

  1. 1.

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

  2. 2.

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

When adversarial examples are defined in the 2\ell_{2}-norm (or for any p2p\geq 2), the bound in equation 4 becomes

1Ucexp(2πϵ2)/(2πϵ).\displaystyle 1-U_{c}\exp(-2\pi\epsilon^{2})/(2\pi\epsilon). (5)

Provided the class distribution is not overly concentrated, equation 5 guarantees adversarial examples with relatively “small” ϵ\epsilon relative to a typical vector. In nn dimensions, the 2\ell_{2} diameter of the cube is n,\sqrt{n}, and so it is reasonable to choose ϵ=O(n)\epsilon=O(\sqrt{n}) in equation 5. In Figure 1, we chose ϵ=10.\epsilon=10. A similarly strong bound of 1Ucnexp(2πϵ2/n)/(2πϵ)1-U_{c}\sqrt{n}\exp(-2\pi\epsilon^{2}/n)/(2\pi\epsilon) holds for the case of the 1\ell_{1}-norm, in which case the diameter is n.n.

Oddly, equation 4 seems particularly weak when the \ell_{\infty} norm is used. In this case, the bound on the right side of equation 4 becomes equation 5 just like in the 2\ell_{2} case. However, \ell_{\infty} adversarial examples are only interesting if we take ϵ<1\epsilon<1, in which case equation 5 becomes vacuous for large nn. This bound can be tightened up in certain situations. If we prove Theorem 2 using the tighter (but messier) bound of equation 2 instead of equation 3, we can replace equation 4 with

1UcΦ^(α+2πϵ)1-U_{c}\hat{\Phi}\left(\alpha+\sqrt{2\pi}\epsilon\right)

for p2,p\geq 2, where Φ^(z)=12πzet2/2dt12πzez2/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(1fc)\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<103f_{c}<10^{-3} for at least one class.

Interestingly, under \ell_{\infty}-norm attacks, guarantees of adversarial examples are much stronger on the sphere (Section 3) than on the cube. One might wonder whether the weakness of Theorem 4 in the \ell_{\infty} case is fundamental, or if this is a failure of our approach. One can construct examples of sets with \ell_{\infty} expansions that nearly match the behavior of equation 5, and so our theorems in this case are actually quite tight. It seems to be inherently more difficult to prove the existence of adversarial examples in the cube using the \ell_{\infty}-norm.

5 What about sparse adversarial examples?

A number of papers have looked at sparse adversarial examples, in which a small number of image pixels, in some cases only one (Su et al. 2017), are changed to manipulate the class label. To study this case, we would like to investigate adversarial examples under the 0\ell_{0} metric. The 0\ell_{0} distance is defined as

d(x,y)=xy0=card{i|xiyi}.d(x,y)=\|x-y\|_{0}=\text{card}\{i\,|\,x_{i}\neq y_{i}\}.

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; Talagrand 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.

Lemma 4 (Isoperimetric inequality on the cube: small pp).

Consider a measurable subset of the cube 𝒜[0,1]n,\mathcal{A}\subset[0,1]^{n}, and a p-norm distance metric d(x,y)=xypd(x,y)=\|x-y\|_{p} for any p0.p\geq 0. We have

vol[𝒜(ϵ,dp)]\displaystyle\vol[\mathcal{A}(\epsilon,d_{p})] 1exp(ϵ2p/n)vol[𝒜], for p>0 and\displaystyle\geq 1-\frac{\exp\left(-\epsilon^{2p}/n\right)}{\vol[\mathcal{A}]}\text{, for }p>0\text{ and} (6)
vol[𝒜(ϵ,d0)]\displaystyle\vol[\mathcal{A}(\epsilon,d_{0})] 1exp(ϵ2/n)vol[𝒜], for p=0.\displaystyle\geq 1-\frac{\exp\left(-\epsilon^{2}/n\right)}{\vol[\mathcal{A}]}\text{, for }p=0. (7)

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.

Theorem 3 (Sparse adversarial examples).

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

12Ucexp(ϵ2/n)\displaystyle 1-2U_{c}\exp(-\epsilon^{2}/n) (8)

one of the following conditions holds:

  1. 1.

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

  2. 2.

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

6 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.

Theorem 4 (Condition for existence of adversarial examples).

Consider the setup of Theorem 2. Choose a class cc that occupies a fraction of the cube fc<12.f_{c}<\frac{1}{2}. Pick an p\ell_{p} norm and set p=min(p,2).p^{*}=\min(p,2).

Let supp(ρc)\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

vol[supp(ρc)]{12exp(πϵ2n12/p), for p>0 orexp(2(ϵnlog22)2/n), for p=0.\displaystyle\vol[\supp(\rho_{c})]\geq\begin{cases}\frac{1}{2}\exp(-\pi\epsilon^{2}n^{1-2/p^{*}}),&\text{ for $p>0$ or}\\ \exp\left(-2\left(\epsilon-\sqrt{\frac{n\log 2}{2}}\right)^{2}/n\right),&\text{ for $p=0.$}\end{cases} (9)

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

It is interesting to consider when Theorem 4 produces non-vacuous bounds. When the 2\ell_{2}-norm is used, the bound becomes vol[supp(ρc)]exp(πϵ2)/2.\vol[\supp(\rho_{c})]\geq\exp(-\pi\epsilon^{2})/2. The diameter of the cube is n,\sqrt{n}, and so the bound becomes active for ϵ=n.\epsilon=\sqrt{n}. Plugging this in, we see that the bound is active whenever the size of the support satisfies vol[supp(ρc)]>12eπn.\vol[\supp(\rho_{c})]>\frac{1}{2e^{\pi n}}. Remarkably, this holds for large nn whenever the support of class cc is larger than (or contains) a hypercube of side length at least eπ0.043.e^{-\pi}\approx 0.043. It is important note, however, that the bound being “active” does not mean there are adversarial examples with a “small” ϵ.\epsilon.

7 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.

Unbounded density functions and low-dimensional data manifolds

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.

8 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 MNIST11 1 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.

Theorem 5.

Suppose ϵ\epsilon and pp are such that, for all MNIST classifiers, a random image from class cc has an ϵ\epsilon-adversarial example (in the 2\ell_{2}-norm) with probability at least pp. Then for all classifiers on bb-MNIST, with integer b1,b\geq 1, a random image from cc has a bϵb\epsilon-adversarial example with probability at least pp.

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

(a) MNIST at different resolutions

(b) MNIST with adver. training

(c) CIFAR-10 vs big MNIST

Figure 4: (a) Robustness of MNIST and “big” MNIST classifiers as a function of ϵ\epsilon. Naturally trained classifiers are less robust with increased dimensionality. (b) With adversarial training, susceptibility curves behave as predicted by Theorems 2 and 5. (c) The susceptibility of CIFAR-10 is compared to big MNIST. Both datasets have similar dimension, but the higher complexity of CIFAR-10 results in far worse susceptibility. Perturbations are measured in the 2\ell_{2}-norm.

Theorem 5 predicts that the perturbation needed to fool all 56×5656\times 56 classifiers is twice that needed to fool all 28×2828\times 28 classifiers. This is reasonable since the 2\ell_{2}-norm of a 56×5656\times 56 image is twice that of its 28×2828\times 28 counterpart. Put simply, fooling big MNIST is just as hard/easy as fooling the original MNIST regardless of resolution. This also shows that for big MNIST, as the expansion factor bb gets larger and ϵ\epsilon is expanded to match, the concentration bound UcU_{c} grows at exactly the same rate as the exponential term in equation 2 shrinks, and there is no net effect on fundamental susceptibility. Also note that an analogous result could be based on any image classification problem (we chose MNIST only for illustration), and any p0.p\geq 0.

We get a better picture of the fundamental limits of MNIST by considering classifiers that are hardened by adversarial training22 2 Adversarial 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.

9 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.

10 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

  • Abramowitz & Stegun (1965) Milton Abramowitz and Irene A Stegun. Handbook of mathematical functions: with formulas, graphs, and mathematical tables, volume 55. Courier Corporation, 1965.
  • Athalye & Sutskever (2017) Anish Athalye and Ilya Sutskever. Synthesizing robust adversarial examples. arXiv preprint arXiv:1707.07397, 2017.
  • Athalye et al. (2018) Anish Athalye, Nicholas Carlini, and David Wagner. Obfuscated gradients give a false sense of security: Circumventing defenses to adversarial examples. arXiv preprint arXiv:1802.00420, 2018.
  • Biggio et al. (2013) Battista Biggio, Igino Corona, Davide Maiorca, Blaine Nelson, Nedim Šrndić, Pavel Laskov, Giorgio Giacinto, and Fabio Roli. Evasion attacks against machine learning at test time. In Joint European conference on machine learning and knowledge discovery in databases, pp. 387–402. Springer, 2013.
  • Bobkov et al. (1997) Sergey G Bobkov et al. An isoperimetric inequality on the discrete cube, and an elementary proof of the isoperimetric inequality in gauss space. The Annals of Probability, 25(1):206–214, 1997.
  • Buckman et al. (2018) Jacob Buckman, Aurko Roy, Colin Raffel, and Ian Goodfellow. Thermometer encoding: One hot way to resist adversarial examples. OpenReview, 2018.
  • Carlini & Wagner (2016) Nicholas Carlini and David Wagner. Defensive distillation is not robust to adversarial examples. arXiv preprint arXiv:1607.04311, 2016.
  • Carlini & Wagner (2017a) Nicholas Carlini and David Wagner. Magnet and “efficient defenses against adversarial attacks” are not robust to adversarial examples. arXiv preprint arXiv:1711.08478, 2017a.
  • Carlini & Wagner (2017b) Nicholas Carlini and David Wagner. Towards evaluating the robustness of neural networks. In Security and Privacy (SP), 2017 IEEE Symposium on, pp. 39–57. IEEE, 2017b.
  • Das et al. (2018) Nilaksh Das, Madhuri Shanbhogue, Shang-Tse Chen, Fred Hohman, Siwei Li, Li Chen, Michael E Kounavis, and Duen Horng Chau. Shield: Fast, practical defense and vaccination for deep learning using JPEG compression. arXiv preprint arXiv:1802.06816, 2018.
  • Dhillon et al. (2018) Guneet S Dhillon, Kamyar Azizzadenesheli, Zachary C Lipton, Jeremy Bernstein, Jean Kossaifi, Aran Khanna, and Anima Anandkumar. Stochastic activation pruning for robust adversarial defense. arXiv preprint arXiv:1803.01442, 2018.
  • Fawzi et al. (2018) Alhussein Fawzi, Hamza Fawzi, and Omar Fawzi. Adversarial vulnerability for any classifier. arXiv preprint arXiv:1802.08686, 2018.
  • Gilmer et al. (2018) Justin Gilmer, Luke Metz, Fartash Faghri, Samuel S Schoenholz, Maithra Raghu, Martin Wattenberg, and Ian Goodfellow. Adversarial spheres. arXiv preprint arXiv:1801.02774, 2018.
  • Goodfellow et al. (2014) Ian J Goodfellow, Jonathon Shlens, and Christian Szegedy. Explaining and harnessing adversarial examples. arXiv preprint arXiv:1412.6572, 2014.
  • Guo et al. (2017) Chuan Guo, Mayank Rana, Moustapha Cisse, and Laurens van der Maaten. Countering adversarial images using input transformations. arXiv preprint arXiv:1711.00117, 2017.
  • Kolter & Wong (2017) J Zico Kolter and Eric Wong. Provable defenses against adversarial examples via the convex outer adversarial polytope. arXiv preprint arXiv:1711.00851, 2017.
  • Kurakin et al. (2016) Alexey Kurakin, Ian Goodfellow, and Samy Bengio. Adversarial examples in the physical world. arXiv preprint arXiv:1607.02533, 2016.
  • Ledoux (2001) Michel Ledoux. The concentration of measure phenomenon. Number 89. American Mathematical Soc., 2001.
  • Lévy & Pellegrino (1951) Paul Lévy and Franco Pellegrino. Problémes concrets d’analyse fonctionnelle. Gauthier-Villars, Paris, 1951.
  • Ma et al. (2018) Xingjun Ma, Bo Li, Yisen Wang, Sarah M Erfani, Sudanthi Wijewickrema, Michael E Houle, Grant Schoenebeck, Dawn Song, and James Bailey. Characterizing adversarial subspaces using local intrinsic dimensionality. arXiv preprint arXiv:1801.02613, 2018.
  • Madry et al. (2017) Aleksander Madry, Aleksandar Makelov, Ludwig Schmidt, Dimitris Tsipras, and Adrian Vladu. Towards deep learning models resistant to adversarial attacks. arXiv preprint arXiv:1706.06083, 2017.
  • Mahloujifar et al. (2018) Saeed Mahloujifar, Dimitrios I Diochnos, and Mohammad Mahmoody. The curse of concentration in robust learning: Evasion and poisoning attacks from concentration of measure. arXiv preprint arXiv:1809.03063, 2018.
  • McDiarmid (1989) C McDiarmid. On the method of bounded differences. London Mathematical Society Lecture Notes, 141:148–188, 1989.
  • Meng & Chen (2017) Dongyu Meng and Hao Chen. MagNet: a two-pronged defense against adversarial examples. In Proceedings of the 2017 ACM SIGSAC Conference on Computer and Communications Security, pp. 135–147. ACM, 2017.
  • Milman & Schechtman (1986) Vitali D Milman and Gideon Schechtman. Asymptotic Theory of Finite Dimensional Normed Spaces. Springer-Verlag, Berlin, Heidelberg, 1986. ISBN 0-387-16769-2.
  • Osserman et al. (1978) Robert Osserman et al. The isoperimetric inequality. Bulletin of the American Mathematical Society, 84(6):1182–1238, 1978.
  • Papernot et al. (2016a) Nicolas Papernot, Patrick McDaniel, and Ian Goodfellow. Transferability in machine learning: from phenomena to black-box attacks using adversarial samples. arXiv preprint arXiv:1605.07277, 2016a.
  • Papernot et al. (2016b) Nicolas Papernot, Patrick McDaniel, Xi Wu, Somesh Jha, and Ananthram Swami. Distillation as a defense to adversarial perturbations against deep neural networks. In Security and Privacy (SP), 2016 IEEE Symposium on, pp. 582–597. IEEE, 2016b.
  • Raghunathan et al. (2018) Aditi Raghunathan, Jacob Steinhardt, and Percy Liang. Certified defenses against adversarial examples. arXiv preprint arXiv:1801.09344, 2018.
  • Ros (2001) Antonio Ros. The isoperimetric problem. Global theory of minimal surfaces, 2:175–209, 2001.
  • Samangouei et al. (2018) Pouya Samangouei, Maya Kabkab, and Rama Chellappa. Defense-GAN: Protecting classifiers against adversarial attacks using generative models. arXiv preprint arXiv:1805.06605, 2018.
  • Shen et al. (2017) Shiwei Shen, Guoqing Jin, Ke Gao, and Yongdong Zhang. APE-GAN: Adversarial perturbation elimination with GAN. ICLR Submission, available on OpenReview, 4, 2017.
  • Simon-Gabriel et al. (2018) Carl-Johann Simon-Gabriel, Yann Ollivier, Bernhard Schölkopf, Léon Bottou, and David Lopez-Paz. Adversarial vulnerability of neural networks increases with input dimension. arXiv preprint arXiv:1802.01421, 2018.
  • Sinha et al. (2018) Aman Sinha, Hongseok Namkoong, and John Duchi. Certifying some distributional robustness with principled adversarial training. OpenReview, 2018.
  • Song et al. (2017) Yang Song, Taesup Kim, Sebastian Nowozin, Stefano Ermon, and Nate Kushman. PixelDefend: Leveraging generative models to understand and defend against adversarial examples. arXiv preprint arXiv:1710.10766, 2017.
  • Su et al. (2017) Jiawei Su, Danilo Vasconcellos Vargas, and Sakurai Kouichi. One pixel attack for fooling deep neural networks. arXiv preprint arXiv:1710.08864, 2017.
  • Sudakov & Tsirelson (1974) Vladimir Sudakov and Boris Tsirelson. Extremal properties of half-spaces for spherically invariant measures. J. Soviet Math., 41, 1974.
  • Szegedy et al. (2013) Christian Szegedy, Wojciech Zaremba, Ilya Sutskever, Joan Bruna, Dumitru Erhan, Ian Goodfellow, and Rob Fergus. Intriguing properties of neural networks. arXiv preprint arXiv:1312.6199, 2013.
  • Talagrand (1995) Michel Talagrand. Concentration of measure and isoperimetric inequalities in product spaces. Publications Mathématiques de l’Institut des Hautes Etudes Scientifiques, 81(1):73–205, 1995.
  • Talagrand (1996) Michel Talagrand. A new look at independence. The Annals of probability, pp. 1–34, 1996.
  • Tramèr et al. (2017a) Florian Tramèr, Alexey Kurakin, Nicolas Papernot, Ian Goodfellow, Dan Boneh, and Patrick McDaniel. Ensemble adversarial training: Attacks and defenses. arXiv preprint arXiv:1705.07204, 2017a.
  • Tramèr et al. (2017b) Florian Tramèr, Nicolas Papernot, Ian Goodfellow, Dan Boneh, and Patrick McDaniel. The space of transferable adversarial examples. arXiv preprint arXiv:1704.03453, 2017b.
  • Vershynin (2017) Roman Vershynin. High-Dimensional Probability. Cambridge University Press, 2017.
  • Wang et al. (2018) Yizhen Wang, Somesh Jha, and Kamalika Chaudhuri. Analyzing the robustness of nearest neighbors to adversarial examples. International Conference of Machine Learning (ICML), 2018.
  • Xie et al. (2017) Cihang Xie, Jianyu Wang, Zhishuai Zhang, Zhou Ren, and Alan Yuille. Mitigating adversarial effects through randomization. arXiv preprint arXiv:1711.01991, 2017.
  • Xu et al. (2017) Weilin Xu, David Evans, and Yanjun Qi. Feature squeezing: Detecting adversarial examples in deep neural networks. arXiv preprint arXiv:1704.01155, 2017.
  • Zantedeschi et al. (2017) Valentina Zantedeschi, Maria-Irina Nicolae, and Ambrish Rawat. Efficient defenses against adversarial attacks. In Proceedings of the 10th ACM Workshop on Artificial Intelligence and Security, pp. 39–49. ACM, 2017.

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.

In the lemma below, we consider the standard Gaussian density in n\mathbb{R}^{n} given by p(x)=1(2π)n/2enx2/2p(x)=\frac{1}{(2\pi)^{n/2}}e^{-nx^{2}/2} and corresponding Gaussian measure μ.\mu. We also define

Φ(z)=12πzet2/2dt,\Phi(z)=\frac{1}{\sqrt{2\pi}}\int_{-\infty}^{z}e^{-t^{2}/2}dt,

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.

Lemma 5 (Gaussian Isoperimetric Inequality).

Of all sets with the same Gaussian measure, the set with 2\ell_{2} ϵ\epsilon-expansion of smallest measure is a half space. Furthermore, for any measurable set 𝒜n,\mathcal{A}\subset\mathbb{R}^{n}, and scalar constant aa such that Φ(a)=μ[𝒜],\Phi(a)=\mu[\mathcal{A}],

μ[𝒜(ϵ,d2)]Φ(a+ϵ).\mu[\mathcal{A}(\epsilon,d_{2})]\geq\Phi(a+\epsilon).

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

This function Φ\Phi maps a random Guassian vector zN(0,I)z\in N(0,I) onto a random uniform vector in the unit cube. To see why, consider a measurable subset Rn.\mathcal{B}\subset R^{n}. If μ\mu is the Gaussian measure on n\mathbb{R}^{n} and σ\sigma is the uniform measure on the cube, then

σ[Φ()]=χΦ()(s)𝑑σ=χΦ()(Φ(z))1det(JΦ)𝑑z=χ(z)𝑑μ=μ[].\sigma[\Phi(\mathcal{B})]=\int\chi_{\Phi(\mathcal{B})}(s)d\sigma=\int\chi_{\Phi(\mathcal{B})}(\Phi(z))\frac{1}{\det(J\Phi)}dz=\int\chi_{\mathcal{B}}(z)d\mu=\mu[\mathcal{B}].

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

Φ(z)Φ(w)p12π(xw)p=12πzwp\|\Phi(z)-\Phi(w)\|_{p}\leq\|\frac{1}{\sqrt{2\pi}}(x-w)\|_{p}=\frac{1}{\sqrt{2\pi}}\|z-w\|_{p}

for any z,wn.z,w\in\mathbb{R}^{n}. From this, we see that for p=min(p,2)p^{*}=\min(p,2)

Φ(z)Φ(w)pn1/p1/2Φ(z)Φ(w)2n1/p2πnzw2\displaystyle\|\Phi(z)-\Phi(w)\|_{p}\leq n^{1/p^{*}-1/2}\|\Phi(z)-\Phi(w)\|_{2}\leq\frac{n^{1/p^{*}}}{\sqrt{2\pi n}}\|z-w\|_{2} (10)

where we have used the identity upn1/min(p,2)1/2u2.\|u\|_{p}\leq n^{1/\min(p,2)-1/2}\|u\|_{2}.

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

Φ(2πnn1/pϵ,d2)𝒜(ϵ,dp).\Phi\mathcal{B}\left(\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon,d_{2}\right)\subset\mathcal{A}(\epsilon,d_{p}).

It follows from equation 10 that

σ[𝒜(ϵ,dp)]μ[(2πnn1/pϵ,d2)].\sigma[\mathcal{A}(\epsilon,d_{p})]\geq\mu\left[\mathcal{B}\left(\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon,d_{2}\right)\right].

Applying Lemma 5, we see that

σ[𝒜(ϵ,dp)]Φ(α+2πnn1/pϵ)\displaystyle\sigma[\mathcal{A}(\epsilon,d_{p})]\geq\Phi\left(\alpha+\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon\right) (11)

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

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

12πxet2𝑑t<ex22πx\frac{1}{\sqrt{2\pi}}\int_{x}^{\infty}e^{-t^{2}}dt<\frac{e^{-x^{2}}}{\sqrt{2\pi}x}

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

Appendix B Proof of Lemma 4

Our proof emulates the method of Talagrand, with modifications that extend the result to other p\ell_{p} norms. We need the following standard inequality. Proof can be found in Talagrand 1995; Talagrand 1996.

Lemma 6 (Talagrand).

Consider a probability space Ω\Omega with measure μ.\mu. For g:Ω[0,1],g:\Omega\to[0,1], we have

Ωmin(et,1gα)𝑑μ×(Ωg𝑑μ)αexp(t2(α+1)8α).\int_{\Omega}\min\left(e^{t},\frac{1}{g^{\alpha}}\right)d\mu\,\times\,\left(\int_{\Omega}gd\mu\right)^{\alpha}\leq\exp\left(\frac{t^{2}(\alpha+1)}{8\alpha}\right).

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

etf(x,𝒜)𝑑x1σα[𝒜]exp(nt2(α+1)8α)\displaystyle\int e^{tf(x,\mathcal{A})}dx\leq\frac{1}{\sigma^{\alpha}[\mathcal{A}]}\exp\left(\frac{nt^{2}(\alpha+1)}{8\alpha}\right) (12)

where f(x,𝒜)=miniy𝒜|xiyi|p=dpp(x,𝒜)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 𝒜\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

etf(x,𝒜)𝑑xσ[𝒜]+𝒜cetf(x,𝒜)𝑑xσ[𝒜]+𝒜c1𝑑xσ[𝒜]+(1σ[𝒜])et1σα[𝒜]exp(t2(α+1)8α).\int e^{tf(x,\mathcal{A})}dx\leq\sigma[\mathcal{A}]+\int_{\mathcal{A}^{c}}e^{{tf(x,\mathcal{A})}}dx\leq\sigma[\mathcal{A}]+\int_{\mathcal{A}^{c}}1dx\leq\sigma[\mathcal{A}]+(1-\sigma[\mathcal{A}])e^{t}\leq\frac{1}{\sigma^{\alpha}[\mathcal{A}]}\exp\left(\frac{t^{2}(\alpha+1)}{8\alpha}\right).

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 𝒜[0,1]n.\mathcal{A}\subset[0,1]^{n}. Define

𝒜ω={zn1|(ω,z)𝒜} and ={zn1|(ω,z)𝒜 for some ω}.\mathcal{A}_{\omega}=\{z\in\mathbb{R}^{n-1}\,|\,(\omega,z)\in\mathcal{A}\}\text{ and }\mathcal{B}=\{z\in\mathbb{R}^{n-1}\,|\,(\omega,z)\in\mathcal{A}\text{ for some }\omega\}.

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

etf(x,𝒜)𝑑xω[0,1]z[0,1]n1etf(z,𝒜ω)𝑑z𝑑xω[0,1]1σα[𝒜ω]exp((n1)t2(α+1)8α).\int e^{tf(x,\mathcal{A})}dx\leq\int_{\omega\in[0,1]}\int_{z\in[0,1]^{n-1}}e^{tf(z,\mathcal{A}_{\omega})}dz\,dx\leq\int_{\omega\in[0,1]}\frac{1}{\sigma^{\alpha}[\mathcal{A}_{\omega}]}\exp\left(\frac{(n-1)t^{2}(\alpha+1)}{8\alpha}\right).

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

etf(x,𝒜)𝑑x(ω,z)[0,1]net(f(x,)+1)et(ω,z)[0,1]netf(x,)etσα[]exp((n1)t2(α+1)8α).\int e^{tf(x,\mathcal{A})}dx\leq\int_{(\omega,z)\in[0,1]^{n}}e^{t(f(x,\mathcal{B})+1)}\leq e^{t}\int_{(\omega,z)\in[0,1]^{n}}e^{tf(x,\mathcal{B})}\leq\frac{e^{t}}{\sigma^{\alpha}[\mathcal{B}]}\exp\left(\frac{(n-1)t^{2}(\alpha+1)}{8\alpha}\right).

Applying equation 6 gives us

etf(x,𝒜)𝑑x\displaystyle\int e^{tf(x,\mathcal{A})}dx ω[0,1]min(etσα[]exp((n1)t2(α+1)8α),1σα[𝒜ω]exp((n1)t2(α+1)8α))\displaystyle\leq\int_{\omega\in[0,1]}\min\left(\frac{e^{t}}{\sigma^{\alpha}[\mathcal{B}]}\exp\left(\frac{(n-1)t^{2}(\alpha+1)}{8\alpha}\right),\frac{1}{\sigma^{\alpha}[\mathcal{A}_{\omega}]}\exp\left(\frac{(n-1)t^{2}(\alpha+1)}{8\alpha}\right)\right)
=exp((n1)t2(α+1)8α)1σα[]ω[0,1]min(et,σα[]σα[𝒜ω]).\displaystyle=\exp\left(\frac{(n-1)t^{2}(\alpha+1)}{8\alpha}\right)\frac{1}{\sigma^{\alpha}[\mathcal{B}]}\int_{\omega\in[0,1]}\min\left(e^{t},\frac{\sigma^{\alpha}[\mathcal{B}]}{\sigma^{\alpha}[\mathcal{A}_{\omega}]}\right). (13)

Now, we apply lemma 6 to equation 13 with g(ω)=α[𝒜ω]/α[]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

1σ[𝒜(ϵ,dp)]=σ{x|f(x)>ϵp}etf(x,A)𝑑xetϵpexp(nt2(α+1)8α)σα[𝒜]etϵp.\displaystyle 1-\sigma[\mathcal{A}(\epsilon,d_{p})]=\sigma\{x\,|\,f(x)>\epsilon^{p}\}\leq\frac{\int e^{tf(x,A)}dx}{e^{t\epsilon^{p}}}\leq\frac{\exp\left(\frac{nt^{2}(\alpha+1)}{8\alpha}\right)}{\sigma^{\alpha}[\mathcal{A}]e^{t\epsilon^{p}}}. (14)

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

1σ[𝒜(ϵ,dp)]exp(2αϵ2pn(α+1))σα[𝒜].\displaystyle 1-\sigma[\mathcal{A}(\epsilon,d_{p})]\leq\frac{\exp\left(-\frac{2\alpha\epsilon^{2p}}{n(\alpha+1)}\right)}{\sigma^{\alpha}[\mathcal{A}]}. (15)

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

1σ[𝒜(ϵ,dp)]exp(ϵ2p/n)σ[𝒜],\displaystyle 1-\sigma[\mathcal{A}(\epsilon,d_{p})]\leq\frac{\exp\left(-\epsilon^{2p}/n\right)}{\sigma[\mathcal{A}]}, (16)

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 ϵ2pn2log(1/σ(𝒜)).\epsilon^{2p}\geq\frac{n}{2}\log(1/\sigma(\mathcal{A})). We arrive at

1σ[𝒜(ϵ,dp)]exp(2n(ϵpnlog(σ1[𝒜])/2)2).\displaystyle 1-\sigma[\mathcal{A}(\epsilon,d_{p})]\leq\exp\left(-\frac{2}{n}\left(\epsilon^{p}-\sqrt{n\log(\sigma^{-1}[\mathcal{A}])/2}\right)^{2}\right). (17)

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.

Choose a class cc with fc12f_{c}\leq\frac{1}{2} and let ={x|𝒞(x)=c}\mathcal{R}=\{x|\mathcal{C}(x)=c\} denote the subset of the cube lying in class cc according to the classifier 𝒞\mathcal{C}. Let ¯\overline{\mathcal{R}} be the complement, who’s p\ell_{p} expansion is denoted ¯(ϵ,dp)\overline{\mathcal{R}}(\epsilon;d_{p}). Because ¯\overline{\mathcal{R}} covers at least half the cube, we can invoke Lemma 3. We have that

vol[¯(ϵ,h)]1δ,\vol[\overline{\mathcal{R}}(\epsilon;h)]\geq 1-\delta,

where

δ={exp(2πn12/pϵ2)2πϵn1/21/p, for Theorem 2 and 2Ucexp(ϵ2/n), for Theorem 3.\displaystyle\delta=\begin{cases}\frac{\exp(-2\pi n^{1-2/p^{*}}\epsilon^{2})}{2\pi\epsilon n^{1/2-1/p^{*}}}\text{, for Theorem \ref{cubeTheorem} and }\\ 2U_{c}\exp(-\epsilon^{2}/n),\text{ for Theorem \ref{sparseTheorem}}.\end{cases} (18)

The set ¯(ϵ,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 ¯c.\overline{\mathcal{R}}_{c}.

Appendix D Proof of Theorem 4

Let 𝒜\mathcal{A} denote the support of pc,p_{c}, and suppose that this support has measure vol[𝒜]=η.\vol[\mathcal{A}]=\eta. We want to show that, for large enough ϵ,\epsilon, the expansion 𝒜(ϵ,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 𝒜(ϵ,dp)\mathcal{A}(\epsilon,d_{p}) overlaps with other classes, and so there must be data points in 𝒜\mathcal{A} with ϵ\epsilon-adversarial examples.

We start with the case p>0,p>0, where we bound 𝒜(ϵ,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

Φ(α)=12παet2/2dt12eα2/2,\Phi(\alpha)=\frac{1}{2\pi}\int_{-\infty}^{\alpha}e^{-t^{2}/2}dt\leq\frac{1}{2}e^{-\alpha^{2}/2},

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

αlog14Φ(α)2.\displaystyle\alpha\geq-\sqrt{\log\frac{1}{4\Phi(\alpha)^{2}}}. (19)

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

vol[𝒜(ϵ,dp)]Φ(α+ϵ)Φ(log14η2+2πnn1/pϵ).\vol[\mathcal{A}(\epsilon,d_{p})]\geq\Phi(\alpha+\epsilon)\geq\Phi\left(-\sqrt{\log\frac{1}{4\eta^{2}}}+\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon\right).

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

2πnn1/pϵ>log14η2.\frac{\sqrt{2\pi n}}{n^{1/p^{*}}}\epsilon>\sqrt{\log\frac{1}{4\eta^{2}}}.

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

vol[𝒜(ϵ,d0)]1exp(2n(ϵnlog(1/η)/2)2).\vol[\mathcal{A}(\epsilon,d_{0})]\geq 1-\exp\left(-\frac{2}{n}\left(\epsilon-\sqrt{n\log(1/\eta)/2}\right)^{2}\right).

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

exp(2n(ϵnlog(1/η)/2)2)<12.\exp\left(-\frac{2}{n}\left(\epsilon-\sqrt{n\log(1/\eta)/2}\right)^{2}\right)<\frac{1}{2}.

which holds if

η>exp(2(ϵnlog2/2)2n).\eta>\exp\left(-\frac{2\left(\epsilon-\sqrt{n\log 2/2}\right)^{2}}{n}\right).

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 b1b\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 b1b\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.