Local density approximation in site-occupation embedding theory

Extended N-centered ensemble density functional theory of double electronic excitations

A recent work (arXiv:2401.04685) has merged N-centered ensembles of neutral and charged electronic ground states with ensembles of neutral ground and excited states, thus providing a general and in-principle exact (so-called extended N-centered) ensemble density functional theory of neutral and charged electronic excitations. This formalism made it possible to revisit the concept of density-functional derivative discontinuity, in the particular case of single excitations from the highest occupied Kohn-Sham (KS) molecular orbital, without invoking the usual "asymptotic behavior of the density" argument. In this work, we address a broader class of excitations, with a particular focus on double excitations. An exact implementation of the theory is presented for the two-electron Hubbard dimer model. A thorough comparison of the true physical ground- and excited-state electronic structures with that of the fictitious ensemble density-functional KS system is also presented. Depending on the choice of the density-functional ensemble as well as the asymmetry of the dimer and the correlation strength, an inversion of states can be observed. In some other cases, the strong mixture of KS states within the true physical system makes the assignment "single excitation" or "double excitation" irrelevant.

Exact Excited-State Functionals of the Asymmetric Hubbard Dimer

The exact functionals associated with the (singlet) ground state and the two singlet excited states of the asymmetric Hubbard dimer at half-filling are calculated using both Levy's constrained search and Lieb's convex formulation. While the ground-state functional is, as is commonly known, a convex function with respect to the density, the functional associated with the doubly excited state is found to be concave. Also, because the density-potential mapping associated with the first excited state is noninvertible, its "functional" is a partial, multivalued function composed of one concave and one convex branch that correspond to two separate domains of the external potential. Remarkably, it is found that, although the one-to-one mapping between density and external potential may not apply (as in the case of the first excited state), each state-specific energy and corresponding universal functional are "functions" whose derivatives are each other's inverse, just as in the ground state formalism.

Comparing correlation components and approximations in Hartree-Fock and Kohn-Sham theories via an analytical test case study

The asymmetric Hubbard dimer is a model that allows for explicit expressions of the Hartree-Fock (HF) and Kohn-Sham (KS) states as analytical functions of the external potential, Δv, and of the interaction strength, U. We use this unique circumstance to establish a rigorous comparison between the individual contributions to the correlation energies stemming from the two theories in the {U, Δv} parameter space. Within this analysis of the Hubbard dimer, we observe a change in the sign of the HF kinetic correlation energy, compare the indirect repulsion energies, and derive an expression for the "traditional" correlation energy, i.e., the one that corrects the HF estimate, in a pure site-occupation function theory spirit [Eq. (45)]. Next, we test the performances of the Liu-Burke and the Seidl-Perdew-Levy functionals, which model the correlation energy based on its weak- and strong-interaction limit expansions and can be used for both the traditional and the KS correlation energies. Our results show that, in the Hubbard dimer setting, they typically work better for the HF reference, despite having been originally devised for KS. These conclusions are somewhat in line with prior assessments of these functionals on various chemical datasets. However, the Hubbard dimer model allows us to show the extent of the error that may occur in using the strong-interaction ingredient for the KS reference in place of the one for the HF reference, as has been carried out in most of the prior assessments.

Local Potential Functional Embedding Theory: A Self-Consistent Flavor of Density Functional Theory for Lattices without Density Functionals

Quantum embedding is a divide and conquer strategy that aims at solving the electronic Schrödinger equation of sizeable molecules or extended systems. We establish in the present work a clearer and in-principle-exact connection between density matrix embedding theory (DMET) and density-functional theory (DFT) within the simple but nontrivial one-dimensional Hubbard model. For that purpose, we use our recent reformulation of single-impurity DMET as a Householder transformed density-matrix functional embedding theory (Ht-DMFET). On the basis of well-identified density-functional approximations, a self-consistent local potential functional embedding theory (LPFET) is formulated and implemented. Combining both LPFET and DMET numerical results with our formally exact density-functional embedding theory reveals that a single statically embedded impurity can in principle describe the density-driven Mott–Hubbard transition, provided that a complementary density-functional correlation potential (which is neglected in both DMET and LPFET) exhibits a derivative discontinuity (DD) at half filling. The extension of LPFET to multiple impurities (which would enable to circumvent the modeling of DDs) and its generalization to quantum chemical Hamiltonians are left for future work.

Ensemble Density Functional Theory of Neutral and Charged Excitations : Exact Formulations, Standard Approximations, and Open Questions

Recent progress in the field of (time-independent) ensemble density-functional theory (DFT) for excited states are reviewed. Both Gross-Oliveira-Kohn (GOK) and N-centered ensemble formalisms, which are mathematically very similar and allow for an in-principle-exact description of neutral and charged electronic excitations, respectively, are discussed. Key exact results, for example, the equivalence between the infamous derivative discontinuity problem and the description of weight dependencies in the ensemble exchange-correlation density functional, are highlighted. The variational evaluation of orbital-dependent ensemble Hartree-exchange (Hx) energies is discussed in detail. We show in passing that state-averaging individual exact Hx energies can lead to severe (although solvable) v-representability issues. Finally, we explore the possibility of using the concept of density-driven correlation, which has been introduced recently and does not exist in regular ground-state DFT, for improving state-of-the-art correlation density-functional approximations for ensembles. The present review reflects the efforts of a growing community to turn ensemble DFT into a rigorous and reliable low-cost computational method for excited states. We hope that, in the near future, this contribution will stimulate new formal and practical developments in the field.

Embedding Methods for Quantum Chemistry: Applications from Materials to Life Sciences

