Title: Affinity-Aware Graph Networks
Venue: NeurIPS 2023
Accepted: yes
Keywords: graph neural networks

message passing

effective resistance

hitting time
Abstract: Graph Neural Networks (GNNs) have emerged as a powerful technique for learning on relational data. Owing to the relatively limited number of message passing steps they perform—and hence a smaller receptive field—there has been significant interest in improving their expressivity by incorporating structural aspects of the underlying graph. In this paper, we explore the use of affinity measures as features in graph neural networks, in particular measures arising from random walks, including effective resistance, hitting and commute times. We propose message passing networks based on these features and evaluate their performance on a variety of node and graph property prediction tasks. Our architecture has low computational complexity, while our features are invariant to the permutations of the underlying graph. The measures we compute allow the network to exploit the connectivity properties of the graph, thereby allowing us to outperform relevant benchmarks for a wide variety of tasks, often with significantly fewer message passing steps. On one of the largest publicly available graph regression datasets, OGB-LSC-PCQM4Mv1, we obtain the best known single-model validation MAE at the time of writing.

Affinity-Aware Graph Networks

Abstract

Graph Neural Networks (GNNs) have emerged as a powerful technique for learning on relational data. Owing to the relatively limited number of message passing steps they perform-and hence a smaller receptive field-there has been significant interest in improving their expressivity by incorporating structural aspects of the underlying graph. In this paper, we explore the use of affinity measures as features in graph neural networks, in particular measures arising from random walks, including effective resistance, hitting and commute times. We propose message passing networks based on these features and evaluate their performance on a variety of node and graph property prediction tasks. Our architecture has low computational complexity, while our features are invariant to the permutations of the underlying graph. The measures we compute allow the network to exploit the connectivity properties of the graph, thereby allowing us to outperform relevant benchmarks for a wide variety of tasks, often with significantly fewer message passing steps. On one of the largest publicly available graph regression datasets, OGB-LSC-PCQM4Mv1, we obtain the best known single-model validation MAE at the time of writing.

1 Introduction

Graph Neural Networks (GNNs) constitute a powerful tool for learning meaningful representations in non-Euclidean domains. GNN models have achieved significant successes in a wide variety of node prediction [18, 31], link prediction [55, 52], and graph prediction [11, 50] tasks. These tasks naturally emerge in a wide range of applications, including autonomous driving [8], neuroimaging [36], combinatorial optimization [14, 34], and recommender systems [51], while they have enabled significant scientific advances in the fields of biomedicine [45], structural biology [23], molecular chemistry [42] and physics [3].

Despite the predictive power of GNNs, it is known that the expressive power of standard GNNs is limited by the 1-Weisfeiler-Lehman (1-WL) test [47]. Intuitively, GNNs possess at most the same power in terms of distinguishing between non-isomorphic (sub-)graphs, while having the added benefit of adapting to the given data distribution. For some architectures, two nodes with different local structures have the same computational graph, thus thwarting distinguishability in a standard GNN. Even though some attempts have been made to address this limitation with higher-order GNNs [32], most traditional GNN architectures fail to distinguish between such nodes.

A common approach to improving the expressive power of GNNs involves encoding richer structural/positional properties. For example, distance-based approaches form the basis for works such as

Position-aware Graph Neural Networks [52], which capture positions/locations of nodes with respect to a set of anchor nodes, as well as Distance Encoding Networks [29], which use the first few powers of the normalized adjacency matrix as node features associated with a set of target nodes.

Here, we take an approach that is inspired by this line of work but departs from it in some crucial ways: we seek to capture both distance and connectivity information using general-purpose node and edge features without the need for specifying any anchor or target nodes.

Contributions : We propose the use of affinity metrics as features in a GNN. Specifically, we consider statistics that arise from random walks in graphs, such as hitting time and commute time between pairs of vertices (see Sections 3.1 and 3.2). We present a means of incorporating these statistics as scalar edge features in a message passing neural network (MPNN) [15] (see Section 3.4). In addition to these scalar features, we present richer vector-valued resistive embeddings (see Section 3.3), which can be incorporated as node or edge feature vectors in the network. Resistive embeddings are a natural way of embedding each node into Euclidean space such that the squared L 2 -distance between nodes recovers the commute time. We show that such embeddings can be incorporated into MPNNs, even for larger graphs , by efficiently approximating them using sketching and dimensionality reduction techniques and prove a novel additive approximation for hitting time (see Section 4).

Moreover, we evaluate our networks on a number of benchmark datasets of diverse scales (see Section 5). First, we show that our networks outperform other baselines on the PNA dataset [9], which includes 6 node and graph algorithmic tasks, showing the ability of affinity measures to exploit structural properties of graphs. We also evaluate the performance on a number of graph and node tasks for datasets in the Open Graph Benchmark (OGB) collection [21], including molecular and citation graphs. In particular, our networks with scalar effective resistance edge features achieve the state of the art on the OGB-LSC PCQM4Mv1 dataset, which was featured in a KDD Cup 2021 competition for large scale graph representation learning.

Finally, we provide intuition for why affinity-based measures are fundamentally different from aforementioned distance-based approaches (see Section 3.5) and bolster it with detailed theoretical and empirical results (see Appendix D) showing favorable results for affinity-based measures.

2 Related Work

Our work builds upon a wealth of graph theoretic and graph representation learning works, while we focus on a supervised, inductive setting.

Even though GNN architectures were originally classified as spectral or spatial, we abstain from this division as recent research has demonstrated some equivalence of the graph convolution process regardless of the choice of convolution kernels [2, 7]. Spectrally-motivated methods are often theoretically founded on the eigendecomposition of the graph Laplacian matrix (or an approximation thereof) and, hence, corresponding convolutions capture different frequencies of the graph signal. Early works in this space include ChebNet [10] and its more efficient 1-hop version by Kipf et al. [24], which offers a linear function on the graph Laplacian spectrum. Levie et al. [28] proposed CayleyNets, an alternative rational filter.

Message passing neural networks (MPNNs) [15] perform a transformation of node and edge representations before and after an arbitrary aggregator (e.g. sum ). Graph attention networks (GATs) [44] aimed to augment the computations of GNNs by allowing graph nodes to 'attend' differently to different edges, inspired by the success of transformers in NLP tasks. One of the most relevant works was proposed by Beaini et al. [4], i.e. directional graph networks (DGN). DGN uses the gradients of the low-frequency eigenvectors of the graph Laplacian, which are known to capture key information about the global structure of the graph and prove that the aggregators they construct using these gradients lead to more discriminative models than standard GNNs according to the 1-WL test. Prior work [32] used higher-order ( k -dimensional) GNNs, based on k -WL, and a hierarchical variant and proved theoretically and experimentally the improved expressivity in comparison to other models.

Other notable works include Graph Isomorphism Networks (GINs) [47], which represent a simple, maximally-powerful GNN over discrete-featured inputs. Hamilton et al. [18] proposed a method to construct node representations by sampling a fixed-size neighborhood of each node and then performing aggregation over it, which led to impressive performance on large-scale inductive benchmarks. Bouritsas et al. [6] use topologically-aware message passing to detect and count graph substructures,

while Bodnar et al. [5] propose a message-passing procedure on cell complexes motivated by a novel color refinement algorithm that proves to be powerful for molecular benchmarks. Meanwhile, Horn et al. [19] propose a layer based on persistent homology to incorporate global topological information.

Expressivity. Techniques that improve a GNN's expressive power largely fall under three broad directions. While we focus on the feature-based direction in this paper, we also acknowledge that it in no way compels the GNN to use the additional provided features. Hence, we briefly survey the other two, as an indication of the future research in affinity-based GNNs we hope this work will inspire.

Another avenue involves modulating the message passing rule to make advantage of the desired computations. Popular recent examples of this include DGN [4] and LSPE [12]. DGNs leverage the graph's Laplacian eigenvectors, but they do not merely use them as input features; instead, they define a directional vector field based on the eigenvectors, and use it explicitly to anisotropically aggregate neighbourhoods. LSPE features a 'bespoke pipeline' for processing positional inputs.

The final direction is to modulate the graph over which messages are passed, usually by adding new nodes that correspond to desired substructures. An early proponent of this is the work of [32], which explicitly performs message passing over k -tuples of nodes at once. Recently, scalable efforts in this direction focus on carefully chosen substructures, e.g., junction trees [13], cellular complexes [5].

3 Affinity Measures and GNNs

3.1 Random Walks, Hitting and Commute Times

Let G = ( V, E ) be a graph of vertices V and edges E ⊆ V × V between them. We define several natural properties of a graph that arise from a random walk. A random walk on G starting from a node u is a Markov chain on the vertex set V such that the initial vertex is u , and at each time step, one moves from the current vertex to a neighbor, chosen with probability proportional to the weight of outgoing edges. We will use π to denote the stationary distribution of this Markov Chain. For random walks on weighted, undirected graphs, we know that π u = d u 2 M , where d u is the weighted degree of node u , and M is the sum of edge weights.

The hitting time H uv from u to v is defined as the expected number of steps for a random walk starting at u to hit v . We can also define the commute time between u and v as K uv = H uv + H vu , the expected round-trip time for a random walk starting at u to reach v and then return to u .

