English

Accelerating inverse Kohn-Sham calculations using reduced density matrices

Chemical Physics 2024-08-06 v1

Abstract

The Ryabinkin-Kohut-Staroverov (RKS) and Kanungo-Zimmerman-Gavini (KZG) methods offer two approaches to find exchange-correlation (XC) potentials from ground state densities. The RKS method utilizes the one- and two-particle reduced density matrices to alleviate any numerical artifacts stemming from a finite basis (e.g., Gaussian- or Slater-type orbitals). The KZG approach relies solely on the density to find the XC potential, by combining a systematically convergent finite-element basis with appropriate asymptotic correction on the target density. The RKS method, being designed for a finite basis, offers computational efficiency. The KZG method, using a complete basis, provides higher accuracy. In this work, we combine both the methods to simultaneously afford accuracy and efficiency. In particular, we use the RKS solution as initial guess to the KZG method to attain a significant 311×3-11\times speedup. This work also presents a direct comparison of the XC potentials from the RKS and the KZG method and their relative accuracy on various weakly and strongly correlated molecules, using their ground state solutions from accurate configuration interaction calculations solved in a Slater orbital basis.

Keywords

Cite

@article{arxiv.2408.02342,
  title  = {Accelerating inverse Kohn-Sham calculations using reduced density matrices},
  author = {Bikash Kanungo and Soumi Tribedi and Paul M. Zimmerman and Vikram Gavini},
  journal= {arXiv preprint arXiv:2408.02342},
  year   = {2024}
}
R2 v1 2026-06-28T18:04:01.723Z