Self-consistent quantum embedding forms a route to bridge local descriptions of electronic structure methods and the electronic structure of very large systems and condensed phases. Existing examples of quantum embedding methods include the density functional embedding of Cortona, Wesolowski and Warshel, based on the single-particle density and Carter, and the dynamical mean-field theory, based on the single-particle Greens function. In this talk, however, I will describe a new formalism for quantum embedding that is based on the single-particle density matrix, and which is based on a self-consistent approximation to the exact entanglement embedding of a subsystem in an environment. Unlike density functional embedding, this formalism can practically describe strongly correlated systems, and unlike dynamical mean-field theory, it is easily combined with traditional ground-state electronic structure methods, such as quantum chemical wavefunction methods and projector and variational Monte Carlo approaches. I will demonstrate the performance of the method on the phase diagram of the 1D and 2D Hubbard models and discuss applications to ab-initio systems.
Self-consistent quantum embedding forms a route to bridge local descriptions of electronic structure methods and the electronic structure of very large systems and condensed phases. Existing examples of quantum embedding methods include the density functional embedding of Cortona, Wesolowski and Warshel, based on the single-particle density and Carter, and the dynamical mean-field theory, based on the single-particle Greens function. In this talk, however, I will describe a new formalism for quantum embedding that is based on the single-particle density matrix, and which is based on a self-consistent approximation to the exact entanglement embedding of a subsystem in an environment. Unlike density functional embedding, this formalism can practically describe strongly correlated systems, and unlike dynamical mean-field theory, it is easily combined with traditional ground-state electronic structure methods, such as quantum chemical wavefunction methods and projector and variational Monte Carlo approaches. I will demonstrate the performance of the method on the phase diagram of the 1D and 2D Hubbard models and discuss applications to ab-initio systems.
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.
We present a multiscale approach to efficiently embed an ab initio correlated chemical fragment described by its energy-weighted density matrices and entangled with a wider mean-field many-electron system. This approach, first presented by Fertitta and Booth [Phys. Rev. B 98, 235132 (2018)], is here extended to account for realistic long-range interactions and broken symmetry states. The scheme allows for a systematically improvable description in the range of correlated fluctuations out of the fragment into the system, via a self-consistent optimization of a coupled auxiliary mean-field system. It is discussed that the method has rigorous limits equivalent to the existing quantum embedding approaches of both dynamical mean-field theory and density matrix embedding theory, to which this method is compared, and the importance of these correlated fluctuations is demonstrated. We derive a self-consistent local energy functional within the scheme and demonstrate the approach for hydrogen rings, where quantitative accuracy is achieved despite only a single atom being explicitly treated.
Conspectus Ab initio modeling of matter has become a pillar of chemical research: with ever-increasing computational power, simulations can be used to accurately predict, for example, chemical reaction rates, electronic and mechanical properties of materials, and dynamical properties of liquids. Many competing quantum mechanical methods have been developed over the years that vary in computational cost, accuracy, and scalability: density functional theory (DFT), the workhorse of solid-state electronic structure calculations, features a good compromise between accuracy and speed. However, approximate exchange-correlation functionals limit DFT's ability to treat certain phenomena or states of matter, such as charge-transfer processes or strongly correlated materials. Furthermore, conventional DFT is purely a ground-state theory: electronic excitations are beyond its scope. Excitations in molecules are routinely calculated using time-dependent DFT linear response; however applications to condensed matter are still limited. By contrast, many-electron wavefunction methods aim for a very accurate treatment of electronic exchange and correlation. Unfortunately, the associated computational cost renders treatment of more than a handful of heavy atoms challenging. On the other side of the accuracy spectrum, parametrized approaches like tight-binding can treat millions of atoms. In view of the different (dis-)advantages of each method, the simulation of complex systems seems to force a compromise: one is limited to the most accurate method that can still handle the problem size. For many interesting problems, however, compromise proves insufficient. A possible solution is to break up the system into manageable subsystems that may be treated by different computational methods. The interaction between subsystems may be handled by an embedding formalism. In this Account, we review embedded correlated wavefunction (CW) approaches and some applications. We first discuss our density functional embedding theory, which is formally exact. We show how to determine the embedding potential, which replaces the interaction between subsystems, at the DFT level. CW calculations are performed using a fixed embedding potential, that is, a non-self-consistent embedding scheme. We demonstrate this embedding theory for two challenging electron transfer phenomena: (1) initial oxidation of an aluminum surface and (2) hot-electron-mediated dissociation of hydrogen molecules on a gold surface. In both cases, the interaction between gas molecules and metal surfaces were treated by sophisticated CW techniques, with the remainder of the extended metal surface being treated by DFT. Our embedding approach overcomes the limitations of conventional Kohn-Sham DFT in describing charge transfer, multiconfigurational character, and excited states. From these embedding simulations, we gained important insights into fundamental processes that are crucial aspects of fuel cell catalysis (i.e., O2 reduction at metal surfaces) and plasmon-mediated photocatalysis by metal nanoparticles. Moreover, our findings agree very well with experimental observations, while offering new views into the chemistry. We finally discuss our recently formulated potential-functional embedding theory that provides a seamless, first-principles way to include back-action onto the environment from the embedded region.
We introduce extended density matrix embedding theory (EDMET), a static quantum embedding theory explicitly self-consistent with respect to local two-body physics. This overcomes the biggest practical and conceptual limitation of more traditional one-body embedding methods, namely the lack of screening and treatment of longer-range interactions. This algebraic zero-temperature embedding augments a local interacting cluster model with a minimal number of bosons from a description of the full system correlations via the random phase approximation, and admits an analytic approach to build a self-consistent Coulomb-exchange-correlation kernel. For extended Hubbard models with nonlocal interactions, this leads to the accurate description of phase transitions, static quantities, and dynamics. We also move towards ab initio systems via the Parriser-Parr-Pople model of conjugated coronene derivatives, finding good agreement with experimental optical gaps.
Strongly correlated electronic systems underpin a wide range of emerging quantum technologies, including solid-state defects for quantum information processing, molecular magnets, spin qubits, and quantum sensing platforms. Predictive computational design of these materials requires electronic structure methods capable of simultaneously describing strong electron correlation, excited-state phenomena, spin-dependent interactions, and environmental effects. While multireference wave function methods provide the necessary accuracy, their computational cost severely limits their application to realistic molecular and condensed-phase systems. This dissertation addresses this challenge through the development and application of quantum embedding and fragmentation methodologies that extend multireference electronic structure theory to increasingly complex quantum materials. The first part of this dissertation develops embedding-based approaches for describing optical properties of strongly correlated defects in extended solids. Complete Active Space Density Matrix Embedding Theory (CAS-DMET) and NEVPT2-DMET are applied to oxygen vacancy defects in magnesium oxide, where a systematic double-extrapolation protocol is introduced to approach both the non-embedding and thermodynamic limits. This framework enables quantitatively reliable calculations of absorption and emission energies while significantly reducing the computational cost of conventional multireference treatments of periodic systems. The second part extends embedding methodologies to the description of spin-dependent electronic structure. Density Matrix Embedding Pair-Density Functional Theory (DME-PDFT) is developed to recover dynamic electron correlation beyond the embedded active space, providing an efficient route for studying transition metal complexes and other strongly correlated molecular systems. Building upon these developments, CAS-DMET is further extended to investigate spin-phonon relaxation in molecular magnets and molecular crystals, demonstrating that quantum embedding can accurately describe spin dynamics and decoherence processes while substantially reducing computational cost. Beyond classical electronic structure calculations, this dissertation explores fragmentation strategies for hybrid quantum-classical algorithms through the development of efficient methodologies within the Localized Active Space Unitary Selective Coupled Cluster (LAS-USCC) framework. These developments improve the scalability of quantum chemistry algorithms for large active spaces and establish connections between embedding concepts and emerging quantum computing approaches. Finally, electronic structure calculations based on time-dependent density functional theory are employed to investigate fluorescent protein chromophores as potential biocompatible quantum sensors. These studies demonstrate how first-principles calculations can provide insight into spin-dependent properties in biologically relevant systems and expand the scope of computational approaches for quantum materials beyond traditional inorganic and molecular platforms. Collectively, the work presented in this dissertation establishes embedding and fragmentation as versatile frameworks for accurately describing the optical and spin properties of strongly correlated quantum systems across multiple length scales. By combining multireference electronic structure methods with physically motivated subsystem decompositions,these approaches bridge the gap between chemical accuracy and computational scalability, providing practical computational tools for the study and design of quantum materials relevant to quantum information science, molecular spintronics, and quantum sensing.
Density matrix embedding theory [G. Knizia and G. K.-L. Chan, Phys. Rev. Lett. 109, 186404 (2012)] and density embedding theory [I. W. Bulik, G. E. Scuseria, and J. Dukelsky, Phys. Rev. B 89, 035140 (2014)] have recently been introduced for model lattice Hamiltonians and molecular systems. In the present work, the formalism is extended to the ab initio description of infinite systems. An appropriate definition of the impurity Hamiltonian for such systems is presented and demonstrated in cases of 1, 2, and 3 dimensions, using coupled cluster theory as the impurity solver. Additionally, we discuss the challenges related to disentanglement of fragment and bath states. The current approach yields results comparable to coupled cluster calculations of infinite systems even when using a single unit cell as the fragment. The theory is formulated in the basis of Wannier functions but it does not require separate localization of unoccupied bands. The embedding scheme presented here is a promising way of employing highly accurate electronic structure methods for extended systems at a fraction of their original computational cost.