Modern quantum chemical calculations commonly use wave function based models or densitry functionaltheory (DFT). Using Green’s functions is an approach that has not been explored much, and is one thathas only been explored relatively recently, with older examples being useful, but limited in scope. Theuse of Green’s functions in literature, historically and contemporary, is outlined and explored, with afocus on the self-consistent second order Green’s function (GF2), as well as the GW-approximation.
Modern quantum chemical calculations commonly use wave function based models or densitry functionaltheory (DFT). Using Green’s functions is an approach that has not been explored much, and is one thathas only been explored relatively recently, with older examples being useful, but limited in scope. Theuse of Green’s functions in literature, historically and contemporary, is outlined and explored, with afocus on the self-consistent second order Green’s function (GF2), as well as the GW-approximation.
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.
Calculations of ground-state and excited-state properties of materials have been one of the major goals of condensed matter physics. Ground-state properties of solids have been extensively investigated for several decades within the standard density functional theory. Excited-state properties, on the other hand, were relatively unexplored in ab initio calculations until a decade ago. The most suitable approach up to now for studying excited-state properties of extended systems is the Green function method. To calculate the Green function one requires the self-energy operator which is non-local and energy dependent. In this article we describe the GW approximation which has turned out to be a fruitful approximation to the self-energy. The Green function theory, numerical methods for carrying out the self-energy calculations, simplified schemes, and applications to various systems are described. Self-consistency issue and new developments beyond the GW approximation are also discussed as well as the success and shortcomings of the GW approximation.
The field of quantum chemistry is concerned with the modelling and simulation of the behaviour of molecular systems on the basis of the fundamental equations of quantum mechanics. Since these equations exhibit an extreme case of the curse of dimensionality (the Schrödinger equation for N electrons being a partial differential equation on \mathbb R^{3N} ), the quantum-chemical simulation of even moderate-size molecules already requires highly sophisticated model-reduction, approximation, and simulation techniques. The workshop brought together selected quantum chemists and physicists, and the growing community of mathematicians working in the area, to report and discuss recent advances on topics such as coupled-cluster theory, direct approximation schemes in full configuration-interaction (FCI) theory, interacting Green’s functions, foundations and computational aspects of density-functional theory (DFT), low-rank tensor methods, quantum chemistry in the presence of a strong magnetic field, and multiscale coupling of quantum simulations.
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.
While the use of Green's function techniques has a long tradition in quantum chemistry, the possibility of propagating the Kadanoff-Baym equations has remained largely unexplored. We have implemented the time-propagation for atoms and diatomic molecules, starting from a system in the groundstate. The initial stage of the calculation requires solving the Dyson equation self-consistently for the equilibrium Green's function. This Green's function contains a huge amount of information, and we have found it particularly interesting to compare the self-consistent total energies to the results of variational energy functionals of the Green's function. We also use time-propagation for calculating linear response functions, as a means for obtaining the excitation energies of the system. We have presently implemented the propagation for the second Born approximation, while the GW approximation has now been implemented for the ground state calculations.
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.