Title: THE LOGICAL EXPRESSIVENESS OF GRAPH NEURAL NETWORKS

Abstract: The ability of graph neural networks (GNNs) for distinguishing nodes in graphs has been recently characterized in terms of the Weisfeiler-Lehman (WL) test for checking graph isomorphism. This characterization, however, does not settle the issue of which Boolean node classifiers (i.e., functions classifying nodes in graphs as true or false) can be expressed by GNNs. We tackle this problem by focusing on Boolean classifiers expressible as formulas in the logic FOC 2 , a well-studied fragment of first order logic. FOC 2 is tightly related to the WL test, and hence to GNNs. We start by studying a popular class of GNNs, which we call AC-GNNs, in which the features of each node in the graph are updated, in successive layers, only in terms of the features of its neighbors. We show that this class of GNNs is too weak to capture all FOC 2 classifiers, and provide a syntactic characterization of the largest subclass of FOC 2 classifiers that can be captured by AC-GNNs. This subclass coincides with a logic heavily used by the knowledge representation community. We then look at what needs to be added to AC-GNNs for capturing all FOC 2 classifiers. We show that it suffices to add readout functions, which allow to update the features of a node not only in terms of its neighbors, but also in terms of a global attribute vector. We call GNNs of this kind ACR-GNNs. We experimentally validate our findings showing that, on synthetic data conforming to FOC 2 formulas, AC-GNNs struggle to fit the training data while ACR-GNNs can generalize even to graphs of sizes not seen during training.

1. INTRODUCTION:

Graph neural networks (GNNs) (Merkwirth & Lengauer, 2005; Scarselli et al., 2009) are a class of neural network architectures that has recently become popular for a wide range of applications dealing with structured data, e.g., molecule classification, knowledge graph completion, and Web page ranking (Battaglia et al., 2018; Gilmer et al., 2017; Kipf & Welling, 2017; Schlichtkrull et al., 2018) . The main idea behind GNNs is that the connections between neurons are not arbitrary but reflect the structure of the input data. This approach is motivated by convolutional and recurrent neural networks and generalize both of them (Battaglia et al., 2018) . Despite the fact that GNNs have recently been proven very efficient in many applications, their theoretical properties are not yet well-understood. In this paper we make a step towards understanding their expressive power by establishing connections between GNNs and well-known logical formalisms. We believe these connections to be conceptually important, as they permit us to understand the inherently procedural behavior of some fragments of GNNs in terms of the more declarative flavor of logical languages.

Two recent papers (Morris et al., 2019; Xu et al., 2019) have started exploring the theoretical properties of GNNs by establishing a close connection between GNNs and the Weisfeiler-Lehman (WL) test for checking graph isomorphism. The WL test works by constructing a labeling of the nodes of the graph, in an incremental fashion, and then decides whether two graphs are isomorphic by comparing the labeling of each graph. To state the connection between GNNs and this test, consider the simple GNN architecture that updates the feature vector of each graph node by combining it with the aggregation of the feature vectors of its neighbors. We call such GNNs aggregate-combine GNNs, or AC-GNNs. The authors of these papers independently observe that the node labeling produced by the WL test always refines the labeling produced by any GNN. More precisely, if two nodes are labeled the same by the algorithm underlying the WL test, then the feature vectors of these nodes produced by any AC-GNN will always be the same. Moreover, there are AC-GNNs that can reproduce the WL labeling, and hence AC-GNNs can be as powerful as the WL test for distinguishing nodes. This does not imply, however, that AC-GNNs can capture every node classifier-that is, a function assigning true or false to every node-that is refined by the WL test. In fact, it is not difficult to see that there are many such classifiers that cannot be captured by AC-GNNs; one simple example is a classifier assigning true to every node if and only if the graph has an isolated node. Our work aims to answer the question of what are the node classifiers that can be captured by GNN architectures such as AC-GNNs.

To start answering this question, we propose to focus on logical classifiers-that is, on unary formulas expressible in first order predicate logic (FO): such a formula classifies each node v according to whether the formula holds for v or not. This focus gives us an opportunity to link GNNs with declarative and well understood formalisms, and to establish conclusions about GNNs drawing upon the vast amount of work on logic. For example, if one proves that two GNN architectures are captured with two logics, then one can immediately transfer all the knowledge about the relationships between those logics, such as equivalence or incomparability of expressiveness, to the GNN setting.

For AC-GNNs, a meaningful starting point to measure their expressive power is the logic FOC 2 , the two variable fragment of first order predicate logic extended with counting quantifiers of the form ∃ ≥N ϕ, which state that there are at least N nodes satisfying formula ϕ (Cai et al., 1992) . Indeed, this choice of FOC 2 is justified by a classical result due to Cai et al. (1992) establishing a tight connection between FOC 2 and WL: two nodes in a graph are classified the same by the WL test if and only if they satisfy exactly the same unary FOC 2 formulas. Moreover, the counting capabilities of FOC 2 can be mimicked in FO (albeit with more than just two variables), hence FOC 2 classifiers are in fact logical classifiers according to our definition.

Given the connection between AC-GNNs and WL on the one hand, and that between WL and FOC 2 on the other hand, one may be tempted to think that the expressivity of AC-GNNs coincides with that of FOC 2 . However, the reality is not as simple, and there are many FOC 2 node classifiers (e.g., the trivial one above) that cannot be expressed by AC-GNNs. This leaves us with the following natural questions. First, what is the largest fragment of FOC 2 classifiers that can be captured by AC-GNNs? Second, is there an extension of AC-GNNs that allows to express all FOC 2 classifiers? In this paper we provide answers to these two questions. The following are our main contributions.

• We characterize exactly the fragment of FOC 2 formulas that can be expressed as AC-GNNs. This fragment corresponds to graded modal logic (de Rijke, 2000), or, equivalently, to the description logic ALCQ, which has received considerable attention in the knowledge representation community (Baader et al., 2003; Baader & Lutz, 2007) . • Next we extend the AC-GNN architecture in a very simple way by allowing global readouts, where in each layer we also compute a feature vector for the whole graph and combine it with local aggregations; we call these aggregate-combine-readout GNNs (ACR-GNNs). These networks are a special case of the ones proposed by Battaglia et al. (2018) for relational reasoning over graph representations. In this setting, we prove that each FOC 2 formula can be captured by an ACR-GNN.

We experimentally validate our findings showing that the theoretical expressiveness of ACR-GNNs, as well as the differences between AC-GNNs and ACR-GNNs, can be observed when we learn from examples. In particular, we show that on synthetic graph data conforming to FOC 2 formulas, AC-GNNs struggle to fit the training data while ACR-GNNs can generalize even to graphs of sizes not seen during training.

2. GRAPH NEURAL NETWORKS:

