Title: ROBUSTNESS VERIFICATION FOR TRANSFORMERS

Abstract: Robustness verification that aims to formally certify the prediction behavior of neural networks has become an important tool for understanding model behavior and obtaining safety guarantees. However, previous methods can usually only handle neural networks with relatively simple architectures. In this paper, we consider the robustness verification problem for Transformers. Transformers have complex self-attention layers that pose many challenges for verification, including cross-nonlinearity and cross-position dependency, which have not been discussed in previous works. We resolve these challenges and develop the first robustness verification algorithm for Transformers. The certified robustness bounds computed by our method are significantly tighter than those by naive Interval Bound Propagation. These bounds also shed light on interpreting Transformers as they consistently reflect the importance of different words in sentiment analysis.

1. INTRODUCTION:

Deep neural networks have been successfully applied to many domains. However, these blackbox models are generally difficult to analyze and their behavior is not guaranteed. Moreover, it has been shown that the predictions of deep networks become unreliable and unstable when tested in unseen situations, e.g., in the presence of small adversarial perturbations to the input (Szegedy et al., 2013; Goodfellow et al., 2014; Lin et al., 2019) . Therefore, neural network verification has become an important tool for analyzing and understanding the behavior of neural networks, with applications in safety-critical applications (Katz et al., 2017; Julian et al., 2019; Lin et al., 2019) , model explanation (Shih et al., 2018) and robustness analysis (Tjeng et al., 2019; Wang et al., 2018c; Gehr et al., 2018; Wong & Kolter, 2018; Singh et al., 2018; Weng et al., 2018; Zhang et al., 2018) .

Formally, a neural network verification algorithm aims to provably characterize the prediction of a network within some input space. For example, given a K-way classification model f : R d → R K , where f i (x) stands for the predicted score of class i, we can verify some linear specification (defined by a vector c) as below:

EQUATION

where S is a predefined input space. In the robustness verification problem, S = {x | x-x 0 p ≤ } is defined as some small p -ball around the original example x 0 , and setting up c = 1 y0 -1 y enables us to verify whether the logit output of class y 0 is always greater than another class y for any input within S. This is a nonconvex optimization problem which makes computing the exact solution challenging, and thus several algorithms are recently proposed to find the lower bounds of Eq. (1) in order to efficiently obtain a safety guarantee (Gehr et al., 2018; Weng et al., 2018; Zhang et al., 2018; Singh et al., 2019) . Moreover, extensions of these algorithms can be used for verifying properties beyond robustness, such as rotation or shift invariant (Singh et al., 2019) , conservation of energy (Qin et al., 2019) and model correctness (Yang & Rinard, 2019) .

However, most of existing verification methods focus on relatively simple neural network architectures, such as feed-forward and recurrent neural networks, while they cannot handle complex structures. In this paper, we develop the first robustness verification algorithm for Transformers (Vaswani et al., 2017) with self-attention layers. Transformers have been widely used in natural language processing (Devlin et al., 2019; Yang et al., 2019; Liu et al., 2019) and many other domains (Parmar et al., 2018; Kang & McAuley, 2018; Li et al., 2019b; Su et al., 2019; Li et al., 2019a) . For frames under perturbation in the input sequence, we aim to compute a lower bound such that when these frames are perturbed within p -balls centered at the original frames respectively and with a radius of , the model prediction is certified to be unchanged. To compute such bounds efficiently, we adopt the linear-relaxation framework (Weng et al., 2018; Zhang et al., 2018) -we recursively propagate and compute linear lower and upper bounds for each neuron w.r.t the input within the perturbation space S.

We resolve several particular challenges in verifying Transformers. First, Transformers with selfattention layers have a complicated architecture. Unlike simpler networks, they cannot be written as multiple layers of affine transformations or element-wise activation functions. Therefore, we need to propagate linear bounds differently for self-attention layers. Second, dot products, softmax, and weighted summation in self-attention layers involve multiplication or division of two variables both under perturbation, namely cross-nonlinearity, which is not present in feed-forward networks. Ko et al. (2019) proposed a gradient descent based approach to find linear bounds, however it is inefficient and poses a computational challenge for Transformer verification since self-attention is the core of Transformers. In contrast, we derive closed-form linear bounds that can be computed in O(1) complexity. Third, in the computation of self-attention, output neurons in each position depend on all input neurons from different positions (namely cross-position dependency), unlike the case in recurrent neural networks where outputs depend on only the hidden features from the previous position and the current input. Previous works (Zhang et al., 2018; Weng et al., 2018; Ko et al., 2019) have to track all such dependency and thus is costly in time and memory. To tackle this, we introduce an efficient bound propagating process in a forward manner specially for self-attention layers, enabling the tighter backward bounding process for other layers to utilize bounds computed by the forward process. In this way, we avoid cross-position dependency in the backward process which is relatively slower but produces tighter bounds. Combined with the forward process, the complexity of the backward process is reduced by O(n) for input length n, while the computed bounds remain comparably tight. Our contributions are summarized below:

• We propose an effective and efficient algorithm for verifying the robustness of Transformers with self-attention layers. To our best knowledge, this is the first method for verifying Transformers.

• We resolve key challenges in verifying Transformers, including cross-nonlinearity and cross-position dependency. Our bounds are significantly tighter than those by adapting Interval Bound Propagation (IBP) (Mirman et al., 2018; Gowal et al., 2018) . • We quantitatively and qualitatively show that the certified bounds computed by our algorithm consistently reflect the importance of input words in sentiment analysis, which justifies that these bounds are meaningful in practice and they shed light on interpreting Transformers.

2. RELATED WORK:

Robustness Verification for Neural Networks. Given an input x 0 and a small region B p (x 0 , ) := {x | xx 0 p ≤ }, the goal of robustness verification is to verify whether the prediction of the neural network is unchanged within this region. This problem can be mathematically formulated as Eq. ( 1). If Eq. ( 1) can be solved optimally, then we can derive the minimum adversarial perturbation of x by conducting binary search on . Equivalently, we obtain the maximum such that any perturbation within B p (x 0 , ) cannot change the predicted label.

Several works focus on solving Eq. (1) exactly and optimally, using mixed integer linear programming (MILP) (Tjeng et al., 2019; Dutta et al., 2018) , branch and bound (BaB) (Bunel et al., 2018) , and satisfiability modulo theory (SMT) (Ehlers, 2017; Katz et al., 2017) . Unfortunately, due to the nonconvexity of model f , solving Eq. ( 1) is NP-hard even for a simple ReLU network (Katz et al., 2017) . Therefore, we can only expect to compute a lower bound of Eq. (1) efficiently by using relaxations. Many algorithms can be seen as using convex relaxations for non-linear activation functions (Salman et al., 2019) , including using duality (Wong & Kolter, 2018; Dvijotham et al., 2018) , abstract domains (Gehr et al., 2018; Singh et al., 2018; Mirman et al., 2018; Singh et al., 2019) , layer-by-layer reachability analysis (Wang et al., 2018b; Weng et al., 2018; Zhang et al., 2018; Gowal et al., 2018) and semi-definite relaxations (Raghunathan et al., 2018; Dvijotham et al., 2019) .

Additionally, robustness verification can rely on analysis on local Lipschitz constants (Hein & Andriushchenko, 2017; Zhang et al., 2019) . However, existing methods are mostly limited to verifying networks with relatively simple architectures, such as feed-forward networks and RNNs (Wang et al., 2018a; Akintunde et al., 2019; Ko et al., 2019) , while none of them are able to handle Transformers.