3.2 Effective Resistance

A closely related quantity is the measure of effective resistances in undirected graphs. This quantity corresponds to the effective resistance if the whole graph was replaced with a circuit where each edge becomes a resistor with resistance equal to the reciprocal of its weight. We will use Res ( u, v ) to denote the effective resistance between nodes u and v . For undirected graphs, it is known that [30] the effective resistance is proportional to the commute time, Res ( u, v ) = 1 2 M K uv .

Our broad goal is to incorporate effective resistances and hitting times as edge features in an MPNN. In Section 3.4, we will show they can provably improve MPNNs' expressivity.

3.3 Resistive Embeddings

Effective resistances allow us to define the resistive embedding , a mapping that associates each node v of a graph G = ( V, E, W ) , where W are the non-negative edge weights, with an embedding vector. Before we specify the resistive embedding, we define a few terms. Let L = D -A be the graph Laplacian of G , where D ∈ R n × n is the diagonal matrix containing the weighted degree of each node and A ∈ R n × n is the adjacency matrix, whose ( i, j ) th entry is equal to the edge weight between i and j , if it exists; and 0 otherwise. Let B be the m × n edge-node incidence matrix, where | V | = n and | E | = m , defined as follows: The i -th row of B corresponds to the i -th edge e i = ( u i , v i ) of G and has a +1 in the u i -th column and a -1 in the v i -th column, while all other entries are zero. Finally we will use C ∈ R m × m to denote the conductance matrix, which is a diagonal matrix with C ii being the weight of i th edge. It is easy to verify that B T CB = L . Even though L is not invertible, its null-space consists of the indicator vectors for every connected component of G . For example, if G is

connected, then L 's nullspace is spanned by the all1 's vector 2 . Hence, for any vector x orthogonal to all1 's, L · L † x = x , where L † is the pseudo-inverse.

We can express effective resistance between any pair of nodes using Laplacian matrices [30] as Res ( u, v ) = ( 1 u -1 v ) T L † ( 1 u -1 v ) , where 1 v is an n -dimensional vector specifying the indicator for node v . We are now ready to define the resistive embedding.

Definition 3.1. (Effective Resistance Embedding) r v = C 1 / 2 BL † 1 v .

A key property is that the effective resistance between two nodes in the graph can be obtained easily from the distance between their corresponding embeddings (see proof in Appendix C):

Lemma 3.2. For any pair of nodes u, v , we have ∥ r u -r v ∥ 2 2 = Res ( u, v ) .

One can easily check that any rotation of r also satisfies Lemma 3.2, since rotations preserve Euclidean distances; more generally, if U is an orthonormal matrix, then U r is also a valid resistive embedding. This poses a challenge if we want to use the resistive embeddings as node or edge features: we want a way to enforce that a (G)NN using them will do so in a way that is invariant or equivariant to any rotations of the embeddings. In our current work, we rely on data augmentation: at every training iteration, we apply random rotations to the input ER embeddings.

Remark. While data augmentation is a popular approach for promoting invariant and equivariant predictions, it is only hinting to the network that such predictions are favourable. It is also possible, in the spirit of the geometric deep learning blueprint [7], to combine ER embeddings with an O ( n ) -equivariant GNN, which rigorously enforces rotational equivariance. A popular approach to building equivariant GNNs has been proposed by [39], though it focuses on the full Euclidean group E ( n ) rather than O ( n ) . We leave this exploration to future work.

Definition 3.3. Let p := ∑ u π u r u be the mean of effective resistance embedding.

We might view p as a 'weighted mean' 3 of r . We will define the hitting time radius, H max , of a given graph as the maximum hitting time between any two nodes:

Definition 3.4 (Hitting Time Radius) . H max := max u,v H u,v .

We will need the following to bound the hitting times we computed:

Lemma 3.5. For any node u , ∥ r u -p ∥ 2 ≤ H max M .

The proof follows easily from the fact that p is a convex combination of all r 's and Jensen's inequality.

3.4 Incorporating Features into MPNNs

We reiterate that our main aim is to demonstrate (theoretically and empirically) that there are good reasons to incorporate affinity-based measures into GNN computations.

In the simplest instance, a method that improves a GNN's expressive power may compute additional features (positional or structural) which would assist the GNN in discriminating between examples it otherwise wouldn't (easily) be able to. These features are then appended to the GNN's inputs for further processing. For example, it has been shown that endowing nodes with a one-hot based identity is already sufficient for improving expressive power [33]; this was then relaxed to any randomly-sampled scalar feature by [38]. It is, of course, possible to create dedicated features that even count substructures of interest [6]. Further, the adjacency information can be factorised [37] or eigendecomposed [12] to provide useful structural embeddings for the GNN.

We will focus our attention on exactly this class of methods, as it is a lightweight and direct way of demonstrating improvements from these computations. Hence, our baselines will all be instances of the MPNN framework [15], which we will attempt to improve by endowing them with affinity-based features. We start by theoretically proving that these features indeed improve expressive power:

Theorem 3.6. MPNNs that make use of any one of (a.) effective resistances, (b.) hitting times, (c.) resistive embeddings are strictly more powerful than the WL-1 test.

Proof. Since the networks in question arise from augmenting standard MPNNs with additional node/edge features, we have that these networks are at least as powerful as the 1-WL test.

In order to show that these networks are strictly more powerful than the 1-WL test, it suffices to show the existence of a graph for which our affinity measure based networks can distinguish between certain nodes that a standard GNN (limited by the 1-WL test) cannot.

We present an example of a 3-regular graph on 8 nodes in Figure 1. It is well-known that a standard GNN that is limited by the 1-WL test cannot distinguish any pair of nodes in a regular graph, as the computation tree rooted at any node in the graph looks identical. However, there are three isomorphism classes of nodes in the above graph (denoted by different colors), namely, V 1 = { 1 , 2 } , V 2 = { 3 , 4 , 7 , 8 } , and V 3 = { 5 , 6 } .

̸

We now show that GNNs with affinity based measures can distinguish between a node in V i and a node in V j , for i = j . We note that the hitting time from a to b depends only on the isomorphism classes of a and b . Thus, we write r i,j as the effective resistance between a node in V i and a node in V j . Note that r i,j = r j,i , and it is easy to verify that:

Hence, it follows that in a message passing step of an MPNN that uses effective resistances, vertices in V 1 , V 2 , and V 3 will aggregate feature multisets { r 1 , 1 , r 1 , 2 , r 1 , 2 } = { 2 / 3 , 185 / 336 , 185 / 336 } , { r 2 , 1 , r 2 , 2 , r 2 , 3 } = { 185 / 336 , 15 / 28 , 209 / 336 } , and { r 2 , 3 , r 2 , 3 , r 3 , 3 } = { 209 / 336 , 209 / 336 , 4 / 7 } , respectively, all of which are all distinct multisets. Hence, such an MPNN can distinguish nodes in V i and V j , i = j for a suitable aggregation function.

̸

If, instead of effective resistance, we use hitting time features or resistive embeddings, our results still hold. This is because, as we showed previously, the effective resistance between nodes is a function of the two hitting times in either direction, as well as of the resistive embeddings of the nodes. In other words, if either hitting time features or resistive embeddings are used as input features for an MPNN, this MPNN would be able to compute the effective resistance features by applying an appropriate function (e.g., Lemma 3.2 for the case of resistive embeddings). Having computed these features, the MPNN can distinguish any two graphs that the MPNN with effective resistance features can.

3.5 Effective Resistance vs. Shortest Path Distance

It is interesting to ask how effective resistances compare with shortest path distances (SPDs) in GNNs, given the plethora of recent works that make use of SPDs (e.g., [49, 52, 29]). The most direct comparison of our effective resistance-based MPNNs would be to use SPDs as edge features in the MPNNs. However, note that SPDs along graph edges are trivial (unlike effective resistances, which incorporate useful information about the global graph structure).

An alternative to edge features would be to use (a) SPDs to a small set of anchor nodes as features in an MPNN (e.g., P-GNN [52]) or (b) a dense featurization incorporating shortest paths between all pairs of nodes (e.g., the dense attention mechanism in Graphormer [49]). We remark that the latter approach typically incurs an O ( n 2 ) overhead, which our MPNN-based approach avoids.

We empirically compare our models to MPNNs that use approaches (a) and (b). Results on the PNA dataset show that our effective resistance-based GNNs outperform these approaches. Furthermore, we complement the empirical results with a theoret-

ical result (see Appendix D) showing that under a limited number of message-passing steps, effective resistance features can allow one to distinguish structures that cannot be done using shortest path features. This provides some insight into why effective resistances can capture structure in GNNs

that SPDs are unable to. We describe the result below. We point the reader to Appendix D for a proof of Theorem D.1.

4 Efficient Computation of Affinity Measures

In order to use our features, it is important that they be computable efficiently. In this section, we show how to compute or approximate the various random walk-based affinity measures. We share the results of our method on large-scale graphs in Section 5.4.

4.1 Reducing Dimensionality of Resistive Embeddings

Given their higher dimensionality, we might find it beneficial to use resistive embeddings as GNN features instead of the effective resistance. Now, the difficulty with using resistive embeddings directly as GNN features is the fact that the embeddings have dimension m , which can be quite large, e.g., up to n 2 for dense graphs. It was shown by [41] that one can reduce the dimensionality of the embed-

