Graph Neural Reaction Diffusion Models Submitted to the editors June 2023. Accepted in March 2024.
Abstract
The integration of Graph Neural Networks (GNNs) and Neural Ordinary and Partial Differential Equations has been extensively studied in recent years. GNN architectures powered by neural differential equations allow us to reason about their behavior, and develop GNNs with desired properties such as controlled smoothing or energy conservation. In this paper we take inspiration from Turing instabilities in a Reaction Diffusion (RD) system of partial differential equations, and propose a novel family of GNNs based on neural RD systems. We blackdemonstrate that our RDGNN is powerful for the modeling of various data types, from homophilic, to heterophilic, and spatio-temporal datasets. We discuss the theoretical properties of our RDGNN, its implementation, and show that it improves or offers competitive performance to state-of-the-art methods.
Benchmark Results
Tested on 8 datasets 32 baseline values reported too low
Tables show presented results in global context (best 3 papers + GCN + MLP).
Classic
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | GraphSAGE+CNA | 0.9418
±0.0033 | |
| 2 | GCN+CNA | 0.9366
±0.0048 | |
| 3 | TransformerConv | 0.9350
±0.0058 | |
| 21 | GCN | 0.9004 best reported | best at pub: Understanding Heterophily for Graph Neural Networks (ICML 2024) |
| 22 | RDGNN-I | 0.8991
±0.0110 | This paper |
| 29 | RDGNN-S | 0.8953
±0.0122 | This paper |
| 39 | SAGE | 0.8911
±0.0012 | |
| 44 | GAT | 0.8897
±0.0014 best reported | best at pub: Learned Low Precision Graph Neural Networks (2020) |
| 45 | GCNII | 0.8893
±0.0137 best reported | best at pub: Language is All a Graph Needs (2023) best today: Large Language Model Meets Graph Neural Network in Knowledge Distillation (AAAI 2024) Difference within reported uncertainty; not counted as underperformance. |
| 51 | H2GCN | 0.8882
±0.0011 best reported | best at pub: Feature Transportation Improves Graph Neural Networks (AAAI 2023) best today: Adaptive Tokenization: On the Hop-Overpriority Problem in Tokenized Graph Learning Models (2025) Difference within reported uncertainty; not counted as underperformance. |
| 59 | PDE-GCNM | 0.8860 | This paper |
| 62 | GREAD | 0.8857
±0.0066 | This paper |
| 79 | ACMII-GCN | 0.8825
±0.0096 | This paper |
| 82 | WRGAT | 0.8820
±0.0226 | This paper |
| 96 | GRAFF | 0.8801
±0.0103 | This paper |
| 103 | GGCN | 0.8795
±0.0105 | This paper |
| 113 | MixHop | 0.8761
±0.0203 | This paper |
| 138 | NSD | 0.8714
±0.0113 | This paper |
| 168 | DMP | 0.8652 | This paper |
| 207 | Geom-GCN | 0.8540
±0.0040 best reported | best at pub: The Expressive Power of Graph Neural Networks: A Survey (2023) Difference within reported uncertainty; not counted as underperformance. |
| 243 | LINKX | 0.8464
±0.0113 | This paper |
| 724 | MLP | 0.7768 |
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | CNA | 0.9575
±0.0058 | |
| 2 | IS-GIB | 0.9390
±0.0187 | |
| 3 | EERM | 0.9112
±0.0145 | |
| 35 | GAT | 0.8080
±0.0026 best reported | best at pub: Feature Expansion for Graph Neural Networks (ICML 2023) |
| 36 | GCNII | 0.8076
±0.0030 best reported | best at pub: Feature Expansion for Graph Neural Networks (ICML 2023) |
| 44 | GCN | 0.8033
±0.0088 best reported | best at pub: Gophormer: Ego-Graph Transformer for Node Classification (2021) |
| 61 | GraphSAGE | 0.7923
±0.0053 | |
| 80 | PDE-GCNM | 0.7848 | This paper |
| 84 | RDGNN-I | 0.7834
±0.0155 | This paper |
| 88 | Geom-GCN | 0.7799
±0.0115 | This paper |
| 94 | RDGNN-S | 0.7782
±0.0148 | This paper |
| 99 | GREAD | 0.7760
±0.0181 | This paper |
| 104 | GRAFF | 0.7730
±0.0185 | This paper |
| 107 | NSD | 0.7714
±0.0157 | This paper |
| 108 | GGCN | 0.7714
±0.0145 | This paper |
| 110 | ACMII-GCN | 0.7712
±0.0158 | This paper |
| 111 | H2GCN | 0.7711
±0.0157 | This paper |
| 117 | DMP | 0.7687 | This paper |
| 119 | WRGAT | 0.7681
±0.0189 | This paper |
| 126 | MixHop | 0.7626
±0.0295 | This paper |
| 312 | LINKX | 0.7319
±0.0099 | This paper |
| 350 | MLP | 0.7289
±0.0010 |
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | LGTL | 0.9518
±0.0011 | |
| 2 | MixHop (LO) | 0.9513
±0.0038 | |
| 3 | GraphSAGE (LO) | 0.9511
±0.0036 | |
| 16 | H2GCN | 0.9361
±0.0013 best reported | best at pub: Optimization-Induced Graph Implicit Nonlinear Diffusion (ICML 2022) |
| 18 | GCN | 0.9296
±0.0015 best reported | best at pub: Node Identifiers: Compact, Discrete Representations for Efficient Graph Learning (ICLR 2024) |
| 20 | GAT | 0.9233
±0.0018 best reported | best at pub: Tackling Over-Smoothing for General Graph Convolutional Networks (2020) |
| 55 | RDGNN-I | 0.9037
±0.0043 | This paper |
| 58 | GCNII | 0.9030
±0.0043 | This paper |
| 62 | GREAD | 0.9023
±0.0055 | This paper |
| 66 | RDGNN-S | 0.9011
±0.0062 | This paper |
| 69 | Geom-GCN | 0.9005
±0.0047 | This paper |
| 70 | MixHop | 0.9004
±0.0141 best reported | best at pub: Language is All a Graph Needs (2023) |
| 71 | GRAFF | 0.9004
±0.0041 | This paper |
| 75 | PDE-GCNM | 0.8993 | This paper |
| 83 | ACMII-GCN | 0.8971
±0.0048 | This paper |
| 88 | NSD | 0.8949
±0.0040 | This paper |
| 96 | DMP | 0.8927 | This paper |
| 98 | GGCN | 0.8915
±0.0037 | This paper |
| 116 | SAGE | 0.8881
±0.0015 | |
| 129 | WRGAT | 0.8852
±0.0092 | This paper |
| 162 | LINKX | 0.8786
±0.0077 | This paper |
| 210 | MLP | 0.8690
±0.0074 |
Heterophilic Graphs
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | P^2GNN | 0.9672
±0.0147 | |
| 2 | DoG | 0.9660
±0.0210 | |
| 3 | GPF-plus | 0.9583
±0.0419 | |
| 9 | RDGNN-S | 0.9459
±0.0597 | This paper |
| 19 | RDGNN-I | 0.9351
±0.0593 | This paper |
| 24 | PDE-GCNM | 0.9324 | This paper |
| 28 | GPRGNN | 0.9295
±0.0131 best reported | best at pub: Node-oriented Spectral Filtering for Graph Neural Networks (2022) best today: GrokFormer: Graph Fourier Kolmogorov-Arnold Transformers (ICML 2024) |
| 34 | MLP | 0.9230
±0.0070 | |
| 39 | GGCN | 0.9213
±0.0305 best reported | best at pub: Global-Local Graph Neural Networks for Node-Classification (2024) best today: Graph Wave Networks (WWW 2025) |
| 40 | H2GCN | 0.9189
±0.0393 best reported | best at pub: Global-Local Graph Neural Networks for Node-Classification (2024) best today: GrokFormer: Graph Fourier Kolmogorov-Arnold Transformers (ICML 2024) Difference within reported uncertainty; not counted as underperformance. |
| 53 | GraphSAGE | 0.9000
±0.0365 | |
| 57 | GREAD | 0.8973
±0.0449 | prev best: Understanding Oversmoothing in Diffusion-Based GNNs From the Perspective of Operator Semigroup Theory (KDD 2024) Improvement within reported uncertainty; not highlighted. |
| 61 | FAGCN | 0.8961
±0.0152 best reported | best at pub: Node-oriented Spectral Filtering for Graph Neural Networks (2022) |
| 68 | DMP | 0.8919 | prev best: Heterophily-Aware Graph Attention Network (2023) Improvement within reported uncertainty; not highlighted. |
| 70 | GCNII | 0.8918
±0.0443 best reported | best at pub: Global-Local Graph Neural Networks for Node-Classification (2024) best today: Piecewise Constant Spectral Graph Neural Network (TMLR 2025) |
| 82 | GRAFF | 0.8838
±0.0453 | This paper |
| 83 | ACMII-GCN | 0.8838
±0.0343 | This paper |
| 98 | G2 | 0.8757
±0.0386 | This paper |
| 104 | WRGAT | 0.8698 best reported | best at pub: Global-Local Graph Neural Networks for Node-Classification (2024) Difference within reported uncertainty; not counted as underperformance. |
| 115 | ACMP-GCN | 0.8620
±0.0300 | This paper |
| 121 | NSD | 0.8595
±0.0551 | This paper |
| 137 | GraphCON | 0.8540
±0.0420 | This paper |
| 181 | MagNet | 0.8330
±0.0610 | This paper |
| 211 | GAT | 0.8162
±0.0645 best reported | best at pub: A Quasi-Wasserstein Loss for Learning Graph Neural Networks (WWW 2023) best today: Equivariance Everywhere All At Once: A Recipe for Graph Foundation Models (2025) |
| 217 | GRAND | 0.8140
±0.0350 best reported | best at pub: Unifying over-smoothing and over-squashing in graph neural networks: A physics informed approach and beyond (2023) best today: When Graph Neural Networks Meet Dynamic Mode Decomposition (ICLR 2024) Difference within reported uncertainty; not counted as underperformance. |
| 257 | GCN | 0.7933
±0.0447 best reported | best at pub: Restructuring Graphs for Higher Homophily via Adaptive Spectral Clustering (AAAI 2022) best today: Graph Wave Networks (WWW 2025) |
| 276 | MixHop | 0.7784
±0.0773 | This paper |
| 325 | LINKX | 0.7460
±0.0837 | This paper |
| 415 | Geom-GCN | 0.6760 best reported | best at pub: Non-Local Graph Neural Networks (2020) Difference within reported uncertainty; not counted as underperformance. |
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | P^2GNN | 0.9541
±0.0272 | |
| 2 | ACM-GCN | 0.9480
±0.0380 | |
| 3 | RF-GCN | 0.9429 | |
| 12 | RDGNN-I | 0.9272
±0.0588 | This paper |
| 14 | RDGNN-S | 0.9243
±0.0551 | This paper |
| 15 | GGCN | 0.9230
±0.0367 best reported | best at pub: Understanding convolution on graphs via energies (TMLR 2022) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) Difference within reported uncertainty; not counted as underperformance. |
| 22 | GPRGNN | 0.9157
±0.0196 best reported | best at pub: Simplified PCNet with Robustness (2024) Difference within reported uncertainty; not counted as underperformance. |
| 24 | MLP | 0.9136
±0.0070 | |
| 29 | PDE-GCNM | 0.8973 | This paper |
| 31 | DMP | 0.8919 | prev best: Heterophily-Aware Graph Attention Network (2023) Improvement within reported uncertainty; not highlighted. |
| 35 | GCNII | 0.8918
±0.0396 best reported | best at pub: Convolutional Neural Networks on Graphs with Chebyshev Approximation, Revisited (NeurIPS 2022) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) |
| 38 | FAGCN | 0.8852
±0.0133 best reported | best at pub: Node-oriented Spectral Filtering for Graph Neural Networks (2022) |
| 49 | G2 | 0.8730
±0.0484 | This paper |
| 52 | GREAD | 0.8703
±0.0495 | |
| 62 | NSD | 0.8649
±0.0471 | This paper |
| 63 | ACMII-GCN | 0.8649
±0.0673 | This paper |
| 68 | H2GCN | 0.8623
±0.0471 best reported | best at pub: Design your own universe: a physics-informed agnostic method for enhancing graph neural networks (2024) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) Difference within reported uncertainty; not counted as underperformance. |
| 83 | ACMP-GCN | 0.8540
±0.0700 | This paper |
| 106 | MagNet | 0.8430
±0.0700 | prev best: Breaking the Entanglement of Homophily and Heterophily in Semi-supervised Node Classification (2023) Improvement within reported uncertainty; not highlighted. |
| 107 | GraphCON | 0.8430
±0.0480 | This paper |
| 112 | GRAFF | 0.8405
±0.0610 | This paper |
| 123 | LINKX | 0.8346
±0.0061 best reported | best at pub: Differentiable Cluster Graph Neural Network (2024) Difference within reported uncertainty; not counted as underperformance. |
| 138 | SAGE | 0.8269 | |
| 142 | GRAND | 0.8246
±0.0709 best reported | best at pub: Dirichlet Energy Enhancement of Graph Neural Networks by Framelet Augmentation (2023) best today: Graph Navier–Stokes Networks (2026) Difference within reported uncertainty; not counted as underperformance. |
| 158 | GAT | 0.8190
±0.0540 best reported | best at pub: Restructuring Graphs for Higher Homophily via Adaptive Spectral Clustering (AAAI 2022) |
| 160 | WRGAT | 0.8162
±0.0390 | This paper |
| 170 | GCN | 0.8082
±0.0360 best reported | best at pub: Differentiable Cluster Graph Neural Network (2024) best today: Bayesian Neighborhood Adaptation for Graph Neural Networks (TMLR 2025) |
| 183 | MixHop | 0.7950
±0.0630 best reported | best at pub: Label-Wise Graph Convolutional Network for Heterophilic Graphs (2021) Difference within reported uncertainty; not counted as underperformance. |
| 342 | Geom-GCN | 0.6412 best reported | best at pub: Steering Graph Neural Networks with Pinning Control (2023) Difference within reported uncertainty; not counted as underperformance. |
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | DiGGR | 0.4535
±0.0353 | |
| 2 | TRIGON | 0.4402
±0.0028 | |
| 3 | DoG | 0.4320
±0.0120 | |
| 38 | GCNII | 0.4082
±0.0179 best reported | best at pub: Feature Expansion for Graph Neural Networks (ICML 2023) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) |
| 42 | MLP | 0.4026
±0.0108 | |
| 52 | GPRGNN | 0.3992
±0.0067 best reported | best at pub: LON-GNN: Spectral GNNs with Learnable Orthonormal Basis (2023) best today: GrokFormer: Graph Fourier Kolmogorov-Arnold Transformers (ICML 2024) |
| 56 | FAGCN | 0.3970 best reported | best at pub: Beyond Low-frequency Information in Graph Convolutional Networks (AAAI 2021) |
| 62 | GGCN | 0.3929
±0.0048 best reported | best at pub: Feature Transportation Improves Graph Neural Networks (AAAI 2023) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) Difference within reported uncertainty; not counted as underperformance. |
| 68 | GraphSAGE | 0.3899
±0.0085 | |
| 71 | H2GCN | 0.3885
±0.0117 best reported | best at pub: ES-GNN: Generalizing Graph Neural Networks Beyond Homophily with Edge Splitting (2022) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) |
| 72 | RDGNN-I | 0.3869
±0.0141 | This paper |
| 75 | WRGAT | 0.3860
±0.0110 best reported | best at pub: ES-GNN: Generalizing Graph Neural Networks Beyond Homophily with Edge Splitting (2022) |
| 78 | RDGNN-S | 0.3855
±0.0134 | This paper |
| 87 | Geom-GCN | 0.3810 best reported | best at pub: Beyond Low-frequency Information in Graph Convolutional Networks (AAAI 2021) |
| 96 | GREAD | 0.3790
±0.0117 | prev best: Understanding Oversmoothing in Diffusion-Based GNNs From the Perspective of Operator Semigroup Theory (KDD 2024) Improvement within reported uncertainty; not highlighted. |
| 100 | NSD | 0.3781
±0.0115 best reported | best at pub: Graph Parsing Networks (ICLR 2024) Difference within reported uncertainty; not counted as underperformance. |
| 120 | G2 | 0.3730
±0.0101 | This paper |
| 122 | DMP | 0.3726
±0.0179 best reported | best at pub: Heterophily-Aware Graph Attention Network (2023) Difference within reported uncertainty; not counted as underperformance. |
| 129 | GRAFF | 0.3711
±0.0108 | This paper |
| 131 | ACMII-GCN | 0.3709
±0.0132 | This paper |
| 135 | GCN | 0.3690 best reported | best at pub: Beyond Low-frequency Information in Graph Convolutional Networks (AAAI 2021) |
| 162 | MixHop | 0.3648
±0.0120 best reported | best at pub: Towards Deep Attention in Graph Neural Networks: Problems and Remedies (ICML 2023) |
| 180 | LINKX | 0.3610
±0.0155 | This paper |
| 182 | GAT | 0.3605
±0.0035 best reported | best at pub: Hierarchical Graph Transformer with Adaptive Node Sampling (NeurIPS 2022) |
| 203 | GRAND | 0.3562
±0.0101 | This paper |
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | IBG-NN | 0.7763
±0.0179 | |
| 2 | Dir-GNN+CNA | 0.7747
±0.0128 | |
| 3 | DHGNN | 0.7684
±0.0161 | |
| 26 | ACMII-GCN | 0.6740
±0.0221 | This paper |
| 27 | RDGNN-S | 0.6596
±0.0211 | This paper |
| 28 | RDGNN-I | 0.6562
±0.0233 | This paper |
| 32 | G2 | 0.6426
±0.0238 | This paper |
| 39 | H2GCN | 0.6190
±0.0140 best reported | best at pub: Optimization-Induced Graph Implicit Nonlinear Diffusion (ICML 2022) best today: Efficient Learning on Large Graphs using a Densifying Regularity Lemma (2025) |
| 40 | LINKX | 0.6181
±0.0180 | This paper |
| 46 | GCN | 0.6067 best reported | best at pub: Understanding Heterophily for Graph Neural Networks (ICML 2024) |
| 51 | GREAD | 0.5922
±0.0144 | This paper |
| 53 | GRAFF | 0.5901
±0.0131 | This paper |
| 65 | NSD | 0.5634
±0.0132 | This paper |
| 67 | GAT | 0.5630
±0.0220 best reported | best at pub: Restructuring Graphs for Higher Homophily via Adaptive Spectral Clustering (AAAI 2022) |
| 69 | GGCN | 0.5567
±0.0187 best reported | best at pub: GPNet: Simplifying Graph Neural Networks via Multi-channel Geometric Polynomials (2022) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) Difference within reported uncertainty; not counted as underperformance. |
| 76 | GPRGNN | 0.5435
±0.0087 best reported | best at pub: LON-GNN: Spectral GNNs with Learnable Orthonormal Basis (2023) best today: Directed Homophily-Aware Graph Neural Network (2025) |
| 96 | WRGAT | 0.4885
±0.0078 | This paper |
| 101 | DMP | 0.4726 | This paper |
| 120 | SAGE | 0.4460
±0.0130 | |
| 127 | MixHop | 0.4380
±0.0148 | This paper |
| 139 | FAGCN | 0.4259
±0.0069 | This paper |
| 164 | GRAND | 0.4005
±0.0150 | This paper |
| 172 | MLP | 0.3930
±0.0079 | |
| 176 | GCNII | 0.3881
±0.0197 best reported | best at pub: GPNet: Simplifying Graph Neural Networks via Multi-channel Geometric Polynomials (2022) best today: $P^2$GNN: Two Prototype Sets to boost GNN Performance (2026) Difference within reported uncertainty; not counted as underperformance. |
| 179 | Geom-GCN | 0.3832
±0.0092 | This paper |
LINKX Benchmarks
| # | Model | Accuracy ↑ | Paper |
|---|---|---|---|
| 1 | LargeScaleNet | 0.6582
±0.0036 | |
| 2 | S² DirGCN | 0.6495
±0.0033 | |
| 3 | FaberNet | 0.6462
±0.0101 | |
| 11 | RDGNN-I | 0.5846
±0.0059 | This paper |
| 13 | LINKX | 0.5600
±0.0134 | This paper |
| 15 | RDGNN-S | 0.5516
±0.0046 | This paper |
| 17 | LINK | 0.5397
±0.0018 | This paper |
| 19 | MixHop | 0.5181
±0.0017 | This paper |
| 22 | GAT | 0.4937
±0.0020 best reported | best at pub: Uplifting Message Passing Neural Network with Graph Original Information (2022) |
| 23 | H2GCN | 0.4910
±0.0010 best reported | best at pub: Label-Wise Graph Convolutional Network for Heterophilic Graphs (2021) Difference within reported uncertainty; not counted as underperformance. |
| 26 | GCNII | 0.4721
±0.0028 | This paper |
| 30 | L Prop (2 hop) | 0.4607
±0.0015 | This paper |
| 31 | GCN | 0.4602
±0.0026 | This paper |
| 32 | GPR-GNN | 0.4597
±0.0026 best reported | best at pub: Convolutional Neural Networks on Graphs with Chebyshev Approximation, Revisited (NeurIPS 2022) |
| 40 | C&S (2 hop) | 0.4217
±0.0027 | This paper |
| 53 | APPNP | 0.3815
±0.0026 | This paper |
| 56 | MLP | 0.3670
±0.0021 | This paper |
| 59 | SGC (2 hop) | 0.3227
±0.0006 | This paper |