In this section we describe the architecture of AC-GNNs and introduce other related notions. We concentrate on the problem of Boolean node classification: given a (simple, undirected) graph G = (V, E) in which each vertex v ∈ V has an associated feature vector x v , we wish to classify each graph node as true or false; in this paper, we assume that these feature vectors are one-hot encodings of node colors in the graph, from a finite set of colors. The neighborhood

N G (v) of a node v ∈ V is the set {u | {v, u} ∈ E}.

The basic architecture for GNNs, and the one studied in recent studies on GNN expressibility (Morris et al., 2019; Xu et al., 2019) , consists of a sequence of layers that combine the feature vectors of every node with the multiset of feature vectors of its neighbors. Formally, let {AGG (i) } L i=1 and {COM (i) } L i=1 be two sets of aggregation and combination functions. An aggregate-combine GNN (AC-GNN) computes vectors x (i) v for every node v of the graph G, via the recursive formula

EQUATION

where each x

(0) v is the initial feature vector x v of v. Finally, each node v of G is classified according to a Boolean classification function CLS applied to x (L) v . Thus, an AC-GNN with L layers is defined as a tuple A = {AGG (i) } L i=1 , {COM (i) } L i=1

, CLS , and we denote by A(G, v) the class (i.e., true or false) assigned by A to each node v in G.1 

There are many possible aggregation, combination, and classification functions, which produce different classes of GNNs (Hamilton et al., 2017; Kipf & Welling, 2017; Morris et al., 2019; Xu et al., 2019) . A simple, yet common choice is to consider the sum of the feature vectors as the aggregation function, and a combination function as

EQUATION

where C (i) and A (i) are matrices of parameters, b (i) is a bias vector, and f is a non-linearity function, such as relu or sigmoid. We call simple an AC-GNN using these functions. Furthermore, we say that an AC-GNN is homogeneous if all AGG (i) are the same and all COM (i) are the same (share the same parameters across layers). In most of our positive results we construct simple and homogeneous GNNs, while our negative results hold in general (i.e., for GNNs with arbitrary aggregation, combining, and classification functions).

The Weisfeiler-Lehman (WL) test is a powerful heuristic used to solve the graph isomorphism problem (Weisfeiler & Leman, 1968) , or, for our purposes, to determine whether the neighborhoods of two nodes in a graph are structurally close or not. Due to space limitations, we refer to (Cai et al., 1992) for a formal definition of the underlying algorithm, giving only its informal description: starting from a colored graph, the algorithm iteratively assigns, for a certain number of rounds, a new color to every node in the graph; this is done in such a way that the color of a node in each round has a one to one correspondence with its own color and the multiset of colors of its neighbors in the previous round. An important observation is that the rounds of the WL algorithm can be seen as the layers of an AC-GNN whose aggregation and combination functions are all injective (Morris et al., 2019; Xu et al., 2019) . Furthermore, as the following proposition states, an AC-GNN classification can never contradict the WL test. Proposition 2.1 (Morris et al., 2019; Xu et al., 2019) . If the WL test assigns the same color to two nodes in a graph, then every AC-GNN classifies either both nodes as true or both nodes as false.

3. CONNECTION BETWEEN GNNS AND LOGIC:



3.1. LOGICAL NODE CLASSIFIERS:

Our study relates the power of GNNs to that of classifiers expressed in first order (FO) predicate logic over (undirected) graphs where each vertex has a unique color (recall that we call these classifiers logical classifiers). To illustrate the idea of logical node classifiers, consider the formula α(x) := Red(x) ∧ ∃y E(x, y) ∧ Blue(y) ∧ ∃z E(x, z) ∧ Green(z) .

(3)

This formula has one free variable, x, which is not bounded by any quantifier of the form ∃ or ∀, and two quantified variables y and z. In general, formulas with one free variable are evaluated over nodes of a given graph. For example, the above formula evaluates to true exactly in those nodes v whose color is Red and that have both a Blue and a Green neighbor. In this case, we say that node v of G satisfies α, and denote this by (G, v) |= α.

Formally, a logical (node) classifier is given by a formula ϕ(x) in FO logic with exactly one free variable. This formula classifies as true those nodes v in G such that (G, v) |= ϕ, while all other nodes (i.e., those with (G, v) |= ϕ) are classified as false. We say that a GNN classifier captures a logical classifier when both classifiers coincide over every node in every possible input graph. Definition 3.1. A GNN classifier A captures a logical classifier ϕ(x) if for every graph G and node v in G, it holds that A(G, v) = true if and only if (G, v) |= ϕ.

3.2. LOGIC FOC 2:

Logical classifiers are useful as a declarative formalism, but as we will see, they are too powerful to compare them to AC-GNNs. Instead, for reasons we explain later we focus on classifiers given by formulas in FOC 2 , the fragment of FO logic that only allows formulas with two variables, but in turn permits to use counting quantifiers.

Let us briefly introduce FOC 2 and explain why it is a restriction of FO logic. The first remark is that reducing the number of variables used in formulas drastically reduces their expressive power.

Consider for example the following FO formula expressing that x is a red node, and there is another node, y, that is not connected to x and that has at least two blue neighbors, z 1 and z 2 :

β(x) := Red(x) ∧ ∃y ¬E(x, y)∧∃z 1 ∃z 2 E(y, z 1 )∧E(y, z 2 )∧z 1 = z 2 ∧Blue(z 1 )∧Blue(z 2 ) .

The formula β(x) uses four variables, but it is possible to find an equivalent one with just three: the trick is to reuse variable x and replace every occurrence of z 2 in β(x) by x. However, this is as far as we can go with this trick: β(x) does not have an equivalent formula with less than three variables.

In the same way, the formula α(x) given in Equation (3) can be expressed using only two variables, x and y, simply by reusing y in place of z.

That being said, it is possible to extend the logic so that some node properties, such as the one defined by β(x), can be expressed with even less variables. To this end, consider the counting quantifier ∃ ≥N for every positive integer N . Analogously to how the quantifier ∃ expresses the existence of a node satisfying a property, the quantifier ∃ ≥N expresses the existence of at least N different nodes satisfying a property. For example, with ∃ ≥2 we can express β(x) by using only two variables by means of the classifier

EQUATION

Based on this idea, the logic FOC 2 allows for formulas using all FO constructs and counting quantifiers, but restricted to only two variables. Note that, in terms of their logical expressiveness, we have that FOC 2 is strictly less expressive than FO (as counting quantifiers can always be mimicked in FO by using more variables and disequalities), but is strictly more expressive than FO 2 , the fragment of FO that allows formulas to use only two variables (as β(x) belongs to FOC 2 but not to FO 2 ).

