Transition metal doping is commonly used for altering the properties of solid-state materials to suit applications in science and technology. Partially filled $d$-shells of transition metal atoms lead to electronic states with diverse spatial and spin symmetries. Chromium(III) cations have shown great potential for designing laser materials and, more recently, for developing spin qubits in quantum applications. They also represent an intriguing class of chemical systems with strongly correlated multi-reference excited states, due to the $d^3$ electron configuration. These states are difficult to describe accurately using single-reference quantum chemical methods such as density functional theory (DFT), the most commonly used method to study the electronic structures of solid-state systems. Recently, the periodic effective Hamiltonian of crystal field (pEHCF) method has been shown to overcome some limitations arising in the calculations of excited $d$-states. In this work, we assess the suitability of DFT and pEHCF to calculate the electronic structure and $d$-$d$ excitations of chromium(III) dopants in wide band gap host materials. The results will aid computational development of novel transition metal-doped materials and provide a deeper understanding of the complex nature of transition metal dopants in solids.
Transition metal doping is commonly used for altering the properties of solid-state materials to suit applications in science and technology. Partially filled d-shells of transition metal atoms lead to electronic states with diverse spatial and spin symmetries. Chromium(III) cations have shown great potential for designing laser materials and, more recently, for developing spin qubits in quantum applications. They also represent an intriguing class of chemical systems with strongly correlated multi-reference excited states, due to the d3 electron configuration. These states are difficult to describe accurately using single-reference quantum chemical methods such as density functional theory (DFT), the most commonly used method to study the electronic structures of solid-state systems. Recently, the periodic effective Hamiltonian of crystal field (pEHCF) method has been shown to overcome some limitations arising in the calculations of excited d-states. In this work, we assess the suitability of DFT and pEHCF to calculate the electronic structure and d-d excitations of chromium(III) dopants in wide band gap host materials. The results will aid computational development of novel transition metal-doped materials and provide a deeper understanding of the complex nature of transition metal dopants in solids.
Understanding emergent many-body phenomena in correlated materials remains one of the grandest challenges in solid-state quantum chemistry calculations. The individuality of electrons diminishes as correlation increases, which leads to the emergence of new features in the system. Seeking to address such problems accurately requires the aid of non-perturbative many-body techniques as well as the corresponding high-quality basis sets. In periodic systems, a given element may be present in different spatial arrangements displaying vastly different physical and chemical properties, and an elemental basis set independent of the physical properties of materials may lead to significant simulation inaccuracies. In fact, with the rapid progress of quantum chemistry methods in condensed-phase simulations, the need for a library of reliable Gaussian basis sets explicitly designed for periodic calculations has become urgent. In the meantime, the development of parameter-free, systematic and reliable quantum chemistry methods for simultaneously treating weak and strong correlations in periodic systems has never stopped. Their ability to treat realistic solid materials needs to be tested, and the obtained results can serve as benchmarks for empirical parameter-guided solid-state calculations. In this thesis, we present numerical studies of various solid-state systems, emphasizing the design and optimization of Gaussian basis sets and assessing the band structure properties of the periodic solid systems. The work begins by introducing the theoretical framework of the electronic structure theory and presenting the fundamentals of a fully self-consistent GW approximation and the self-energy embedding theory (SEET) in Chapters 2 and 3. It continues with a description of the Gaussian basis set optimization scheme we devised for solid-state quantum chemistry calculations. The scheme we present is designed to avoid a lack of material specificity within a given basis set by simultaneously minimizing the total energy of the system and optimizing the band energies when compared to the reference plane wave calculation while accounting for the overlap matrix condition number. We compare the quality of the Gaussian basis sets generated via our method against the existing basis sets. And the optimization scheme is tested to yield improved results. We finally present a quantitative study of electronic properties of realistic materials with the treatment of Green's function based weakly correlated self-consistent GW approximation and the strongly correlated self-energy embedding theory. In the transition metal oxide BiVO3, we found the inaccurate illustration of orbitals including t_{2g} in transition metal Vanadium by DFT. Our calculation shows a systematic trend of DFT failure in illustrating both insulating and metallic solutions, demonstrating the dangers of using DFT when the experimental data is scarce.
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.
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.
Accurate modeling of transition metal-containing compounds is of great interest due to their wide-ranging and significant applications. These systems present several challenges from an electronic structure perspective, including significant multi-reference characters and many chemically relevant orbitals. A further complication arises from the so-called double d-shell effect, which is known to cause a myriad of issues in the treatment of first-row transition metals with both single- and multi-reference methods. While this effect has been well documented for several decades, a comprehensive understanding of its consequences and underlying causes is still evolving. Here, we characterize the second d-shell effect by analyzing the information entropy of correlated wavefunctions in a periodic series of 3d and 4d transition metal molecular hydrides and oxides. These quantum information techniques provide unique insight into the nuanced electronic structure of these species and are powerful tools for the study of weak and strong correlations in the transition metal d manifold.
In the previous chapters DMFT has been developed as a methodology for calculation of dynamical and thermodynamic properties of many-body systems, with applications to model systems. In order for DMFT to be a general method for quantitative calculations, there must be systematic procedures to derive all the ingredients from first principles. In this chapter we consider the most prominent examples, materials with partially filled 3 d , 4 f , or 5 f localized, atomic-like states, and set up the characteristic form of the hamiltonian. This is the basis for DFT+DMFT methods discussed here and used in the examples in Ch. 20, and for the combination of DMFT with many-body perturbation theory in Ch. 21. Simpler techniques, including self-interaction correction, DFT+U, and DFT+Gutzwiller, are also useful and are summarized in the last sections of this chapter. Up to this point, dynamical mean-field theory has been developed for model systems with Hubbard-type hamiltonians, where the hopping matrix elements and interactions are considered as parameters. For problems with strong interactions, DMFT leads to renormalized bands, satellites in the one-particle spectra, magnetic phase transitions and local moments, metal–insulator transitions, and many other phenomena observed in real materials, as brought out in Ch. 2. The purpose of this chapter is to provide the framework for a quantitative theory of these materials. These are difficult problems, experimentally, theoretically, and computationally, and the first step is to identify the aspects of the materials that must be treated explicitly as the local strongly interacting system, and the “rest” that are treated by a different method. This division of the problem is useful because a real material involves many degrees of freedom, i.e., the electronic states must be described in a basis sufficiently large to describe all the relevant states. However, the part that can be treated by non-perturbative solvers (Ch. 18) is limited by the computational cost that scales rapidly with the number of states included in the embedded site or cluster. Thus one must choose a small number of localized orbitals that inevitably are non-unique. This may seem like a step backward for quantitative methods; however, it opens a window for methods to treat these interesting problems and it presents a challenge to overcome the limitations.
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.