Many-Body Permutationally Invariant Polynomial Neural Network Potential Energy Surface for N4

Jun Li, Zoltan Varga, Donald G. Truhlar, Hua Guo

Research output: Contribution to journalArticlepeer-review

3 Scopus citations


A potential energy surface (PES) for high-energy collisions between nitrogen molecules is useful for modeling chemical dynamics in shock waves and plasmas. In the present work, we fit the many-body (MB) component of the ground-state PES of N4 to an analytic function using neural networks (NNs) with permutationally invariant polynomials (PIPs). The fitting dataset of the N4 system is an extension of one used previously, extended with 4859 new CASPT2 points and 13 new CCSD(T) points to reach a total of 21 406 points. The MB-PIP-NN fit required a very complete coverage of the geometry domain in order to get a physical fit, and we devised several tactical steps to achieve this, including trajectory calculations, the comparison of the NN fit with PIP fits with mixed-exponential Gaussian bond order variables, and searching geometry regions with sparse data coverage. With these efforts, the final dataset is more suitable for a NN fit. The energy range of the dataset is much wider than those used in other systems previously fitted by NNs, and there are more rugged surface regions than usual due to locally avoided crossings. The performance of the new MB-PIP-NN fit is compared to that of another new fit to the same data by least-squares methods not employing NNs, and the mean unsigned deviation from the electronic structure calculations is reduced by a factor of 3, although the computing time to calculate a force for a trajectory calculation is larger by an order of magnitude. The accuracy achievable with the NN fit for this difficult system is very impressive, and we anticipate that the NN method can give useful fits to difficult cases that cannot be achieved by more conventional methods.

Original languageEnglish (US)
Pages (from-to)4822-4832
Number of pages11
JournalJournal of Chemical Theory and Computation
Issue number8
StatePublished - Aug 11 2020

Bibliographical note

Funding Information:
The authors are grateful to Rubén Meana-Pañeda and Yuliya Paukku for helpful contributions at the beginning of this project and to Linyao Zhang for his help with the diatom–diatom trajectory calculations. J.L. is supported by the National Natural Science Foundation of China (grant nos. 21573027 and 21973009), the Chongqing Municipal Natural Science Foundation (grant no. cstc2019jcyj-msxmX0087), and Chongqing Talents Plan for Young Talents (CQYC201905047). The work of Z.V. is supported by the Air Force Office of Scientific Research (grant no. FA9550-19-1-0219). The work of D.G.T. and H.G. is supported in part by the U.S. Department of Energy, Office of Science, Office of Basic Energy Sciences under award DE-SC0015997. The Alexander von Humboldt Foundation is thanked for a Humboldt Fellowship for Experienced Researchers (to J.L.) and a Humboldt Research Award (to H.G.).

Publisher Copyright:
© 2020 American Chemical Society.

PubMed: MeSH publication types

  • Journal Article


Dive into the research topics of 'Many-Body Permutationally Invariant Polynomial Neural Network Potential Energy Surface for N<sub>4</sub>'. Together they form a unique fingerprint.

Cite this