The works included in this thesis aim to extend the applications of pair-coupled cluster doubles (pCCD) theory. pCCD produces reliable results for strongly correlated chemical systems and performs better than other traditional quantum chemical methods. Still, achieving greater accuracy requires post-pCCD corrections, which steeply increase computational cost, limiting its widespread application. Quantum embedding approaches, which reserve accurate but high-cost methods for the active site of chemical change and approximated methods for other regions of a chemical system, can be used to circumvent this problem. In this context, two new pCCD-based embedding formulations have been implemented and tested. Both the pCCD-in-DFT and pCCD-in-pCCD methods produce reliable results across a wide range of intermolecular interactions. In particular, the latter method shows significantly better performance, as it does not depend on approximated density functionals. Additionally, estimation of molecular properties is also pertinent for practical applications with pCCD. As the traditional response-theory-based formulation for pCCD analyzed here did not meet the required level of accuracy, an alternate approach based on the expectation-value method was implemented for pCCD and post-pCCD methods. This approach produces excellent agreement with the references at a much lower computational cost, as we bypass the need to solve the response equations. Thus, the pCCD-in-pCCD embedding method and the expectation-value-pCCD approach provide routes for large-scale simulations of challenging, complex systems.
The works included in this thesis aim to extend the applications of pair-coupled cluster doubles (pCCD) theory. pCCD produces reliable results for strongly correlated chemical systems and performs better than other traditional quantum chemical methods. Still, achieving greater accuracy requires post-pCCD corrections, which steeply increase computational cost, limiting its widespread application. Quantum embedding approaches, which reserve accurate but high-cost methods for the active site of chemical change and approximated methods for other regions of a chemical system, can be used to circumvent this problem. In this context, two new pCCD-based embedding formulations have been implemented and tested. Both the pCCD-in-DFT and pCCD-in-pCCD methods produce reliable results across a wide range of intermolecular interactions. In particular, the latter method shows significantly better performance, as it does not depend on approximated density functionals. Additionally, estimation of molecular properties is also pertinent for practical applications with pCCD. As the traditional response-theory-based formulation for pCCD analyzed here did not meet the required level of accuracy, an alternate approach based on the expectation-value method was implemented for pCCD and post-pCCD methods. This approach produces excellent agreement with the references at a much lower computational cost, as we bypass the need to solve the response equations. Thus, the pCCD-in-pCCD embedding method and the expectation-value-pCCD approach provide routes for large-scale simulations of challenging, complex systems.
The accurate and efficient description of strongly correlated systems remains an important challenge for computational methods. Doubly occupied configuration interaction (DOCI), in which all electrons are paired and no correlations which break these pairs are permitted, can in many cases provide an accurate account of strong correlations, albeit at combinatorial computational cost. Recently, there has been significant interest in a method we refer to as pair coupled cluster doubles (pCCD), a variant of coupled cluster doubles in which the electrons are paired. This is simply because pCCD provides energies nearly identical to those of DOCI, but at mean-field computational cost (disregarding the cost of the two-electron integral transformation). Here, we introduce the more complete pair extended coupled cluster doubles (pECCD) approach which, like pCCD, has mean-field cost and reproduces DOCI energetically. We show that unlike pCCD, pECCD also reproduces the DOCI wave function with high accuracy. Moreover, pECCD yields sensible albeit inexact results even for attractive interactions where pCCD breaks down.
Perhaps the most important approximations to the electronic structure problem in quantum chemistry are those based on coupled cluster and density functional theories. Coupled cluster theory has been called the ``gold standard'' of quantum chemistry due to the high accuracy that it achieves for weakly correlated systems. Kohn--Sham density functionals based on semilocal approximations are, without a doubt, the most widely used methods in chemistry and material science because of their high accuracy/cost ratio. The root of the success of coupled cluster and density functionals is their ability to efficiently describe the dynamic part of the electron correlation. However, both traditional coupled cluster and density functional approximations may fail catastrophically when substantial static correlation is present. This severely limits the applicability of these methods to a plethora of important chemical and physical problems such as, e.g., the description of bond breaking, transition states, transition metal-, lanthanide- and actinide-containing compounds, and superconductivity. In an attempt to tackle this problem, nonstandard (single-reference) coupled cluster-based techniques that aim to describe static correlation have been recently developed: pair coupled cluster doubles (pCCD) and singlet-paired coupled cluster doubles (CCD0). The ability to describe static correlation in pCCD and CCD0 comes, however, at the expense of important amounts of dynamic correlation so that the high accuracy of standard coupled cluster becomes unattainable. Thus, the reliable and efficient description of static and dynamic correlation in a simultaneous manner remains an open problem for quantum chemistry and many-body theory in general. In this thesis, different ways to combine pCCD and CCD0 with density functionals in order to describe static and dynamic correlation simultaneously (and efficiently) are explored. The combination of wavefunction and density functional methods has a long history in quantum chemistry (practical implementations have appeared in the literature since the 1970s). However, this kind of techniques have not achieved widespread use due to problems such as double counting of correlation and the symmetry dilemma---the fact that wavefunction methods respect the symmetries of Hamiltonian, while modern functionals are designed to work with broken symmetry densities. Here, particular mathematical features of pCCD and CCD0 are exploited to avoid these problems in an efficient manner. The two resulting families of approximations, denoted as pCCD+DFT and CCD0+DFT, are shown to be able to describe static and dynamic correlation in standard benchmark calculations. Furthermore, it is also shown that CCD0+DFT lends itself to combination with correlation from the direct random phase approximation (dRPA). Inclusion of dRPA in the long-range via the technique of range-separation allows for the description of dispersion correlation, the remaining part of the correlation. Thus, when combined with the dRPA, CCD0+DFT can account for all three-types of electron correlation that are necessary to accurately describe molecular systems. Lastly, applications of CCD0+DFT to actinide chemistry are considered in this work. The accuracy of CCD0+DFT for predicting equilibrium geometries and vibrational frequencies of actinide molecules and ions is assessed and compared to that of well-established quantum chemical methods. For this purpose, the f0 actinyl series (UO22+, NpO23+, PuO24+), the isoelectronic NUN, and Thorium (ThO, ThO2+) and Nobelium (NoO, NoO2) oxides are studied. It is shown that the CCD0+DFT description of these species agrees with available experimental data and is comparable with the results given by the highest-level calculations that are possible for such heavy compounds while being, at least, an order of magnitude lower in computational cost.
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 present a multilevel embedding framework for quantum chemistry calculations on a quantum computer. In our framework, a quantum algorithm treats the strongly correlated active space, while a high-level wave-function method such as coupled cluster theory or multireference perturbation theory recovers the remaining correlation in the surrounding region. A sampling-based quantum algorithm, quantum-selected configuration interaction, bridges the quantum and classical treatments. The entire calculation is embedded in a low-cost density functional theory description of the surrounding environment using Manby's projection technique. We apply the framework to organic, metal-organic, and metallic systems, computing bond dissociation energies, adsorption energies, and reaction barriers using only the subset of qubits of a 144-qubit superconducting quantum computer at the University of Osaka and achieving $\sim$1 kcal/mol agreement with classical references for a Menshutkin SN2 reaction inside a carbon nanotube. Our results may open the way to quantitatively reliable quantum-classical hybrid calculations for large-scale chemical systems.
Tailored coupled cluster theory represents a computationally inexpensive way to describe static and dynamical electron correlation effects. In this work, we scrutinize the performance of various coupled cluster methods tailored by electronic wave functions of polynomial cost. Specifically, we focus on frozen-pair coupled cluster (fpCC) methods, which are tailored by pair-coupled cluster doubles (pCCD), and coupled cluster theory tailored by matrix product state wave functions optimized by the density matrix renormalization group (DMRG) algorithm. As test system, we selected a set of various small- and medium-sized molecules containing diatomics (N 2, F 2, C 2, CN +, CO, BN, BO +, and Cr 2 ) and molecules (ammonia, ethylene, cyclobutadiene, benzene, hydrogen chains, rings, and cuboids) for which the conventional single-reference coupled cluster singles and doubles (CCSD) method is not able to produce accurate results for spectroscopic constants, potential energy surfaces, and barrier heights. Most importantly, DMRG-tailored and pCCD-tailored approaches yield similar errors in spectroscopic constants and potential energy surfaces compared to accurate theoretical and/or experimental reference data. Although fpCC methods provide a reliable description for the dissociation pathway of molecules featuring single and quadruple bonds, they fail in the description of triple or hextuple bond-breaking processes or avoided crossing regions.
Doubly occupied configuration interaction (DOCI) with optimized orbitals often accurately describes strong correlations while working in a Hilbert space much smaller than that needed for full configuration interaction. However, the scaling of such calculations remains combinatorial with system size. Pair coupled cluster doubles (pCCD) is very successful in reproducing DOCI energetically, but can do so with low polynomial scaling (N(3), disregarding the two-electron integral transformation from atomic to molecular orbitals). We show here several examples illustrating the success of pCCD in reproducing both the DOCI energy and wave function and show how this success frequently comes about. What DOCI and pCCD lack are an effective treatment of dynamic correlations, which we here add by including higher-seniority cluster amplitudes which are excluded from pCCD. This frozen pair coupled cluster approach is comparable in cost to traditional closed-shell coupled cluster methods with results that are competitive for weakly correlated systems and often superior for the description of strongly correlated systems.