The following result establishes a classical connection between FOC 2 and the WL test. Together with Proposition 2.1, this provides a justification for our choice of logic FOC 2 for measuring the expressiveness of AC-GNNs. Proposition 3.2 (Cai et al., 1992) . For any graph G and nodes u, v in G, the WL test colors v and u the same after any number of rounds iff u and v are classified the same by all FOC 2 classifiers.

3.3. FOC 2 AND AC-GNN CLASSIFIERS:

Having Propositions 2.1 and 3.2, one may be tempted to combine them and claim that every FOC 2 classifier can be captured by an AC-GNN. Yet, this is not the case as shown in Proposition 3.3 below.

In fact, while it is true that two nodes are declared indistinguishable by the WL test if and only if they are indistinguishable by all FOC 2 classifiers (Proposition 3.2), and if the former holds then such nodes cannot be distinguished by AC-GNNs (Proposition 2.1), this by no means tells us that every FOC 2 classifier can be expressed as an AC-GNN.

Proposition 3.3. There is an FOC 2 classifier that is not captured by any AC-GNN.

One such FOC 2 classifier is γ(x) in Equation ( 4), but there are infinitely many and even simpler FOC 2 formulas that cannot be captured by AC-GNNs. Intuitively, the main problem is that an AC-GNN has only a fixed number L of layers and hence the information of local aggregations cannot travel further than at distance L of every node along edges in the graph. For instance, the red node in γ(x) may be farther away than the node with the blue neighbours, which means that AC-GNNs would never be able to connect this information. Actually, both nodes may even be in different connected components of a graph, in which case no number of layers would suffice.

The negative result of Proposition 3.3 opens up the following important questions.

1. What kind of FOC 2 classifiers can be captured by AC-GNNs? 2. Can we capture FOC 2 classifiers with GNNs using a simple extension of AC-GNNs?

We provide answers to these questions in the next two sections.

4. THE EXPRESSIVE POWER OF AC-GNNS:

Towards answering our first question, we recall that the problem with AC-GNN classifiers is that they are local, in the sense that they cannot see across a distance greater than their number of layers. Thus, if we want to understand which logical classifiers this architecture is capable of expressing, we must consider logics built with similar limitations in mind. And indeed, in this section we show that AC-GNNs capture any FOC 2 classifier as long as we further restrict the formulas so that they satisfy such a locality property. This happens to be a well-known restriction of FOC 2 , and corresponds to graded modal logic (de Rijke, 2000) or, equivalently, to description logic ALCQ (Baader et al., 2003) , which is fundamental for knowledge representation: for instance, the OWL 2 Web Ontology Language (Motik et al., 2012; W3C OWL Working Group, 2012 ) relies on ALCQ.

The idea of graded modal logic is to force all subformulas to be guarded by the edge predicate E. This means that one cannot express in graded modal logic arbitrary formulas of the form ∃yϕ(y), i.e., whether there is some node that satisfies property ϕ. Instead, one is allowed to check whether some neighbor y of the node x where the formula is being evaluated satisfies ϕ. That is, we are allowed to express the formula ∃y (E(x, y) ∧ ϕ(y)) in the logic as in this case ϕ(y) is guarded by E(x, y). We can define this fragment of FO logic using FO syntax as follows. A graded modal logic formula is either Col(x), for Col a node color, or one of the following, where ϕ and ψ are graded modal logic formulas and N is a positive integer:

¬ϕ(x), ϕ(x) ∧ ψ(x), ∃ ≥N y (E(x, y) ∧ ϕ(y)).

Notice then that the formula δ(x) := Red(x) ∧ ∃y E(x, y) ∧ Blue(y) is in graded modal logic, but the logical classifier γ(x) in Equation ( 4) is not, because the use of ¬E(x, y) as a guard is disallowed. As required, we can now show that AC-GNNs can indeed capture all graded modal logic classifiers.

Proposition 4.1. Each graded modal logic classifier is captured by a simple homogeneous AC-GNN.

The key idea of the construction is that the vectors' dimensions used by the AC-GNN to label nodes, represent the sub-formulas of the captured classifier. Thus, if a feature in a node is 1 then the node satisfies the corresponding sub-formula, and the opposite holds after evaluating L layers, where L is the "quantifier depth" of the classifier (which does not depend on the graph). The construction uses simple, homogeneous AC-GNNs with the truncated relu non-linearity max(0, min(x, 1)). The formal proof of Proposition 4.1, as well as other formal statements, can be found in the Appendix.

An interesting question that we leave as future work is to investigate whether the same kind of construction can be done with AC-GNNs using different aggregate and combine operators than the ones we consider here; for instance, using max instead of sum to aggregate the feature vectors of the neighbors, or using other non-linearity such as sigmoid, etc.

The relationship between AC-GNNs and graded modal logic goes further: we can show that graded modal logic is the "largest" class of logical classifiers captured by AC-GNNs. This means that the only FO formulas that AC-GNNs are able to learn accurately are those in graded modal logic.

Theorem 4.2. A logical classifier is captured by AC-GNNs if and only if it can be expressed in graded modal logic.

The backward direction of this theorem is Proposition 4.1, while the proof of the forward direction is based on a recently communicated extension of deep results in finite model theory (Otto, 2019) . We point out that the forward direction holds no matter which aggregate and combine operators are considered, i.e., this is a limitation of the architecture for AC-GNNs, not of the specific functions that one chooses to update the features.

5 GNNS FOR CAPTURING FOC 2

5.1. GNNS WITH GLOBAL READOUTS:

In this section we tackle our second question: which kind of GNN architecture we need to capture all FOC 2 classifiers? Recall that the main shortcoming of AC-GNNs for expressing such classifiers is their local behavior. A natural way to break such a behavior is to allow for a global feature computation on each layer of the GNN. This is called a global attribute computation in the framework of Battaglia et al. (2018) . Following the recent GNN literature (Gilmer et al., 2017; Morris et al., 2019; Xu et al., 2019) , we refer to this global operation as a readout.

Formally, an aggregate-combine-readout GNN (ACR-GNN) extends AC-GNNs by specifying readout functions {READ (i) } L i=1 , which aggregate the current feature vectors of all the nodes in a graph. Then, the vector x (i) v of each node v in G on each layer i, is computed by the following formula, generalizing Equation (1):

x (i) v = COM (i) x (i-1) v , AGG (i) { {x (i-1) u | u ∈ N G (v)} } , READ (i) { {x (i-1) u | u ∈ G} } . (5)

Intuitively, every layer in an ACR-GNN first computes (i.e., "reads out") the aggregation over all the nodes in G; then, for every node v, it computes the aggregation over the neighbors of v; and finally it combines the features of v with the two aggregation vectors. All the notions about AC-GNNs extend to ACR-GNNs in a straightforward way; for example, a simple ACR-GNN uses the sum as the function READ (i) in each layer, and the combination function i) with a matrix R (i) , generalizing Equation (2).

