We establish the theoretical foundations for embedding a correlated wave function in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation which equates two, in principle distinct, kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-DFT formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based DMRG-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (PDFT) instead of a semilocal GGA does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.
We establish the theoretical foundations for embedding a correlated wave function in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation which equates two, in principle distinct, kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-DFT formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based DMRG-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (PDFT) instead of a semilocal GGA does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.
We establish theoretical foundations for embedding a correlated wavefunction in an environment formed by Kohn-Sham orbitals. We show that introducing an approximation that equates two, in principle, distinct kinetic-energy functionals yields an embedding functional identical to the projection-based wavefunction-in-density functional theory (DFT) formulation of Miller and co-workers. We demonstrate that this functional is inherently nonvariational: its minimum is not guaranteed to coincide with the exact ground-state energy and remains bounded from above by it. Building on this formal framework, we analyze the dominant sources of error in projection-based density matrix renormalization group-in-DFT embedding with approximate exchange-correlation (xc) functionals. Using molecules with dissociating covalent bonds as a diagnostic example, we demonstrate that the primary source of error is the nonadditive exchange-correlation energy describing the nonclassical coupling between the active subsystem and its environment. Eliminating the fractional-spin error by employing a pair-density xc functional (pair-density functional theory) instead of a semilocal generalized gradient approximation (GGA) does not remedy this deficiency, because the inaccuracy stems from self-interaction effects at the subsystem-environment interface.
The projection-based wave function in density functional theory (WF-in-DFT) embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [ Beran, P. J. Phys. Chem. Lett. 2023, 14, 716–722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H 20 chain and the cleavage of a triple CN bond in propionitrile.
The projection-based wave function (WF)-in-DFT embedding enables an efficient description of both the energetics and properties of large and complex chemical systems, with accuracy exceeding that of pure DFT. Recently, we have proposed using the density matrix renormalization group (DMRG) as the WF method for molecules containing strongly correlated fragments [Beran, P. et al. J. Phys. Chem. Lett. 2023, 14, 3, 716-722]. In this work, we demonstrate that the accuracy of the DMRG-in-DFT approach is primarily limited by the approximate treatment of the coupling between the active component and its environment through nonadditive exchange-correlation functionals. To address this issue, we combine exact exchange to reduce the nonadditive exchange error with a multireference adiabatic connection (AC) scheme to recover nonadditive correlation. The performance of the improved DMRG-in-DFT embedding is illustrated on two prototypical strongly correlated systems: the dissociation of the H20 chain and the cleavage of a triple CN bond in propionitrile.
Aiming to combine density functional theory (DFT) and wave-function theory, we study a mapping from the many-body interacting system to an effectively interacting Kohn-Sham system instead of a noninteracting Kohn-Sham system. Because a ground state of effectively interacting systems requires having a solution for the correlated many-body wave functions, this provides a natural framework to many-body wave-function theories such as the configuration interaction and the coupled-cluster method in the formal theoretical framework of DFT. Employing simple one-dimensional two-electron systems---namely, the one-dimensional helium atom, the hydrogen molecule, and the heteronuclear diatomic molecule---we investigate properties of many-body wave functions and exact exchange-correlation potentials of effectively interacting Kohn-Sham systems. As a result, we find that the asymptotic behavior of the exact exchange-correlation potential can be controlled by optimizing that of the effective interaction. Furthermore, the typical features of the exact noninteracting Kohn-Sham system, namely, a spiky feature and a step feature in the exchange-correlation potential for the molecular dissociation limit, can be suppressed by a proper choice of the effective interaction. These findings open a possibility to construct numerically robust and efficient exchange-correlation potentials and functionals based on the effectively interacting Kohn-Sham scheme.
Providing sufficient accuracy at low-computational cost, not surprisingly, density-functional theory ( DFT ) is a standard tool in theoretical material science. Similarly, quantum-electrodynamical density-functional theory ( QEDFT ) provides a promising route to describe experiments at the interface of quantum optics and quantum chemistry. Formally, in the Kohn-Sham ( KS ) formulation of DFT and QEDFT the exact many-particle Schrödinger equation is replaced by a set of coupled non-linear one-particle differential equations. The electron-electron, the electron-ion and the electron-photon interactions of the interacting system are mimicked by their corresponding effective KS potential with an, in principle unknown exchange-correlation ( xc ) term. In practice, DFT and QEDFT require approximations for this xc potential. The quality of which fully determines how much the DFT and QEDFT observables deviate from the exact ones. Understanding the shortcomings and improving the quality of approximate functionals makes the study of the exact functional of model systems indispensable. Here, we employ exact diagonalization for the Hamiltonian of different types of real-space lattice models in one spatial dimension to access the ground-state wave function, and thus to explicitly construct the exact functionals of DFT and QEDFT . To study the functional behavior of purely electronic systems, we consider two singlet-state electrons moving in an external potential with two connected atomic wells of variable depth and distance. To study coupled electron-photon systems, we consider the Rabi Hamiltonian with variable electron-photon coupling strength. By varying the model parameters, we systematically construct the functionals for systems from weak to strong coupling between the building blocks of the global system. Using the concept of the effective density-to- potential map, we discuss the intrasystem steepening which is an exact condition on the xc functional that captures the static correlation present in the system and carve out how it relates to the well-known intersystem derivative discontinuity. We further construct arbitrary observables as functionals of the internal variables which are directly affected by such features of the effective density-to-potential map. By constructing the reduced dipole map in the four-particle case of coupled electron-photon system, we show that the density-to-potential map helps to identify and understand recently discovered physical effects in the field of strongly coupled electron-photon systems. This work provides a contribution to the development of approximate function- als required in practical density-functional calculations of arbitrary experimentally accessible observables. Both the intrasystem steepening and the intersystem derivative discontinuity play a crucial role for systems with highly-localized electrons. Including both features in approximate functionals will improve the accu- racy of DFT evaluation of observables for such systems. Those results provide novel insight into the properties of the exact xc potentials in QEDFT . This work opens a new route to construct new functionals using coarse-grained or embed- ding approaches within the general and exact frameworks of DFT and QEDFT.
Standard approximations for the exchange-correlation (XC) functional in Kohn-Sham density functional theory (KS-DFT) typically lead to unacceptably large errors when applied to strongly-correlated electronic systems. Partition-DFT (PDFT) is a formally exact reformulation of KS-DFT in which the ground-state density and energy of a system are obtained through self-consistent calculations on isolated fragments, with a partition energy representing the \textit{inter}-fragment interactions. Here we show how typical errors of the local density approximation (LDA) in KS-DFT can be largely suppressed through a simple approximation, the generalized overlap approximation (GOA), for the partition energy in PDFT. Our method is illustrated on simple models of one-dimensional strongly-correlated linear hydrogen chains. The GOA, when used in combination with the LDA for the fragments, improves the LDA dissociation curves of hydrogen chains and produces results that are comparable to those of spin-unrestricted LDA, but without breaking the spin symmetry. GOA also induces a correction to the LDA electron density that partially captures the correct density dimerization in strongly-correlated hydrogen chains. Moreover, with an additional correction to the partition energy, the approximation is shown to produce dissociation energies in quantitative agreement to calculations based on the Density Matrix Renormalization Group method.