ICML 2025

Efficient and Scalable Density Functional Theory Hamiltonian Prediction through Adaptive Sparsity

Erpai Luo, Xinran Wei, Lin Huang, Yunyang Li, Hang Yang, Zun Wang, Chang Liu, Zaishuo Xia, Jia Zhang, Bin Shao

Main figure for Efficient and Scalable Density Functional Theory Hamiltonian Prediction through Adaptive Sparsity

Abstract

Hamiltonian matrix prediction is pivotal in computational chemistry, serving as the foundation for determining a wide range of molecular properties. While SE(3) equivariant graph neural networks have achieved remarkable success in this domain, their substantial computational cost--driven by high-order tensor product (TP) operations--restricts their scalability to large molecular systems with extensive basis sets. To address this challenge, we introduce SPHNet, an efficient and scalable equivariant network, that incorporates adaptive SParsity into Hamiltonian prediction. SPHNet employs two innovative sparse gates to selectively constrain non-critical interaction combinations, significantly reducing tensor product computations while maintaining accuracy. To optimize the sparse representation, we develop a Three-phase Sparsity Scheduler, ensuring stable convergence and achieving high performance at sparsity rates of up to 70%. Extensive evaluations on QH9 and PubchemQH datasets demonstrate that SPHNet achieves state-of-the-art accuracy while providing up to a 7x speedup over existing models. Beyond Hamiltonian prediction, the proposed sparsification techniques also hold significant potential for improving the efficiency and scalability of other SE(3) equivariant networks, further broadening their applicability and impact. Our code can be found at https://github.com/microsoft/SPHNet.

Proposed method: SPHNet

Benchmark Results

Tested on 1 dataset Non-standard metric on 1 dataset

Tables show presented results in global context (best 3 papers + GCN + MLP).

Quantum Chemistry

MD17

Global leaderboard →

Default metric is MAE. The paper reports 0 values on MAE , 4 values on H , 4 values on epsilon and 4 values on psi .

Sorted by H. Highlighted rows are from other papers.

#ModelMAE ↓ defaultH ↓epsilon ↓psi ↓Paper
1
PhiSNet (ori)
1.07e-05
—
——
2
QHNet
—
1.08e-05
3.38e-05
0.9999
MAE —
epsilon 3.38e-05
psi 0.9999
This paper
3
DEQHNet
1.51e-05
—
——
4
PhiSNet
—
1.57e-05
8.55e-05
1.0000
MAE —
epsilon 8.55e-05
psi 1.0000
This paper
6
SPHNet
—
1.94e-05
8.23e-05
1.0000
MAE —
epsilon 8.23e-05
psi 1.0000
This paper
7
SchNOrb
—
0.0002
0.0003
1.0000
MAE —
epsilon 0.0003
psi 1.0000
This paper
SOTA performance (not significantly different from the best) Proposed method (this paper) Baseline / comparison (this paper) Global context (other papers)