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 work 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 work 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.