Transferable Dynamic Molecular Charge Assignment Using Deep Neural Networks

• Published in: Journal of chemical theory and computation 2018-09, Vol.14 (9), p.4687-4698
• Main Authors: Nebgen, Benjamin, Lubbers, Nicholas, Smith, Justin S, Sifain, Andrew E, Lokhov, Andrey, Isayev, Olexandr, Roitberg, Adrian E, Barros, Kipton, Tretiak, Sergei
• Format: Article
• Language: English
• Subjects: Artificial neural networks; Benchmarks; Charge transfer; Chemical bonds; Computer simulation; Density functional theory; Extensibility; Infrared spectroscopy; INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CHEMISTRY; Machine learning; MATHEMATICS AND COMPUTING; Neural networks; Organic chemistry; Parameterization; Proteins; Spectra
• ISSN: 1549-9618; 1549-9626
• DOI: 10.1021/acs.jctc.8b00524

The ability to accurately and efficiently compute quantum-mechanical partial atomistic charges has many practical applications, such as calculations of IR spectra, analysis of chemical bonding, and classical force field parametrization. Machine learning (ML) techniques provide a possible avenue for the efficient prediction of atomic partial charges. Modern ML advances in the prediction of molecular energies [i.e., the hierarchical interacting particle neural network (HIP-NN)] has provided the necessary model framework and architecture to predict transferable, extensible, and conformationally dynamic atomic partial charges based on reference density functional theory (DFT) simulations. Utilizing HIP-NN, we show that ML charge prediction can be highly accurate over a wide range of molecules (both small and large) across a variety of charge partitioning schemes such as the Hirshfeld, CM5, MSK, and NBO methods. To demonstrate transferability and size extensibility, we compare ML results with reference DFT calculations on the COMP6 benchmark, achieving errors of 0.004e– (elementary charge). This is remarkable since this benchmark contains two proteins that are multiple times larger than the largest molecules in the training set. An application of our atomic charge predictions on nonequilibrium geometries is the generation of IR spectra for organic molecules from dynamical trajectories on a variety of organic molecules, which show good agreement with calculated IR spectra with reference method. Critically, HIP-NN charge predictions are many orders of magnitude faster than direct DFT calculations. These combined results provide further evidence that ML (specifically HIP-NN) provides a pathway to greatly increase the range of feasible simulations while retaining quantum-level accuracy.