PhysNet: A Neural Network for Predicting Energies, Forces, Dipole Moments and Partial Charges

Oliver T. Unke, M. Meuwly

Main figure for PhysNet: A Neural Network for Predicting Energies, Forces, Dipole Moments and Partial Charges

Abstract

In recent years, machine learning (ML) methods have become increasingly popular in computational chemistry. After being trained on appropriate ab initio reference data, these methods allow to accurately predict the properties of chemical systems, circumventing the need for explicitly solving the electronic Schr\"odinger equation. Because of their computational efficiency and scalability to large datasets, deep neural networks (DNNs) are a particularly promising ML algorithm for chemical applications. This work introduces, a DNN architecture designed for predicting energies, forces and dipole moments of chemical systems. \ achieves state-of-the-art performance on the QM9, MD17 and ISO17 benchmarks. Further, two new datasets are generated in order to probe the performance of ML models for describing chemical reactions, long-range interactions, and condensed phase systems. It is shown that explicitly including electrostatics in energy predictions is crucial for a qualitatively correct description of the asymptotic regions of a potential energy surface (PES). \ models trained on a systematically constructed set of small peptide fragments (at most eight heavy atoms) are able to generalize to considerably larger proteins like deca-alanine (Ala$_10$): The optimized geometry of helical Ala$_10$ predicted by \ is virtually identical to ab initio results ($RMSD=0.21$). By running unbiased molecular dynamics (MD) simulations of Ala$_10$ on the -PES in gas phase, it is found that instead of a helical structure, Ala$_10$ folds into a ``wreath-shaped'' configuration, which is more stable than the helical form by $0.46$ kcal mol$^-1$ according to the reference ab initio calculations.

Benchmark Results

Tested on 1 dataset

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

Quantum Chemistry

#ModelMAE ↓Paper
1
QHNet
0.0000
2
QHNetV2
0.0000
3
PhiSNet
0.0000
25
PhysNet-ens5
0.0200
This paper
24
GDML
actual #24 vs paper #24
0.0200
actual #24 vs paper #24
paper 0.0200
paper 0.0200 Δ at pub 0.0500
29
PhysNet
0.0300
This paper
39
DTNN
0.0400
This paper
46
SchNet
actual #46 vs paper #46
0.0500
actual #46 vs paper #46
paper 0.0500
paper 0.0500 Δ at pub 0.0300
62
HIP-NN
0.0600
This paper
SOTA performance (not significantly different from the best) Proposed method (this paper) Baseline / comparison (this paper) Global context (other papers) Best at pub — reported value improved over the best known when published