COM (i) (x 1 , x 2 , x 3 ) = f x 1 C (i) + x 2 A (i) + x 3 R (i) + b (

5.2. ACR-GNNS AND FOC 2:

To see how a readout function could help in capturing non-local properties, consider again the logical classifier γ(x) in Equation ( 4), that assigns true to every red node v as long as there is another node not connected with v having two blue neighbors. We have seen that AC-GNNs cannot capture this classifier. However, using a single readout plus local aggregations one can implement this classifier as follows. First, define by B the property "having at least 2 blue neighbors". Then an ACR-GNN that implements γ(x) can (1) use one aggregation to store in the local feature of every node if the node satisfies B, then (2) use a readout function to count how many nodes satisfying B exist in the whole graph, and (3) use another local aggregation to count how many neighbors of every node satisfiy B. Then γ is obtained by classifying as true every red node having less neighbors satisfying B than the total number of nodes satisfying B in the whole graph. It turns out that the usage of readout functions is enough to capture all non-local properties of FOC 2 classifiers.

Theorem 5.1. Each FOC 2 classifier can be captured by a simple homogeneous ACR-GNN.

The construction is similar to that of Proposition 4.1 and uses simple, homogeneous ACR-GNNsthat is, the readout function is just the sum of all the local node feature vectors. Moreover, the readout functions are only used to deal with subformulas asserting the existence of a node that is not connected to the current node in the graph, just as we have done for classifier γ(x). As an intermediate step in the proof, we use a characterization of FOC 2 using an extended version of graded modal logic, which was obtained by Lutz et al. (2001) . We leave as a challenging open problem whether FOC 2 classifiers are exactly the logical classifiers captured by ACR-GNNs.

5.3. COMPARING THE NUMBER OF READOUT LAYERS:

The proof of Theorem 5.1 constructs GNNs whose number of layers depends on the formula being captured-that is, readout functions are used unboundedly many times in ACR-GNNs for capturing different FOC 2 classifiers. Given that a global computation can be costly, one might wonder whether this is really needed, or if it is possible to cope with all the complexity of such classifiers by performing only few readouts. We next show that actually just one readout is enough. However, this reduction in the number of readouts comes at the cost of severely complicating the resulting GNN.

Formally, an aggregate-combine GNN with final readout (AC-FR-GNN) results out of using any number of layers as in the AC-GNN definition, together with a final layer that uses a readout function, according to Equation (5). Theorem 5.2. Each FOC 2 classifier is captured by an AC-FR-GNN.

The AC-FR-GNN in the proof of this theorem is not based on the idea of evaluating the formula incrementally along layers, as in the proofs of Proposition 4.1 and Theorem 5.1, and it is not simple (note that AC-FR-GNNs are never homogeneous). Instead, it is based on a refinement of the GIN architecture proposed by Xu et al. (2019) to obtain as much information as possible about the local neighborhood in graphs, followed by a readout and combine functions that use this information to deal with non-local constructs in formulas. The first component we build is an AC-GNN that computes an invertible function mapping each node to a number representing its neighborhood (how big is this neighborhood depends on the classifier to be captured). This information is aggregated so that we know for each different type of a neighborhood how many times it appears in the graph. We then use the combine function to evaluate FOC 2 formulas by decoding back the neighborhoods.

6. EXPERIMENTAL RESULTS:

We perform experiments with synthetic data to empirically validate our results. The motivation of this section is to show that the theoretical expressiveness of ACR-GNNs, as well as the differences between AC-and ACR-GNNs, can actually be observed when we learn from examples. We perform two sets of experiments: experiments to show that ACR-GNNs can learn a very simple FOC 2 node classifier that AC-GNNs cannot learn, and experiments involving complex FOC 2 classifiers that need more intermediate readouts to be learned. We implemented our experiments in the PyTorch Geometric library (Fey & Lenssen, 2019) . Besides testing simple AC-GNNs, we also tested the GIN network proposed by Xu et al. (2019) (we consider the implementation by Fey & Lenssen (2019) and adapted it to classify nodes). Our experiments use synthetic graphs, with five initial colors encoded as one-hot features, divided in three sets: train set with 5k graphs of size up to 50-100 nodes, test set with 500 graphs of size similar to the train set, and another test set with 500 graphs of size bigger than the train set. We tried several configurations for the aggregation, combination and readout functions, and report the accuracy on the best configuration. Accuracy in our experiments is computed as the total number of nodes correctly classified among all nodes in all the graphs in the dataset. In every case we run up to 20 epochs with the Adam optimizer. More details on the experimental setting, data, and code can be found in the Appendix. We finally report results on a real benchmark (PPI) where we did not observe an improvement of ACR-GNNs over AC-GNNs.

Separating AC-GNNs and ACR-GNNs We consider a very simple FOC 2 formula defined by α(x) := Red(x) ∧ ∃y Blue(y), which is satisfied by every red node in a graph provided that the graph contains at least one blue node. We tested with line-shaped graphs and Erdös-Renyi (E-R) random graphs with different connectivities. In every set (train and test) we consider 50% of graphs not containing any blue node, and 50% containing at least one blue node (around 20% of nodes are in the true class in every set). For both types of graphs, already single-layer ACR-GNNs showed perfect performance (ACR-1 in Table 1 ). This was what we expected given the simplicity of the property being checked. In contrast, AC-GNNs and GINs (shown in Table 1 as AC-L and GIN-L, representing AC-GNNs and GINs with L layers) struggle to fit the data. For the case of the line-shaped graph, they were not able to fit the train data even by allowing 7 layers. For the case of random graphs, the performance with 7 layers was considerably better. In a closer look at the performance for different connectivities of E-R graphs, we found an improvement for AC-GNNs when we train them with more dense graphs (details in the Appendix). This is consistent with the fact that AC-GNNs are able to move information of local aggregations to distances up to their 6) number of layers. This combined with the fact that random graphs that are more dense make the maximum distances between nodes shorter, may explain the boost in performance for AC-GNNs.

Complex FOC 2 properties In the second experiment we consider classifiers α i (x) constructed as

EQUATION

where ∃ [N,M ] stands for "there exist between N and M nodes" satisfying a given property. Observe that each α i (x) is in FOC 2 , as ∃ [N,M ] can be expressed by combining ∃ ≥N and ¬∃ ≥M +1 . We created datasets with E-R dense graphs and labeled them according to α 1 (x), α 2 (x), and α 3 (x), ensuring in each case that approximately half of all nodes in our dataset satisfy every property. Our experiments show that when increasing the depth of the formula (existential quantifiers with negations inside other existential quantifiers) more layers are needed to increase train and test accuracy (see Table 2 ). We report ACR-GNNs performance up to 3 layers (ACR-L in Table 2 ) as beyond that we did not see any significant improvement. We also note that for the bigger test set, AC-GNNs and GINs are unable to substantially depart from a trivial baseline of 50%. We tested these networks with up to 10 layers but only report the best results on the bigger test set. We also test AC-FR-GNNs with two and three layers (AC-FR-L in Table 2 ). As we expected, although theoretically using a single readout gives the same expressive power as using several of them (Theorem 5.2), in practice more than a single readout can actually help the learning process of complex properties.

