One of the primary challenges in quantum chemistry is the accurate modeling of strong electron correlation. While multireference methods effectively capture such correlation, their steep scaling with system size prohibits their application to large molecules and extended materials. Quantum embedding offers a promising solution by partitioning complex systems into manageable subsystems. In this Review, we highlight recent advances in multireference density matrix embedding and localized active space self-consistent field approaches for complex molecules and extended materials. We discuss both classical implementations and the emerging potential of these methods on quantum computers. By extending classical embedding concepts to the quantum landscape, these algorithms have the potential to expand the reach of multireference methods in quantum chemistry and materials.
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 major challenge in molecular electronic structure theory arises from chemical systems where the mean-field picture based upon a single electronic configuration fails to properly describe mechanics of many electrons. Such complex electronic states should be accounted for using a correlated superposition of multiple determinants and are referred to as multireference in quantum chemistry. Application domain of traditional multireference methods is severely limited to systems containing a small number of correlated electrons because of their exponential complexity. In recent years, much attention has been drawn to the density matrix renormalization group (DMRG) theory as a promising alternative. It has been introduced in ab initio quantum chemistry calculations and shown to be a highly-scalable multireference method that can handle full quantum degrees of freedom for the strongly-correlated wavefunction with unprecedented large configuration space. In this article, we provide an explicative introduction to the underlying theory and formalism built into the DMRG method and give a brief overview of our methodological extensions that combine DMRG and active-space correlation models.
Quantum chemistry models atoms and molecules using quantum mechanics, with the electronic structure problem at its core. This problem, which determines the energy and properties of electrons in a molecule, is exponentially hard to solve exactly on classical computers. To approximate solutions, effective models introduce layers of simplifications. However, these approximations can fail for strongly correlated systems, necessitating alternative approaches. Quantum computers, capable of processing quantum information natively, offer a promising path to simulate molecular systems beyond classical reach, with potential applications in catalysis, materials science, and photochemistry. Despite rapid advancements in quantum hardware, practical quantum computations remain challenging due to limited resources. This thesis optimizes quantum resource usage in computational chemistry, introducing novel embedding techniques and exploring applications beyond ground-state energies. Chapter 1 establishes the foundation by reviewing quantum chemistry methods, from Hartree-Fock and Density Functional Theory (DFT) to active space approaches. It also provides an introduction to quantum computing for chemistry, focusing on near-term algorithms. The following chapters present novel methods aimed at reducing quantum resource requirements and tackling problems classical methods struggle with. Chapter 2 introduces a technique to reduce the
λ-norm of the electronic structure Hamiltonian, a key measure of quantum resource demand. By transforming molecular orbitals, this method enhances sampling efficiency in Noisy Intermediate-Scale Quantum (NISQ) algorithms and reduces runtime in Fault-Tolerant Quantum Computing (FTQC). Localized orbitals outperform canonical ones, and an orbital-optimization scheme is proposed to minimize the
λ-norm. Benchmarks on hydrogen chains, FeMoco, and Ruthenium complexes demonstrate its effectiveness. Chapter 3 presents the FragPT2 framework for simulating large systems with multiple interacting fragments. This method embeds each fragment in the mean field of others and applies multi-reference perturbation theory to capture inter-fragment correlations such as dispersion and charge transfer efficiently. The approach reduces computational overhead and enables scalable simulations, as shown in applications to \ce{N2} dimers, aromatic dimers, and butadiene. Chapter 4 focuses on photochemical systems, where multiple electronic states must be treated equally. The state-averaged orbital-optimized variational quantum eigensolver (SAOOVQE) is extended with analytical energy gradients and non-adiabatic couplings, enabling geometry optimizations and simulations of non-adiabatic dynamics on quantum hardware. This is demonstrated by locating the conical intersection in formaldimine, highlighting SAOOVQE’s potential for studying photoisomerization in biomolecules. Building on this, Chapter 5 develops a quantum algorithm to detect conical intersections, where the Born-Oppenheimer approximation breaks down. These intersections are characterized by a discrete Berry phase (0 or
π), representing their topological nature. The algorithm extracts this phase as a single bit of information, ensuring resilience to noise. Analytical proof and numerical validation using a toy molecular model illustrate its effectiveness for studying light-perception biochemistry and related photochemical processes. Finally, Chapter 6 explores machine learning density functionals trained on noisy quantum data. While DFT is computationally efficient, its accuracy depends on approximated exchange-correlation functionals. This work proposes a neural network-based approach to learn the exact functional from quantum data, showing that training on noisy outputs from NISQ devices can mitigate sampling errors. Applications to the Fermi-Hubbard model and other systems demonstrate its potential to generalize and provide accurate results under noise. This thesis advances quantum computational chemistry by optimizing quantum resource usage, developing new embedding strategies, and expanding the scope of quantum simulations beyond ground-state properties. These contributions pave the way for practical quantum algorithms that could revolutionize molecular simulations in the near future.
Chemistry is considered as one of the more promising applications to science of near-term quantum computing. Recent work in transitioning classical algorithms to a quantum computer has led to great strides in improving quantum algorithms and illustrating their quantum advantage. Because of the limitations of near-term quantum computers, the most effective strategies split the work over classical and quantum computers. There is a proven set of methods in computational chemistry and materials physics that has used this same idea of splitting a complex physical system into parts that are treated at different levels of theory to obtain solutions for the complete physical system for which a brute force solution with a single method is not feasible. These methods are variously known as embedding, multi-scale, and fragment techniques and methods. We review these methods and then propose the embedding approach as a method for describing complex biochemical systems, with the parts not only treated with different levels of theory, but computed with hybrid classical and quantum algorithms. Such strategies are critical if one wants to expand the focus to biochemical molecules that contain active regions that cannot be properly explained with traditional algorithms on classical computers. While we do not solve this problem here, we provide an overview of where the field is going to enable such problems to be tackled in the future.
Abstract Quantum computing has shown great potential in various quantum chemical applications such as drug discovery, material design, and catalyst optimization. Although significant progress has been made in quantum simulation of simple molecules, ab initio simulation of solid-state materials on quantum computers is still in its early stage, mostly owing to the fact that the system size quickly becomes prohibitively large when approaching the thermodynamic limit. In this work, we introduce an orbital-based multi-fragment approach on top of the periodic density matrix embedding theory, resulting in a significantly smaller problem size for the current near-term quantum computer. We demonstrate the accuracy and efficiency of our method compared with the conventional methodologies and experiments on solid-state systems with complex electronic structures. These include spin polarized states of a hydrogen chain (1D-H), the equation of state of a boron nitride layer (h-BN) as well as the magnetic ordering in nickel oxide (NiO), a prototypical strongly correlated solid. Our results suggest that quantum embedding combined with a chemically intuitive fragmentation can greatly advance quantum simulation of realistic materials, thereby paving the way for solving important yet classically hard industrial problems on near-term quantum devices.