We describe a recent implementation of the combined
GW and dynamical mean field method (
GW+DMFT) for the two-dimensional Hubbard model with onsite and nearest-neighbor repulsion. We clarify the relation of the
GW+DMFT scheme to alternative approaches in the literature, and discuss the corresponding approximations to the free-energy functional of the model. Furthermore, we describe a numerically exact technique for the solution of the
GW+DMFT equations, namely, the hybridization expansion continuous-time algorithm for impurity models with retarded interactions. We compute the low-temperature phase diagram of the half-filled extended Hubbard model, addressing the metal-insulator transition at small intersite interactions and the transition to a charge-ordered state for stronger intersite repulsions.
GW+DMFT introduces a nontrivial momentum dependence into the many-body self-energy and polarization. We find that the charge fluctuations included in the present approach have a larger impact on the latter than on the former. Finally, within the
GW+DMFT framework, as in extended DMFT, the intersite repulsion translates into a frequency dependence of the local effective interaction. We analyze this dependence and show how it affects the local spectral function.
The
GW approximation is the standard tool for quasiparticle predictions in materials, yet its regime of validity in correlated systems remains poorly quantified, because \textit{ab initio} vertex corrections are computationally prohibitive. Using exact diagonalization of the half-filled extended Hubbard model on finite rings as a numerically exact reference, we construct the corresponding model-space
GW theory on the identical Hilbert space and quantify its error as a function of local (
U) and non-local (
V) interaction strength. We find that the required vertex correction changes character across the phase diagram: in the weak-coupling regime the effective vertex
Γeff<1, reflecting the suppression of RPA charge fluctuations by exact short-range correlations, whereas in the Mott regime
Γeff grows monotonically (to
∼3 at
U=8t for
N=6), reflecting the local dynamical self-energy structure required to open the Hubbard gap. Vertex corrections in the electron-hole (polarizability) channel are shown to \emph{worsen} the gap error, indicating that reproducing the Mott gap requires dynamical self-energy structure rather than improved screening. For
V=0, static COHSEX is accidentally exact at a single crossover
U∗≈3.5t; finite
V, through non-local Fock exchange, splits this point into a double-crossover window that collapses toward weak coupling. Even at the crossover, however, the exact spectral function retains Hubbard-band structure that no static functional reproduces, so gap agreement does not imply functional accuracy. These results yield quantitative diagnostics for the reliability of
GW in correlated materials.
In this thesis, the GW approximation (GWA, Green's function G times screened interaction W) and the Hubbard I approximation (HIA) are combined in a non-self-consistent one-shot calculation to determine the electronic structure of a one-dimensional strongly-correlated model. The scheme was chosen to incorporate both screening effects through the GWA and strong on-site correlations through the atomic HIA. The resulting self-energies of both methods are summed together and another term is subtracted to correct double-counting. The examined double-counting terms were the local self-energy of the GWA result, the GWA self-energy of the atomic model used in the HIA and lastly the impurity self-energy within the GWA. The self-energy of the atomic model in the HIA was solved analytically. It was found that all three approaches yield non-causal features, which increase with the on-site Coulomb repulsion, in the resulting spectral function. The local GWA self-energy was observed to perform best in terms of causality and computational effort. Changes in the resulting quasi-particle structure in the spectral function showed that screening effects and local correlations were included successfully.
The time evolution in quantum many-body systems after external excitations is attracting high interest in many fields, including dense plasmas, correlated solids, laser-excited materials, or fermionic and bosonic atoms in optical lattices. The theoretical modeling of these processes is challenging, and the only rigorous quantum-dynamics approach that can treat correlated fermions in two and three dimensions is nonequilibrium Green functions (NEGF). However, NEGF simulations are computationally expensive due to their
T3 scaling with the simulation duration
T. Recently,
T2 scaling was achieved with the generalized Kadanoff-Baym ansatz (GKBA), for the second-order Born (SOA) self energy, which has substantially extended the scope of NEGF simulations. In a recent Letter [Schl\"unzen et al., Phys. Rev. Lett. 124, 076601 (2020).] we demonstrated that GKBA-NEGF simulations can be efficiently mapped onto coupled time-local equations for the single-particle and two-particle Green functions on the time diagonal, hence the method has been called the G1-G2 scheme. This allows one to perform the same simulations with order
T1 scaling, both for SOA and
GW self energies giving rise to a dramatic speedup. Here we present more details on the G1-G2 scheme, including derivations of the basic equations including results for a general basis, for Hubbard systems, and for jellium. Also, we demonstrate how to incorporate initial correlations into the G1-G2 scheme. Further, the derivations are extended to a broader class of self energies, including the
T matrix in the particle-particle and particle-hole channels and the dynamically-screened-ladder approximation. Finally, we demonstrate that, for all self energies, the CPU-time scaling of the G1-G2 scheme with the basis dimension
Nb can be improved compared to our first report: The overhead compared to the original GKBA is not more than an additional factor
Nb, even for Hubbard systems.
We present a fully self-consistent combined GW and dynamical mean field (DMFT) study of the extended two-dimensional Hubbard model. The inclusion of the local dynamical vertex stemming from the DMFT self-energy and polarization is shown to cure the known problems of self-consistent GW. We calculate momentum-resolved spectral functions, two-particle polarizations, and electron-loss spectra, as well as the effective dynamical interaction induced by nonlocal screening. The momentum-dependence introduced by GW into the extended DMFT description leads to a narrowing of the quasiparticle width and more pronounced Hubbard bands in the metallic regime as one approaches the charge-ordering transition. It further affects the shape of collective modes, giving rise to dispersive plasmon-like long-wavelength and stripe modes.
We explore the two-dimensional extended Hubbard model on the triangular lattice in the high doping regime. On-site and nearest-neighbour repulsive interactions are treated in a non-perturbative way by means of Extended Dynamical Mean Field Theory. We compute the low-temperature phase diagram, displaying a metallic phase and a symmetry-broken phase for strong intersite repulsions. We describe the correlation effects on both single-particle and two-particle observables in the metallic phase. Whereas single-particle spectra feature a Hubbard satellite typical of strongly correlated systems, local susceptibilities remain close to their non-interacting limit, even for large on-site repulsions. We argue that this behaviour is typical of the strongly doped case. We also report a region in parameter space with negative static local screening.