Predicting The One-Particle Density Matrix With Machine Learning

S. Hazra,U. Patil,S. Sanvito
DOI: https://doi.org/10.1021/acs.jctc.4c00042
2024-06-05
Abstract:Two of the most widely used electronic structure theory methods, namely Hartree-Fock and Kohn-Sham density functional theory, both requires the iterative solution of a set of Schrödinger-like equations. The speed of convergence of such self-consistent field process depends on the complexity of the system under investigation and on the initial guess for the density matrix. An initial density matrix close to the ground-state one will effectively allow one to cut out many of the self-consistent steps necessary to achieve convergence. Here, we predict the density matrix of Kohn-Sham density functional theory by constructing a neural network, which uses the atomic positions as only information. Such neural network provides an initial guess for the density matrix far superior to any other recipes available, in particular for molecules with metallic bonds. Furthermore, the quality of such neural-network density matrix is good enough for the evaluation of interatomic forces. This allows us to run accelerated {\it ab-initio} molecular dynamics with little to no self-consistent steps.
Chemical Physics,Computational Physics
What problem does this paper attempt to address?
The purpose of this paper is to predict the one-particle density matrix in Kohn-Sham density functional theory (DFT) using machine learning. In electronic structure theory methods such as Hartree-Fock and DFT, solving the iterative process similar to the Schrödinger equation is a key step, and its convergence speed depends on the system complexity, the algorithm used, and the quality of the initial density matrix. The paper proposes to construct a neural network that uses only atomic position information to predict the density matrix, in order to provide a better initial guess than existing methods. This method can accelerate the calculation of ab-initio molecular dynamics and reduce the required self-consistent steps. In the study, they compared the initial density matrices generated by different strategies and used three molecules (H2O, S2O, and [Fe(H2O)6]2+) to validate the quality of the density matrix predicted by the neural network, which can be used to evaluate interatomic forces and non-self-consistent energy calculations. Finally, the paper demonstrates the results of non-self-consistent molecular dynamics simulations using machine learning predicted density matrices, compared to fully converged simulations.