PPI We also tested AC-and ACR-GNNs on the Protein-Protein Interaction (PPI) benchmark (Zitnik & Leskovec, 2017) . We chose PPI since it is a node classification benchmark with different graphs in the train set (as opposed to other popular benchmarks for node classification such as Core or Citeseer that have a single graph). Although the best results for both classes of GNNs on PPI were quite high (AC: 97.5 F1, ACR: 95.4 F1 in the test set), we did not observe an improvement when using ACR-GNNs. Chen et al. (2019) recently observed that commonly used benchmarks are inadequate for testing advanced GNN variants, and ACR-GNNs might be suffering from this fact.

7. FINAL REMARKS:

Our results show the theoretical advantages of mixing local and global information when classifying nodes in a graph. Recent works have also observed these advantages in practice, e.g., Deng et al. (2018) 2019) introduced the idea of a "star node" that stores global information of the graph. As mentioned before, our work is close in spirit to that of Xu et al. (2019) and Morris et al. (2019) establishing the correspondence between the WL test and GNNs. In contrast to our work, they focus on graph classification and do not consider the relationship with logical classifiers.

Regarding our results on the links between AC-GNNs and graded modal logic (Theorem 4.2), we point out that very recent work of Sato et al. (2019) establishes close relationships between GNNs and certain classes of distributed local algorithms. These in turn have been shown to have strong correspondences with modal logics (Hella et al., 2015) . Hence, variants of our Proposition 4.1 could be obtained by combining these two lines of work (but it is not clear if this combination would yield AC-GNNs that are simple). However, these works do not investigate the impact of having non-local computations (such as the readouts that we consider), hence our results on the relationships between FO an ACR-GNNs (Theorem 5.1 and 5.2) do not follow from these.

Morris et al. ( 2019) also studied k-GNNs, which are inspired by the k-dimensional WL test. In k-GNNs, graphs are considered as structures connecting k-tuples of nodes instead of just pairs of them. We plan to study how our results on logical classifiers relate to k-GNNs, in particular, with respect to the logic FOC k that extends FOC 2 by allowing formulas with k variables, for each fixed k > 1. Recent work has also explored the extraction of finite state representations from recurrent neural networks as a way of explaining them (Weiss et al., 2018; Koul et al., 2019; Oliva & Lago-Fernández, 2019) . We would like to study how our results can be applied for extracting logical formulas from GNNs as possible explanations for their computations.

We note that ϕ = ϕ L and thus, the last component of each feature vector after evaluating L layers in every node gets a value 1 if and only if the node satisfies ϕ. We will then be able to use a final classification function CLS that simply extracts that particular component.

Formally, the simple homogeneous AC-GNN A ϕ has L layers and uses the aggregation and combine functions

AGG(X) = x∈X x, COM(x, y) = σ xC + yA + b ,

where A, C ∈ R L×L , and b ∈ R L are defined next, and σ is the truncated ReLU activation defined by σ(x) = min(max(0, x), 1). The entries of the -th columns of A, C, and b depend on the sub-formulas of ϕ as follows:

Case 0. if ϕ (x) = Col(x) with Col one of the (base) colors, then C = 1, Case 1. if ϕ (x) = ϕ j (x) ∧ ϕ k (x) then C j = C k = 1 and b = -1, Case 2. if ϕ (x) = ¬ϕ k (x) then C k = -1 and b = 1, Case 3. if ϕ (x) = ∃ ≥N (E(x, y) ∧ ϕ k (y)) then A k = 1 and b = -N + 1,

and all other values in the -th columns of A, C, and b are 0.

We now prove that A ϕ indeed captures ϕ. Let G = (V, E) be a colored graph. For every node v in G we consider the initial feature vector x (0) v = (x 1 , . . . , x L ) such that x = 1 if sub-formula ϕ is the initial color assigned to v, and x = 0 otherwise. By definition, AC-GNN A ϕ will iterate the aggregation and combine functions defined above for L rounds (L layers) to produce feature vectors x (i) v for every node v ∈ G and = 1, . . . , L as follows:

EQUATION

We next prove that for every ϕ ∈ sub(ϕ), every i ∈ { , . . . , L}, and every node v in G it holds that

EQUATION

(i) v ) is the -th component of x (i) v -that is, the -th component of x (i)