Transformers and Self-Attentive Models. Transformers (Vaswani et al., 2017) based on selfattention mechanism, further with pre-training on large-scale corpora, such as BERT (Devlin et al., 2019) , XLNet (Yang et al., 2019) , RoBERTa (Liu et al., 2019) , achieved state-of-the-art performance on many NLP tasks. Self-attentive models are also useful beyond NLP, including VisualBERT on vision and language applications (Li et al., 2019b; Su et al., 2019) , image transformer for image generation (Parmar et al., 2018) , acoustic models for speech recognition (Zhou et al., 2018) , sequential recommendation (Kang & McAuley, 2018) and graph embedding (Li et al., 2019a) .

The robustness of NLP models has been studied, especially many methods have been proposed to generate adversarial examples (Papernot et al., 2016; Jia & Liang, 2017; Zhao et al., 2017; Alzantot et al., 2018; Cheng et al., 2018; Ebrahimi et al., 2018; Shi et al., 2019) . In particular, Hsieh et al. (2019) showed that Transformers are more robust than LSTMs. However, there is not much work on robustness verification for NLP models. Ko et al. (2019) 

3. METHODOLOGY:

We aim to verify the robustness of a Transformer whose input is a sequence of frames

X = [x (1) , x (2) , • • • , x (n) ].

We take binary text classification as a running example, where x (i) is a word embedding and the model outputs a score y c (X) for each class c (c ∈ {0, 1}). Nevertheless, our method for verifying Transformers is general and can also be applied in other applications.

For a clean input sequence

EQUATION

0 ] correctly classified by the model, let P = {r 1 , r 2 , • • • , r t }(1 ≤ r k ≤ n) be the set of perturbed positions, where t is the number of perturbed positions. Thus the perturbed input belongs to

S := {X = [x (1) , x (2) , • • • , x (n) ] : x (r k ) - x (r k ) 0 p ≤ , 1 ≤ k ≤ t, x (i) = x (i) 0 , ∀i / ∈ P }.

Assuming that c is the gold class, the goal of robustness verification is to compute

min X∈S y c (X) -y 1-c (X) := δ .

If δ > 0, the output score of the correct class is always larger than the incorrect one for any input within S . As mentioned previously, computing the exact values of δ is NP-hard, and thus our goal is to efficiently compute a lower bound δ L ≤ δ .

3.1. BASE FRAMEWORK:

We obtain δ L (X) by computing the bounds of each neuron when X is perturbed within S (δ L can be regarded as a final neuron). A Transformer layer can be decomposed into a number of sub-layers, where each sub-layer contains neurons after some operations. These operations can be categorized into three categories: 1) linear transformations, 2) unary nonlinear functions, and 3) operations in self-attention. Each sub-layer contains n positions in the sequence and each position contains a group of neurons. We assume that the Transformer we verify has m sub-layers in total, and the value of the j-th neuron at the i-th position in the l-th sub-layer is Φ (l,i) j (X), where Φ (l,i) (X) is a vector for the specified sub-layer and position. Specially, Φ (0,i) = x (i) taking l = 0. We aim to compute a global lower bound f (l,i),L j and a global upper bound f

(l,i),U j of Φ (l,i) j (X) for X ∈ S .

We compute bounds from the first sub-layer to the last sub-layer. For neurons in the l-th layer, we aim to represent their bounds as linear functions of neurons in a previous layer, the l -th layer:

EQUATION

where Λ (l,i,l ,k),L , ∆ (l,i,l ),L and Λ (l,i,l ,k),U , ∆ (l,i,l ),U are parameters of linear lower and upper bounds respectively. Using linear bounds enables us to efficiently compute bounds with a reasonable tightness. We initially have Λ (l,i,l,i),L = Λ (l,i,l,i),U = I and ∆ (l,i,l),L = ∆ (l,i,l),U = 0. Thereby the right-hand-side of Eq. ( 2) equals to Φ (l,i) j (X) when l = l. Generally, we use a backward process to propagate the bounds to previous sub-layers, by substituting Φ (l ,i) with linear functions of previous neurons. It can be recursively conducted until the input layer l = 0. Since Φ (0,k) = x (k) = x (k) 0 (∀k / ∈ P ) is constant, we can regard the bounds as linear functions of the perturbed embeddings Φ (0,r k ) = x (r k ) (1 ≤ k ≤ t), and take the global bounds for x (r k ) ∈ B p (x

EQUATION

f (l,i),U j = t k=1 Λ (l,i,0,r k ),U j,: q + n k=1 Λ (l,i,0,k),U j,:

x

EQUATION

where 1/p + 1/q = 1 with p, q ≥ 1. These steps resemble to CROWN (Zhang et al., 2018) which is proposed to verify feed-forward networks. We further support verifying self-attentive Transformers which are more complex than feed-forward networks. Moreover, unlike CROWN that conducts a fully backward process, we combine the backward process with a forward process (see Sec. 3.3) to reduce the computational complexity of verifying Transformers.

3.2. LINEAR TRANSFORMATIONS AND UNARY NONLINEAR FUNCTIONS:

Linear transformations and unary nonlinear functions are basic operations in neural networks. We show how bounds Eq. ( 2) at the l -th sub-layer are propagated to the (l -1)-th layer.

Linear Transformations If the l -th sub-layer is connected with the (l -1)-th sub-layer with a linear transformation l ) where W (l ) , b (l ) are parameters of the linear transformation, we propagate the bounds to the (l -1)-th layer by substituting Φ (l ,k) (X):

Φ (l ,k) (X) = W (l ) Φ (l -1,k) (X) + b (

Λ (l,i,l -1,k),L/U = Λ (l,i,l ,k),L/U W (l ) , ∆ (l,i,l -1),L/U = ∆ (l,i,l ),L/U + n k=1 Λ (l,i,l ,k),L/U b (l ) ,

where "L/U " means that the equations hold for both lower bounds and upper bounds respectively.

Unary Nonlinear Functions If the l -th layer is obtained from the (l -1)-th layer with an unary nonlinear function

Φ (l ,k) j (X) = σ (l ) (Φ (l -1,k) j (X)

), to propagate linear bounds over the nonlinear function, we first bound σ (l ) (Φ

(l -1,k) j (X)) with two linear functions of Φ (l -1,k) j (X): α (l ,k),L j Φ (l -1,k) j (X) + β (l ,k),L j ≤ σ (l ) (Φ (l -1,k) j (X)) ≤ α (l ,k),U j Φ (l -1,k) j (X) + β (l ,k),U j

, where α

(l ,k),L/U j , β (l ,k),L/U j

are parameters such that the inequation holds true for all Φ (l -1,k) j (X) within its bounds computed previously. Such linear relaxations can be done for different functions, respectively. We provide detailed bounds for functions involved in Transformers in Appendix B.

We then back propagate the bounds:

Λ (l,i,l -1,k),L/U :,j = α (l ,k),L/U j Λ (l,i,l ,k),L/U :,j,+ + α (l ,k),U/L j Λ (l,i,l ,k),L/U :,j,- , ∆ (l,i,l -1),L/U j = ∆ (l,i,l ),L/U j + n k=1 β (l ,k),L/U j Λ (l,i,l ,k),L/U :,j,+ + β (l ,k),U/L j Λ (l,i,l ,k),L/U :,j,-

, where Λ (l,i,l ,k),L/U :,j,+ and Λ (l,i,l ,k),L/U :,j,-mean to retain positive and negative elements in vector

Λ (l,i,l ,k),L/U :,j

respectively and set other elements to 0.

3.3. SELF-ATTENTION MECHANISM:

Self-attention layers are the most challenging parts for verifying Transformers. We assume that Φ (l-1,i) (X) is the input to a self-attention layer. We describe our method for computing bounds for one attention head, and bounds for different heads of the multi-head attention in Transformers can be easily concatenated. Φ (l-1,i) (X) is first linearly projected to queries q (l,i) (X), keys k (l,i) (X), and values v (l,i) (X) with different linear projections, and their bounds can be obtained as described in Sec. 3.2. We also keep their linear bounds that are linear functions of the perturbed embeddings. For convenience, let rt) , where ⊕ indicates vector concatenation, and thereby we represent the linear bounds as linear functions of x (r) :

