Exchange Correlation Potentials from Full Configuration Interaction in a Slater Orbital Basis

Soumi Tribedi,Duy-Khoi Dang,Bikash Kanungo,Vikram Gavini,Paul M. Zimmerman
DOI: https://doi.org/10.1063/5.0157942
2023-07-08
Abstract:Ryabinkin-Kohut-Staroverov (RKS) theory builds a bridge between wave function theory and density functional theory by using quantities from the former to produce accurate exchange-correlation potentials needed by the latter. In this work, the RKS method is developed and tested alongside Slater atomic orbital basis functions for the first time. To evaluate this approach, Full Configuration Interaction computations in the Slater orbitals are employed to give quality input to RKS method, allowing full correlation to be present along with correct nuclei cusps and asymptotic decay of the wavefunction. The RKS method will be shown to be an efficient algorithm to arrive at exchange correlation potentials without unphysical artifacts in moderately-sized basis sets. Furthermore, enforcement of the nuclear cusp conditions will be shown to be vital for the success of the Slater-basis RKS method. Examples of weakly and strongly correlated molecular systems will demonstrate the main features of Slater RKS.
Chemical Physics
What problem does this paper attempt to address?
The problem that this paper attempts to solve is to develop and test a new method for extracting accurate exchange - correlation potentials from full configuration interaction (FCI) calculations of Slater atomic orbital basis sets. Specifically, the authors combine wave - function theory with density functional theory (DFT) through the Ryabinkin - Kohut - Staroverov (RKS) theory to generate more accurate exchange - correlation potentials. To achieve this goal, the authors implemented the RKS method in Slater orbital basis sets for the first time and used FCI calculations as input, ensuring complete correlation as well as correct nuclear cusp conditions and the asymptotic decay behavior of the wave function. ### Main problems: 1. **Improving the accuracy of exchange - correlation potentials**: Traditional DFT methods have deficiencies in dealing with self - interaction, delocalization errors, and strong - correlation problems, resulting in inaccurate exchange - correlation potentials. This paper aims to improve the accuracy of exchange - correlation potentials by improving the wave function and electron density. 2. **Satisfying nuclear cusp conditions**: Correctly describing the singular behavior of electron density near the nucleus (i.e., Kato cusp conditions) is crucial for obtaining accurate exchange - correlation potentials. This paper shows how to enforce these conditions in Slater orbital basis sets by modifying the SCF program. 3. **Efficient calculation**: Using GPU - accelerated algorithms and the resolution - of - the - identity (RI) to efficiently calculate Slater - type orbital integrals, making the RKS method more feasible in practical applications. ### Solutions: - **SlaterRKS method**: Combines the advantages of wave - function theory and DFT, and calculates the exchange - correlation potential by comparing the Kohn - Sham equation and the local energy balance equation in wave - function theory. - **FCI calculation**: Uses high - precision FCI calculations as a reference to ensure that the wave function contains all dynamic and static correlation effects. - **Slater - type orbitals**: Compared with commonly used Gaussian - type orbitals, Slater - type orbitals can better describe the cusp conditions and long - range asymptotic behavior of electron density near the nucleus. - **Nuclear cusp conditions**: Enforce Kato cusp conditions by Handy's method to ensure that the behavior of molecular orbitals near the nucleus is in line with physical expectations. ### Results: Through a series of test cases of atoms and molecules, such as the hydrogen molecule, lithium - hydrogen molecule, water molecule, and methylene (CH₂), the SlaterRKS method shows its performance in different basis sets. In particular, when dealing with strongly correlated systems, the SlaterRKS method can produce exchange - correlation potentials that are closer to the exact solution than traditional DFT methods. ### Formula summary: - Definition of exchange - correlation potential: \[ v_{XC}(\mathbf{r})=\frac{\delta E_{XC}[\rho(\mathbf{r})]}{\delta \rho(\mathbf{r})} \] - RKS working equation: \[ v_{XC}(\mathbf{r}) = v_{XC,Slater}^{\text{WF}}(\mathbf{r})+\frac{\tau^{\text{WF}}(\mathbf{r})}{\rho^{\text{WF}}(\mathbf{r})}-\frac{\tau^{\text{KS}}(\mathbf{r})}{\rho^{\text{KS}}(\mathbf{r})}+\epsilon^{\text{KS}}(\mathbf{r})-\epsilon^{\text{WF}}(\mathbf{r}) \] - Nuclear cusp condition: \[ \left.\frac{\partial \phi_i}{\partial r}\right|_{r = R_B}=-Z_B\phi_i(R_B) \]