Density functional theory based neural network force fields from energy decompositions

In order to develop force fields (FF) for molecular dynamics simulations that retain the accuracy of ab initio density functional theory (DFT), we developed a machine learning protocol based on an energy decomposition scheme that extracts atomic energies from DFT calculations. Our DFT to FF (DFT2FF) approach provides almost hundreds of times more data for the DFT energies, which dramatically improves accuracy with less DFT calculations. In addition, we use piecewise cosine basis functions to systematically construct symmetry invariant features into the neural network model. We illustrate this DFT2FF approach for amorphous silicon where only 800 DFT configurations are sufficient to achieve an accuracy of 1 meV/atom for energy and 0.1 eV/A for forces. We then use the resulting FF model to calculate the thermal conductivity of amorphous Si based on long molecular dynamics simulations. The dramatic speedup in training in our DFT2FF protocol allows the adoption of a simulation paradigm where an accurate and problem specific FF for a given physics phenomenon is trained on-the-spot through a quick DFT precalculation and FF training.