ding while approximately preserving Euclidean distances. The idea is to use a random projection via a constructive version of the Johnson-Lindenstrauss Lemma:

Lemma 4.1 (Constructive Johnson-Lindenstrauss) . Let x 1 , x 2 , . . . , x n ∈ R d be a set of n points; and let α 1 , α 2 , . . . , α m ∈ R n be fixed linear combinations. Suppose Π is a k × d matrix whose entries are chosen i.i.d. from a Gaussian N (0 , 1) and consider ̂ x i := 1 √ k Π x i .

Then, it follows that for k ≥ C log( mn ) /ϵ 2 , with probability 1 -o (1) , we have ∥ ∑ j α i,j ̂ x j ∥ 2 = (1 ± ϵ ) ∥ ∑ j α i,j x j ∥ 2 2 for every 1 ≤ i ≤ m .

Using triangle inequality, we can see that the inner products between fixed linear combinations are preserved up to some additive error:

Corollary 4.2. For any fixed vectors α, β ∈ R n , if we let X := ∑ i α i x i , ̂ X := ∑ i α i ̂ x i and similarly Y := ∑ i β i x i , ̂ Y := ∑ i β i ̂ x i ; then:

Therefore, we can choose a desired ϵ > 0 and k = O (log( n ) /ϵ 2 ) and instead use ̂ r : V → R k as the embedding, where

for a randomly chosen k × d matrix Π whose entries are i.i.d. Gaussians from N (0 , 1) . Then, by lemma 3.2 and lemma 4.1, we have that for every edge ( u, v ) ∈ E ,

with probability at least 1 -1 n 2 .

So the computation of random embeddings, ̂ r , requires solving O (( n + m ) log n/ϵ 2 ) many Laplacian linear systems. By using one of the nearly linear time Laplacian solvers [26], we can compute the random embeddings in the near-linear time. Hence the total running time becomes O ( ( n + m ) log 3 / 2 n poly log log n/ϵ 2 ) .

4.2 Fast Computation of Hitting Times

Note that it is not clear whether there is a fast method for computing hitting times similar to commute times / effective resistances. The naive approach involves solving a linear system for each edge, resulting in a running time of at least Ω( nm ) , which is prohibitive. One of our technical contributions in this paper is a novel method for fast computation of hitting times. In particular, we will show how to use the approximate effective resistance embeddings, ̂ r , to obtain an estimate for hitting times with additive error. Proofs can be found in Appendix C.

Let ̂ p := ∑ u π u ̂ r u . Just like ̂ r being an approximation of r , ̂ p is an approximation of p . Consider the following quantity:

We will use this quantity as an approximation of H u,v . In the following part, we will bound the difference between H u,v and ̂ H u,v . Our starting point will be expressing H u,v in terms of the effective resistance embeddings.

Given Lemma 4.3, we can establish the desired additive approximation property of ̂ H u,v .

5 Experiments

As previously discussed, our empirical evaluation seeks to show benefits from endowing standard expressive GNNs with additional affinity-based features. All architectures we experiment with will therefore conform to the MPNN blueprint [15], which we describe in Appendix B. When relevant, we may also note results that a particular strong baseline (e.g., DGN [4], Graphormer [48]) achieves on a dataset of interest. Note that these baselines modulate the message passing procedure rather than appending features and are hence a different category to our method-their performance is provided for indicative reasons only. Where appropriate, we use ' DGN (features) ' to refer to an MPNN that uses eigenvector flows as additional edge features, without modulating the mechanism. We also use 'ER (GNN)' to refer to the use of scalar effective resistance features.

5.1 PNA dataset

The first dataset we explore is the PNA dataset [9], which captures a multimodal setting. This consists of a collection of node tasks , i.e. (1) Single-source shortest paths, (2) Eccentricity and (3) Laplacian features, as well as graph tasks , i.e. (4) Connectivity, (5) Diameter and (6) Spectral radius. PNA dataset is a set of structured graph tasks, and it complements our other datasets. As we can see in Table 1, even adding a single feature, effective resistance (ER GNN), yields the best average score compared to other models. Using hitting times as edge features improve upon effective resistances. However once we combine all ER features, which include effective resistances, hitting times as well as node and edge embeddings, we get the best scores. On these structured tasks, we can see that the affinity based measures provide a significant advantage.

We also note that incorporating affinity measures as additional features to the 'DGN (features)' MPNN model also provides improvement over the standard 'DGN (features)' baseline. We have accordingly included the results of a model that uses (both node and edge) ER embeddings and hitting time features along with DGN features.

5.2 Small molecule classification: ogbg-molhiv

The ogbg-molhiv dataset is a molecular property prediction dataset comprised of molecular graphs without spatial information (such as atom coordinates). Each graph corresponds to a molecule, with nodes representing atoms and edges representing chemical bonds. Each node has an associated 9 -dimensional feature, containing atomic number and chirality, as well as other additional atom features such as formal charge and whether the atom is in the ring or not. The goal is to predict whether a molecule inhibits HIV virus replication or not. Our results are given in Table 2.

On this dataset, effective resistances provide an improvement over the standard MPNN. We achieve the best performance using ER node embeddings and hitting times with random rotations. With these features, our network achieves 79 . 13 % ± 0 . 358 test accuracy, which is close to DGN.

