Quantum defect embedding theory (QDET) is a many-body embedding method designed to describe condensed systems with correlated electrons localized within a given region of space, for example spin defects in semiconductors and insulators. Although the QDET approach has been successful in predicting the electronic properties of several point defects, several limitations of the method remain. In this work, we propose multiple advances to the QDET formalism. We derive a double-counting correction that consistently treats the frequency dependence of the screened Coulomb interaction, and we illustrate the effect of including unoccupied orbitals in the active space. In addition, we propose a method to describe hybridization effects between the active space and the environment, and we compare the results of several impurity solvers, providing further insights into improving the reliability and applicability of the method. We present results for defects in diamond and for molecular qubits, including a detailed comparison with experiments.
Quantum defect embedding theory (QDET) is a many-body embedding method designed to describe condensed systems with correlated electrons localized within a given region of space, for example spin defects in semiconductors and insulators. Although the QDET approach has been successful in predicting the electronic properties of several point defects, several limitations of the method remain. In this work, we propose multiple advances to the QDET formalism. We derive a double-counting correction that consistently treats the frequency dependence of the screened Coulomb interaction, and we illustrate the effect of including unoccupied orbitals in the active space. In addition, we propose a method to describe hybridization effects between the active space and the environment, and we compare the results of several impurity solvers, providing further insights into improving the reliability and applicability of the method. We present results for defects in diamond and for molecular qubits, including a detailed comparison with experiments.
Quantum embedding theories are promising approaches to investigate strongly correlated electronic states of active regions of large-scale molecular or condensed systems. Notable examples are spin defects in semiconductors and insulators. We present a detailed derivation of a quantum embedding theory recently introduced, which is based on the definition of effective Hamiltonians. The effect of the environment on a chosen active space is accounted for through screened Coulomb interactions evaluated using density functional theory. Importantly, the random phase approximation is not required, and the evaluation of virtual electronic orbitals is circumvented with algorithms previously developed in the context of calculations based on many-body perturbation theory. In addition, we generalize the quantum embedding theory to active spaces composed of orbitals that are not eigenstates of Kohn-Sham Hamiltonians. Finally, we report results for spin defects in semiconductors.
A quantitative description of the excited electronic states of point defects and impurities is crucial for understanding materials properties, and possible applications of defects in quantum technologies. This is a considerable challenge for computational methods, since Kohn-Sham density functional theory (DFT) is inherently a ground-state theory, while higher-level methods are often too computationally expensive for defect systems. Recently, embedding approaches have been applied that treat defect states with many-body methods, while using DFT to describe the bulk host material. We implement such an embedding method, based on Wannierization of defect orbitals and the constrained random-phase approximation approach, and perform systematic characterization of the method for three distinct systems with current technological relevance: a carbon dimer replacing a B and N pair in bulk hexagonal BN (CBCN), the negatively charged nitrogen-vacancy center in diamond (${\mathrm{NV}}^{\ensuremath{-}}$), and an Fe impurity on the Al site in wurtzite AlN (FeAl). In the context of these test-case defects, we demonstrate that crucial considerations of the methodology include convergence of the bulk screening of the active-space Coulomb interaction, the choice of exchange-correlation functional for the initial DFT calculation, and the treatment of the ``double-counting'' correction. For CBCN we show that the embedding approach gives many-body states in agreement with analytical results on the Hubbard dimer model, which allows us to elucidate the effects of the DFT functional and double-counting correction. For the ${\mathrm{NV}}^{\ensuremath{-}}$ center, our method demonstrates good quantitative agreement with experiments for the zero-phonon line of the triplet-triplet transition. Finally, we illustrate challenges associated with this method for determining the energies and orderings of the complex spin multiplets in FeAl.
Point defects are of interest for many applications, from quantum sensing to modifying bulk properties of materials. Because of their localized orbitals, the electronic states are often strongly correlated, which has led to a proliferation of quantum embedding techniques to treat this correlation. In these techniques, most of the one-body states are treated with a weakly correlated theory such as density functional theory, and certain one-body states are singled out as an active space to be treated using an effective interaction. We assess these techniques for iron and chromium defects in aluminum nitride using quantum Monte Carlo (QMC) calculations on identical Hamiltonians. For these systems, we find the dominant errors in the embedding arise from the one-body crystal-field splitting in the d orbitals inherited from density functional theory (DFT), rather than double counting corrections, with the screened interactions also affected by the DFT orbitals. Strikingly, the best double counting recipe is opposite in these two cases. Because excitation energies can agree while the underlying wave functions do not, diagnosing these errors requires detailed information about the many-body wave functions, which QMC provides.
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.
We investigate the electronic properties of an exemplar transition metal impurity in an insulator, with the goal of accurately describing strongly correlated defect states. We consider iron in aluminum nitride, a material of interest for hybrid quantum technologies, and we carry out calculations with quantum embedding methods, density matrix embedding theory (DMET) and quantum defect embedding theory (QDET), and with spin-flip time-dependent density functional theory (TDDFT). We show that both DMET and QDET accurately describe the ground state and low-lying excited states of the defect and that TDDFT yields photoluminescence spectra in agreement with experiments. In addition, we provide a detailed discussion of the convergence of our results as a function of the active space used in the embedding methods, thus defining a protocol to obtain converged data directly comparable with experiments.
The self-energy embedding theory (SEET), in which the active space self-energy is embedded in the self-energy obtained from a perturbative method treating the non-local correlation effects, was recently developed in our group. In SEET, the double counting problem does not appear and the accuracy can be improved either by increasing the perturbation order or by enlarging the active space. This method was first calibrated for the 2D Hubbard lattice showing promising results. In this paper, we report an extension of SEET to quantum chemical ab initio Hamiltonians for applications to molecular systems. The self-consistent second-order Green's function method is used to describe the non-local correlations, while the full configuration interaction method is carried out to capture strong correlation within the active space. Using few proof-of-concept examples, we show that SEET yields results of comparable quality to n-electron valence state second-order perturbation theory with the same active space, and furthermore, the full active space can be split into smaller active spaces without further implementation. Moreover, SEET avoids intruder states and does not require any high-order reduced density matrices. These advantages show that SEET is a promising method to describe physical and chemical properties of challenging molecules requiring large active spaces.