x (r) = x (r1) ⊕ x (r2) ⊕ • • • x (

Ω (l,i),q/k/v,L j,: x (r) + Θ (l,i),q/k/v,L j ≤ (q/k/v) (l,i) j (X) ≤ Ω (l,i),q/k/v,U j,: x (r) + Θ (l,i),q/k/v,U j

, where q/k/v and q/k/v mean that the inequation holds true for queries, keys and values respectively. We then bound the output of the self-attention layer starting from

q (l,i) (X), k (l,i) (X), v (l,i) (X).

Bounds of Multiplications and Divisions:

We bound multiplications and divisions in the selfattention mechanism with linear functions. We aim to bound bivariate function z = xy or z = x y (y > 0) with two linear functions

z L = α L x + β L y + γ L and z U = α U x + β U y + γ U , where x ∈ [l x , u x ], y ∈ [l y , u y ] are bounds of x, y obtained previously. For z = xy, we derive optimal parameters: α L = l y , α U = u y , β L = β U = l x , γ L = -l x l y , γ U = -l x u y .

We provide a proof in Appendix C. However, directly bounding z = x y is tricky; fortunately, we can bound it indirectly by first bounding a unary function y = 1 y and then bounding the multiplication z = xy.

A Forward Process For the self-attention mechanism, instead of using the backward process like CROWN (Zhang et al., 2018) , we compute bounds with a forward process which we will show later that it can reduce the computational complexity. Attention scores are computed from q (l,i) (X) and k (l,i) (X):

S (l) i,j = (q (l,i) (X)) T k (l,j) (X) = d qk k=1 q (l,i) k (X)k (l,j) k (X)

, where d qk is the dimension of q (l,i) (X) and k (l,j) (X). For each multiplication q

(l,i) k (X)k (l,j) k (X)

, it is bounded by:

q (l,i) k (X)k (l,j) k (X) ≥ α (l,i,j),L k q (l,i) k (X) + β (l,i,j),L k k (l,j) k (X) + γ (l,i,j),L k , q (l,i) k (X)k (l,j) k (X) ≤ α (l,i,j),U k q (l,i) k (X) + β (l,i,j),U k k (l,j) k (X) + γ (l,i,j),U k .

We then obtain the bounds of S (l) i,j : Ω (l,i),s,L j,:

x (r) + Θ (l,i),s,L j ≤ S (l) i,j ≤ Ω (l,i),s,U j,: x (r) + Θ (l,i),s,U j , Ω (l,i),s,L/U j,: = α (l,i,j),L/U k ( α (l,i,j),L/U k >0 Ω (l,i),q,L/U k,: + α (l,i,j),L/U k <0 Ω (l,i),q,U/L k,: )+ β (l,i,j),L/U k ( β (l,i,j),L/U k >0 Ω (l,j),k,L/U k,: + β (l,i,j),L/U k <0 Ω (l,j),k,U/L k,:

),

Θ (l,i),s,L/U j = α (l,i,j),L/U k ( α (l,i,j),L/U k >0 Θ (l,i),q,L/U k + α (l,i,j),L/U k <0 Θ (l,i),q,U/L k )+ β (l,i,j),L/U k ( β (l,i,j),L/U k >0 Θ (l,j),k,L/U k + β (l,i,j),L/U k <0 Θ (l,j),k,U/L k ) + d qk k=1 γ (l,i,j),L/U k .

In this way, linear bounds of q (l,i) (X) and k (l,i) (X) are forward propagated to S (l) i,j . Attention scores are normalized into attention probabilities with a softmax, i.e.

S(l):

i,j = exp(S

(l) i,j )/( n k=1 exp(S (l) i,k )), where S(l) i,j is a normalized attention probability. exp(S (l) i,j

) is an unary nonlinear function and can be bounded by α

(l),L/U i,j S (l) i,j + β (l),L/U i,j

. So we forward propagate bounds of S (l) i,j to bound exp(S (l) i,j ) with Ω (l,i),e,L/U j,:

x (r) + Θ (l,i),e,L/U j , where:

Ω (l,i),e,L/U j,: = α (l),L/U i,j Ω (l,i),s,L/U j,: Θ (l,i),e,L/U j = α (l),L/U i,j Θ (l,i),s,L/U j + β (l),L/U i,j α (l),L/U i,j ≥ 0, Ω (l,i),e,L/U j,: = α (l),L/U i,j Ω (l,i),s,U/L j,: Θ (l,i),e,L/U j = α (l),L/U i,j Θ (l,i),s,U/L j + β (l),L/U i,j α (l),L/U i,j < 0.

By summing up bounds of each exp(S i,k ) ready, we forward propagate the bounds to S(l) i,j with a division similarly to bounding q

(l,i) k (X)k (l,j) k (X).

The output of the self-attention Φ (l,i) (X) is obtained with a summation of v (l,j) (X) weighted by attention probability

S(l) i,k : Φ (l,i) j (X) = n k=1 S(l) i,k v (l,k) j

(X), which can be regarded as a dot product of S(l) i and ṽ(l,j) k (X), where ṽ(l,j)

k (X) = v (l,k) j

(X) whose bounds can be obtained from those of v

(l,k) j (X)

with a transposing. Therefore, bounds of S(l) i,k and ṽ(l,j)

k (X) can be forward propagated to Φ (l,i) (X)

similarly to bounding S (l) i,j . In this way, we obtain the output bounds of the self-attention: Ω (l ,i),Φ,L j,:

x (r) + Θ (l ,i),Φ,L j ≤ Φ (l ,i) (X) ≤ Ω (l ,i),Φ,U j,: x (r) + Θ (l ,i),Φ,U j .

(5) rt) . We can split Ω (l ,i),Φ,L/U j,:

EQUATION

into t vectors with equal dimensions, Ω (l ,i,1),Φ,L/U j,:

, Ω (l ,i,2),Φ,L/U j,:

, • • • , Ω (l ,i,t),Φ,L/U j,:

, such that Eq. ( 5) becomes

EQUATION

)

Backward Process to Self-Attention Layers When computing bounds for a later sub-layer, the lth sub-layer, using the backward process, we directly propagate the bounds at the the closest previous self-attention layer assumed to be the l -th layer, to the input layer, and we skip other previous sublayers. The bounds propagated to the l -th layer are as Eq. ( 2). We substitute Φ (l ,k) (X) with linear bounds in Eq. ( 6):

Λ (l,i,0,r k ),L/U j,: = n k =1 (Λ (l,i,l ,k ),L/U j,:,+ Ω (l ,k ,k),Φ,L/U j,: +Λ (l,i,l ,k ),L/U j,:,- Ω (l ,k ,k),Φ,U/L j,: )(1 ≤ k ≤ t), Λ (l,i,0,k),L/U j,: = 0 (∀k / ∈ P ), ∆ (l,i,0),L/U j = ∆ (l,i,l ,L/U ) j + n k=1 Λ (l,i,l ,k),L/U j,:,+ Θ (l ,k),Φ,L/U j + Λ (l,i,l ,k),L/U j,:,- Θ (l ,k),Φ,U/L j .

