The accurate ab initio simulation of molecules and periodic solids with diagrammatic perturbation theory is an important task in quantum chemistry, condensed matter physics, and materials science. In this article, we present the WeakCoupling module of the open-source software package Green, which implements fully self-consistent diagrammatic weak coupling simulations, capable of dealing with real materials in the finite-temperature formalism. The code is licensed under the permissive MIT license. We provide self-consistent GW (scGW) and self-consistent second-order Green's function perturbation theory (GF2) solvers, analysis tools, and post-processing methods. This paper summarizes the theoretical methods implemented and provides background, tutorials and practical instructions for running simulations.
The accurate ab initio simulation of molecules and periodic solids with diagrammatic perturbation theory is an important task in quantum chemistry, condensed matter physics, and materials science. In this article, we present the WeakCoupling module of the open-source software package Green, which implements fully self-consistent diagrammatic weak coupling simulations, capable of dealing with real materials in the finite-temperature formalism. The code is licensed under the permissive MIT license. We provide self-consistent GW (scGW) and self-consistent second-order Green's function perturbation theory (GF2) solvers, analysis tools, and post-processing methods. This paper summarizes the theoretical methods implemented and provides background, tutorials and practical instructions for running simulations.
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.
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.
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.
We present an implementation of a fully self-consistent finite temperature second order Green's function perturbation theory (GF2) within the diagrammatic Monte Carlo framework. In contrast to the previous implementations of stochastic GF2 ({\it J. Chem. Phys.},{\bf 151}, 044144 (2019)), the current self-consistent stochastic GF2 does not introduce a systematic bias of the resulting electronic energies. Instead, the introduced implementation accounts for the stochastic errors appearing during the solution of the Dyson equation. We present an extensive discussion of the error handling necessary in a self-consistent procedure resulting in dressed Green's function lines. We test our method on a series of simple molecular examples.
Despite recent advances, systematic quantitative treatment of the electron correlation problem in extended systems remains a formidable task. Systematically improvable Green's function methods capable of quantitatively describing weak and at least qualitatively strong correlations appear as promising candidates for computational treatment of periodic systems. We present a periodic implementation of temperature-dependent self-consistent 2nd-order Green's function (GF2) method, where the self-energy is evaluated in the basis of atomic orbitals. Evaluating the real-space self-energy in atomic orbitals and solving the Dyson equation in k-space are the key components of a computationally feasible algorithm. We apply this technique to the one-dimensional hydrogen lattice--a prototypical crystalline system with a realistic Hamiltonian. By analyzing the behavior of the spectral functions, natural occupations, and self-energies, we claim that GF2 is able to recover metallic, band insulating, and at least qualitatively Mott regimes. We observe that the iterative nature of GF2 is essential to the emergence of the metallic and Mott phases.