v has a 1 if and only if v satisfies ϕ in G. In the rest of the proof we will be continuously using the value of (x

EQUATION

We proceed to prove (8) by induction on the number of sub-formulas of every ϕ . If ϕ has one sub-formula, then ϕ (x) = Col(x) with Col a base color. We next prove that (x

v ) = 1 if and only if v has Col as its initial color. Since ϕ (x) = Col(x) we know that C = 1 and C k = 0 for every k = (see Case 0 above). Moreover, we know that b = 0 and A k = 0 for every k. Then, from Equation ( 9) we obtain that

(x (1) v ) = σ L k=1 (x (0) v ) k C k + {v,u}∈E L k=1 (x (0) u ) k A k + b = σ (x (0) v ) .

Then, given that (x

EQUATION

v ) = 0 otherwise, we have that (x

EQUATION

v ) = 0 otherwise. From this it is easy to prove that for every i ≥ 1 the vector (x (i) v ) satisfies the same property. Now assume that ϕ has more than one Proposition C.1. If a logical classifier α is not equivalent to any graded modal logic formula, then there is no AC-GNN that captures α.

To prove this proposition, we will need the following definition, which is standard in modal logics theory.

Definition C.2. Let G be a graph (simple, undirected and node-colored), v be a node in G, and L ∈ N. The unravelling of v in G at depth L, denoted by Unr L G (v), is the (simple undirected nodecolored) graph that is the tree having -a node (v, u 1 , . . . , u i ) for each path (v, u 1 , . . . , u i ) in G with i ≤ L, -an edge between (v, u 1 , . . . , u i-1 ) and (v, u 1 , . . . , u i ) when {u i-1 , u i } is an edge in G (assuming that u 0 is v), and -each node (v, u 1 , . . . , u i ) colored the same as u i in G.

We then observe the following.

Observation C.3. Let G and G be two graphs, and v and v be two nodes in G and G , respectively. Then for every L ∈ N, the WL test assigns the same color to v and v at round L if and only if there is an isomorphism between Unr L G (v) and Unr L G (v ) sending v to v .

We will write Unr L G (v) Unr L G (v ) to denote the existence of the isomorphism as in this observation. To prove Proposition C.1, we first rephrase Proposition 2.1 in terms of unravellings. Proof. Follows directly from Proposition 2.1 and Observation C.3.

The crucial part of the proof of Proposition C.1 is the following non-trivial result, intuitively establishing that the fragment of unary FO formulas that only depend on the unravelling of a node is exactly the graded modal logic.

Theorem C.5 (Otto, 2019) . Let α be a unary FO formula. If α is not equivalent to a graded modal logic formula then there exist two graphs G, G and two nodes v in G and u in G such that Unr L G (v) Unr L G (v ) for every L ∈ N and such that u |= α in G but u |= α in G .

Proof. This directly follows from the van Benthem & Rosen characterization obtained in (Otto, 2019, Theorem 2. 2) for finite structures (graphs), by noticing that for the notion of graded bisimulation ∼ # introduced in this note, we have that G, u ∼ # G , u if and only if we have that Unr L G (v) Unr L G (v ) for every L ∈ N. We point out here that the fact that the edge relation in G is undirected in our setting (as opposed to E being directed in (Otto, 2019) ), and the fact that every node can only have one color in our setting (as opposed to being able to satisfy multiple "unary predicates" in (Otto, 2019) ) are inessential, and that the proof of (Otto, 2019, Theorem 2. 2) carries over to this setting.

We can now gather all of these to prove Proposition C.1.

Proof of Proposition C.1. Let α be a logical classifier (i.e., a unary FO formula) that is not equivalent to any graded modal logic formula. Assume for a contradiction that there exists an AC-GNN A α that captures α. Since α is not equivalent to any graded modal logic formula, by Theorem C.5 there exist two graphs G, G and two nodes v in G and u in G such that Unr L G (v) Unr L G (v ) for every L ∈ N and such that ( ) u |= α in G but u |= α in G . Since we have that Unr L G (v) Unr L G (v ) for every L ∈ N, by Proposition C.4 we should have that A α (G, u) = A α (G , u ). But this contradicts ( ) and the fact that A α is supposed to capture α.

D PROOF OF THEOREM 5.1:

We first recall the theorem. Theorem 5.1. Each FOC 2 classifier can be captured by a simple homogeneous ACR-GNN.

To prove the theorem, we will use a characterization of the unary FOC 2 formulas provided by (Lutz et al., 2001) that uses a specific modal logic. That logic is defined via what are called modal parameters. We adapt the definitions of (Lutz et al., 2001) to deal with simple undirected node-colored graphs. Definition D.1. A modal parameter is an expression built from the following grammar:

S ::= id | e | S ∪ S | S ∩ S | ¬S.

Given an undirected colored graph G = (V, E) and a node v of G, the interpretation of S on v is the set ε S (v) ⊆ V defined inductively as follows:

-if S = id then ε S (v) := {v};

-if S = e then ε S (v) := {u | {u, v} ∈ E}; -if S = S 1 ∪ S 2 then ε S (v) := ε S1 (v) ∪ ε S2 (v); -if S = S 1 ∩ S 2 then ε S (v) := ε S1 (v) ∩ ε S2 (v); -if S = ¬S then ε S (v) := V \ ε S (v).

The modal logic EMLC consists of all the unary formulas that are built with the following grammar:

ϕ ::= C | ϕ ∧ ϕ | ¬ϕ | S ≥N ϕ,

where C ranges over node colors, S over modal parameters, and N over N. The semantics of the first four constructs is defined as expected, and for an undirected colored graph G = (V, E) and node v ∈ V , we have (G, v) |= S ≥N ϕ if and only if there exist at least N nodes u in ε S (v) such that (G, u) |= ϕ. Example D.2. On an undirected graph G = (V, E), the EMLC formula ¬e ≥2 ( e ≥3 Green) holds on a node v ∈ V if v has at least two nonadjacent nodes u (and since our graphs have no self-loops, v could be u) such that u has at least three green neighbors.

The following theorem is essentially a reformulation of (Lutz et al., 2001 , Theorem 1) to our context (Lutz et al. (2001) show this for FO 2 without counting quantifiers and for EMLC without counting, but an inspection of the proofs reveals that the result extends to counting quantifiers). Theorem D.3 (Lutz et al., 2001, Theorem 1) . For every EMLC formula, there exists an equivalent FOC 2 unary formula. Conversely, for every unary FOC 2 formula, there exists an equivalent EMLC formula.

In order to simplify the proof, we will use the following lemma. Lemma D.4. Let ϕ be an EMLC formula. Then there exists an EMLC formula ϕ equivalent to ϕ such that each modal parameter appearing in ϕ is one of the following: a) id, thus representing the current node; b) e, thus representing the neighbours of the current node; c) ¬e ∩ ¬id, thus representing the nodes distinct from the current node and that are not neighbours of the current node; d) id ∪ e, thus representing the current node and its neighbors; e) ¬id, thus representing all the nodes distinct from the current node: f) ¬e, thus representing the nodes that are not neighbours of the current node (note that this includes the current node);

In the following proof we will use the machinery introduced in Appendices C and D. We will also make use of a particular AC-GNN with L layers, which we call A L primes , that maps every node v in a graph G to a natural number representing the complete unravelling of v of depth L in G (note that we do not claim that this AC-GNN can be realized in practice, this construction is mostly for theoretical purposes). Let primes : N → N be the function such that primes(i) is the i-th prime number indexed from 0. For instance, we have that primes(0) = 2, primes(1) = 3, etc. Now consider the function f(•, •) that has as input a pair (c, X) where c ∈ N and X is a multiset of numbers in N, and produces a number in N as output, defined as follows

f(c, { {x 1 , x 2 , . . . , x k } }) = 2 c × k i=1 primes(x i + 1).

It is not difficult to prove that, as defined above, f(•, •) is an injective function. Thus using the results by Xu et al. (2019) (see the proof of their Theorem 3) we know that f can be used to implement the combine and aggregate operators of an AC-GNN such that for every graph G, after L layers, the color (natural number) assigned to every node in G has a one to one correspondence with the color assigned to that node in the L-th iteration of the WL test over G. We call this AC-GNN A L primes . Observation E.1. We note that Xu et al. (2019) also constructed an injective function that has (c, X) as inputs where c ∈ N and X is a multiset of elements in N (see their Lemma 5 and Corollary 6). Nevertheless we cannot directly use that construction as it assumes the existence of a fixed N such that the size of all multisets are bounded by N . This would put also a bound of N on the maximum number of neighbors in the input graphs. Thus we developed a new function (using an encoding based on prime numbers) to be able to deal with general graphs of unbounded degree.