We take global bounds as Eq. (3) and Eq. ( 4) to obtain the bounds of the l-th layer.

Advantageous of Combining the Backward Process with a Forward Process Introducing a forward process can significantly reduce the complexity of verifying Transformers. With the backward process only, we need to compute Λ (l,i,l ,k) and ∆ (l,i,l ) (l ≤ l), where the major cost is on Λ (l,i,l ,k) and there are O(m 2 n 2 ) such matrices to compute. The O(n 2 ) factor is from the dependency between all pairs of positions in the input and output respectively, which makes the algorithm inefficient especially when the input sequence is long. In contrast, the forward process represents the bounds as linear functions of the perturbed positions only instead of all positions by computing Ω (l,i) and Θ (l,i) . Imperceptible adversarial examples may not have many perturbed positions (Gao et al., 2018; Ko et al., 2019) , and thus we may assume that the number of perturbed positions, t, is small. The major cost is on Ω (l,i) while there are only O(mn) such matrices and the sizes of Λ (l,i,l ,k) and Ω (l,i) are relatively comparable for a small t. We combine the backward process and the forward process. The number of matrices Ω in the forward process is O(mn), and for the backward process, since we do not propagate bounds over self-attention layers and there is no crossposition dependency in other sub-layers, we only compute Λ (l,i,l ,k) such that i = k, and thus the number of matrices Λ is reduced to O(m 2 n). Therefore, the total number of matrices Λ and Ω we compute is O(m 2 n) and is O(n) times smaller than O(m 2 n 2 ) when only the backward process is used. Moreover, the backward process makes bounds tighter compared to solely the forward one, as we explain in Appendix D.

4. EXPERIMENTS:

To demonstrate the effectiveness of our algorithm, we compute certified bounds for several sentiment classification models and perform an ablation study to show the advantage of combining the backward and forward processes. We also demonstrate the meaningfulness of our certified bounds with an application on identifying important words.

4.1. DATASETS AND MODELS:

We use two datasets: Yelp (Zhang et al., 2015) and SST (Socher et al., 2013) Table 1 : Clean accuracies and computed bounds for 1-position perturbation. Bounds include upper bounds (obtained by an enumeration based method), certified lower bounds by IBP and our method respectively. We also report the gap between upper bounds and our lower bounds (represented as the percentage of lower bounds relative to upper bounds). We compute bounds for each possible option of perturbed positions and report the minimum ("Min") and average ("Avg") among them. We compute certified lower bounds for different models on different datasets. We include 1-position perturbation constrained by 1 / 2 / ∞ -norms and 2-position perturbation constrained by 2 -norm. We compare our lower bounds with those computed by the Interval Bound Propagation (IBP) (Gowal et al., 2018) baseline. For 1-position perturbation, we also compare with upper bounds computed by enumerating all the words in the vocabulary and finding the word closest to the original one such that the word substitution alters the predicted label. This method has an exponential complexity with respect to the vocabulary size and can hardly be extended to perturbations on 2 or more positions; thus we do not include upper bounds for 2-position perturbation. For each example, we enumerate possible options of perturbed positions (there are n t options), and we integrate results from different options by taking the minimum or average respectively. We report the average results on 10 correctly classified random test examples with sentence lengths no more than 32 for 1-position perturbation and 16 for 2-position perturbation. Table 1 and Table 2 present the results for 1-position and 2position perturbation respectively. Our certified lower bounds are significantly larger and thus tighter than those by IBP. For 1-position perturbation, the lower bounds are consistently smaller than the upper bounds, and the gap between the upper bounds and our lower bounds is reasonable compared with that in previous work on verification of feed-forward networks, e.g., in (Weng et al., 2018; Zhang et al., 2018) the upper bounds are in the order of 10 times larger than lower bounds. This demonstrates that our proposed method can compute robustness bounds for Transformers in a similar quality to the bounds of simpler neural networks. In the following, we show the effectiveness of combining the backward process with a forward process. We compare our proposed method (Backward & Forward) with two variations: 1) Fully-Forward propagates bounds in a forward manner for all sub-layers besides self-attention layers; 2) Fully-Backward computes bounds for all sub-layers including self-attention layers using the backward bound propagation and without the forward process. We compare the tightness of bounds and computation time of the three methods. We use smaller models with the hidden sizes reduced by 75%, and we use 1-position perturbation only, to accommodate Fully-Backward with large computational cost. Experiments are conducted on an NVIDIA TITAN X GPU. Table 3 presents the results. Bounds by Fully-Forward are significantly looser while those by Fully-Backward and Backward & Forward are comparable. Meanwhile, the computation time of Backward & Forward is significantly shorter than that of Fully-Backward. This demonstrates that our method of combining the backward and forward processes can compute comparably tight bounds much more efficiently.

4.3. EFFECTIVENESS OF COMBINING THE BACKWARD PROCESS WITH:



4.4. IDENTIFYING WORDS IMPORTANT TO PREDICTION:

The certified lower bounds can reflect how sensitive a model is to the perturbation of each input word. Intuitively, if a word is more important to the prediction, the model is more sensitive to its perturbation. Therefore, the certified lower bounds can be used to identify important words. In the following, we conduct an experiment to verify whether important words can be identified by our certified lower bounds. We use a 1-layer Transformer classifier under 1-position perturbation constrained by 2 -norm. The certified lower bounds are normalized by the norm of the unperturbed word embeddings respectively, when they are used for identifying the most/least important words. We compare our method with two baselines that also estimate local vulnerability: 1) Upper uses upper bounds; 2) Gradient identifies the word whose embedding has the largest 2 -norm of gradients as the most important and vice versa. Quantitative Analysis on SST SST contains sentiment labels for all phrases on parse trees, where the labels range from very negative (0) to very positive (4), and 2 for neutral. For each word, assuming its label is x, we take |x -2|, i.e., the distance to the neutral label, as the importance score, since less neutral words tend to be more important for the sentiment polarity of the sentence. We evaluate on 100 random test input sentences and compute the average importance scores of the most or least important words identified from the examples. In Table 4 , compared to the baselines ("Upper" and "Grad"), the average importance score of the most important words identified by our lower bounds are the largest, while the least important words identified by our method have the smallest average score. This demonstrates that our method identifies the most and least important words more accurately compared to baseline methods.

Qualitative Analysis on Yelp:

We further analyze the results on a larger dataset, Yelp. Since Yelp does not provide per-word sentiment labels, importance scores cannot be computed as on SST. Thus, we demonstrate a qualitative analysis. We use 10 random test examples and collect the words identified as the most and least important word in each example. In Table 4 , most words identified as the most important by certified lower bounds are exactly the words reflecting sentiment polarities (boldfaced words), while those identified as the least important words are mostly stopwords. Baseline methods mistakenly identify more words containing no sentiment polarity as the most important. This again demonstrates that our certified lower bounds identify word importance better than baselines and our bounds provide meaningful interpretations in practice. While gradients evaluate the sensitivity of each input word, this evaluation only holds true within a very small neighborhood (where the classifier can be approximated by a first-order Taylor expansion) around the input sentence. Our certified method gives valid lower bounds that hold true within a large neighborhood specified by a perturbation set S, and thus it provides more accurate results.

5. CONCLUSION:

We propose the first robustness verification method for Transformers, and tackle key challenges in verifying Transformers, including cross-nonlinearity and cross-position dependency. Our method computes certified lower bounds that are significantly tighter than those by IBP. Quantitative and qualitative analyses further show that our bounds are meaningful and can reflect the importance of different words in sentiment analysis.

