Quantitative simulation of electronic structure of solids requires treating local and nonlocal electron correlations on an equal footing. We present a new ab initio formulation of Green's function embedding which, unlike dynamical mean-field theory that uses noninteracting bath, derives bath representation with general two-particle interactions in a systematically improvable manner. The resulting interacting-bath dynamical embedding theory (ibDET) utilizes an efficient real-axis coupled-cluster solver to compute the self-energy, approaching the full system limit at much reduced cost. When combined with the GW theory, GW+ibDET achieves good agreement with experimental spectral properties across a range of semiconducting, insulating, and metallic materials. Our approach also enables quantifying the role of nonlocal electron correlation in determining material properties and addressing the long-standing debate on the bandwidth narrowing of metallic sodium.
Quantitative simulation of electronic structure of solids requires treating local and nonlocal electron correlations on an equal footing. We present a new ab initio formulation of Green's function embedding which, unlike dynamical mean-field theory that uses noninteracting bath, derives bath representation with general two-particle interactions in a systematically improvable manner. The resulting interacting-bath dynamical embedding theory (ibDET) utilizes an efficient real-axis coupled-cluster solver to compute the self-energy, approaching the full system limit at much reduced cost. When combined with the GW theory, GW+ibDET achieves good agreement with experimental spectral properties across a range of semiconducting, insulating, and metallic materials. Our approach also enables quantifying the role of nonlocal electron correlation in determining material properties and addressing the long-standing debate on the bandwidth narrowing of metallic sodium.
We present a molecular extension of our recently proposed Green’s function embedding method, interacting-bath dynamical embedding theory (ibDET), for computing charged excitation energies at the GW and EOM-CCSD levels. Starting from atom-centered impurities, we construct bath representations that capture the frequency-dependent entanglement between the impurity and its environment and can be systematically improved via the construction of cluster-specific natural orbitals. Utilizing a GW or coupled-cluster Green’s function solver, the self-energy of the full system is assembled from all embedding problems to obtain the interacting Green’s function. We show that ibDET provides accurate spectral properties with much reduced cost for a broad range of systems, including conjugated molecules and nanoclusters. Compared with full-system results, the errors in the predicted ionization potentials and electron affinities are around 0.1 eV or smaller, while each embedding problem includes only a small fraction of the total orbital space. This work provides an efficient and scalable framework for computing spectral properties of molecular systems.
We present an efficient ab initio dynamical mean-field theory (DMFT) implementation for quantitative simulations in solids. Our DMFT scheme employs ab initio Hamiltonians defined for impurities comprising the full unit cell or a supercell of atoms and for realistic quantum chemical basis sets. We avoid double counting errors by using Hartree-Fock as the low-level theory. Intrinsic and projected atomic orbitals (IAO + PAO) are chosen as the local embedding basis, facilitating numerical bath truncation. Using an efficient integral transformation and coupled-cluster Green's function impurity solvers, we are able to handle embedded impurity problems with several hundred orbitals. We apply our ab initio DMFT approach to study a hexagonal boron nitride monolayer, crystalline silicon, and nickel oxide in the antiferromagnetic phase, with up to 104 and 78 impurity orbitals in the spin-restricted and unrestricted cluster DMFT calculations and over 100 bath orbitals. We show that our scheme produces accurate spectral functions compared to both benchmark periodic coupled-cluster computations and experimental spectra.
Periodic equation-of-motion coupled-cluster theory with single and double excitations (EOM-CCSD) has shown promise for quantitative calculations of band structures in solids. However, its steep computational scaling has limited calculations to relatively coarse k-point meshes, leading to sizable finite-size errors and discrepant estimates of thermodynamic-limit band gaps in recent benchmarks. In this work, we revisit EOM-CCSD band gaps for ten semiconductors and insulators using interacting-bath dynamical embedding theory (ibDET), a systematically improvable Green's function embedding framework that enables dense Brillouin-zone sampling at modest computational cost. By pushing the k-point sampling up to 10 × 10 × 10, well beyond the system sizes accessible in canonical periodic EOM-CCSD calculations, we significantly reduce finite-size errors and obtain stable thermodynamic-limit extrapolations. We further compare G0W0@PBE, G0W0@HF, and EOM-CCSD on an equal footing using the same numerical settings in PySCF. We find that EOM-CCSD yields a mean absolute error of 0.32 eV relative to experimental band gaps for a test set of ten semiconductors and insulators, which is lower than that of G0W0@PBE. For ZnO, EOM-CCSD also accurately describes the Zn 3d-band binding energy, despite overestimating the band gap. These results demonstrate that ibDET offers a practical route to high-accuracy many-body electronic structure calculations in periodic systems.
We present the concept, derivation, and implementation of dynamical configuration interaction, a quantum embedding theory that combines Green's function methodology with the many-body wave function. In a strongly correlated active space, we use full configuration interaction (CI) to describe static correlation exactly. We add energy-dependent corrections to the CI Hamiltonian which, in principle, include all remaining correlations derived from the bath space surrounding the active space. Next, we replace the exact Hamiltonian in the bath with one of excitations defined over a correlated ground state. This transformation is naturally suited to the methodology of many-body Green's functions. In this space, we use a modified GW/Bethe-Salpeter equation procedure to calculate excitation energies. Combined with an estimate of the ground-state energy in the bath, we can efficiently compute the energy-dependent corrections, which correlate the full set of orbitals, for very low computational cost. We present dimer dissociation curves for H2 and N2 in good agreement with exact results. Additionally, excited states of N2 and C2 are in excellent agreement with benchmark theory and experiment. By combining the strengths of two disciplines, we achieve a balanced description of static and dynamic correlation in a fully ab initio, systematically improvable framework.
In strongly correlated materials, quantum fluctuations taking place at the level of atomic sites play a central role. These local processes are captured nonperturbatively in dynamical mean field theory (DMFT). Nonlocal mechanisms, however, are also suspected to be essential for the proper description of some of these materials, including the well-known high-temperature cuprate superconductors and the more exotic systems of atoms adsorbed on semiconducting surfaces.In this thesis, the effect of nonlocal, potentially long-range fluctuations beyond the DMFT framework is investigated from two angles.First, long-ranged interactions in adatom systems, computed from first principles within the constrained random-phase approximation (cRPA), are treated within the extended DMFT (EDMFT) framework and its combination with the diagrammatic GW method (GW +EDMFT). They are shown to explain some material trends observed in these compounds. In addition, thanks to the spatial and temporal information delivered by GW+EDMFT, they solve existing puzzles raised by different experimental conclusions drawn from different probes.Second, a new approach to the Hubbard model including nonlocal correlation effects is proposed. Formally based on the local expansion of the triply-irreducible functional, this approach, dubbed TRILEX, approximates the renormalized three-leg electron-boson vertex by a local but frequency-dependent vertex computed from an impurity model with retarded density-density and spin-spin interactions. This local vertex is used to construct frequency- and momentum-dependent self-energies capturing nonlocal correlation effects. By construction, this method reconciles two theoretical pictures of cuprate superconductors, namely the spin fluctuation theory at weak interaction strengths and the doped-Mott-insulator picture at strong coupling.This novel approach is applied to the single-band Hubbard model. It is shown to reduce to DMFT at strong coupling and high temperatures, and to fluctuation-exchange theories at weak coupling, while substantial deviations from both limits are observed at intermediate coupling. There, upon doping, antiferromagnetic spin fluctuations are observed to lead to a momentum-differentiated spectral intensity on the Fermi surface, which is reminiscent of the formation of Fermi arcs in photoemission measurements on cuprate superconductors.These two methodological developments are intrinsically linked to algorithmic aspects which are reviewed in the last part of this work. In particular, an open-source library dedicated to massively parallel computations in the field of quantum many-body physics, TRIQS, is presented, with a focus on error estimation through the binning and jackknife methods. The hybridization-expansion continuous-time quantum Monte-Carlo algorithm in the segment picture -- central to the TRILEX method -- and its state-of-the-art implementation are explained. Finally, as a further illustration, a new lightweight method relying on this algorithm and designed for realistic systems, the combined screened-exchange and dynamical DMFT method (SEx+DDMFT), is presented together with an application to the BaCo2As2 compound.
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.