Proof of Theorem 5.2. Let α be an FOC 2 unary formula, and let ϕ be an equivalent EMLC formula that uses only modal parameters of the form given by Lemma D.4. We construct an ACR-FR-GNN A ϕ capturing ϕ and hence α.

Let L be the quantifier depth of ϕ (i.e., the deepest nesting of S ≥N quantifiers). For a subformula ϕ of ϕ, we also define the nesting depth nd ϕ (ϕ ) of ϕ in ϕ to be the number of modal parameters under which ϕ is in ϕ. The first L -1 layers of A ϕ are the same as those of A L-1 primes , which do not use readouts. With Observation C.3 at hand and using the fact that the inverses of the aggregation and combination functions of A L-1 primes are computable, this ensures that, after L -1 layers, for any graph G and node v in G, we can compute from A L-1 primes (G, v) the unravelling Unr L-1 G (v). Thus, we can assume without loss of generality (by modifying the last combination function for instance), that after L -1 layers A ϕ computes Unr L-1 G (v) in every node v of G. We then use a readout whose output is a natural number representing the multiset { {Unr L-1 G (v) | v node in G} }; for instance, we can encode this multiset using the same technique that we use for A primes . Again, since this technique uses functions with computable inverses, we can assume without loss of generality that the output of this readout is actually the multiset { {Unr L-1 G (v) | v node in G} }. Finally, we use a final combination function COM (L) , that uses only the feature of the current node and the output of the readout-that is, the final feature of a node v is COM (L) (Unr L-1 G (v), { {Unr L-1 G (u) | u node in G} }). We now explain how we define COM (L) . By induction on the structure of ϕ, for every subformula ϕ of ϕ, we do the following: for every node v in G and every node u in Unr L-1 G (v) that is at depth (i.e., the distance from v) at most nd ϕ (ϕ ) in the tree Unr L-1 G (v), we will label u by either ϕ or by ¬ϕ . We do so to ensure that ( ) for every node v in G and every node u = (v, u 1 , . . . , u i ) in Unr L-1 G (v), we label u by ϕ if and only if (G, u i ) |= ϕ . We explain our labeling process by induction on the structure of ϕ, and one can easily check in each case that ( ) will hold by induction. Let v be a node in G and u be a node in Unr L-1 G (v) that is at depth at most nd ϕ (ϕ ) in the unravelling. Case 1. If ϕ is a color Col, we label u by ϕ if u is of that color, and by ¬ϕ otherwise.

Case 2. If ϕ is ϕ 1 ∧ ϕ 2 , then observe that we have nd ϕ (ϕ ) = nd ϕ (ϕ 1 ) = nd ϕ (ϕ 2 ), so that u is at depth at most both nd ϕ (ϕ 1 ) and nd ϕ (ϕ 2 ) in the unravelling Unr L-1 (v). Thus, we know that we have already labeled u by either ϕ 1 or ¬ϕ 1 , and also by either ϕ 2 or ¬ϕ 2 . We then label u by ϕ if u is already labeled by ϕ 1 and ϕ 2 , and we label it by ¬ϕ otherwise. Erdös-Renyi graphs These are random graphs in which one specifies the number N of nodes and the number M of edges. For this experiment we consider as extreme cases the case in which graphs contain the same number of nodes and edges and graphs in which the number of edges is twice the number of nodes.

Some statistics of the datasets are shown in Table 3 .

EXPERIMENTS FOR DENSE ERD ÖS-RENYI GRAPHS:

We also took a closer look at the performance for different connectivities of random graphs (Table 4 ). We define the set "Erdös-Renyi + k%" as a set of graphs in which the number of edges is k% larger than the number of nodes. For example, "Erdös-Renyi + 100%" contains random graphs in which the number of egdes doubles the number of nodes. We see a consistent improvement in the performance of AC-GNNs and GINs when we train and test them with more dense graphs and more layers (Table 4 ).

DATA FOR THE EXPERIMENT WITH CLASSIFIER α i (x) IN EQUATION ( 6)

For this case we only consider dense Erdös-Renyi synthetic graphs. For the train set we consider graphs with nodes varying from 40 to 50 nodes and edges from 280 to 350 and similarly for the first test set. For the bigger test set, we consider graphs with nodes from 51 to 60 with edges ranging from 360 and 480. For labeling we consider the following formulas (starting from α 0 (x) := Blue(x)):

α 1 (x) := ∃ [8,10] y α 0 (y) ∧ ¬E(x, y) , α 2 (x) := ∃ [10,20] y α 1 (y) ∧ ¬E(x, y) , α 3 (x) := ∃ [10,30] y α 2 (y) ∧ ¬E(x, y) .

The choices of the intervals for every classifier were for the pourpose of having approximately half of the nodes in the random graphs marked as true. Statistics of the datasets are shown in Table 5 .

:

For graph classification, which we do not consider in this paper, the classification function CLS inputs the multiset { {x(L) v | v ∈ V } }and outputs a class for the whole graph. Such a function is often called readout in previous work(Morris et al.,

2019;Xu et al., 2019). In this paper, however, we use the term readout to refer to intermediate global operations performed while computing features for nodes (see Section 5).

https://imfd.cl/en/

Figures and Tables:

Proposition C.4. Let G and G be two graphs with nodes v in G and v in G such that Unr L G (v) Unr L G (v ) for every L ∈ N. Then for any AC-GNN A, we have A(G, u) = A(G , u ).

Results on synthetic data for nodes labeled by classifier α(x) := Red(x) ∧ ∃y Blue(y)