Zhengli Zhao, Dheeru Dua, and Sameer Singh. Generating natural adversarial examples. arXiv preprint arXiv:1710.11342, 2017.

Shiyu Zhou, Linhao Dong, Shuang Xu, and Bo Xu. Syllable-based sequence-to-sequence speech recognition with the transformer in mandarin chinese. arXiv preprint arXiv:1804.10752, 2018.

B LINEAR BOUNDS OF UNARY NONLINEAR FUNCTIONS:

We show in Sec. 3.2 that linear bounds can be propagated over unary nonlinear functions as long as the unary nonlinear functions can be bounded with linear functions. Such bounds are determined for each neuron respectively, according to the bounds of the input for the function. Specifically, for a unary nonlinear function σ(x), with the bounds of x obtained previously as x ∈ [l, u], we aim to derive a linear lower bound α L x + β L and a linear upper bound α U x + β U , such that

α L x + β L ≤ σ(x) ≤ α U x + β U (∀x ∈ [l, u]),

where parameters α L , β L , α U , β U are dependent on l, u and designed for different functions σ(x) respectively. We introduce how the parameters are determined for different unary nonlinear functions involved in Transformers such that the linear bounds are valid and as tight as possible. Bounds of ReLU and tanh has been discussed by Zhang et al. (2018) , and we further derive bounds of e x , 1

x , x 2 , √ x. x 2 and √ x are only used when the layer normalization is not modified for experiments to study the impact of our modification. For the following description, we define the endpoints of the function to be bounded within range (l, r) as (l, σ(l)) and (u, σ(u)). We describe how the lines corresponding to the linear bounds of different functions can be determined, and thereby parameters α L , β L , α U , β U can be determined accordingly.

ReLU For ReLU activation, σ(x) = max(x, 0). ReLU is inherently linear on segments (-∞, 0] and [0, ∞) respectively, so we make the linear bounds exactly σ(x) for u ≤ 0 or l ≥ 0; and for l < 0 < u, we take the line passing the two endpoints as the upper bound; and we take σ L (x) = 0 when u < |l| and σ L (x) = x when u ≥ |l| as the lower bound, to minimize the gap between the lower bound and the original function.

Tanh For tanh activation, σ(x) = 1-e -2x