Note that we provide GSN [6], HIMP [13], and GRWNN [35] as baselines. HIMP and GSN exploit structural/toplogical properties; the former is specifically designed for learning on molecular graphs, while the latter incorporates graph substructures. In addition to the comparison to GSN as a baseline, we additionally train an MPNN with substructure counts as additional features, namely counts of k -cycles ( k = 3 , 4 ) and k -paths ( k = 3 ) as motifs (in the same vein as the approach of [6], obtaining a test ROC-AUC of 76.41% ± 0.42 (note that this improves on the base MPNN test ROC-AUC of 74.67% ± 0.19). Furthermore, if we also add affinity measures (in addition to the aforementioned substructure counts), we obtain a test ROC-AUC of 78.50% ± 0.51. This suggests that affinity measures provide a performance or expressivity lift beyond that of substructure counts.

5.3 Multi-task molecular classification: ogbg-molpcba

The ogbg-molpcba dataset comprises molecular graphs without spatial information (such as atom coordinates). The aim is to classify them across 128 different biological activities. We follow the baseline MPNN architecture from [17], including the use of Noisy Nodes.

Mirroring the evaluation protocol of [17], Table 3 compares the performance of incorporating ER and hitting time (HT) features into the baseline MPNN models with Noisy Nodes, at various depths. What can be noticed is that models utilising affinity-based features are capable of reaching as well as exceeding peak test performance (in terms of mean average precision). However, what's important is the effect of these features at lower depths: it is possible to achieve comparable or better levels of performance with half the layers , when utilising ER or HT features. This result illustrates the potential benefit affinity-based computations can have on molecular benchmarks, especially when no spatial geometry is provided as input.

5.4 Scaling to larger graphs

Next, we present results on ogbn-arxiv and ogbn-mag, transductive datasets with large graphs.

5.4.1 Large citation network: ogbn-arxiv

Most expressive GNNs that rely on computation of structural features have not been scaled beyond small molecular datasets (such as the ones discussed in prior sections). This is due to the fact that computing them requires (time or storage) complexity which is at least quadratic in the graph size-making them inapplicable even for modest-sized graphs. This is, however, not the case for our proposed affinity-based metrics. We demonstrate this by scalably computing them on a larger-scale node classification benchmark, ogbn-arxiv (a citation network with the goal of predicting the arXiv category of each paper). ogbn-arxiv has 169,343 nodes and 1,166,243 edges, making quadratic approaches infeasible. Using the combinatorial multigrid preconditioner [27, 25], we constructed the effective resistances on this graph in an hour on a standard MacBook Pro 2019 laptop.

Table 4: ogbg-molpcba performance compared to other baselines.

Model

Test Mean AP

GCN (w/ virtual node)

GIN (w/ virtual node)

DeeperGCN

HIMP

MPNN + NN + HT (ours)

24.24

±

0.34

27.03

27.81

27.39

±

±

±

28.32

%

0.23

0.38

0.17

±

0.13

Table 2: Test % AUC-ROC on the MolHIV dataset. Our results are averaged over five seeds.

As MPNN models overfit this transductive dataset quite easily, the dominant approach to tackling it are graph attention networks (GATs) [44]. Accordingly, we trained a simple four-layer GAT on this dataset, achieving 72.02% ± 0.05 test accuracy. This compares with 71.97% ± 0.24 reported for a related attentional baseline on the leaderboard [54], indicating that our baseline performance is relevant.

ER embeddings on ogbn-arxiv need to be exceptionally high-dimensional to achieve accurate ER estimates ( ∼ 11,000 dimensions), hence we were unable to use them here. However, incorporating ER scalar features into our GAT model yielded a statistically-significant improvement of 72.14% ± 0.03 test accuracy. Hitting time features improve this result further to 72.25 % ± 0.04 test accuracy. This demonstrates that our affinity-based metrics can yield useful improvements even on larger scale graphs, which are traditionally out of reach for methods like DGN [4] due to computational complexity limitations.

Reliable global leaderboarding with respect to ogbn-arxiv is difficult, as state-of-the art approaches rely either on privileged information (such as raw text of the paper abstracts), incorporating node labels as features [46], post-processing the predictions [22], or various related tricks [46]. With that in mind, we report for convenience that the current state-of-the-art performance for ogbn-arxiv without using raw text is 76.11% ± 0.09 test accuracy, achieved by GIANT-XRT+DRGAT.

5.4.2 Heterogeneous citation network: ogbn-mag

We additionally present experiments on an even larger-scale network, ogbn-mag, a heterogeneous network derived from the Microsoft Academic Graph whose nodes consist of four types of entities (papers, authors, institutions, fields of study) and whose edges capture directed relations (e.g., a paper cites another paper, an author is affiliated with an institution, etc.). The task is to predict the venues of papers. The ogbn-mag dataset consists of 111,059,956 nodes and 1,615,685,872 edges.

As in the case of ogbn-arxiv, ER embeddings require an exceptionally high number of dimensions, given the large size of the network. Hence, we incorporated ER scalar features. As a baseline, we used a GraphSAINT [53] model that uses R-GCN aggregation [40], which reported a test accuracy of 47.51% ± 0.22 in a relevant leaderboard entry. After incorporating ER scalar features into the same model, we obtained a statistically-signficant improvement of 47.99 % ± 0.23 in test accuracy.

5.5 Large scale graph regression: OGB-LSC PCQM4Mv1

We finally include experimental results for one of the largest-scale publicly available graph regression tasks: the PCQM4Mv1 dataset from the OGB Large Scale Challenge [20]. PCQM4M is a quantum chemistry dataset spanning 4 million small molecules, with a task to predict the HOMO-LUMO gap, an important quantum-chemical property. It is anticipated that structural features such as ER could be of great help on this task, as the v1 version of it is provided without any structural information, and the molecule's geometry is assumed critical for predicting the gap. We report the single-model validation performance on this dataset, in line with previous works [17, 48, 1].

PCQM4Mv1 comprises molecular graphs which consist of bonds and atom types, and no 3D or 2D coordinates. We reuse the experimental setup and architecture from [17], with only one difference: appending the effective resistance to the edge features. Additionally, we compare against an equivalent model

which uses molecular conformations estimated by RDKit as an additional feature. This gives us a baseline which leverages an explicit estimate of the molecular geometry.

Our results are summarised in Table 5. We once again see a powerful synergy of effective resistanceendowed GNNs and Noisy Nodes [17], allowing us to significantly reduce the number of layers (to 32) and outperform the 50-layer MPNN result in [17]. Further, we improve on the singlemodel performance of both the Graphormer [48] (which won the original PCQM4M contest after ensembling), and an equivalent model to ours which uses molecular conformers from RDKit. This illustrates how ER features can be competitive in geometry-relevant tasks even against features that inherently encode an estimate of the molecule's spatial geometry.

Lastly, we remark that, to the best of our knowledge, our result is the best published single-model result on the large-scale PCQM4M-v1 benchmark to date, and the only single model result with validation MAE under 0.120. We hope this will inspire future investigation on affinity-related GNNs for molecular tasks, especially in settings where spatial geometry is not reliably available.

6 Conclusions

In this paper, we proposed a message passing network based on random walk based affinity measures.

We believe that the comprehensive theoretical and practical results presented in our paper have solidified affinity-based computations as a strong component of a graph representation learner's toolbox. Our proposal carefully balances theoretical expressive power, empirical performance, and scalability to large graphs. While adding affinity measures as node/edge features provides the most direct route to incorporating them in a graph network, an interesting future direction would be to explore variants of GNN message functions that explicitly make use of affinity-based computations .

Acknowledgments and Disclosure of Funding

We would like to specially thank Jonathan Godwin for his extensive support in setting up the codebase for large-scale GNNs, and Gabriele Corso for his support in setting up the PNA dataset. Additionally, we would like to thank Beatrice Bevilacqua, Wilfried Bounsi, Larisa Markeeva and Pete Battaglia for reviewing the paper prior to submission.

This research was funded by Google Research and Google DeepMind.

References

Ravichandra Addanki, Peter W Battaglia, David Budden, Andreea Deac, Jonathan Godwin, Thomas Keck, Wai Lok Sibon Li, Alvaro Sanchez-Gonzalez, Jacklynn Stott, Shantanu Thakoor, et al. Large-scale graph representation learning with very deep gnns and self-supervision. arXiv preprint arXiv:2107.09422 , 2021.

Muhammet Balcilar, Guillaume Renton, Pierre Héroux, Benoit Gaüzère, Sébastien Adam, and Paul Honeine. Analyzing the expressive power of graph neural networks in a spectral perspective. In International Conference on Learning Representations , 2021.

Victor Bapst, Thomas Keck, A Grabska-Barwi´ nska, Craig Donner, Ekin Dogus Cubuk, Samuel S Schoenholz, Annette Obika, Alexander WR Nelson, Trevor Back, Demis Hassabis, et al. Unveiling the predictive power of static structure in glassy systems. Nature Physics , 16(4):448454, 2020.

Dominique Beaini, Saro Passaro, Vincent Létourneau, Will Hamilton, Gabriele Corso, and Pietro Liò. Directional graph networks. In International Conference on Machine Learning (ICML) , pages 748-758. PMLR, 2021.

Cristian Bodnar, Fabrizio Frasca, Nina Otter, Yu Guang Wang, Pietro Liò, Guido F Montufar, and Michael Bronstein. Weisfeiler and lehman go cellular: Cw networks. Advances in Neural Information Processing Systems , 34, 2021.

Giorgos Bouritsas, Fabrizio Frasca, Stefanos Zafeiriou, and Michael M Bronstein. Improving graph neural network expressivity via subgraph isomorphism counting. arXiv preprint arXiv:2006.09252 , 2020.

Michael M Bronstein, Joan Bruna, Taco Cohen, and Petar Veliˇ ckovi´ c. Geometric deep learning: Grids, groups, graphs, geodesics, and gauges. arXiv preprint arXiv:2104.13478 , 2021.

Siheng Chen, Sufeng Niu, Tian Lan, and Baoan Liu. Pct: Large-scale 3d point cloud representations via graph inception networks with applications to autonomous driving. In 2019 IEEE international conference on image processing (ICIP) , pages 4395-4399. IEEE, 2019.

Gabriele Corso, Luca Cavalleri, Dominique Beaini, Pietro Liò, and Petar Veliˇ ckovi´ c. Principal neighbourhood aggregation for graph nets. Advances in Neural Information Processing Systems (NeurIPS) , 33, 2020.

Michaël Defferrard, Xavier Bresson, and Pierre Vandergheynst. Convolutional neural networks on graphs with fast localized spectral filtering. Advances in neural information processing systems , 29:3844-3852, 2016.

David K Duvenaud, Dougal Maclaurin, Jorge Iparraguirre, Rafael Bombarell, Timothy Hirzel, Alan Aspuru-Guzik, and Ryan P Adams. Convolutional networks on graphs for learning molecular fingerprints. Advances in Neural Information Processing Systems , 28:2224-2232, 2015.

Vijay Prakash Dwivedi, Anh Tuan Luu, Thomas Laurent, Yoshua Bengio, and Xavier Bresson. Graph neural networks with learnable structural and positional representations. arXiv preprint arXiv:2110.07875 , 2021.

Matthias Fey, Jan-Gin Yuen, and Frank Weichert. Hierarchical inter-message passing for learning on molecular graphs. arXiv preprint arXiv:2006.12179 , 2020.

Maxime Gasse, Didier Chételat, Nicola Ferroni, Laurent Charlin, and Andrea Lodi. Exact combinatorial optimization with graph convolutional neural networks. arXiv preprint arXiv:1906.01629 , 2019.

Justin Gilmer, Samuel S Schoenholz, Patrick F Riley, Oriol Vinyals, and George E Dahl. Neural message passing for quantum chemistry. In International Conference on Machine Learning (ICML) , pages 1263-1272. PMLR, 2017.

Jonathan Godwin, Thomas Keck, Peter Battaglia, Victor Bapst, Thomas Kipf, Yujia Li, Kimberly Stachenfeld, Petar Veliˇ ckovi´ c, and Alvaro Sanchez-Gonzalez. Jraph: A library for graph neural networks in jax., 2020.

Jonathan Godwin, Michael Schaarschmidt, Alexander L Gaunt, Alvaro Sanchez-Gonzalez, Yulia Rubanova, Petar Veliˇ ckovi´ c, James Kirkpatrick, and Peter Battaglia. Simple GNN regularisation for 3d molecular property prediction and beyond. In The Tenth International Conference on Learning Representations , 2022.

William L Hamilton, Rex Ying, and Jure Leskovec. Inductive representation learning on large graphs. In Proceedings of the 31st International Conference on Neural Information Processing Systems , pages 1025-1035, 2017.

Max Horn, Edward De Brouwer, Michael Moor, Yves Moreau, Bastian Rieck, and Karsten M. Borgwardt. Topological graph neural networks. In The Tenth International Conference on Learning Representations, ICLR 2022, Virtual Event, April 25-29, 2022 . OpenReview.net, 2022.

Weihua Hu, Matthias Fey, Hongyu Ren, Maho Nakata, Yuxiao Dong, and Jure Leskovec. Ogblsc: A large-scale challenge for machine learning on graphs. arXiv preprint arXiv:2103.09430 , 2021.

Weihua Hu, Matthias Fey, Marinka Zitnik, Yuxiao Dong, Hongyu Ren, Bowen Liu, Michele Catasta, and Jure Leskovec. Open graph benchmark: Datasets for machine learning on graphs. arXiv preprint arXiv:2005.00687 , 2020.

Qian Huang, Horace He, Abhay Singh, Ser-Nam Lim, and Austin R Benson. Combining label propagation and simple models out-performs graph neural networks. arXiv preprint arXiv:2010.13993 , 2020.

John Jumper, Richard Evans, Alexander Pritzel, Tim Green, Michael Figurnov, Olaf Ronneberger, Kathryn Tunyasuvunakool, Russ Bates, Augustin Žídek, Anna Potapenko, et al. Highly accurate protein structure prediction with alphafold. Nature , 596(7873):583-589, 2021.

Thomas N. Kipf and Max Welling. Semi-supervised classification with graph convolutional networks. In International Conference on Learning Representations (ICLR) , 2017.

Ioannis Koutis and Gary L. Miller. CMG: Combinatorial Multigrid. https://github.com/ ikoutis/cmg-solver .

Ioannis Koutis, Gary L. Miller, and Richard Peng. Approaching optimality for solving SDD linear systems. SIAM J. Comput. , 43(1):337-354, 2014.

Ioannis Koutis, Gary L. Miller, and David Tolliver. Combinatorial preconditioners and multilevel solvers for problems in computer vision and image processing. Computer Vision and Image Understanding , 115(12):1638-1646, 2011. Special issue on Optimization for Vision, Graphics and Medical Imaging: Theory and Applications.

Ron Levie, Federico Monti, Xavier Bresson, and Michael M Bronstein. Cayleynets: Graph convolutional neural networks with complex rational spectral filters. IEEE Transactions on Signal Processing , 67(1):97-109, 2018.

Pan Li, Yanbang Wang, Hongwei Wang, and Jure Leskovec. Distance encoding: Design provably more powerful neural networks for graph representation learning. In Hugo Larochelle, Marc'Aurelio Ranzato, Raia Hadsell, Maria-Florina Balcan, and Hsuan-Tien Lin, editors, Advances in Neural Information Processing Systems 33: Annual Conference on Neural Information Processing Systems 2020, NeurIPS 2020, December 6-12, 2020, virtual , 2020.

László Lovász. Random walks on graphs: A survey. Combinatorics , 2:1-46, 1993.

Sitao Luan, Mingde Zhao, Xiao-Wen Chang, and Doina Precup. Break the ceiling: Stronger multi-scale deep graph convolutional networks. Advances in Neural Information Processing Systems , 32:10945-10955, 2019.

Christopher Morris, Martin Ritzert, Matthias Fey, William L Hamilton, Jan Eric Lenssen, Gaurav Rattan, and Martin Grohe. Weisfeiler and leman go neural: Higher-order graph neural networks. In Proceedings of the AAAI Conference on Artificial Intelligence , volume 33, pages 4602-4609, 2019.

Ryan Murphy, Balasubramaniam Srinivasan, Vinayak Rao, and Bruno Ribeiro. Relational pooling for graph representations. In International Conference on Machine Learning , pages 4663-4673. PMLR, 2019.

Vinod Nair, Sergey Bartunov, Felix Gimeno, Ingrid von Glehn, Pawel Lichocki, Ivan Lobov, Brendan O'Donoghue, Nicolas Sonnerat, Christian Tjandraatmadja, Pengming Wang, et al. Solving mixed integer programs using neural networks. arXiv preprint arXiv:2012.13349 , 2020.

Giannis Nikolentzos and Michalis Vazirgiannis. Geometric random walk graph neural networks via implicit layers. 2022.

Sarah Parisot, Sofia Ira Ktena, Enzo Ferrante, Matthew Lee, Ricardo Guerrero, Ben Glocker, and Daniel Rueckert. Disease prediction using graph convolutional networks: application to autism spectrum disorder and alzheimer's disease. Medical image analysis , 48:117-130, 2018.

Jiezhong Qiu, Yuxiao Dong, Hao Ma, Jian Li, Kuansan Wang, and Jie Tang. Network embedding as matrix factorization: Unifying deepwalk, line, pte, and node2vec. In Proceedings of the eleventh ACM international conference on web search and data mining , pages 459-467, 2018.

Ryoma Sato, Makoto Yamada, and Hisashi Kashima. Random features strengthen graph neural networks. In Proceedings of the 2021 SIAM International Conference on Data Mining (SDM) , pages 333-341. SIAM, 2021.

Victor Garcia Satorras, Emiel Hoogeboom, and Max Welling. E(n) equivariant graph neural networks. In Marina Meila and Tong Zhang, editors, Proceedings of the 38th International Conference on Machine Learning, ICML 2021, 18-24 July 2021, Virtual Event , volume 139 of Proceedings of Machine Learning Research , pages 9323-9332. PMLR, 2021.

Michael Sejr Schlichtkrull, Thomas N. Kipf, Peter Bloem, Rianne van den Berg, Ivan Titov, and Max Welling. Modeling relational data with graph convolutional networks. In Aldo Gangemi, Roberto Navigli, Maria-Esther Vidal, Pascal Hitzler, Raphaël Troncy, Laura Hollink, Anna

Tordai, and Mehwish Alam, editors, The Semantic Web - 15th International Conference, ESWC 2018, Heraklion, Crete, Greece, June 3-7, 2018, Proceedings , volume 10843 of Lecture Notes in Computer Science , pages 593-607. Springer, 2018.

Daniel A. Spielman and Nikhil Srivastava. Graph sparsification by effective resistances. SIAM J. Comput. , 40(6):1913-1926, 2011.

Jonathan M Stokes, Kevin Yang, Kyle Swanson, Wengong Jin, Andres Cubillos-Ruiz, Nina M Donghia, Craig R MacNair, Shawn French, Lindsey A Carfrae, Zohar Bloom-Ackermann, et al. A deep learning approach to antibiotic discovery. Cell , 180(4):688-702, 2020.

Prasad Tetali. Random walks and effective resistance of networks. J. Theoretical Probability , 1:101-109, 1991.

Petar Veliˇ ckovi´ c, Guillem Cucurull, Arantxa Casanova, Adriana Romero, Pietro Lio, and Yoshua Bengio. Graph attention networks. In International Conference on Learning Representations (ICLR) , 2018.

Juexin Wang, Anjun Ma, Yuzhou Chang, Jianting Gong, Yuexu Jiang, Ren Qi, Cankun Wang, Hongjun Fu, Qin Ma, and Dong Xu. scgnn is a novel graph neural network framework for single-cell rna-seq analyses. Nature communications , 12(1):1-11, 2021.

Yangkun Wang, Jiarui Jin, Weinan Zhang, Yong Yu, Zheng Zhang, and David Wipf. Bag of tricks for node classification with graph neural networks. arXiv preprint arXiv:2103.13355 , 2021.

Keyulu Xu, Weihua Hu, Jure Leskovec, and Stefanie Jegelka. How powerful are graph neural networks? In International Conference on Learning Representations , 2018.

Chengxuan Ying, Tianle Cai, Shengjie Luo, Shuxin Zheng, Guolin Ke, Di He, Yanming Shen, and Tie-Yan Liu. Do transformers really perform bad for graph representation? ArXiv , abs/2106.05234, 2021.

Chengxuan Ying, Tianle Cai, Shengjie Luo, Shuxin Zheng, Guolin Ke, Di He, Yanming Shen, and Tie-Yan Liu. Do transformers really perform badly for graph representation? In Marc'Aurelio Ranzato, Alina Beygelzimer, Yann N. Dauphin, Percy Liang, and Jennifer Wortman Vaughan, editors, Advances in Neural Information Processing Systems 34: Annual Conference on Neural Information Processing Systems 2021, NeurIPS 2021, December 6-14, 2021, virtual , pages 28877-28888, 2021.

Rex Ying, Dylan Bourgeois, Jiaxuan You, Marinka Zitnik, and Jure Leskovec. GNNexplainer: Generating explanations for graph neural networks. Advances in neural information processing systems , 32:9240, 2019.

Rex Ying, Ruining He, Kaifeng Chen, Pong Eksombatchai, William L Hamilton, and Jure Leskovec. Graph convolutional neural networks for web-scale recommender systems. In Proceedings of the 24th ACM SIGKDD International Conference on Knowledge Discovery & Data Mining , pages 974-983, 2018.

Jiaxuan You, Rex Ying, and Jure Leskovec. Position-aware graph neural networks. In International Conference on Machine Learning , pages 7134-7143. PMLR, 2019.

Hanqing Zeng, Hongkuan Zhou, Ajitesh Srivastava, Rajgopal Kannan, and Viktor K. Prasanna. Graphsaint: Graph sampling based inductive learning method. In 8th International Conference on Learning Representations, ICLR 2020, Addis Ababa, Ethiopia, April 26-30, 2020 . OpenReview.net, 2020.

Jiani Zhang, Xingjian Shi, Junyuan Xie, Hao Ma, Irwin King, and Dit-Yan Yeung. Gaan: Gated attention networks for learning on large and spatiotemporal graphs. arXiv preprint arXiv:1803.07294 , 2018.

Muhan Zhang and Yixin Chen. Link prediction based on graph neural networks. Advances in Neural Information Processing Systems , 31:5165-5175, 2018.

A Hyperparameters for PNA dataset.

In this section we provide the hyperparameters used for the different models on the PNA multitask benchmark. We train all models for 2000 steps and with 3 layers. The remaining hyperparameters for hidden size of each layer, learning rate, number of message passing steps (only valid for MPNN models), number of rotation matrices and same example frequency (when relevant) are provided in Table 6.

B Details of MPNN Framework

As discussed previously, the architectures used in the experiments conform to the MPNN framework, which allows affinity measures to be added as additional node and edge features. We describe the details here for completeness.

Assume that our input graph, G = ( V , E ) , has node features x u ∈ R n , edge features x uv ∈ R m and graph-level features x G ∈ R l , for nodes u, v ∈ V and edges ( u, v ) ∈ E . We provide encoders f n : R n → R k , f e : R m → R k and f g : R l → R k that transform these inputs into a latent space:

Our MPNN then performs several message passing steps:

where H ( t ) = ( { h ( t ) u } u ∈V , { h ( t ) uv } ( u,v ) ∈E , h ( t ) G ) contains all of the latents at a particular processing step t ≥ 0 .

This process is iterated for T steps, recovering final latents H ( T ) . These can then be decoded into node-, edge-, and graph-level predictions (as required), using analogous decoder functions g n , g e and g g :

Generally, f and g are simple MLPs, whereas we use the MPNN update rule for P . It computes message vectors, m ( t ) uv , to be sent across the edge ( u, v ) , and then aggregates them in the receiver nodes as follows:

The message function ψ t +1 and the update function ϕ t +1 are both MLPs. All of our models have been implemented using the jraph library [16].

We incorporate edge-based affinity features (e.g., effective resistances and hitting times) in f e and node-based affinity features (e.g., resistive embeddings) in f n . Note that node-based affinity

features may also naturally be incorporated as edge features by concatenating the node features at the endpoints.

Occasionally, the dataset in question will be easy to overfit with the most general form of message function (see (6)). In these cases, we resort to assuming that ψ factorises into an attention mechanism :

where the attention function a is scalar-valued. We will refer to this particular MPNN baseline as a graph attention network ( GAT ) [44].

C Omitted Proofs

Lemma 3.2. For any pair of nodes u, v , we have ∥ r u -r v ∥ 2 2 = Res ( u, v ) .

Proof.

Corollary 4.2. For any fixed vectors α, β ∈ R n , if we let X := ∑ i α i x i , ̂ X := ∑ i α i ̂ x i and similarly Y := ∑ i β i x i , ̂ Y := ∑ i β i ̂ x i ; then:

Proof. Since ⟨ X,Y ⟩ = 1 4 ( ∥ X + Y ∥ 2 -∥ X -Y ∥ 2 ) , we can bound A = ∣ ∣ ⟨ X,Y ⟩ - ⟨ ̂ X, ̂ Y ⟩ ∣ ∣ from above as:

where (8) follows from Lemma 4.1 with probability 1 -o (1) , by our choice of k (as Lemma 4.1 guarantees that each of ∥ ̂ X + ̂ Y ∥ 2 = (1 ± ϵ ) ∥ X + Y ∥ 2 and ∥ ̂ X -̂ Y ∥ 2 = (1 ± ϵ ) ∥ X -Y ∥ 2 holds with probability 1 -o (1) , and one can take a union bound over the two events).

Proof. Consider the following expression of hitting times in terms of commute times by [43].

Dividing both sides of eq. (9) and using the relation K u,v = 2 M Res ( u, v ) , we see that:

Let's focus on the inner summation. After expanding out the squared norms, we see that:

Substituting this back into eq. (10), we can express 1 2 M H u,v as:

Lemma 4.4. | ̂ H u,v -H u,v | ≤ 3 ϵH max .

Proof. Using Lemma 4.3, we see that

where we used Corollary 4.2 in the first inequality and Definition 3.4 in the last inequality.

D Comparison: Effective Resistances vs. Shortest Path Distances

Given that effective resistance (ER) captures times associated with random walks in a graph, it is tempting to ask how effective resistances compare to shortest path distances (SPDs) between nodes in a graph. Indeed, for some simple graphs, e.g., trees, shortest path distances and effective resistances turn out to be identical. However, in general, effective resistances and shortest path distances behave quite differently.

Nevertheless, it is tempting to ask how effective resistance features compare to SPD features in GNNs, especially as there have been a number of recent model architectures that make use of SPD features (e.g., Graphormer [49], Position-Aware GNNs [52], DE-GNN [29]). We first note that the most natural direct comparison of our ER-based MPNNs with SPD-based networks does not quite make sense. The reason is that the analogous comparison would be to determine the effect of replace ERs with SPDs as features in our MPNNs. However, since our networks only use ER features along edges of the given graph, the corresponding SPD features would then be trivial (as the SPD between two nodes directly connected by an edge in the graph is 1, resulting in a constant feature on every edge)!

As a result, graph learning architectures that use SPDs typically either (a.) use a densely-connected network (e.g., Graphormer [49], which uses a densely-connected attention mechanism) that incurs O ( n 2 ) overhead, or (b.) pick a small set of anchor nodes or landmark nodes to which SPDs from all other nodes are computed and incorporated as node features (e.g., Position-Aware GNNs [52], DE-GNN [29]). We stress that the former approach generally modifies the graph (by connecting all pairs of nodes) and therefore does not fall within the standard MPNN approach, while the latter includes architectures that fall within the MPNN paradigm.

Furthermore, we note that DE-GNNs are arguably one of the closest proposals to ours, as they compute distance-encoded features. These features can be at least as powerful as our proposed affinity-based features if polynomially many powers of the adjacency matrix are used. However, for all but the smallest graphs, using this many powers will be impractical-in fact, [29] only use powers of A up to 3, which would not be able to reliably approximate affinity-based features. We also

observe that the DE-GNN paper is concerned with learning representations of small sets of nodes (e.g., node-, link-, and triangle-prediction) and does not show how to handle graph prediction tasks, which the authors mention as possible future work. This makes a direct comparison of our methods with DE-GNNs difficult.

D.1 Empirical Results

In an effort to empirically compare the expressivity of ER features with that of SPD features, we once again perform experiments on the PNA dataset, picking the following baselines that make use of SPD features:

The first baseline is roughly an MPNN with Graphormer-based features . More precisely, it is a densely-connected MPNN with SPDs from the original graph as edge features. In order to retain the structure of the original graph, we also use additional edge features to indicate whether or not an edge in the dense (complete) graph is a true edge of the original graph. We also explore the use of the centrality encoding (in-degree and out-degree embeddings) from Graphormer as additional node features.

The second baseline is the Position-Aware GNN (P-GNN), which makes use of 'anchor sets' of nodes and encodes distances to these nodes.

The results of these baselines are shown in Table 7. In particular, we note that our ER-based MPNNs outperform all aforementioned baselines.

Table 7: Results on the PNA dataset for MPNNs with Graphormer-based features (yellow) as well as SPD-based P-GNNs (orange). Here, CE refers to the centrality encoding , which is incorporated in the relevant MPNNs as additional node features. Similarly, SPD refers to shortest path distance features - in the relevant MPNNs, shortest path distances between all pairs of nodes in the graph are incorporated as edge features, along with an additional edge feature indicating whether an edge exists in the input graph. Therefore, the MPNN baselines are all variants of the same model with additional node/edge features. Similarly, P-GNN [52] uses SPD features with respect to a set of chosen anchor nodes . The average score metric is, as before, the average of the log( MSE ) metric over all six tasks, as in Table 1.

D.2 Theory: ER vs. SPD

In addition to experimental results, we would like to provide some theory for why effective resistances can capture structure in GNNs that SPDs are unable to.

̸

We will call an initialization function u ↦→ h (0) u on nodes of a graph node-based if it assigns values that are independent of the edges of the graph. Such an initialization is, however, allowed to depend on node identities (e.g., for the single-source shortest path problem from a source s , one might find it natural to define h (0) s = 0 and h (0) u = + ∞ for all u = s ).

Consider the task of computing 'single-source effective resistances,' i.e., the effective resistance from a particular node to every other node. We show that a GNN with a limited number of message passing steps cannot possibly learn single-source effective resistances, even to nearby nodes.

Theorem D.1. Suppose we fix k > 0 . Then, given any node-based initialization function h (0) u , it is impossible for a GNN to compute single-source effective resistances from a given node w to any nodes within a k -hop neighborhood.

More specifically, for any update rule

̸

there exists a graph G = ( V, E ) and u ∈ V such that after k rounds of message passing, h ( k ) v = Res ( u, v ) for some v = u within a k -hop neighborhood of u .

̸

On the other hand, there exists an initialization with respect to which k rounds of message passing will compute the correct shortest path distances to all nodes within k -hop neighborhood.

Note that the assumption on the initialization function in the above theorem is reasonable because enabling the use of arbitrary, unrestricted functions would allow for the possibility of precomputing effective resistances in the graph and trivially incorporating them as node features, which would defeat the purpose of computing them using message-passing.

Proof. Consider the following set of graphs, each on 4 k +1 nodes:

Figure 2: Both of the above graphs are on 4 k +1 vertices, labeled v 0 , v 1 , . . . , v 4 k . The only difference is a single edge, i.e., the graph on the left has an edge between v 2 k and v 2 k +1 , while the one on the right does not have this edge.

v

0

v

1

v

2

v

2

k

v

2

k

+1

v

4

k

v

0

v

1

v

2

v

2

k

v

2

k

+1

v

4

k

Let w = v 0 be the source and let { h (0) v : v ∈ V } be a 'local' node feature initialization. Note that for any GNN (i.e., update and aggregation rules in (11), add the formal update rule somewhere), the computation tree after k rounds of message passing is identical for nodes v 0 , v 1 , . . . , v k , v 3 k +1 , v 3 k +2 , . . . , v 4 k (i.e., the nodes within the k -hop neighborhood of v 0 ) in both G and G ′ . This is because the only difference between G and G ′ is the existence of the edge between v 2 k and v 2 k +1 , and this edge is beyond a k -hop neighborhood centered at any one of the aforementioned nodes. Therefore, we will necessarily have that h ( k ) v i is identical in both G and G ′ for i = 1 , . . . , k, 3 k +1 , 3 k +2 , . . . , 4 k .

̸

However, it is easy to calculate the effective resistances in both graphs. In G , we have Res G ( v 0 , v i ) = i (4 k +1 -i ) 4 k +1 , while in G ′ , we have Res G ′ ( v 0 , v i ) = min { i, 4 k +1 -i } . Therefore, Res G ( v 0 , v i ) = Res G ′ ( v 0 , v i ) for all i = 1 , 2 , . . . , k, 3 k +1 , 3 k +2 , . . . , 4 k .

It follows that for any i = 1 , 2 , . . . , k, 3 k +1 , 3 k +2 , . . . , 4 k , the execution of k message passing steps of a GNN cannot result in h ( k ) v i = Res ( v 0 , v i ) for both G and G ′ , which proves the first claim of the theorem.

̸

For the second part (regarding single-source shortest paths), observe that single-source shortest path distances can, indeed, be realized via aggregation and update rules for a message passing network. In particular, for k rounds of message passing, it is possible to learn shortest path distances of all nodes within a k -hop neighborhood. Specifically, for a source w , we can use the following setup: Take h w = 0 and h u = ∞ for all u = w . Moreover, for any edge ( u, v ) , let the edge feature x uv ∈ R simply be the weight of ( u, v ) in the graph. Then, take the update rule (11) with f e , ψ t +1 as identity functions and

It is clear that the above update rule simply simulates the execution of an iteration of the Bellman-Ford algorithm. Therefore, k message passing steps will simulate k iterations of Bellman-Ford, resulting in correct shortest path distances from the source w for every node within a k -hop neighborhood.

Peer reviews:
{"IS_META_REVIEW": false, "comments": "SUMMARY:\nThis work proposes to incorporate affinity measures as features into message-passing networks (MPNNs) in order to enhance the expressivity without enlarging the computational cost notably.\nThe authors introduce three examples of random-walk-based affinity measures, e.g., effective resistance, hitting, and commute times, and provide efficient computation and low-rank approximation for them.\nOn multiple graph benchmarks, ranging from graph-level tasks to transductive node-classification on large-scale graphs, the proposed affinity measures can remarkably improve the performance of MPNNs and reach state-of-the-art performance among models without privileged information.\n\nSTRENGTHS:\n1. According to my best knowledge, this work is the first to propose that incorporating affinity-based measures into MPNNs can improve the expressive power of MPNNs, going beyond the 1-WL algorithm.\n2. This work also provides several techniques to improve the efficiency of the introduced affinity-based measures.\n3. The authors show that MPNNs+affinity measures are strictly more powerful than the 1-WL algorithm via the example of regular graphs.\n4. The empirical experiments demonstrate the effectiveness of the proposed methods on graph-level tasks and transductive node-level tasks. The experiment on large-scale graphs also shows the efficiency of the proposed technique.\n\nWEAKNESSES:\n1. The paper does not have more throughout comparisons with other positional encoding enhanced graph neural networks, theoretically and empirically. For details, please see the question section.\n2. The baselines in empirical comparisons seem to have different numbers of parameters compared to the proposed method. It will be better to have a more detailed explanation on the choices of hyperparameters.\n\nQUESTIONS:\n(Loukas, 2020; Dwivedi et al., 2021) mention that encoding can enhance the ability of MPNNs to uniquely distinguish nodes and hence improve the expressive power of MPNNs.\nI am curious whether the affinity measures enhance MPNNs in a similar way. Or they can provide extra valuable information that non-affinity-based encoding can not provide. \nIn other words, if there exists a non-affinity-based encoding able to distinguish isomorphic graphs as well as the proposed affinity measures, such as subgraph counting (Bouritsas et al., 2022), will incorporating the proposed affinity measures into MPNNs enjoy extra advantages over such an encoding? Whether graph-level tasks and node-level tasks have different conclusions regarding this?\n\nIf possible, it would be perfect if you could provide more comparisons theoretically/empirically against several PE-enhanced MPNNs (Dwivedi et al., 2021; Zhang et al., 2021; Wang et al., 2022; Lim et al., 2022; Li et al., 2020; Bouritsas et al., 2022) and potentially the PE for Graph Transformers (Zhang et al., 2023; Ma et al., 2023).\n\n\n\n\n- Loukas, A. (2020). How hard is to distinguish graphs with graph neural networks? Adv. Neural Inf. Process. Syst.\n- Dwivedi, V. P., Luu, A. T., Laurent, T., Bengio, Y., & Bresson, X. (2021). Graph Neural Networks with Learnable Structural and Positional Representations. Proc. Int. Conf. Learn. Representations. \n- Bouritsas, G., Frasca, F., Zafeiriou, S. P., & Bronstein, M. (2022). Improving Graph Neural Network Expressivity via Subgraph Isomorphism Counting. IEEE Transactions on Pattern Analysis and Machine Intelligence, 1–1.\n- Zhang, Z., Cui, P., Pei, J., Wang, X., & Zhu, W. (2021). Eigen-GNN: A Graph Structure Preserving Plug-in for GNNs. IEEE Transactions on Knowledge and Data Engineering.\n- Wang, H., Yin, H., Zhang, M., & Li, P. (2022). Equivariant and Stable Positional Encoding for More Powerful Graph Neural Networks. Proc. Int. Conf. Learn. Representations.\n- Lim, D., Robinson, J. D., Zhao, L., Smidt, T., Sra, S., Maron, H., & Jegelka, S. (2022). Sign and Basis Invariant Networks for Spectral Graph Representation Learning. ICLR 2022 Workshop on Geometrical and Topological Representation Learning. \n- Li, P., Wang, Y., Wang, H., & Leskovec, J. (2020). Distance Encoding: Design Provably More Powerful Neural Networks for Graph Representation Learning. Adv. Neural Inf. Process. Syst. \n- Zhang, B., Luo, S., Wang, L., & He, D. (2023). Rethinking the Expressive Power of GNNs via Graph Biconnectivity. Proc. Int. Conf. Learn. Representations. \n- Ma, L., Lin, C., Lim, D., Romero-Soriano, A., K. Dokania, Coates, M., H.S. Torr, P., & Lim, S.-N. (2023). Graph Inductive Biases in Transformers without Message Passing. Proc. Int. Conf. Mach. Learn.\n\nLIMITATIONS:\nNo limitations are provided.", "TITLE": "Official Review", "RECOMMENDATION": 7, "REVIEWER_CONFIDENCE": 3, "SOUNDNESS_CORRECTNESS": 3, "DATE": "2023-06-26"}

{"IS_META_REVIEW": false, "comments": "SUMMARY:\nThe paper proposes a strategy to strengthen the node and edge features in order to enhance the expressiveness of a message passing neural network. In particular, the paper introduces a set of effective resistance (ER) features, including node-level resistive embedding, which further derives two edge-level affinity measures: hitting time and commuting time/effective resistance. Additionally, the paper employs Constructive Johnson-Lindenstrauss to randomly project resistive embeddings to a lower-dimensional space, which is utilized to expedite the computation of hitting time and commuting time. Finally, the paper present the empirical evaluations on the proposed ER features.\n\nSTRENGTHS:\n1. The paper is well-written: the motivation is clearly presented, and most sections are straightforward to follow.\n2. The fast approximation of computing hitting time and commuting time is novel. The dimensionality of resistive embedding is reduced without the requirement of specifying a set of anchor nodes.\n3. The experiment on large graphs demonstrates that the computation cost of adding these features to MPNN is manageable.\n\nWEAKNESSES:\n1. The paper would benefit from providing stronger evidence to support certain claims made throughout the text. Please refer to the \"Questions\" section for details.\n\nQUESTIONS:\n1. I would like to seek clarification from the authors regarding the connection between \"1-WL\" and \"standard GNN\" in the context of node classification, as demonstrated in the paper to support the claim made by Theorem 3.6. The provided example gives me the impression that a \"standard GNN\" represents each node using a spanning tree rooted at the node of interest (referred to as a computation tree) and distinguishes nodes based on the isomorphism of these spanning trees.\n\n   To facilitate readers' understanding of this concept, I kindly request the authors to provide a specific example of the computation tree that a standard GNN employs to represent a node. It would also be helpful if the authors could provide some justification or reference regarding the connection between the node representation learned by standard GNNs and the utilization of spanning trees.\n\n2. By utilizing the initial node feature as the one-hot encoding of the node's index, denoted as $\\textbf{h}_v^{(0)} = \\textbf{e}_v$, and employing the forward function $\\mathbf{H}^{(t+1)} = \\mathbf{A}\\mathbf{H}^{(t)}$, we can observe that after 4 layers, the node embeddings generated by the standard GNN successfully differentiate the 3 color categories. Specifically, if we define a score $s_v = <\\textbf{h}_v^{(4)}, \\textbf{e}_v>$, the value for nodes {1,2} is 15, for nodes {3,4,7,8} it is 17, and for nodes {5,6} it is 19. This observation contradicts the claim that \"a standard GNN that is limited by the 1-WL test cannot distinguish any pair of nodes in a regular graph.\" Consequently, I have some follow-up questions based on this discrepancy.\n  \n   (1) What is the specific constraint on the input features to a standard GNN?\n   \n   (2)  Will the benefit from resistive embedding become less evident as we increase the depth of an MPNN?\n\n3. In Appendix D.2, a comparison is conducted between resistive embedding and shortest-path-distance embedding in terms of their ability to differentiate nodes within a cycle and the corresponding nodes in the path obtained by removing one edge from the cycle. \nThe example effectively illustrates that node embeddings based on shortest-path-distance are insufficient in capturing the divergence that can be successfully captured by resistive embeddings. However, this raises a valid concern regarding real-life networks where edges may be lost due to noise. In such cases, will the strong ability of resistive embeddings to distinguish between graphs that are close in terms of edit distance potentially make them susceptible to overfitting the noise?\n\nLIMITATIONS:\nN/A", "TITLE": "Official Review", "RECOMMENDATION": 7, "REVIEWER_CONFIDENCE": 4, "SOUNDNESS_CORRECTNESS": 3, "DATE": "2023-07-03"}

{"IS_META_REVIEW": false, "comments": "SUMMARY:\nThis paper proposes to use affinity measures as additional features that can incorporate in common standard MPNNs and theoretically and empirically shows the improvement in expressiveness and performance in some datasets without loss of much scalability. The main contributions are listed below.\n1. Generality. The proposed approach can be adopted in standard MPNNs.\n\n2. Expressiveness. Incorporating the proposed features can make standard MPNNs more powerful than the 1-WL test.\n\n3. Scalability. The affinity measures used in this paper are based on the random walk. By approximation, the features derived from the proposed measures can be incorporated with high scalability.\n\n4. Performance. MPNNs equipped with the proposed additional features can make progress or even rank 1 in some datasets.\n\nSTRENGTHS:\n1. This paper was inspired by effective resistances and proposed HT, ER and resistive embeddings all help common MPNNs in both performance and expressiveness. As a feature-based approach, it can adapt to many standard MPNNs and does not need to specify any anchor or target nodes to capture both distance and connectivity information using general-purpose node and edge features.\n\n2. By Approximation and fast computation of Resistive Embeddings and HT, the proposed features can scale up to large graphs in time and space complexity.\n\n3. The experimental results are impressive. With relatively low computation cost, the proposed method can have a competitive performance on small datasets and higher performance on homogeneous/heterogeneous citation networks. Further, the best performance achieved on PCQM4M-v1 is impressive.\n\n4. Both theoretical proof and empirical evaluations are provided.\n\nWEAKNESSES:\n1. The motivation needs to be clarified more clearly. What the main drawbacks of the existing approaches are and which the proposed one surmounts? It’s said that the proposed approach “without the need for specifying any anchor or target nodes” and the scalability problem is mentioned in the experiments section. However, the advantages should be clarified further and organized better since there are other feature-based works and other categories of works aiming at a similar goal. Some practical evaluation or comparison may help if possible.\n\n2. To demonstrate the superiority and necessity of the proposed affinity feature-based approach, other similar approaches need to be compared. Similar to the proposed affinity-based features, approaches incorporating other additional features (e.g. Random features or other dedicated designed features), and other methods that can improve GNNs’ expressiveness (e.g. other 2 directions mentioned in Expressivity para in Sec. 2) should also be compared in scalability and performance.\n\n3. Although by approximation, the computation of random embeddings can be done in near-linear time and can be scaled to a large dataset, it’s better to provide the actual runtime of the proposed and baseline methods.\n\n4. Theorem 3.6 claims that MPNNs that make use of any one of the proposed features are strictly more powerful than the WL-1 test. However, the experiments for it are not enough. Experiments on ogbg-molhiv show that it outperforms substructure counts, but the reasons are not clear. Experiments on some (artificial) datasets (e.g. [1][2]) which are designed for testing expressiveness would be better.\n\n5. More well-organized writing should be considered. The notations scattering in Sec. 3.1-3.3 bring some difficulty in reading. The different baselines, settings, and experimental results in each subsection in Sec. 5 are scattered, either in the tables or in the text, which may weaken the empirical demonstration effects of experiments and cause some confusion. Some nouns or notations are introduced but not explained clearly, such as ER scalar features, node or edge features used in experiments, e_v in Eq. (1), and random rotations. Introducing them in a centralized and unified manner may help.\n\nQUESTIONS:\n1. It’s said in Sec. 5.4 that ER embedding dimension become impractically high and authors alter to use ER scalar features. What are the theoretical and empirical (experimental) differences between the two kinds of features? If ER scalar features in small graphs can make the same progress with ER embedding? If not, why?\n\n2. Incorporating hitting time features, resistive embeddings, and effective resistance are all proven to be more powerful than 1-WL in Sec. 3.4. Why does using ER/HT in experiments have different results? Is there any suggestion about selecting which feature?\n\n3. The performance of the proposed method does not surpass the DGN in Table 2 and does not show advantages when more layers are in Table 3. More explanations are expected although the proposed method shows advantages in other aspects.\n\n4. How the proposed features help in heterogeneous networks may be interesting. Possible explanations are expected.\n\nLIMITATIONS:\n1. As said in W3, the actual runtime in large graphs may be much higher than (scalable) baseline methods.\n\n2. As said in Q2, experiments use different features (ER embedding/HT …), thus suggestion about selecting feature is needed.\n\n3. When demonstrating the expressiveness of the proposed approach, experiments on some datasets tailored for evaluating expressive power are expected. For instance, one can synthesize a dataset containing instances indistinguishable by 1-WL(e.g. [1][2]) to see if the proposed approach can distinguish.\n\n[1] Sato R, Yamada M, Kashima H. Random features strengthen graph neural networks[C]//Proceedings of the 2021 SIAM international conference on data mining (SDM). Society for Industrial and Applied Mathematics, 2021: 333-341.\n[2] Abboud R, Ceylan I I, Grohe M, et al. The surprising power of graph neural networks with random node initialization[J]. arXiv preprint arXiv:2010.01179, 2020.", "TITLE": "Official Review", "RECOMMENDATION": 4, "REVIEWER_CONFIDENCE": 3, "SOUNDNESS_CORRECTNESS": 3, "DATE": "2023-07-06"}

comments: SUMMARY:
In this paper, the authors present MPNN using affinity measures as node and edge features. As the affinity measures, the authors propose effective resistance, hitting times, and resistive embeddings. The authors demonstrate the effectiveness of using the affinity measures through experiments on various datasets. Experimental results on large graph datasets prove the good scalability of the scheme.

STRENGTHS:
1. I considered the need for a comparison with shortest path distances (SPDs), however, the related description in Section D (in the supplementary material) is convincing.

2. The proposed scheme has low computational complexity and good scalability.

3. The paper is well-written and easy to understand.

WEAKNESSES:
It appears that the setting for MPNN models used in the experiments is not fair. According to Table 6 in the supplementary material, the number of layers is set to 3 for all models. Do the authors consider this setting to be fair? Furthermore, the hidden size of each layer is set differently for each method. What is the basis for setting the hidden size for each method? For example, the hidden size for GAT is set to 64 while that for HT + ER (rand rot) is set to 512.

minor comment
- fom -> from (Line 112)

QUESTIONS:
Q1. Do the authors think that the experiments with the current hyper-parameter setting are fair?

Q2. Could the authors explain on what basis they set the hyper-parameters for each method?

Q3. Is it possible to improve performance by adding the proposed features to DGN?

LIMITATIONS:
The authors do not describe the potential negative social impact of their work while present future work to develop the work.
TITLE: Official Review
SOUNDNESS_CORRECTNESS: 3
DATE: 2023-07-07