Equivariant Graph Neural Networks (GNNs) are essential for physically consistent molecular simulations but suffer from high computational costs and memory bottlenecks, especially with high-order representations. While low-bit quantization offers a solution, applying it naively to rotation-sensitive features destroys the SO(3)-equivariant structure, leading to significant errors and violations of conservation laws. To address this issue, in this work, we propose a Geometric-Aware Quantization (GAQ) framework that compresses and accelerates equivariant models while rigorously preserving continuous symmetry in discrete spaces. Our approach introduces three key contributions: (1) a Magnitude-Direction Decoupled Quantization (MDDQ) scheme that separates invariant lengths from equivariant orientations to maintain geometric fidelity; (2) a symmetry-aware training strategy that treats scalar and vector features with distinct quantization schedules; and (3) a robust attention normalization mechanism to stabilize gradients in low-bit regimes. Experiments on the rMD17 benchmark demonstrate that our W4A8 models match the accuracy of FP32 baselines (9.31 meV vs. 23.20 meV) while reducing Local Equivariance Error (LEE) by over 30x compared to naive quantization. On consumer hardware, GAQ achieves 2.39x inference speedup and 4x memory reduction, enabling stable, energy-conserving molecular dynamics simulations for nanosecond timescales.
Figures & tables
Fig. 1: Overview of the proposed equivariant quantization framework. (a) MDDQ : Decouples equivariant vectors into magnitude ∥v∥ and direction u to preserve geometric orientation under low precision. (b) Branch-separated QAT : Treats invariant and equivariant features differently with a staged training schedule. (c) Robust Attention : Stabilizes dot-products via ℓ2 normalization and temperature scaling τ . (d) Results Summary : Our method achieves 2.37 – 2.73× faster inference and ∼4× memory reduction with ultra-low equivariance error ( LEE≈0.15 ).
Fig. 2: Detailed architecture and quantization pipeline. Inputs Z and r are processed through decoupled invariant and equivariant paths, employing specialized quantization strategies (Linear vs. MDDQ) for each branch.
Architecture
Cfull (FP32)
ℓmax
Cquant (k-bit)
Gain Cquant/Cfull
PaiNN [ 28 ]
O(n⟨N⟩4F)
1
ρk
O(n⟨N⟩4F)ρk
SpookyNet [ 29 ]
O(n⟨N⟩(ℓmax+1)2F)
2
ρk
O(n⟨N⟩(ℓmax+1)2F)ρk
NequIP [ 1 ]
O(n⟨N⟩(ℓmax+1)6F)
3
ρk
O(n⟨N⟩(ℓmax+1)6F)ρk
So3krates [ 2 ]
O(n⟨N⟩((ℓmax+1)2+F))
1
ρk
O(n⟨N⟩((ℓmax+1)2+F))ρk
TABLE I: Complexity with and without quantization. Cfull is per-layer asymptotic cost in full precision; Cquant is the cost when both weights and activations are k -bit (here ρk≡k/32 ). Quantization yields a constant-factor speedup ≈ρk without changing scaling in n , ⟨N⟩ , F , or ℓmax . For So3krates, we utilize ℓmax=1 as discussed in [ 2 ] for efficiency and compatibility with our vector quantization scheme, though the architecture supports higher orders.
Method
Bits (W/A)
E-MAE
F-MAE
Stability
FP32 Baseline
32 / 32
23.20
21.20
Stable
Naive INT8
8 / 8
118.20
102.39
Degraded
SVQ-KMeans
8 / 8
226.20 ∗
N/A
Diverged
Degree-Quant [ 22 ]
8 / 8
63.20
58.90
Stable
Ours (GAQ)
4 / 8
9.31
22.60
Stable
∗ Validation loss stagnated due to gradient fracture.
TABLE II: Performance comparison on Azobenzene (rMD17). Energy MAE in meV and Force MAE in meV/Å. ”Stability” indicates whether the training converged.
Fig. 3: Energy conservation in NVE dynamics (1 ns). The Naive INT8 model ( red ) exhibits catastrophic energy divergence and explosion within 100 ps due to symmetry breaking. In contrast, our W4A8 model ( orange ) maintains excellent stability comparable to the FP32 baseline ( gray ). Note that the realistic thermal fluctuations observed in our model confirm that the simulation retains correct physical dynamics without long-term drift ( <0.15 meV/atom/ps).
Method
LEE (meV/Å)
Remark
FP32 Baseline
≈0.0
Exact equivariance
Naive INT8
5.23
Broken symmetry
Degree-Quant
2.10
Partially preserved
Ours (W4A8)
0.15
Preserved
TABLE III: Symmetry analysis. LEE measures the sensitivity of predictions to input rotation (lower is better). Our method reduces the symmetry error by >30× compared to naive quantization.
Operation
FP32
Ours (W4A8)
Speedup
Memory I/O (Weights)
120.5
30.1
4.0 ×
Compute (GEMM)
45.0
25.0
1.8 ×
Quant Overhead
0.0
5.2
-
Attention
15.2
15.2
1.0 ×
Total Latency
180.7
75.5
2.39 ×
TABLE IV: Latency breakdown. Time in μ s. Our method achieves speedup primarily by reducing memory I/O.