Quantum mechanical embedding methods hold the promise to transform not just the way calculations are performed, but to significantly reduce computational costs and improve scaling for macro-molecular systems containing hundreds if not thousands of atoms. The field of embedding has grown increasingly broad with many approaches of different intersecting flavors. In this perspective, we lay out the methods into two streams: QM:MM and QM:QM, showcasing the advantages and disadvantages of both. We provide a review of the literature, the underpinning theories including our contributions, and we highlight current applications with select examples spanning both materials and life sciences. We conclude with prospects and future outlook on embedding, and our view on the use of universal test case scenarios for cross-comparisons of the many available (and future) embedding theories.

Electron transport through a spin crossover junction. Perspectives from a wavefunction-based approach

The present paper reports the application of a computational framework, based on the quantum master equation, the Fermi's golden Rule, and conventional wavefunction-based methods, to describe electron transport through a spin crossover molecular junction (Fe(bapbpy) (NCS)₂, 1, bapbpy = N-(6-(6-(Pyridin-2-ylamino)pyridin-2-yl)pyridin-2-yl)-pyridin-2-amine). This scheme is an alternative to the standard approaches based on the relative position and nature of the frontier orbitals, as it evaluates the junction's Green's function by means of accurate state energies and wavefunctions. In the present work, those elements are calculated for the relevant states of the high- and low-spin species of 1, and they are used to evaluate the output conductance within a given range of bias- and gate-voltages. The contribution of the ground and low-lying excited states to the current is analyzed, and inspected in terms of their 2S + 1 Ms-states. In doing so, it is shown the relevance of treating not only the ground state in its maximum-Ms projection, as usually done in most computational-chemistry packages, but the whole spectrum of low-energy states of the molecule. Such improved representation of the junction has a notable impact on the total conductivity and, more importantly, it restores the equivalence between alpha and beta transport, which means that no spin polarization is observed in the absence of Zeeman splitting. Finally, this work inspects the strong- and weak-points of the suggested theoretical framework to understand electron transport through molecular switchable materials, identifies a pathway for future improvement, and offers a new insight into concepts that play a key role in spintronics.

Quantum Embedding Theories

In complex systems, it is often the case that the region of interest forms only one part of a much larger system. The idea of joining two different quantum simulations-a high level calculation on the active region of interest, and a low level calculation on its environment-formally defines a quantum embedding. While any combination of techniques constitutes an embedding, several rigorous formalisms have emerged that provide for exact feedback between the embedded system and its environment. These three formulations: density functional embedding, Green's function embedding, and density matrix embedding, respectively, use the single-particle density, single-particle Green's function, and single-particle density matrix as the quantum variables of interest. Many excellent reviews exist covering these methods individually. However, a unified presentation of the different formalisms is so far lacking. Indeed, the various languages commonly used, functional equations for density functional embedding, diagrammatics for Green's function embedding, and entanglement arguments for density matrix embedding, make the three formulations appear vastly different. In this Account, we introduce the basic equations of all three formulations in such a way as to highlight their many common intellectual strands. While we focus primarily on a straightforward theoretical perspective, we also give a brief overview of recent applications and possible future developments. The first section starts with density functional embedding, where we introduce the key embedding potential via the Euler equation. We then discuss recent work concerning the treatment of the nonadditive kinetic potential, before describing mean-field density functional embedding and wave function in density functional embedding. We finish the section with extensions to time-dependence and excited states. The second section is devoted to Green's function embedding. Here, we use the Dyson equation to obtain equations that parallel as closely as possible the density functional embedding equations, with the hybridization playing the role of the embedding potential. Embedding a high-level self-energy within a low-level self-energy is treated analogously to wave function in density functional embedding. The numerical computation of the high-level self-energy allows us to briefly introduce the bath representation in the quantum impurity problem. We then consider translationally invariant systems to bring in the important dynamical mean-field theory. Recent developments to incorporate screening and long-range interactions are discussed. The third section concerns density matrix embedding. Here, we first highlight some mathematical complications associated with a simple Euler equation derivation, arising from the open nature of fragments. This motivates the density matrix embedding theory, where we use the Schmidt decomposition to represent the entanglement through bath orbitals. The resulting impurity plus bath formulation resembles that of dynamical mean-field theory. We discuss the numerical self-consistency associated with using a high-level correlated wave function with a mean-field low-level treatment, and connect the resulting numerical inversion to that used in density functional embedding. We finish with perspectives on the future of all three methods.

Polarizable Embedding Density Matrix Renormalization Group

The polarizable embedding (PE) approach is a flexible embedding model where a preselected region out of a larger system is described quantum mechanically, while the interaction with the surrounding environment is modeled through an effective operator. This effective operator represents the environment by atom-centered multipoles and polarizabilities derived from quantum mechanical calculations on (fragments of) the environment. Thereby, the polarization of the environment is explicitly accounted for. Here, we present the coupling of the PE approach with the density matrix renormalization group (DMRG). This PE-DMRG method is particularly suitable for embedded subsystems that feature a dense manifold of frontier orbitals which requires large active spaces. Recovering such static electron-correlation effects in multiconfigurational electronic structure problems, while accounting for both electrostatics and polarization of a surrounding environment, allows us to describe strongly correlated electronic structures in complex molecular environments. We investigate various embedding potentials for the well-studied first excited state of water with active spaces that correspond to a full configuration-interaction treatment. Moreover, we study the environment effect on the first excited state of a retinylidene Schiff base within a channelrhodopsin protein. For this system, we also investigate the effect of dynamical correlation included through short-range density functional theory.