Description of electron correlation is crucial to chemical accuracy in quantum chemical calculations. However, the interaction of electrons in a system is an insoluble many-body problem. Constructing approximations to describe electron correlation is a challenging task, but there has been success in the quantum chemistry community using wavefunction methods and density functional theory. However, there are still challenges to be overcome in areas such as theoretical solid state chemistry, which requires description of large systems and the use of finite temperature. Large system sizes and finite temperature can be difficult to treat solely with the currently available methods. Therefore, a new class of methods, based on the temperature-dependent Green's function, is implemented and explored. This thesis is toward investigating the use of temperature-dependent Green's functions for ab-initio quantum chemistry. That is, we are treating a quantum chemical Hamiltonian with realistic electron interactions. While this formalism has been applied to model systems in the condensed matter community, its has been used much less by the quantum chemistry community. Therefore, the numerical behavior and accuracy of Green's function methods for quantum chemical calculations is relatively unknown. This work investigates the ability of the temperature-dependent Green's function, which is an ensemble formalism, to give access to temperature dependent thermodynamic quantities such as free energy and entropy when calculated in a second-order and perturbative manner (GF2). We find that this method is able to give good accuracy for lower temperatures and excellent accuracy for higher temperatures for a molecular case and is able to qualitatively describe a simple model of a solid. The results of this work are presented in chapter 3. Although Green's functions have a clear connection with spectra at zero temperature, it is not straight forward to obtain spectra from the finite temperature Green's function, which is calculated on the imaginary frequency axis. Therefore, we must investigate methods to obtain spectral quantities in a consistent and reliable manner from the imaginary axis. Chapter 4 investigates several methods to do so and we compare our results with experimental and highly accurate benchmark data. We find that it is possible to obtain spectral quantities that can differ by several electron volts, even if the same level of theory is used to obtain the Green's function. This reiterates that finding a spectrum from the imaginary axis is nontrivial and that one must exercise caution when comparing spectral quantities that were calculated using different techniques, even if they are treated with the same theoretical accuracy. Accessing larger systems with the Green's function requires the use of quantum embedding. Quantum embedding describes nontrivial electron interaction between a highly accurate ``active space'' or ``impurity'' and a larger, lower level ``environment''. It is challenging to construct an impurity solver that is reliable at low temperatures. In chapter 5 of this work, we implement and test a temperature dependent configuration interaction impurity solver for quantum embedding. This solver can be used in quantum embedding schemes such as dynamical mean field theory and self-energy embedding theory for larger systems. Overall, this thesis has made progress toward using Green's functions for ab-initio quantum chemistry at finite temperature. Groundwork has been laid for using this formalism to calculate thermodynamics and spectra using a Green's function with realistic electron interactions and to explore quantum embedding using an impurity solver at low temperatures.
Description of electron correlation is crucial to chemical accuracy in quantum chemical calculations. However, the interaction of electrons in a system is an insoluble many-body problem. Constructing approximations to describe electron correlation is a challenging task, but there has been success in the quantum chemistry community using wavefunction methods and density functional theory. However, there are still challenges to be overcome in areas such as theoretical solid state chemistry, which requires description of large systems and the use of finite temperature. Large system sizes and finite temperature can be difficult to treat solely with the currently available methods. Therefore, a new class of methods, based on the temperature-dependent Green's function, is implemented and explored. This thesis is toward investigating the use of temperature-dependent Green's functions for ab-initio quantum chemistry. That is, we are treating a quantum chemical Hamiltonian with realistic electron interactions. While this formalism has been applied to model systems in the condensed matter community, its has been used much less by the quantum chemistry community. Therefore, the numerical behavior and accuracy of Green's function methods for quantum chemical calculations is relatively unknown. This work investigates the ability of the temperature-dependent Green's function, which is an ensemble formalism, to give access to temperature dependent thermodynamic quantities such as free energy and entropy when calculated in a second-order and perturbative manner (GF2). We find that this method is able to give good accuracy for lower temperatures and excellent accuracy for higher temperatures for a molecular case and is able to qualitatively describe a simple model of a solid. The results of this work are presented in chapter 3. Although Green's functions have a clear connection with spectra at zero temperature, it is not straight forward to obtain spectra from the finite temperature Green's function, which is calculated on the imaginary frequency axis. Therefore, we must investigate methods to obtain spectral quantities in a consistent and reliable manner from the imaginary axis. Chapter 4 investigates several methods to do so and we compare our results with experimental and highly accurate benchmark data. We find that it is possible to obtain spectral quantities that can differ by several electron volts, even if the same level of theory is used to obtain the Green's function. This reiterates that finding a spectrum from the imaginary axis is nontrivial and that one must exercise caution when comparing spectral quantities that were calculated using different techniques, even if they are treated with the same theoretical accuracy. Accessing larger systems with the Green's function requires the use of quantum embedding. Quantum embedding describes nontrivial electron interaction between a highly accurate active space'' or impurity'' and a larger, lower level ``environment''. It is challenging to construct an impurity solver that is reliable at low temperatures. In chapter 5 of this work, we implement and test a temperature dependent configuration interaction impurity solver for quantum embedding. This solver can be used in quantum embedding schemes such as dynamical mean field theory and self-energy embedding theory for larger systems. Overall, this thesis has made progress toward using Green's functions for ab-initio quantum chemistry at finite temperature. Groundwork has been laid for using this formalism to calculate thermodynamics and spectra using a Green's function with realistic electron interactions and to explore quantum embedding using an impurity solver at low temperatures.
Electronic structure of materials from first principles is essential for one-to-one comparisons between theories and experiments. However, the many-body effects resulting from electron correlation and their interplay with other degrees of freedom pose a great challenge to design an ab initio approximate method with good accuracy and efficiency at the same time. In this thesis, we present our recent developments of ab initio many-body theory for solids based on the finite-temperature Green’s function formalism. We begin by reviewing the basic concepts and definitions of realistic electronic Hamiltonians (Ch. 2) as well as Green’s functions at finite temperature (Ch. 3). Overviews of self-consistent many-body perturbation theories and quantum embedding theories are then presented in Chs. 4 and 5 as complementary materials for the following chapters. In Ch. 6, we present algorithms and implementation details for the finite-temperature self-consistent GW (scGW) method in Gaussian-type orbitals for solids. Our implementation is based on the finite-temperature Green’s function formalism in which all equations are solved on the imaginary axis, without resorting to analytical continuation during the self-consistency. No quasiparticle approximation is employed and all matrix elements of self-energies are explicitly evaluated. We show agreement with other, differently formulated finite-temperature scGW implementations. By migrating computationally intensive calculations to GPUs, we obtain scalable results on large supercomputers with nearly optimal performance. Inclusion of relativistic effects in scGW is discussed in Ch. 7. Our formulation of relativistic scGW is based on the exact two-component formalism with one-electron approximation (X2C1e) and non-relativistic Coulomb interactions. Our theory allows us to study scalar relativistic effects, spin-orbit coupling, and their interplay with electron correlation without adjustable parameters. We examine the effect of the X2C1e by comparison to the established four-component formalism and reach excellent agreement. The simplicity of X2C1e enables construction of higher order theories, such as embedding theories, on top of perturbative calculations. In Ch. 8, we use self-energy embedding theory (SEET) to study the spectra of the prototypical undistorted cubic perovskites SrVO3 and SrMnO3. In the strongly correlated metallic SrVO3, we find that the usual attribution of the satellite peaks at -1.8 eV to Hund or Hubbard physics in the t2g orbitals is inconsistent with our calculations. In the strongly correlated insulator SrMnO3 we recover insulating behavior due to a feedback effect between the strongly correlated orbitals and the weakly correlated environment. Our calculation shows a systematic convergence of spectral features as the space of strongly correlated orbitals is enlarged, paving the way to a systematic parameter-free study of correlated perovskites. Lastly, in Ch. 9, we apply the Green’s function coupled cluster singles and doubles (GFCCSD) impurity solver to realistic impurity problems in strongly correlated solids within the framework of SEET. We describe the implementation details of our Green’s function coupled cluster (GFCC) solver, investigate its performance, and highlight potential advantages and problems on examples of impurities created for antiferromagnetic MnO and paramagnetic SrMnO3. GFCCSD provides satisfactory descriptions for weakly and moderately correlated impurities with sizes that are intractable by existing accurate impurity solvers such as exact diagonalization. However, when correlations become strong, the singles and doubles approximation used in GFCC could lead to instabilities in searching for the particle number present in impurity problems.
Recent progress in materials synthesis combined with advances in computational science have opened up an opportunity to design new materials with specific properties. Semiconductors, heterostructures, and thin layered films are used for energy harvesting and solar cells. Transition metal oxides and metallic surfaces with deposited magnetic molecules are used for signal conversion, nonvolatile memory modules and spintronics. Due to growing complexity of the new materials the most efficient approach to their design requires combined experimental and theoretical effort. In order to advance our understanding of the complicated physics of the modern materials new theoretical methods based on controlled, reliable, computationally efficient and systematically improvable approximations describing correlation and finite-temperature effects are needed. Density functional theory (DFT) remains a method of choice for materials science calculations due to its low computational cost. However, despite all recent efforts, the majority DFT applications are limited to weakly correlated systems and zero temperature. The lack of systematic improvability, inaccurate treatment of strongly correlated systems will continue to restrict the domain of DFT applications. The multi-reference wave-function methods can provide an accurate description of the strongly correlated systems. Unfortunately, these methods are computationally expensive. In spite of recent improvement the computational complexity will likely to hinder the applications of wave-function based methods to materials science in foreseeable future. Finite-temperature Green's function methods offer several advantages such as systematic improvability, access to the excitation spectrum, thermodynamic properties and straightforward implementation of embedding frameworks. It has been realized long time ago and has been successfully used in condensed matter physics to study model systems. The applications of such methods to realistic systems are impeded due to lack of theoretical approaches providing chemical accuracy and efficient computational algorithms. In this dissertation, these shortcomings are addressed by developing self-energy embedding theory and its extensive comparison to established quantum chemistry methods. New numerical algorithms for efficient application of the finite-temperature Green's function methods to molecular and extended systems were developed. Properties of the second-order Green's function theory (GF2) in the context of quantum chemistry were investigated. Additionally, a new range-separated hybrid functional combining GF2 and DFT was developed and benchmarked.
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.
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.
This dissertation concerns the quantum many-body problem, which is the problem of predicting the properties of systems of several quantum particles from the first principles of quantum mechanics. Included under this umbrella are various problems of fundamental importance in quantum chemistry, condensed matter physics, and materials science. Of particular interest is the electronic structure problem, the problem of determining the state of the electrons in a system with fixed atomic nuclei. Since direct numerical solution of the many-body Schrödinger equation is intractable even for systems of moderate size, a diverse array of approximate methods has been developed. The broad goals of this dissertation are to improve the mathematical understanding of certain widely-used approximations, as well as to propose new methods. Roughly speaking, we consider three (overlapping) categories of methods: Green's function methods, embedding methods, and variational methods.One can understand Green's function methods in terms of many-body perturbation theory, which computes series expansions of physical quantities about a non-interacting reference system. These expansions can be expressed graphically in terms of Feynman diagrams, which can in turn be reorganized, in some cases, into an expansion in terms of so-called bold diagrams. Green's function methods can be specified by choosing a subset of bold diagrams to approximate the sum. At the same time, such methods can be understood in terms of an object known as the Luttinger-Ward (LW) functional, which admits a representation in terms of the bold diagrams. Many aspects of these constructions are purely formal, and indeed the existence of the fermionic LW functional as a single-valued functional has recently been called into question. To contribute to the understanding of these issues, we provide rigorous proofs of the combinatorial construction and analytic interpretation of the bold diagrams in the simplified setting of a classical field theory. In this setting we also provide a rigorous non-perturbative construction of the LW functional via convex duality and prove several key properties, including continuity up to the boundary of its domain and asymptotics in the limit of large interaction.Quantum embedding methods, meanwhile, view a large system as being composed of smaller fragments that are treated with high accuracy and embedded in the larger system in a mutually consistent way. Inspired by a connection between the boundary analysis of the LW functional and embedding, we perform similar analysis for the 1-RDM theory for fermionic systems, which is also developed via convex duality, illustrating a relation to fermionic embedding methods such as the density matrix embedding theory (DMET).Another embedding method of note is the dynamical mean-field theory (DMFT), which is at the same time a Green's function method that can be understood in terms of the LW functional. DMFT relies on the solution of impurity problems, which specify the embedding of an interacting system into a non-interacting bath. Underlying DMFT is a result about the sparsity pattern of the self-energy matrix for impurity problems, which to our knowledge has not been proved in the literature. We provide a rigorous proof of this result in various classical and quantum settings. We go on to investigate the fermionic DMFT in depth, identifying the key mathematical structures that appear in the algorithmic loop for solving it and using these to prove the well-posedness of this loop, in a certain sense.Finally, we introduce a suite of new approaches to the quantum many-body problem that provide variational lower bounds to the ground-state energy. These methods, which combine the themes of convexity and embedding, are based on novel convex relaxations of the variational principles for the ground-state energies of many-body systems. To begin, we recover a second-quantized version of the formalism of strictly correlated electrons (SCE), which yields an exact expression for the exchange-correlation functional in Kohn-Sham density functional theory in the limit of infinite Coulomb repulsion in terms of the solution of a multi-marginal optimal transport problem. We introduce a semidefinite relaxation method for approximately solving this problem and obtaining a lower bound for the ground-state energy. The ideas underlying this relaxation are generalized considerably, outside the context of SCE, to yield much tighter lower bounds, which we validate numerically for both quantum spin systems and fermionic systems. We also describe how these relaxation methods can be interpreted as embedding methods via convex duality.
The description of electron correlation has been a critical problem among theoretical and computational chemistry researchers. To describe the physics of many new materials, this interaction is crucial. Typically, chemists and physicists have constructed approximations to classic electronic structure methods - wavefunction methods and density functional theory - to attempt to solve this problem. Both have had many successes in the computational chemistry field, but are hindered by either their computational cost or ability to rigorously describe the electron correlation. More recently, researchers have come back to investigating how Green's functions can be of use to study this correlation. This thesis focuses on the second-order Green's function which has shown successes in its moderate computational cost and ability to describe electron correlation. However, in solid-state chemistry where the systems to be studied are quite large, the method still must be implemented by other means. Instead of exploring this method deterministically, this thesis moves towards investigating the method via a stochastic method. While stochastic methods such as diagrammatic Monte Carlo are quite common in the physics community, they are far lesser explored among theoretical chemists. This work attempts to develop a computationally feasible way to implement diagrammatic Monte Carlo within the second-order Green's function scheme. This method is preferred to be implemented in a fully self-consistent matter. Thus, chapter 4 of this work will explore the algorithmics required to sample for the second-order self-energy component of the method. Further in chapter 5, the correct statistical exploration required to obtain important expectation values such as the one-body and two-body energies are discussed. It is found that due to the non-linearity of the data, non-parametric statistical resampling methods are required. In this work, it is shown that via a jackknife algorithm built into the second-order Green's function scheme, the stochastic error from a Monte Carlo evaluation of the self-energy can be controlled. To show the power of this analysis, in chapter 6 self-consistency is shown to be possible via calculations of a few model systems as well as larger more realistic chemical systems. Finally, chapter 7 of this thesis segues into a problem within the chemistry community. Coding skills are at the core of developing the methods described in this thesis; however, coding is not required in the chemistry curriculum in the United States. This chapter describes a curriculum that encouraged chemistry students to develop coding skills in a low stakes environment. While the overall hypothesis going into the study was that students who use coding to study a quantitative problem in chemistry will increase their understanding of that problem, the study left with further results. Via surveys, it was discovered that many chemistry students desire to learn coding and feel that developing the skill positively impacts their studies and future goals. Overall, this thesis has made great progress in developing a diagrammatic Monte Carlo technique that could be of interest to the chemistry community. The method has proven to be computationally feasibly and quantitatively correct in comparison to the analogous deterministic method. This work has formulated a scheme that can provide an accessible approach to solving other Green's function methods of interest and hopefully bring us closer to finding a method that is computationally approachable to quantitatively describing electronic correlation.