The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
The workhorse method of computational materials science is undeniably density functional theory in the Kohn-Sham framework of approximate exchange and correlation energy functionals. However, the need for highly accurate predictions of ground and excited state properties in materials science motivates the further development and exploration of alternative as well as complementary techniques. Among these alternative approaches, quantum chemical wavefunction based theories and in particular coupled cluster theory hold the promise to fill a gap in the toolbox of computational materials scientists. Coupled cluster (CC) theory provides a compelling framework of approximate infinite-order perturbation theory in the form of an exponential of cluster operators describing the true quantum many-body effects of the electronic wave function at a computational cost that, despite being significantly more expensive than DFT, scales polynomially with system size. The hierarchy of size-extensive approximate methods established in the framework of CC theory achieves systematic improvability for many materials properties. This is in contrast to currently available density functionals that often suffer from uncontrolled approximations that limit the accuracy in the prediction of materials properties. In this tutorial-style review we will introduce basic concepts of coupled cluster theory and recent developments that increase its computational efficiency for calculations of molecules, solids and materials in general. We will touch upon the connection between coupled cluster theory and the random-phase approximation to bridge the gap between traditional quantum chemistry and many-body Green’s function theories that are widely-used in the field of solid state physics. We will discuss various approaches to improve the computational performance without compromising accuracy. These approaches include large-scale parallel design as well as techniques that reduce the prefactor of the computational complexity. A central part of this article discusses the convergence of calculated properties to the thermodynamic limit which is of significant importance for reliable predictions of materials properties and constitutes an additional challenge compared to calculations of large molecules. Furthermore we mention technical aspects of computer code implementations of periodic coupled cluster theories in different numerical frameworks of the one-electron orbital basis; the projector-augmented-wave formalism using a plane wave basis set and the numeric atom-centered-orbital with resolution-of-identity.
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.
Many-body perturbation theories (MBPT) (both single and multireference), by design, are capable of representing electron correlation very well. However, the practical implementation of these methods for many-electron systems beyond a certain order is extremely difficult, if not impossible. For instance, there are hardly any calculations beyond third order at the multireference and fifth order at the single reference of perturbation theories. Thus, these methods are typically limited to systems where the higher-order electron correlation effect is not significant. A practical solution to this problem is the coupled cluster theory, which was first introduced by Coester [20] and Coester and Kümmel [21] in nuclear physics and was later adopted in the electron structure theory by Čížek [23, 106] and Paldus and Čížek [107]. The coupled-cluster (CC) theory is a non-perturbative all-order theory. This stateof-the-art theory has now been recognized as the basis of modern quantitative quantum chemistry and atomic physics for high-precision calculations. Its superior accuracy is rooted in the exponential ansatz for the wave function, which provides an energy expression that scales correctly with the number of electrons and under certain provisions leads to a faster convergence of the many-body expansion.
Connection between GW and Extended Coupled Cluster
Coupled-cluster (CC) theory and Green's function many-body perturbation theory (MBPT) have long evolved as distinct yet complementary frameworks for describing electronic correlation. While CC methods employ exponential wavefunction parametrizations that guarantee size extensivity and systematic improvability, Green's function approaches such as the GW approximation describe quasiparticle and optical excitations through diagrammatic resummations. Recent analyses have established a formal correspondence between these frameworks: the GW approximation is equivalent to an equation-of-motion (EOM) treatment of the direct-ring coupled-cluster doubles (drCCD) method. Within this context, the extended CC (ECC) ansatz offers a unified framework connecting CC and MBPT. This formulation bridges CC-based and Green's function-based methods, providing novel avenues for incorporating vertex corrections within a CC framework that preserves a sum-over-state representation of the self-energy and lead to potentially systematically improvable Green's function approaches.
Coupled-cluster (CC) theory and Green’s function many-body perturbation theory (MBPT) have long evolved as distinct yet complementary frameworks for describing electronic correlation. While CC methods employ exponential wave function parametrizations that guarantee size extensivity and systematic improvability, Green’s function approaches, such as the GW approximation, describe quasiparticle and optical excitations through diagrammatic resummations. Recent analyses have established a formal correspondence between these frameworks: the GW approximation is equivalent to an equation-of-motion (EOM) treatment of the direct-ring coupled-cluster doubles (drCCD) method. Within this context, the extended CC (ECC) ansatz offers a unified framework connecting CC and MBPT. This formulation bridges CC-based and Green’s function-based methods, providing novel avenues for incorporating vertex corrections within a CC framework that preserves a sum-over-state representation of the self-energy and leads to potentially systematically improvable Green’s function approaches.
Achieving materials-specific predictions for large, realistic molecules and solids with strong electron correlations remains a long-standing challenge in quantum chemistry. This dissertation develops scalable ab initio quantum many-body methods for crystalline materials, demonstrating that rigorous optimization of theoretical approximations and numerical algorithms enables predictive computation even in strongly correlated systems. The work tackles three interconnected challenges, namely describing material-specific properties of strongly correlated superconductors, achieving linear scaling with k-points for large-scale periodic calculations, and establishing reliable methods for electron-phonon interactions. The following chapters address each challenge in turn. Chapter 2 demonstrates that ab initio density matrix embedding theory, combined with symmetry-breaking strategies for superconducting order parameters, reproduces experimental trends in cuprate superconductivity, including pressure and layer effects on transition temperatures. By directly solving the electronic Schrödinger equation from material structures, these calculations capture multi-orbital covalency effects that simplified model Hamiltonians neglect, revealing their essential role in the spin fluctuations driving pairing. This predictive capability establishes a foundation for efficient computational screening of superconducting materials. Extending these calculations to larger k-point meshes and broader phase diagram explorations, however, demands more efficient computational algorithms. To overcome this computational bottleneck, Chapter 3 introduces interpolative separable density fitting framework. This development enables thermodynamic limit calculations with up to 1000 k-points, as validated on diamond, carbon dioxide, nickel monoxide, and cuprate systems. Integration with density matrix embedding theory and local natural orbital methods further yields converged ground-state energies from accurate correlated wavefunction methods such as CCSD(T) for crystalline materials in the thermodynamic limit. With computational efficiency established, the remaining challenge lies in extending the theoretical framework beyond purely electronic degrees of freedom. Chapter 4 addresses this need by establishing exponential ansatz wavefunctions as competitive approaches for electron-phonon systems. Unlike previous treatments based on perturbation theory, systematic benchmarking of coupled cluster theory and variational Lang-Firsov methods on the Holstein model demonstrates, for the first time, that these approaches accurately describe polaron formation across coupling regimes when combined with appropriate reference state optimization. These results provide a foundation for incorporating phonon effects into ab initio studies of superconducting materials. In addition to these methodological advances, GPU acceleration and MPI parallelization enhance computational throughput, enabling efficient large-scale calculations. Together, these developments constitute a computational framework for predictive ab initio quantum many-body calculations in crystalline materials, opening new avenues for understanding and designing strongly correlated quantum materials.
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.