Search papers, labs, and topics across Lattice.
This paper theoretically analyzes projection-based Density Matrix Renormalization Group in Density Functional Theory (DMRG-in-DFT) embedding methods, demonstrating that they are inherently non-variational, even with the exact exchange-correlation functional. The authors show that the method's functional minimum is not guaranteed to coincide with the exact ground-state energy. Through analysis of molecules with dissociating covalent bonds, they identify the nonadditive exchange-correlation energy describing the coupling between the active subsystem and its environment as the primary source of error, which cannot be remedied by simply eliminating fractional-spin error.
Projection-based DMRG-in-DFT, a popular quantum embedding method, is fundamentally flawed and cannot achieve exactness, even with perfect exchange-correlation functionals.
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.