REVIEW 4 major objections 6 minor 1 cited by
Demystifying MPNNs: Message Passing as Merely Efficient Matrix Multiplication
T0 review · 4 major / 6 minor · reviewed 2026-08-09 · deepseek-v4-flash
Pith's one-line read Message passing is just efficient matrix multiplication
desk verdict The central A^k equivalence is unsupported (invalid UAT step, concrete counterexample), but the empirical depth-vs-neighborhood decomposition and loop analysis deserve referee time. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The load-bearing object is the $k$-th power of the adjacency matrix, $A^k$, whose nonzero entries count directed paths of length exactly $k$. The identity $H^{(k)} = \sigma((W \odot A)^k X W^{(k)})$ carries the argument: it says depth in an MPNN computes powers of the normalized adjacency matrix. The proof's enabling step is the binomial expansion $(A+I)^k = \sum_{i=0}^k \binom{k}{i} A^i$, which converts the self-loop formula into a sum of $k$-hop features. The path-counting interpretation of matrix powers, proved by induction in Lemma 2.5, connects algebraic powers to graph-theoretic neighborhoods.
What would settle it
Train a $k$-layer GCN without self-loops and, in parallel, a single-layer GCN whose adjacency matrix is $A^k$, using the same weight initialization and training budget on a small graph; if the two models' outputs or test accuracies differ substantially (beyond what function-approximation error would allow), the claimed equivalence is falsified, because the closed form predicts they should be interchangeable.
Extended reading notes
Core claim
The paper's central discovery is the closed-form equivalence expressed in Lemma 2.7: for a $k$-layer GCN without self-loops, $H^{(k)} = \sigma((W \odot A)^k X W^{(k)})$, where $W$ is the normalization weight matrix. The proof proceeds by induction and uses the Universal Approximation Theorem to justify dropping the activation inside each layer. For a GCN with self-loops, the analogous formula gives $\sigma((W \odot (A+I))^k X W^{(k)})$, which by the binomial theorem is a linear combination of all powers $A^0$ through $A^k$; GraphSAGE produces a similar combination through its separate self-feature path. The paper also proves loop-driven inheritance rules (self-loops add one hop, undirected edges add two, $m$-cycles add $m$) and uses them to explain over-smoothing and to identify gradient-related degradation in sparse graphs. For uniform features, it shows a $(k+1)$-layer network equals a $k$-layer network using node degrees as features.
Load-bearing premise
The proof that a $k$-layer network equals a single layer on $A^k$ assumes that removing the nonlinear activation inside each layer, justified by the Universal Approximation Theorem, does not change the network's behavior; the theorem only guarantees that some network can approximate a function, not that deleting this specific activation is harmless.
Editorial extensions
If this is right
- A $k$-layer GCN without self-loops uses exactly the $k$-hop neighbors of each node, not a cumulative aggregation of all shorter hops.
- Because $(A+I)^k$ expands into all powers of $A$, adding self-loops (or working on undirected graphs) makes $k$-hop neighborhoods contain shorter-hop neighbors, which the paper identifies as a mechanism for over-smoothing.
- Deep GCNs can maintain accuracy up to 50 layers on sparse directed graphs when self-loops and undirected symmetrization are removed, indicating gradient problems rather than over-smoothing are the main failure mode there.
- With uniform node features, node degree becomes the effective feature: a $(k+1)$-layer MPNN with constant features is equivalent to a $k$-layer MPNN using degree as the only feature.
- Message passing is a memory-efficient alternative to explicitly forming $A^k$ and multiplying it with $X$, since iterative multiplication avoids storing the dense power matrix.
Reading between the lines
- If the central equivalence is accepted, then tuning depth is effectively tuning the power of the adjacency matrix; this suggests that layer count and neighborhood radius are interchangeable up to memory constraints, and that residual or normalization tricks that preserve the linear form may be more effective than architectural modifications aimed at over-smoothing.
- The paper's distinction between density-driven over-smoothing (loops, undirected edges) and gradient-driven degradation (sparse directed graphs) could be turned into a diagnostic: compare the $k$-layer model's accuracy with the single-layer $A^k$ model's accuracy; a large gap with comparable neighborhoods points to gradient or optimization issues rather than over-smoothing.
- A testable extension: on graphs where adding self-loops hurts, the theory predicts the harm scales with the density of $(A+I)^k$ versus $A^k$; this can be checked by ablating self-loops on a range of homophily and heterophily datasets.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The paper aims to provide a theoretical account of message-passing neural networks (MPNNs), claiming three main results: (1) a k-layer MPNN operating on adjacency matrix A is approximately equivalent to a single-layer MPNN operating on A^k, so that message passing is essentially a memory-efficient way to compute powers of the adjacency matrix; (2) different loop structures (self-loops, two-node loops, longer cycles) cause hop coexistence and thereby contribute to over-smoothing; and (3) node classification tasks can be split into structure-feature hybrid and structure-only types, with the latter reducible to degree-based feature prediction. The theoretical core is Lemma 2.7 (and Lemmas 2.8, 2.9), whose proofs in Appendix A rely on removing the inner nonlinearity via the Universal Approximation Theorem. The empirical sections compare deep GCNs, single-layer GCNs on k-hop neighborhoods, and hybrid architectures on citation and web graphs, and they study normalization schemes with uniform features.
Significance. If the central equivalence were valid, the paper would provide a clean and useful interpretation of deep GNNs as iterative matrix-power computation, and it would connect loop structure to over-smoothing through hop coexistence. The empirical comparison of depth growth versus k-hop neighborhood growth on directed sparse graphs is a useful experiment, and the observation that row normalization with uniform features destroys degree information is sound. The paper also provides a code link and reports standard deviations over multiple splits. However, the main theoretical claim is not established: the proof of the key lemmas uses an invalid application of the Universal Approximation Theorem, and a concrete ReLU counterexample shows that the claimed A^k equivalence is false as stated. The paper's central theoretical contribution is therefore unsupported, even though some of its empirical observations may be of independent interest.
major comments (4)
- [Appendix A.2, Eqs. (5)-(6)] The induction step is not a proof. After substituting the induction hypothesis, the expression σ(W⊙(A+I) σ((W⊙(A+I))^k X W^(k)) W^(k+1)) is replaced by σ((W⊙(A+I))^(k+1) X W^(k+1)) with the sole justification 'Thanks to the Universal Approximation Theorem'. UAT asserts that some sufficiently wide network can approximate a continuous function; it does not license deleting a fixed nonlinearity from a fixed compositional expression, nor does it imply that the composition is exactly representable as σ(M X W') for a single matrix M. The same unjustified step is used in Appendix A.3 for Lemma 2.7 and in Appendix A.4, Eqs. (9)-(10), for Lemma 2.9. Since this step is the load-bearing part of the paper's central claim, the proof is invalid.
- [Section 2.2, Lemma 2.7] The claimed identity is false for standard ReLU GCNs. Consider the directed graph with edges 1→4, 2→4, 3→5, 4→6, and 5→6, using the paper's convention A_{ij}=1 for an edge i→j, node feature vector X=[1,1,-1,0,0,0]^T, identity weight matrices, and ReLU activation. The 2-layer GCN without self-loops gives output 2 at node 6, while the single-layer GCN on A^2 gives ReLU(1)=1 at node 6. As functions, the two-layer output is ReLU(ReLU(x1+x2)+ReLU(x3)) and the single-layer output is ReLU(x1+x2+x3), which differ on an open set. Thus the A^k equivalence, and with it the 'efficient matrix multiplication' interpretation, is unsupported for nonlinear MPNNs.
- [Introduction and Section 2.3] The paper's advertised result is that a k-layer MPNN on adjacency matrix A is approximately equivalent to a single-layer MPNN on A^k, but Lemma 2.7 actually states H^(k)=σ((W⊙A)^k X W^(k)), where W⊙A is the normalization-mask form of the adjacency matrix introduced in Appendix E.1, Eq. (12). For a symmetrically normalized GCN, (D^{-1/2}AD^{-1/2})^k is not A^k, and the distinction affects both the neighborhood interpretation and the experiments in Section 4, which use powers of the original adjacency matrix. The formal statement of the paper's own lemma therefore does not match the headline claim.
- [Appendix A.4, Lemma 2.9] The proof of Lemma 2.9 repeats the same invalid UAT step: Eq. (9) contains inner σ terms applied to sums of powers of (W⊙A), and Eq. (10) removes those σ terms with no justification other than the same UAT sentence. In addition, the conclusion line of the proof states H^(n)=σ((W⊙(A+I))^n X W^(n)), which is the GCN-with-self-loop formula rather than the GraphSAGE expression claimed in Lemma 2.9. Since the derivation is central to the paper's treatment of GraphSAGE, this lemma is also unsupported.
minor comments (6)
- [Lemma 2.3] The lemma uses 'kth order node feature' but the statement is about p-hop neighbors; the variables p and k should be made consistent.
- [Equation (6)] The expression cW^(k+1) contains an unexplained symbol c; this appears to be a typo.
- [Notation throughout] The symbol W is overloaded: W⊙A denotes a fixed normalization mask, while W^(l) denotes learnable feature-transformation matrices. These should be distinguished by different symbols, especially in Lemmas 2.7-2.9.
- [Section 4.2.1] The statement that row normalization with homogeneous features makes 'all nodes get identical representations' assumes zero biases and identical initial features; the conditions should be stated explicitly.
- [Appendix C.2.1] The citation 'Rossi et al. (?)' is incomplete and should be filled in.
- [Table 1 and Figures 4-5] The density of the k-hop adjacency matrix is used to support the over-smoothing discussion, but the definition is not precise about whether density counts nonzero entries or path counts, nor about how normalization affects the reported values.
Circularity Check
No circular derivation: the A^k equivalence is not fit-derived or self-referential; the UAT-based deletion of inner nonlinearities is an unsupported proof step, and the only self-citations are non-load-bearing pointers.
full rationale
The derivation chain contains no step where a claimed prediction is fitted from data or defined in terms of its own target. Lemma 2.7 and 2.8 (Appendix A.2/A.3) are inductive proofs built on the layer update H^(k+1)=σ(W⊙(A+I)H^(k)W^(k)); the closed forms are not definitions of the output, and no fitted constant is used to force them. The one non-algebraic move is Eq. 5 to Eq. 6 in Appendix A.2 (and its parallel in A.3), where UAT is invoked to justify removing the inner σ from the composition. That step is mathematically unsupported—UAT does not license deleting a fixed activation from a fixed layerwise network—and a two-layer ReLU example contradicts the exact A^k equality, so the proof is not reliable as written. Unsupportedness is a correctness flaw, however, not a circular reduction: the conclusion is not an input to the proof. The only self-citations (Jiang et al., 2024) point to hyperparameter settings in the footnote of Section 4.2.1 and to background on higher-order models in E.1/E.2; they do not carry the central A^k argument. The score of 2 reflects that minor non-load-bearing self-citation; no circular step is present.
Assumptions & free parameters
assumptions (3)
- ad hoc to paper Universal Approximation Theorem justifies removing nonlinear activations inside each layer of the unrolled network without changing the function (Eq. 5-6, Appendix A.2).
- domain assumption Graphs are directed with adjacency entries in non-negative integers; undirected graphs are treated as bidirected directed graphs.
- domain assumption The conclusion that gradient issues dominate in sparse graphs generalizes from the tested citation and web datasets.
Cite this review
Pith. "Pith review of Demystifying MPNNs: Message Passing as Merely Efficient Matrix Multiplication." pith.science (2026). https://pith.science/paper/J6JZGOYB
@misc{pith2026250200140,
author = {Pith},
title = {Pith review of: Demystifying MPNNs: Message Passing as Merely Efficient Matrix Multiplication},
year = {2026},
howpublished = {\url{https://pith.science/paper/J6JZGOYB}},
note = {Machine review of arXiv:2502.00140}
}
abstract
While Graph Neural Networks (GNNs) have achieved remarkable success, their design largely relies on empirical intuition rather than theoretical understanding. In this paper, we present a comprehensive analysis of GNN behavior through three fundamental aspects: (1) we establish that \textbf{$k$-layer} Message Passing Neural Networks efficiently aggregate \textbf{$k$-hop} neighborhood information through iterative computation, (2) analyze how different loop structures influence neighborhood computation, and (3) examine behavior across structure-feature hybrid and structure-only tasks. For deeper GNNs, we demonstrate that gradient-related issues, rather than just over-smoothing, can significantly impact performance in sparse graphs. We also analyze how different normalization schemes affect model performance and how GNNs make predictions with uniform node features, providing a theoretical framework that bridges the gap between empirical success and theoretical understanding.
Figures
Figures from the paper (5 more)
Forward citations
Cited by 1 Pith paper
-
Machine learning of electronic structure and atomistic properties from the external potential
Representing the external nuclear potential as an atomic-orbital matrix gives a symmetry-preserving ML input whose powers implement equivariant message passing and capture long-range interactions in property and opera...
Reference graph
Works this paper leans on
-
[1]
Graph Neural Networks Use Graphs When They Shouldn 't
Bechler-Speicher, M., Amos, I., Gilad-Bachrach, R., and Globerson, A. Graph Neural Networks Use Graphs When They Shouldn 't. June 2024
work page 2024
-
[2]
I., Bronstein, M., Webb, S., and Rossi, E
Chamberlain, B., Rowbottom, J., Gorinova, M. I., Bronstein, M., Webb, S., and Rossi, E. GRAND : Graph Neural Diffusion . In Proceedings of the 38th International Conference on Machine Learning , pp.\ 1407--1418. PMLR, July 2021. ISSN: 2640-3498
work page 2021
-
[3]
Understanding the Representation Power of Graph Neural Networks in Learning Graph Topology
Dehmamy, N., Barabasi, A.-L., and Yu, R. Understanding the Representation Power of Graph Neural Networks in Learning Graph Topology . In Advances in Neural Information Processing Systems , volume 32. Curran Associates, Inc., 2019
work page 2019
-
[4]
Gilmer, J., Schoenholz, S. S., Riley, P. F., Vinyals, O., and Dahl, G. E. Neural Message Passing for Quantum Chemistry . In Proceedings of the 34th International Conference on Machine Learning , pp.\ 1263--1272. PMLR, July 2017. ISSN: 2640-3498
work page 2017
-
[5]
Inductive representation learning on large graphs
Hamilton, W., Ying, Z., and Leskovec, J. Inductive representation learning on large graphs. Advances in neural information processing systems, 30, 2017
2017
-
[6]
Approximation Capabilities of Muitilayer Feedforward Networks
Hornik, K. Approximation Capabilities of Muitilayer Feedforward Networks
-
[7]
Multilayer feedforward networks are universal approximators
Hornik, K., Stinchcombe, M., and White, H. Multilayer feedforward networks are universal approximators. Neural networks, 2 0 (5): 0 359--366, 1989. ISSN 0893-6080
work page 1989
-
[8]
Scale-aware Message Passing For Graph Node Classification
Jiang, Q., Wang, C., Lones, M., and Pang, W. Scale Invariance of Graph Neural Networks , December 2024. arXiv:2411.19392 [cs]
work page Pith review arXiv 2024
Show all 25 references
-
[9]
Kipf, T. N. and Welling, M. Semi-supervised classification with graph convolutional networks. arXiv preprint arXiv:1609.02907, 2016
2016 arXiv
-
[10]
Deeper Insights Into Graph Convolutional Networks for Semi - Supervised Learning
Li, Q., Han, Z., and Wu, X.-m. Deeper Insights Into Graph Convolutional Networks for Semi - Supervised Learning . Proceedings of the AAAI Conference on Artificial Intelligence, 32 0 (1), April 2018. ISSN 2374-3468. doi:10.1609/aaai.v32i1.11604. Number: 1
2018 doi
-
[11]
Gated Graph Sequence Neural Networks , September 2017
Li, Y., Tarlow, D., Brockschmidt, M., and Zemel, R. Gated Graph Sequence Neural Networks , September 2017. URL http://arxiv.org/abs/1511.05493. arXiv:1511.05493 [cs]
2017 arXiv
-
[12]
K., Liu, X., and Murata, T
Maurya, S. K., Liu, X., and Murata, T. Simplifying approach to node classification in Graph Neural Networks . Journal of Computational Science, 62: 0 101695, July 2022. ISSN 1877-7503. doi:10.1016/j.jocs.2022.101695
2022
-
[13]
Mitigating Oversmoothing Through Reverse Process of GNNs for Heterophilic Graphs
Park, M., Heo, J., and Kim, D. Mitigating Oversmoothing Through Reverse Process of GNNs for Heterophilic Graphs . In Proceedings of the 41st International Conference on Machine Learning , June 2024
2024
-
[14]
C.-C., Lei, Y., and Yang, B
Pei, H., Wei, B., Chang, K. C.-C., Lei, Y., and Yang, B. Geom-gcn: Geometric graph convolutional networks. arXiv preprint arXiv:2002.05287, 2020
2002 arXiv
-
[15]
Rossi, E., Charpentier, B., Di Giovanni, F., Frasca, F., Günnemann, S., and Bronstein, M. M. Edge directionality improves learning on heterophilic graphs. In Learning on Graphs Conference , pp.\ 25--1. PMLR, 2024
2024
-
[16]
Multi-scale attributed node embedding
Rozemberczki, B., Allen, C., and Sarkar, R. Multi-scale attributed node embedding. Journal of Complex Networks, 9 0 (2): 0 cnab014, 2021
2021
-
[17]
K., Bronstein, M
Rusch, T. K., Bronstein, M. M., and Mishra, S. A Survey on Oversmoothing in Graph Neural Networks , March 2023. arXiv:2303.10993 [cs]
2023 arXiv
-
[18]
Graph Neural Networks with maximal independent set-based pooling: Mitigating over-smoothing and over-squashing
Stanovic, S., Gaüzère, B., and Brun, L. Graph Neural Networks with maximal independent set-based pooling: Mitigating over-smoothing and over-squashing. Pattern Recognition Letters, 187: 0 14--20, January 2025. ISSN 01678655. doi:10.1016/j.patrec.2024.11.004
2025 doi
-
[19]
S., and Lim, A
Tong, Z., Liang, Y., Sun, C., Rosenblum, D. S., and Lim, A. Directed graph convolutional network. arXiv preprint arXiv:2004.13970, 2020
2004 arXiv
-
[20]
C.-W., and Han, J
Xie, Y., Li, S., Yang, C., Wong, R. C.-W., and Han, J. When Do GNNs Work : Understanding and Improving Neighborhood Aggregation . IJCAI'20: Proceedings of the Twenty-Ninth International Joint Conference on Artificial Intelligence, \ IJCAI\ 2020 , 2020 0 (1), July 2020. doi:10....
2020 doi
-
[21]
How Powerful are Graph Neural Networks ?, February 2019
Xu, K., Hu, W., Leskovec, J., and Jegelka, S. How Powerful are Graph Neural Networks ?, February 2019. arXiv:1810.00826 [cs]
2019 arXiv
-
[22]
Magnet: A neural network for directed graphs
Zhang, X., He, Y., Brugnone, N., Perlmutter, M., and Hirn, M. Magnet: A neural network for directed graphs. Advances in neural information processing systems, 34: 0 27003--27015, 2021
2021
-
[23]
W., Rao, N., Katariya, S., Wang, Z., and Subbian, K
Zheng, W., Huang, E. W., Rao, N., Katariya, S., Wang, Z., and Subbian, K. Cold Brew : Distilling Graph Node Representations with Incomplete or Missing Neighborhoods , March 2022. URL http://arxiv.org/abs/2111.04840. arXiv:2111.04840 [cs]
2022 arXiv
-
[24]
Beyond homophily in graph neural networks: Current limitations and effective designs
Zhu, J., Yan, Y., Zhao, L., Heimann, M., Akoglu, L., and Koutra, D. Beyond homophily in graph neural networks: Current limitations and effective designs. Advances in neural information processing systems, 33: 0 7793--7804, 2020
2020
-
[25]
write newline
" write newline "" before.all 'output.state := FUNCTION n.dashify 't := "" t empty not t #1 #1 substring "-" = t #1 #2 substring "--" = not "--" * t #2 global.max substring 't := t #1 #1 substring "-" = "-" * t #2 global.max substring 't := while if t #1 #1 substring * t #2 gl...
Reviewed August 9, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.