1+e -2x . tanh is concave for l ≥ 0, and thus we take the line passing the two endpoints as the lower bound and take a tangent line passing ((l+u)/2, σ((l+u)/2) as the upper bound. For u ≤ 0, tanh is convex, and thus we take the line passing the two endpoints as the upper bound and take a tangent line passing ((l + u)/2, σ((l + u)/2) as the lower bound. For l < 0 < u, we take a tangent line passing the right endpoint and (d L , σ(d L ))(d L ≤ 0) as the lower bound, and take a tangent line passing the left endpoint and (d U , σ(d U ))(d U ≥ 0) as the upper bound. d L and d U can be found with a binary search.

Exp σ(x) = exp(x) = e x is convex, and thus we take the line passing the two endpoints as the upper bound and take a tangent line passing (d, σ(d)) as the lower bound. Preferably, we take d = (l + u)/2. However, e x is always positive and used in the softmax for computing normalized attention probabilities in self-attention layers, i.e., exp(S (l) i,j ) and

n k=1 exp(S (l) i,k ). n k=1 exp(S (l) i,k )

appears in the denominator of the softmax, and to make reciprocal function 1

x finitely bounded, the range of x should not pass 0. Therefore, we impose a constraint to force the lower bound function to be always positive, i.e., σ

L (l) > 0, since σ L (l) is monotonously increasing. σ L d (x) = e d (x-d)+e d is the tangent line passing (d, σ(d)). So the constraint σ L d (l) > 0 yields d < l + 1. Hence we take d = min((l + u)/2, l + 1 -∆ d ) where ∆ d is a small real value to ensure that d < l + 1 such as ∆ d = 10 -2 .

Reciprocal For the reciprocal function, σ(x) = 1

x . It is used in the softmax and layer normalization and its input is limited to have l > 0 by the lower bounds of exp(x), and √

x. With l > 0, σ(x) is convex. Therefore, we take the line passing the two endpoints as the upper bound. And we take the tangent line passing ((l + u)/2, σ((l + u)/2)) as the lower bound.

Square For the square function, σ(x) = x 2 . It is convex and we take the line passing the two endpoints as the upper bound. And we take a tangent line passing (d, σ(d)) (d ∈ [l, u] ) as the lower bound. We still prefer to take d = (l + u)/2. x 2 appears in the variance term of layer normalization and is later passed to a square root function to compute a standard derivation. To make the input to the square root function valid, i.e., non-negative, we impose a constraint σ L (x) ≥ 0(∀x ∈ [l, u]). σ L d (x) = 2d(x -d) + d 2 is the tangent line passing (d, σ(d) ). For u ≤ 0, x 2 is monotonously decreasing, the constraint we impose is equivalent to σ L (u) = 2du -d 2 ≥ 0, and with d ≤ 0, we have d ≥ 2u. So we take d = max((l + u)/2, 2u). For l ≥ 0, x 2 is monotonously increasing, and thus the constraint we impose is equivalent to σ L (l) = 2dl -d 2 ≥ 0, and with d ≥ 0, we have d ≤ 2l. So we take d = max((l + u)/2, 2l). And for l < 0 < u, since σ L d (0) = -d 2 is negative for d = 0 while d = 0 yields a valid lower bound, we take d = 0.

Square root For the square root function, σ(x) = √ x. It is used the to compute a standard derivation in layer normalization and its input is limited to be non-negative by the lower bounds of x 2 , and thus l ≥ 0. σ(x) is concave, and thus we take the line passing the two endpoints as the lower bound and take the tangent line passing ((l + u)/2, σ((l + u)/2)) as the upper bound.

C LINEAR BOUNDS OF MULTIPLICATIONS AND DIVISIONS:

We provide a mathematical proof of optimal parameters for linear bounds of multiplications used in Sec. 3.3. We also show that linear bounds of division can be indirectly obtained from bounds of multiplications and the reciprocal function.

For each multiplication, we aim to bound z = xy with two linear bounding planes z L = α L x + β L y+γ L and z U = α U x+β U y+γ U , where x and y are both variables and x ∈ [l x , u x ], y ∈ [l y , u y ] are concrete bounds of x, y obtained from previous layers, such that:

z L = α L x + β L y + γ L ≤ z = xy ≤ z U = α U x + β U y + γ U ∀(x, y) ∈ [l x , u x ] × [l y , u y ].

Our goal is to determine optimal parameters of bounding planes, i.e., α L , β L , γ L , α U , β U , γ U , such that the bounds are as tight as possible.

C.1 LOWER BOUND OF MULTIPLICATIONS:

We define a difference function F L (x, y) which is the difference between the original function z = xy and the lower bound z L = α L x + β L y + γ L :

F L (x, y) = xy -(α L x + β L y + γ L ).

To make the bound as tight as possible, we aim to minimize the integral of the difference function F L (x, y) on our concerned area (x, y) ∈ [l x , u x ] × [l y , u y ], which is equivalent to maximizing

EQUATION

while F L (x, y) ≥ 0 (∀(x, y) ∈ [l x , u x ] × [l y , u y ]). For an optimal bounding plane, there must exist a point (x 0 , y 0 ) ∈ [l x , u x ] × [l y , u y ] such that F L (x 0 , y 0 ) = 0 (otherwise we can validly increase γ L to make V L larger). To ensure that F L (x, y) ≥ 0 within the concerned area, we need to ensure that the minimum value of F L (x, y) is non-negative. We show that we only need to check cases when (x, y) is any of (l x , l y ), (l x , u y ), (u x , l y ), (u x , u y ), i.e., points at the corner of the considered area. The partial derivatives of F L are:

∂F L ∂x = y -α L , ∂F L ∂y = x -β L .

If there is (x 1 , y 1 ) ∈ (l x , u x ) × (l y , u y ) such that F L (x 1 , y 1 ) ≤ F (x, y) (∀(x, y) ∈ [l x , u x ] × [l y , u y ]), ∂F L ∂x (x 1 , y 1 ) = ∂F L ∂y (x 1 , y 1 ) = 0 should hold true. Thereby ∂F L ∂x (x, y), ∂F L ∂y (x, y) < 0 (∀(x, y) ∈ [l x , x 1 ) × [l y , y 1 )), and thus F L (l x , l y ) < F L (x 1 , y 1 ) and (x 1 , y 1 ) cannot be the point with the minimum value of F L (x, y). On the other hand, if there is (x 1 , y 1 )(x 1 = l x , y 1 ∈ (l y , u y )), i.e., on one border of the concerned area but not on any corner, ∂F L ∂y (x 1 , y 1 ) = 0 should hold true. Thereby, ∂F L ∂y (x, y) = ∂F L ∂y (x 1 , y) = 0 (∀(x, y), x = x 1 = l x ), and F L (x 1 , y 1 ) = F L (x 1 , l y ) = F L (l x , l y ). This property holds true for the other three borders of the concerned area. Therefore, other points within the concerned area cannot have smaller function value F L (x, y), so we only We compare the tightness of bounds by the forward process and the backward process respectively. To illustrate the difference, for simplicity, we consider a m-layer feed-forward network Φ (0) = x, y (l) = W (l) Φ (l-1) (x) + b (l) , Φ (l) (x) = σ(y (l) (x))(0 < l ≤ m), where x is the input vector, W (l) and b (l) are the weight matrix and the bias vector for the l-th layer respectively, y (l) (x) is the pre-activation vector of the l-th layer, Φ (l) (x) is the vector of neurons in the l-th layer, and σ(•) is an activation function. Before taking global bounds, both the backward process and the forward process bound Φ (l) j (x) with linear functions of x. When taking global bounds as Eq. (3) and Eq. ( 4), only the norm of weight matrix is directly related to the in binary search for certified lower bounds. Therefore, we try to measure the tightness of the computed bounds using the difference between weight matrices for lower bounds and upper bounds respectively. We show how it is computed for the forward process and the backward process respectively.

D.1 THE FORWARD PROCESS:

For the forward process, we bound each neuron Φ (l) j (x) with linear functions:

Ω (l),L j,: x + Θ (l),L j ≤ Φ (l) j (x) ≤ Ω (l),U j,: x + Θ (l),U j .

To measure the tightness of the bounds, we are interested in Ω (l),L , Ω (l),U , and also Ω (l),U -Ω (l),L . Initially, Ω (0),L/U = I, Θ (0),L/U = 0, Ω (0),U -Ω (0),L = 0. We can forward propagate the bounds of Φ (l-1) (x) to y (l) (x):

Ω (l),y,L j,: x + Θ (l),y,L j ≤ y (l) j (x) ≤ Ω (l),y,U j,: x + Θ (l),y,U j ,

where Ω (l),y,L/U j,:

= W (l) j,i >0 W (l) j,i Ω (l-1),L/U i,: + W (l) j,i <0 W (l) j,i Ω (l-1),U/L i,:

,

Θ (l),y,L/U = W (l) j,i >0 W (l) j,i Θ (l-1),L/U i + W (l) j,i <0 W (l) j,i Θ (l-1),U/L i + b (l) .

With the global bounds of y (l) (x) that can be obtained with Eq. (3) and Eq. ( 4), we bound the activation function:

α (l),L j y (l) (x) + β (l),L j ≤ σ(y (l) j (x)) ≤ α (l),U j y (l) (x) + β (l),U j .

And then bounds can be propagated from Φ (l-1) (x) to Φ (l) (x):

Ω (l),L/U j,: = α (l),L/U j Ω (l),y,L/U j,: α (l),L/U j ≥ 0 α (l),L/U j Ω (l),y,U/L j,: α (l),L/U j < 0 , Θ (l),L/U j = α (l),L/U j Θ (l),y,L/U j + β (l),L/U j α (l),L/U j ≥ 0 α (l),L/U j Θ (l),y,U/L j + β (l),L/U j α (l),L/U j < 0 .

Therefore, Ω (l),U j,:

-Ω (l),L j,:

= (α (l),U j -α (l),L j )|W (l) j |(Ω (l-1),U j,:

-Ω (l-1),L j,:

) (9) illustrates how the tightness of the bounds is changed from earlier layers to later layers.

D.2 THE BACKWARD PROCESS AND DISCUSSIONS:

For the backward process, we bound the neurons in the l-th layer with linear functions of neurons in a previous layer, the l -th layer:

Φ (l,l ),L = Λ (l,l ),L j,: Φ (l ) (x) + ∆ (l,l ),L j ≤ Φ (l) (x) ≤ Λ (l,l ),U j,: Φ (l ) (x) + ∆ (l,l ),U j = Φ (l,l ),U .

We have shown in Sec. 3.1 how such bounds can be propagated to l = 0, for the case when the input is sequential. For the nonsequential case we consider here, it can be regarded as a special case when the input length is 1. So we can adopt the method in Sec. 3.1 to propagate bounds for the feed-forward network we consider here. We are interested in Λ (l,l ),L , Λ (l,l ),U and also Λ (l,l ),U -Λ (l,l ),L . Weight matrices of linear bounds before taking global bounds are Λ (l,0),L and Λ (l,0),U which are obtained by propagating the bounds starting from Λ (l,l),L = Λ (l,l),U = I. According to bound propagation described in Sec. 3.2, Λ (l,l -1),U :,j -Λ (l,l -1),L :,j = (α (l ),U j (Λ (l,l ),U :,j,+ -Λ (l,l ),L :,j,-)-α (l ),L j (Λ (l,l ),L :,j,+ -Λ (l,l ),U :,j,-))W (l ) (10) illustrates how the tightness bounds can be measured during the backward bound propagation until l = 0.

There is a W (l ) in Eq. ( 10) instead of |W (l ) | in Eq. ( 9). The norm of (Ω (l),U j,:

-Ω (l),L j,: ) in Eq. ( 9) can quickly grow large as l increases during the forward propagation when W (l) j is greater than 1, while this generally holds true for neural networks to have W (l) j greater than 1 in feedforward layers. While in Eq. ( 10), W (l ) j can have both positive and negative elements and tends to allow cancellations for different W (l ) j,i , and thus the norm of (Λ (l,l -1),U :,j -Λ (l,l -1),L :,j ) tends to be smaller. Therefore, the bounds computed by the backward process tend to be tighter than those by the forward framework, which is consistent with our experiment results in Table 3 .

E IMPACT OF MODIFYING THE LAYER NORMALIZATION:

The original Transformers have a layer normalization after the embedding layer, and two layer normalization before and after the feed-forward part respectively in each Transformer layer. We modify the layer normalization, f (x) = w(x -µ)/σ + b, where x is d-dimensional a vector to be normalized, µ and σ are the mean and standard derivation of {x i } respectively, and w and b are gain and bias parameters respectively. σ = (1/d) d i=1 (x i -µ) 2 + s where s is a smoothing constant. It involves (x i -µ) 2 whose linear lower bound is loose and exactly 0 when the range of the x i -µ crosses 0. When the p -norm of the perturbation is relatively larger, there can be many x i -µ with ranges crossing 0, which can cause the lower bound of σ to be small and thereby the upper bound of f i (x) to be large. This can make the certified bounds loose. To tackle this, we modify the layer normalization into f (x) = w(x -µ) + b by removing the standard derivation term. We use an experiment to study the impact of this modification. We compare the clean accuracies and certified bounds of the models with modified layer normalization to models with standard layer normalization and with no layer normalization respectively. Table 5 presents the results. Certified lower bounds of models with no layer normalization or our modification are significantly tighter than those of corresponding models with the standard layer normalization. Meanwhile, the clean accuracies of the models with our modification are comparable with those of the models with the standard layer normalization (slightly lower on Yelp and slightly higher on SST). This demonstrates that it appears to be worthwhile to modify the layer normalization in Transformers for easier verification. 

:



Figures and Tables:

), linear bounds can be further propagated to n k=1 exp(S (l) i,k ). With bounds of exp(S

. Yelp consists of 560,000/38,000 examples in the training/test set and SST consists of 67,349/872/1,821 examples in the training/development/test set. Each example is a sentence or a sentence segment (for the training data of SST only) labeled with a binary sentiment polarity.We verify the robustness of Transformers trained from scratch. For the main experiments, we consider N -layer models (N ≤ 3), with 4 attention heads, hidden sizes of 256 and 512 for self-attention and feed-forward layers respectively, and we use ReLU activations for feed-forward layers. We remove the variance related terms in layer normalization, making Transformers verification bounds tighter while the clean accuracies remain comparable (see Appendix E for discussions). Although our method can be in principal applied to Transformers with any number of layers, we do not use large-scale pre-trained models such as BERT because they are too challenging to be tightly verified with the current technologies.

Table content: <table><tr><td colspan="4">4.2 CERTIFIED BOUNDS</td><td/><td/><td/><td/><td/><td/></tr><tr><td/><td/><td/><td>p</td><td colspan="2">Upper Min Avg</td><td colspan="2">Lower (IBP) Min Avg</td><td>Lower (Ours) Min Avg</td><td colspan="2">Ours vs Upper Min Avg</td></tr><tr><td/><td/><td/><td>1</td><td>9.085</td><td>13.917</td><td>1.4E-4</td><td>3.1E-4</td><td colspan="2">1.423 1.809 16%</td><td>13%</td></tr><tr><td/><td>1</td><td>91.5</td><td>2</td><td>0.695</td><td>1.005</td><td>1.4E-4</td><td>3.1E-4</td><td colspan="2">0.384 0.483 55%</td><td>48%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.117</td><td>0.155</td><td>1.4E-4</td><td>3.1E-4</td><td colspan="2">0.034 0.043 29%</td><td>27%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">10.228 15.452</td><td>1.4E-7</td><td>2.2E-7</td><td>0.389 0.512</td><td>4%</td><td>3%</td></tr><tr><td>Yelp</td><td>2</td><td>91.5</td><td>2</td><td>0.773</td><td>1.103</td><td>1.4E-7</td><td>2.2E-7</td><td colspan="2">0.116 0.149 15%</td><td>14%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.122</td><td>0.161</td><td>1.4E-7</td><td>2.2E-7</td><td>0.010 0.013</td><td>9%</td><td>8%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="5">11.137 15.041 4.3E-10 7.1E-10 0.152 0.284</td><td>1%</td><td>2%</td></tr><tr><td/><td>3</td><td>91.6</td><td>2</td><td>0.826</td><td>1.090</td><td colspan="3">4.3E-10 7.1E-10 0.042 0.072</td><td>5%</td><td>7%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.136</td><td>0.187</td><td colspan="3">4.3E-10 7.1E-10 0.004 0.006</td><td>3%</td><td>3%</td></tr><tr><td/><td/><td/><td>1</td><td>7.418</td><td>8.849</td><td>2.4E-4</td><td>2.7E-4</td><td colspan="2">2.503 2.689 34%</td><td>30%</td></tr><tr><td/><td>1</td><td>83.2</td><td>2</td><td>0.560</td><td>0.658</td><td>2.4E-4</td><td>2.7E-4</td><td colspan="2">0.418 0.454 75%</td><td>69%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.091</td><td>0.111</td><td>2.4E-4</td><td>2.7E-4</td><td colspan="2">0.033 0.036 36%</td><td>32%</td></tr><tr><td/><td/><td/><td>1</td><td>6.781</td><td>8.367</td><td>3.6E-7</td><td>3.8E-7</td><td colspan="2">1.919 1.969 28%</td><td>24%</td></tr><tr><td>SST</td><td>2</td><td>83.5</td><td>2</td><td>0.520</td><td>0.628</td><td>3.6E-7</td><td>3.8E-7</td><td colspan="2">0.305 0.315 59%</td><td>50%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.085</td><td>0.105</td><td>3.6E-7</td><td>3.8E-7</td><td colspan="2">0.024 0.024 28%</td><td>23%</td></tr><tr><td/><td/><td/><td>1</td><td>6.475</td><td>7.877</td><td colspan="4">5.7E-10 6.7E-10 1.007 1.031 16%</td><td>13%</td></tr><tr><td/><td>3</td><td>83.9</td><td>2</td><td>0.497</td><td>0.590</td><td colspan="4">5.7E-10 6.7E-10 0.169 0.173 34%</td><td>29%</td></tr><tr><td/><td/><td/><td>∞</td><td>0.084</td><td>0.101</td><td colspan="4">5.7E-10 6.7E-10 0.013 0.014 16%</td><td>13%</td></tr></table>

A FORWARD PROCESS Comparison of certified lower bounds and computation time (sec) by different methods.

Table content: <table><tr><td colspan="2">Dataset Acc.</td><td/><td>Fully-Forward Min Avg Time</td><td>Fully-Backward Min Avg Time</td><td>Backward &amp; Forward Min Avg Time</td></tr><tr><td/><td/><td>1</td><td>2.122 2.173 12.6</td><td colspan="2">3.485 3.737 141.4 3.479 3.729 24.0</td></tr><tr><td>Yelp</td><td>91.3</td><td>2</td><td>0.576 0.599 12.4</td><td colspan="2">0.867 0.947 140.4 0.866 0.946 26.0</td></tr><tr><td/><td/><td>∞</td><td>0.081 0.084 12.6</td><td colspan="2">0.123 0.136 143.9 0.123 0.136 26.4</td></tr><tr><td/><td/><td>1</td><td>1.545 1.592 13.7</td><td colspan="2">1.891 1.961 177.6 1.891 1.961 26.5</td></tr><tr><td>SST</td><td>83.3</td><td>2</td><td>0.352 0.366 12.6</td><td colspan="2">0.419 0.439 178.8 0.419 0.439 24.3</td></tr><tr><td/><td/><td>∞</td><td>0.048 0.050 14.6</td><td colspan="2">0.058 0.061 181.3 0.058 0.061 24.3</td></tr></table>

Average importance scores of the most/least important words identified from 100 examples respectively on SST by different methods. For the most important words identified, larger important scores are better, and vice versa. Additionally, we show most/least important words identified from 10 examples on the Yelp dataset. Boldfaced words are considered to have strong sentiment polarities, and they should appear as most important words rather than least important ones.

Table content: <table><tr><td>Method</td><td>Importance Score (on SST)</td><td>Words Identified from 10 Examples on the Yelp Dataset (split by " / ")</td></tr><tr><td/><td/><td>Most Important Words or Symbols</td></tr><tr><td>Grad</td><td>0.47</td><td/></tr><tr><td/><td/><td>Least Important Words or Symbols</td></tr><tr><td>Grad</td><td>0.40</td><td/></tr></table>

Dataset N LayerNorm Acc.

Table content: <table><tr><td/><td/><td/><td>p</td><td>Min</td><td colspan="2">Upper</td><td>Avg</td><td>Lower Min Avg</td><td>Ours vs Upper Min Avg</td></tr><tr><td/><td/><td/><td>1</td><td colspan="6">189.934 199.265 0.010 0.022 5.3E-5 1.1E-4</td></tr><tr><td/><td>Standard</td><td>91.7</td><td>2</td><td colspan="2">15.125</td><td colspan="2">15.384</td><td colspan="2">0.008 0.019 5.5E-4 1.3E-3</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">2.001</td><td colspan="2">3.066</td><td colspan="2">0.002 0.005 8.7E-4 1.7E-3</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">8.044</td><td colspan="2">12.948</td><td>1.360 1.684</td><td>17%</td><td>13%</td></tr><tr><td>1</td><td>None</td><td>91.4</td><td>2</td><td colspan="2">0.580</td><td colspan="2">0.905</td><td>0.363 0.447</td><td>63%</td><td>49%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.086</td><td colspan="2">0.127</td><td>0.033 0.040</td><td>38%</td><td>32%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">9.085</td><td colspan="2">13.917</td><td>1.423 1.809</td><td>16%</td><td>13%</td></tr><tr><td/><td>Ours</td><td>91.5</td><td>2</td><td colspan="2">0.695</td><td colspan="2">1.005</td><td>0.384 0.483</td><td>55%</td><td>48%</td></tr><tr><td>Yelp</td><td/><td/><td>∞ 1</td><td colspan="6">0.117 190.476 201.092 0.002 0.004 1.2E-5 1.8E-5 0.155 0.034 0.043 29% 27%</td></tr><tr><td/><td>Standard</td><td>92.0</td><td>2</td><td colspan="2">15.277</td><td colspan="2">15.507</td><td colspan="2">0.001 0.002 9.0E-5 1.6E-4</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">2.022</td><td colspan="2">2.901</td><td colspan="2">0.000 0.000 9.5E-5 1.3E-4</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">8.112</td><td colspan="2">15.225</td><td>0.512 0.631</td><td>6%</td><td>4%</td></tr><tr><td>2</td><td>None</td><td>91.5</td><td>2</td><td colspan="2">0.587</td><td colspan="2">1.042</td><td>0.123 0.154</td><td>21%</td><td>15%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.081</td><td colspan="2">0.140</td><td>0.010 0.013</td><td>13%</td><td>9%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">10.228</td><td colspan="2">15.452</td><td>0.389 0.512</td><td>4%</td><td>3%</td></tr><tr><td/><td>Ours</td><td>91.5</td><td>2</td><td colspan="2">0.773</td><td colspan="2">1.103</td><td>0.116 0.149</td><td>15%</td><td>14%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.122</td><td colspan="2">0.161</td><td>0.010 0.013</td><td>9%</td><td>8%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="6">190.777 194.961 0.008 0.015 4.2E-5 7.8E-5</td></tr><tr><td/><td>Standard</td><td>83.0</td><td>2</td><td colspan="2">15.549</td><td colspan="2">15.630</td><td colspan="2">0.006 0.013 4.1E-4 8.2E-4</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">2.241</td><td colspan="2">2.504</td><td colspan="2">0.001 0.003 5.2E-4 1.2E-3</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">6.921</td><td colspan="2">8.417</td><td>2.480 2.659</td><td>36%</td><td>32%</td></tr><tr><td>1</td><td>None</td><td>83.0</td><td>2</td><td colspan="2">0.527</td><td colspan="2">0.628</td><td>0.411 0.447</td><td>78%</td><td>71%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.089</td><td colspan="2">0.109</td><td>0.032 0.035</td><td>36%</td><td>32%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">7.418</td><td colspan="2">8.849</td><td>2.503 2.689</td><td>34%</td><td>30%</td></tr><tr><td/><td>Ours</td><td>83.2</td><td>2</td><td colspan="2">0.560</td><td colspan="2">0.658</td><td>0.418 0.454</td><td>75%</td><td>69%</td></tr><tr><td>SST</td><td/><td/><td>∞ 1</td><td colspan="6">0.091 191.742 196.365 0.002 0.004 9.6E-6 1.9E-5 0.111 0.033 0.036 36% 32%</td></tr><tr><td/><td>Standard</td><td>82.5</td><td>2</td><td colspan="2">15.554</td><td colspan="2">15.649</td><td colspan="2">0.001 0.003 7.0E-5 1.7E-4</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">2.252</td><td colspan="2">2.513</td><td colspan="2">0.000 0.000 6.6E-5 1.9E-4</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">6.742</td><td colspan="2">8.118</td><td>1.821 1.861</td><td>27%</td><td>23%</td></tr><tr><td>2</td><td>None</td><td>83.5</td><td>2</td><td colspan="2">0.515</td><td colspan="2">0.610</td><td>0.298 0.306</td><td>58%</td><td>50%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.085</td><td colspan="2">0.103</td><td>0.023 0.024</td><td>28%</td><td>23%</td></tr><tr><td/><td/><td/><td>1</td><td colspan="2">6.781</td><td colspan="2">8.367</td><td>1.919 1.969</td><td>28%</td><td>24%</td></tr><tr><td/><td>Ours</td><td>83.5</td><td>2</td><td colspan="2">0.520</td><td colspan="2">0.628</td><td>0.305 0.315</td><td>59%</td><td>50%</td></tr><tr><td/><td/><td/><td>∞</td><td colspan="2">0.085</td><td colspan="2">0.105</td><td>0.024 0.024</td><td>28%</td><td>23%</td></tr></table>

Clean accuracies, upper bounds, certified lower bounds by our method of models with different layer normalization settings.

Table content: <table/>

Peer review (iclr2020):
Contributions: This paper develops an algorithm for verifying the robustness of transformers with self-attention layers when the inputs for one input word embedding are perturbed. Unlike previous work the present work can deal with cross nonlinearity and cross position dependency and the lower bounds derived in the paper are much tighter than the Interval Boundary Propagation (IBP) method which uses backward propagation. The core contribution is expounded by developing bounds for multiplication (xy) and division (x/y) and using this to compute tight bounds on self-attention layer computations. Further by introducing a forward process then combining it with a backward process, they can substantially reduce computation time as compared to the IBP method (in experiments demonstrated to be over an order of magnitude). Experimental results: Very reasonable and well thought out experiments demonstrate that a) the lower bounds produced by the method are very tight (much better than IBP for instance), and the backward&forward approach is an order of magnitude faster than the fully backward approach though they are both equally tight (the fully forward method is much looser by contrast) The experiments on using the bound to detect the word whose perturbation will cause biggest impact (ie the most important word for explaining the decision) is much less convincing. In particular in Table 3, the second example shows that the method identifiews *a*, *the* *our* etc, which at least do not appear to be most salient to this reviewer. [Question to authors: perhaps I misunderstood something here -- if so please explain and help me understand it better]. Evaluation / Suggestions. Overall I liked the paper quite a bit. I may not have fully understood something about table 3 and would welcome some explanation about why the results shown there indicate that their method works better int he second example.