Approximate Functionals for Multistate Density Functional Theory

Alexander Humeniuk
DOI: https://doi.org/10.26434/chemrxiv-2024-xb9jr
2024-03-15
Abstract:Recently Lu and Gao [J. Phys. Chem. Lett. 13, 33, 7762-7769 (2022)] published a new, rigorous density functional theory for excited states and proved that the projection of the kinetic and electron repulsion operators into the subspace of the lowest electronic states are a universal functional of the matrix density D(r). This is the first attempt to find an approximation to the multistate universal functional F[D(r)]. It is shown that F (i) does not explicitly depend on the number of electronic states and (ii) is an analytic matrix functional. The Thomas-Fermi-Dirac-von Weizsäcker model and the correlation energy of the homogeneous electron gas are turned into matrix functionals guided by two principles: That each matrix functional should transform properly under basis set transformations, and that the ground state functional should be recovered for a single electronic state. Lieb-Oxford-like bounds on the average kinetic and electron repulsion energies in the subspace are given. When evaluated on the numerically exact matrix density of LiF, this simple approximation reproduces the matrix elements of the Hamiltonian in the basis of the exact eigenstates accurately for all bond lengths. In particular the off-diagonal elements of the effective Hamiltonian that come from the interactions of different electronic states can be calculated with the same or better accuracy than the diagonal elements. Unsurprisingly, the largest error comes from the kinetic energy functional. More exact conditions that constrain the functional form of F are needed to go beyond the local density approximation.
Chemistry
What problem does this paper attempt to address?
This paper discusses a key issue in Multistate Density Functional Theory (MSDFT) - how to handle static correlation with DFT in the presence of electronic state interactions. Existing methods often struggle to accurately calculate the coupling terms between different electronic states. The authors propose a new approach that applies the form of the ground state functional to matrix density in order to approximate the expression of the multistate functional, especially considering the analytical properties of the matrix density. This method can be used to calculate the diagonal and off-diagonal elements of the effective Hamiltonian, and it performs well in numerical tests on LiF molecules. The study also points out the need for more accurate conditions to constrain the functional form beyond local density approximation.