Methods based on the many-body Green's function are generally accepted as the path forward beyond Kohn-Sham based density functional theory, in order to compute from first principles electronic structure of materials with strong correlations and excited-state properties in nano- and materials science. Here we present an efficient method to compute the screened Coulomb interactionW, the crucial and computationally most demanding ingredient in the GW method, within the framework of the all-electron Linearized Augmented Plane Wave method. We use the method to compute from first principles, within the constrained random phase approximation (c-RPA), the frequency-dependent screened Hubbard U-matrix defined for a Wannier basis in which we downfold the many-body Hamiltonian for La2CuO4, the canonical parent compound of several cuprate high-temperature superconductors. These results were computed at scale on the Cray XT5 at ORNL, sustaining 1.30 petaflop. We discuss the details of the algorithm and its implementation that allowed us to reach high efficiency and short time to solution on today's petaflop supercomputers.
Methods based on the many-body Green's function are generally accepted as the path forward beyond Kohn-Sham based density functional theory, in order to compute from first principles electronic structure of materials with strong correlations and excited-state properties in nano- and materials science. Here we present an efficient method to compute the screened Coulomb interactionW, the crucial and computationally most demanding ingredient in the GW method, within the framework of the all-electron Linearized Augmented Plane Wave method. We use the method to compute from first principles, within the constrained random phase approximation (c-RPA), the frequency-dependent screened Hubbard U-matrix defined for a Wannier basis in which we downfold the many-body Hamiltonian for La2CuO4, the canonical parent compound of several cuprate high-temperature superconductors. These results were computed at scale on the Cray XT5 at ORNL, sustaining 1.30 petaflop. We discuss the details of the algorithm and its implementation that allowed us to reach high efficiency and short time to solution on today's petaflop supercomputers.
In this lecture we present many-body perturbation theory as a method to determine quasiparticle excitations in solids, especially electronic band structures, accurately from first principles.The main ingredient is the electronic self-energy that, in principle, contains all many-body exchange and correlation effects beyond the Hartree potential.As its exact mathematical expression is unknown, approximations must be used in practical calculations.The approximation is obtained using a systematic algebraic approach on the basis of Green function techniques.It constitutes an expansion of the self-energy up to linear order in the screened Coulomb potential, which describes the interaction between the quasiparticles and includes dynamic screening through the creation of exchange-correlation holes around the bare particles.The implementation of the approximation relies on a perturbative treatment starting from density functional theory.Besides a detailed mathematical discussion we focus on the underlying physical concepts and show some illustrative applications.
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.
The well-known and well-used Density Functional Theory (DFT) is a variational theory for the ground-state. There are no guarantees that it delivers good results for excited state properties such as quasiparticle energies or wave function. All that is guaranteed is that the total energy is exact, the chemical potential is correct, and that the eective electronic states (Kohn-Sham states) generate the right electron density at the variational minimum. Luttinger-Ward approaches to electronic structure are attractive as beyond-DFT schemes since, in principle, they oer a variational functional of the one-electron Green’s function that yields both the correct total energy and all one-particle properties at the variational extremum. However, compared to DFT, these methods face a number of additional challenges beyond the familiar problem of choosing an exchange-correlation functional: (a) they have a much greater computational cost compared to DFT, (b) one must choose some set of trial Green’s functions, and (c) one must locate an extremum (not necessarily a minimum). This talk aims both to be pedagogical and to provide a summary of our most recent work in this area. We focus mainly on the random-phase approximation (RPA) for the correlation functional which yields the GW self-energy for the Green’s function. We will describe ecient approaches for computing the RPA correlation energy; how various types of approximations to the correlation energy create a ladder of approximate self-energy operators and what type of physics is included in these approximations (exchange, screening, Hubbard U, etc.); and how the extremization is in fact highly problematic. We will end with some ideas for making practical progress with the extremization and their relation to the successful Quasiparticle Self-consistent GW (QSGW) scheme [1].
In principle, the Luttinger-Ward Green's-function formalism allows one to compute simultaneously the total energy and the quasiparticle band structure of a many-body electronic system from first principles. We present approximate and exact expressions for the correlation energy within the GW-random-phase approximation that are more amenable to computation and allow for developing efficient approximations to the self-energy operator and correlation energy. The exact form is a sum over differences between plasmon and interband energies. The approximate forms are based on summing over screened interband transitions. We also demonstrate that blind extremization of such functionals leads to unphysical results: imposing physical constraints on the allowed solutions (Green's functions) is necessary. Finally, we present some relevant numerical results for atomic systems.
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.
A new scheme of first-principles computation for strongly correlated electron systems is proposed. This scheme starts from the local-density approximation (LDA) at high-energy band structure, while the low-energy effective Hamiltonian is constructed by a downfolding procedure using combinations of the constrained-LDA and the GW method. The obtained low-energy Hamiltonian is solved by the path-integral renormalization-group method, where spatial and dynamical fluctuations are fully considered. An application to Sr2VO4 shows that the scheme is powerful in agreement with experimental results. It further predicts a nontrivial orbital-stripe order.