Table content: <table><tr><td/><td/><td colspan="2">Line Train</td><td colspan="2">Line Test</td><td colspan="2">E-R Train</td><td colspan="2">E-R Test</td></tr><tr><td/><td/><td/><td colspan="3">same-size bigger</td><td/><td/><td colspan="2">same-size bigger</td></tr><tr><td/><td>AC-5</td><td>0.887</td><td colspan="2">0.886</td><td>0.892</td><td/><td>0.951</td><td>0.949</td><td>0.929</td></tr><tr><td/><td>AC-7</td><td>0.892</td><td colspan="2">0.892</td><td>0.897</td><td/><td>0.967</td><td>0.965</td><td>0.958</td></tr><tr><td/><td>GIN-5</td><td>0.861</td><td colspan="2">0.861</td><td>0.867</td><td/><td>0.830</td><td>0.831</td><td>0.817</td></tr><tr><td/><td>GIN-7</td><td>0.863</td><td colspan="2">0.864</td><td>0.870</td><td/><td>0.818</td><td>0.819</td><td>0.813</td></tr><tr><td/><td>ACR-1</td><td>1.000</td><td colspan="2">1.000</td><td>1.000</td><td/><td>1.000</td><td>1.000</td><td>1.000</td></tr><tr><td/><td>α1 Train</td><td colspan="2">α1 Test</td><td colspan="2">α2 Train</td><td colspan="2">α2 Test</td><td>α3 Train</td><td colspan="2">α3 Test</td></tr><tr><td/><td/><td colspan="2">same-size bigger</td><td/><td colspan="3">same-size bigger</td><td/><td colspan="2">same-size bigger</td></tr><tr><td>AC</td><td>0.839</td><td>0.826</td><td>0.671</td><td>0.694</td><td colspan="2">0.695</td><td>0.667</td><td>0.657</td><td>0.636</td><td>0.632</td></tr><tr><td>GIN</td><td>0.567</td><td>0.566</td><td>0.536</td><td>0.689</td><td colspan="2">0.693</td><td>0.672</td><td>0.656</td><td>0.643</td><td>0.580</td></tr><tr><td colspan="2">AC-FR-2 1.000</td><td>1.000</td><td>1.000</td><td>0.863</td><td colspan="2">0.860</td><td>0.694</td><td>0.788</td><td>0.775</td><td>0.770</td></tr><tr><td colspan="2">AC-FR-3 1.000</td><td>1.000</td><td>0.825</td><td>0.840</td><td colspan="2">0.823</td><td>0.604</td><td>0.787</td><td>0.767</td><td>0.771</td></tr><tr><td>ACR-1</td><td>1.000</td><td>1.000</td><td>1.000</td><td>0.827</td><td colspan="2">0.834</td><td>0.726</td><td>0.760</td><td>0.762</td><td>0.773</td></tr><tr><td>ACR-2</td><td>1.000</td><td>1.000</td><td>1.000</td><td>0.895</td><td colspan="2">0.897</td><td>0.770</td><td>0.800</td><td>0.799</td><td>0.771</td></tr><tr><td>ACR-3</td><td>1.000</td><td>1.000</td><td>1.000</td><td>0.903</td><td colspan="2">0.902</td><td>0.836</td><td>0.817</td><td>0.802</td><td>0.748</td></tr></table>



Table content: <table/>

use global-context aware local descriptors to classify objects in 3D point clouds, You et al. (2019) construct node features by computing shortest-path distances to a set of distant anchor nodes, and Haonan et al. (

Table content: <table/>

# Graphs Avg. # Nodes Avg. # Edges Avg. # Positive Synthetic data for the experiment with classifier α(x) := Red(x) ∧ ∃y Blue(y)

Table content: <table><tr><td/><td>Line train</td><td/><td colspan="2">5,000</td><td>75</td><td>74</td><td/><td>18</td><td/></tr><tr><td/><td>Line test</td><td/><td>500</td><td/><td>75</td><td>74</td><td/><td>18</td><td/></tr><tr><td/><td colspan="2">Line test bigger</td><td>500</td><td colspan="2">148</td><td>147</td><td/><td>36</td><td/></tr><tr><td/><td colspan="2">Erdös-Renyi train</td><td colspan="2">5,000</td><td>75</td><td>115</td><td/><td>18</td><td/></tr><tr><td/><td colspan="2">Erdös-Renyi test</td><td>500</td><td/><td>75</td><td>115</td><td/><td>18</td><td/></tr><tr><td/><td colspan="2">Erdös-Renyi test bigger</td><td>500</td><td colspan="2">148</td><td>226</td><td/><td>36</td><td/></tr><tr><td/><td colspan="3">Erdös-Renyi + 20%</td><td colspan="3">Erdös-Renyi + 50%</td><td colspan="3">Erdös-Renyi + 100%</td></tr><tr><td/><td>Train Acc.</td><td colspan="2">Test Acc.</td><td>Train Acc.</td><td colspan="2">Test Acc.</td><td>Train Acc.</td><td colspan="2">Test Acc.</td></tr><tr><td/><td/><td colspan="2">same-size bigger</td><td/><td colspan="2">same-size bigger</td><td/><td colspan="2">same-size bigger</td></tr><tr><td>AC-2</td><td>0.810</td><td>0.807</td><td>0.778</td><td>0.829</td><td>0.835</td><td>0.791</td><td>0.861</td><td>0.864</td><td>0.817</td></tr><tr><td>AC-5</td><td>0.940</td><td>0.937</td><td>0.901</td><td>0.975</td><td>0.971</td><td>0.958</td><td>0.994</td><td>0.994</td><td>0.993</td></tr><tr><td>AC-7</td><td>0.963</td><td>0.961</td><td>0.946</td><td>0.983</td><td>0.978</td><td>0.981</td><td>0.995</td><td>0.995</td><td>0.995</td></tr><tr><td>GIN-2</td><td>0.797</td><td>0.795</td><td>0.771</td><td>0.813</td><td>0.818</td><td>0.784</td><td>0.838</td><td>0.840</td><td>0.803</td></tr><tr><td>GIN-5</td><td>0.838</td><td>0.836</td><td>0.819</td><td>0.846</td><td>0.847</td><td>0.833</td><td>0.841</td><td>0.844</td><td>0.838</td></tr><tr><td>GIN-7</td><td>0.838</td><td>0.840</td><td>0.803</td><td>0.841</td><td>0.844</td><td>0.838</td><td>0.784</td><td>0.788</td><td>0.773</td></tr><tr><td>ACR-1</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td><td>1.000</td></tr></table>

Detailed results for Erdös-Renyi synthetic graphs with different connectivities

Table content: <table/>

Peer review (iclr2020):
The paper utilizes recent insights into the relationship between the Weisfeiler-Lehman (WL) test for checking graph isomorphism and Graph Neural Networks (GNNs) in order to characterize the class of node classifiers that can be captured by a specific GNN architecture, called aggregate-combine GNN (AC-GNN). The primary contribution of this work is the identification of the logical classifiers that can be represented within an AC-GNN, a fragment of first order logic called graded modal logic, as well as an extention of AC-GNNs with the ability to capture a strictly more expressive fragment of first order logic, called ACR-GNN. Both of these results are supported by formal proofs and derivations, adding not only theoretical value to the presented work, but also intuitive insights behind the reasons of the AC-GNN limitations. In addition, the presented experiments demonstrate the practical implications of the results, with the exception of some datasets where no significance difference in performance between AC-GNNs and ACR-GNNs was found. I believe that the contributions of this paper are significant. On the one hand, studying theoretical properties of GNNs facilitates their transparency and highlights the range of applications they can be used in. On the other hand, although the GNN variant introduced by the author is a special case of an existing class of GNNs, the motivation behind its introduction is different and follows naturally from the discussion within the paper. The fact that no actual difference in performance between AC-GNNs and ACR-GNNs was noticed in the only non-synthetic dataset used in the experiment should prompt the author to run experiments with more real life datasets, in order to empirically verify the results, but this is a minor point.