In this paper, we investigate how nonlocal correlations affect, selectively, the physics of correlated electrons over different energy scales, from the Fermi level to the band edges. This goal is achieved by applying a diagrammatic extension of dynamical mean field theory (DMFT), the dynamical vertex approximation $(\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A})$, to study several spectral and thermodynamic properties of the unfrustrated Hubbard model in two and three dimensions. Specifically, we focus first on the low-energy regime by computing the electronic scattering rate and the quasiparticle mass renormalization for decreasing temperatures at a fixed interaction strength. This way, we obtain a precise characterization of the several steps through which the Fermi-liquid physics is progressively destroyed by nonlocal correlations. Our study is then extended to a broader energy range, by analyzing the temperature behavior of the kinetic and potential energy, as well as of the corresponding energy distribution functions. Our findings allow us to identify a smooth but definite evolution of the nature of nonlocal correlations by increasing interaction: They either increase or decrease the kinetic energy w.r.t. DMFT depending on the interaction strength being weak or strong, respectively. This reflects the corresponding evolution of the ground state from a nesting-driven (Slater) to a superexchange-driven (Heisenberg) antiferromagnet (AF), whose fingerprints are, thus, recognizable in the spatial correlations of the paramagnetic phase. Finally, a critical analysis of our numerical results of the potential energy at the largest interaction allows us to identify possible procedures to improve the ladder-based algorithms adopted in the dynamical vertex approximation.
In this paper, we investigate how nonlocal correlations affect, selectively, the physics of correlated electrons over different energy scales, from the Fermi level to the band edges. This goal is achieved by applying a diagrammatic extension of dynamical mean field theory (DMFT), the dynamical vertex approximation (D\ensuremathΓA), to study several spectral and thermodynamic properties of the unfrustrated Hubbard model in two and three dimensions. Specifically, we focus first on the low-energy regime by computing the electronic scattering rate and the quasiparticle mass renormalization for decreasing temperatures at a fixed interaction strength. This way, we obtain a precise characterization of the several steps through which the Fermi-liquid physics is progressively destroyed by nonlocal correlations. Our study is then extended to a broader energy range, by analyzing the temperature behavior of the kinetic and potential energy, as well as of the corresponding energy distribution functions. Our findings allow us to identify a smooth but definite evolution of the nature of nonlocal correlations by increasing interaction: They either increase or decrease the kinetic energy w.r.t. DMFT depending on the interaction strength being weak or strong, respectively. This reflects the corresponding evolution of the ground state from a nesting-driven (Slater) to a superexchange-driven (Heisenberg) antiferromagnet (AF), whose fingerprints are, thus, recognizable in the spatial correlations of the paramagnetic phase. Finally, a critical analysis of our numerical results of the potential energy at the largest interaction allows us to identify possible procedures to improve the ladder-based algorithms adopted in the dynamical vertex approximation.
In this paper, we investigate the impact of nonlocal correlations on charge fluctuations in the two-dimensional single-band Hubbard model close to the Mott metal-to-insulator transition, employing the ladder dynamical vertex approximation. At half filling and for interaction strengths and temperatures where the system is in the Mott insulating phase, charge fluctuations are strongly suppressed. Under these conditions, dynamical mean-field theory (DMFT) calculations predict a strong enhancement of the charge susceptibility at small (electron or hole) doping. However, these DMFT results include only the effects of purely local correlations despite the importance of nonlocal correlations in two-dimensional systems. We have, hence, carried out ladder dynamical vertex approximation $(\mathrm{lD}\mathrm{\ensuremath{\Gamma}}\mathrm{A})$ simulations which allow for the inclusion of such nonlocal correlation effects while retaining the local ones of DMFT. Our $\mathrm{lD}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ numerical data show that close to half filling the large uniform charge susceptibility of DMFT is strongly suppressed by nonlocal fluctuations but gradually increases with (electron) doping. At a certain doping value, charge fluctuations eventually become larger in $\mathrm{lD}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ with respect to DMFT, indicating that the absence of nonlocal correlations underestimates the mobility of the charge carriers in this parameter regime. This metallization effect is also reflected in an enhancement of the $\mathrm{lD}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ kinetic and potential energies and a corresponding reduction of the (absolute value of the) $\mathrm{lD}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ Matsubara self-energy with respect to DMFT.
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.
In the last decades, dynamical mean-field theory (DMFT) and its diagrammatic extensions have been successfully applied to describe local and nonlocal correlation effects in correlated electron systems. Unfortunately, except for the exact solution, it is impossible to fulfill both the Pauli principle and conservation laws at the same time. Consequently, fundamental observables such as the kinetic and potential energies are ambiguously defined. In this work, we propose an approach to overcome the ambiguity in the calculation of the potential energy within the ladder dynamical vertex approximation ($\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$) by introducing an effective mass renormalization parameter in both the charge and the spin susceptibility of the system. We then apply our method to the half-filled single-band Hubbard model on a three-dimensional bipartite cubic lattice. We find (i) at weak-to-intermediate coupling, a reasonable modification of the transition temperature TN to the antiferromagnetically ordered state with respect to previous ladder $\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ calculations without charge renormalization. This is in good agreement with dual fermion and Monte Carlo results; (ii) the renormalization of charge fluctuations in our new approach leads to a unique value for the potential energy which is substantially lower than corresponding ones from DMFT and non-self-consistent ladder $\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$; and (iii) the hierarchy of the kinetic energies between the DMFT and the ladder $\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ in the weak coupling regime is restored by the consideration of charge renormalization.
We have implemented the dynamical vertex approximation $(\mathrm{D}\ensuremath{\Gamma}\mathrm{A})$ in its full parquet-based version to include spatial correlations on all length scales and in all scattering channels. The algorithm is applied to study the electronic self-energies and the spectral properties of finite-size one-dimensional Hubbard models with periodic boundary conditions (nanoscopic Hubbard rings). From a methodological point of view, our calculations and their comparison to the results obtained within dynamical mean-field theory, plain parquet approximation, and the exact numerical solution allow us to evaluate the performance of the $\mathrm{D}\ensuremath{\Gamma}\mathrm{A}$ algorithm in the most challenging situation of low dimensions. From a physical perspective, our results unveil how nonlocal correlations affect the spectral properties of nanoscopic systems of various sizes in different regimes of interaction strength.
In this work, we adapt the formalism of the dynamical vertex approximation ($\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$), a diagrammatic approach including many-body correlations beyond the dynamical mean-field theory, to the case of attractive on-site interactions. We start by exploiting the ladder approximation of the $\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ scheme, in order to derive the corresponding equations for the nonlocal self-energy and vertex functions of the attractive Hubbard model. Second, we prove the validity of our derivation by showing that the results obtained in the particle-hole symmetric case fully preserve the exact mapping between the attractive and the repulsive models. It will be shown how this property can be related to the structure of the ladders, which makes our derivation applicable for any approximation scheme based on ladder diagrams. Finally, we apply our $\mathrm{D}\mathrm{\ensuremath{\Gamma}}\mathrm{A}$ algorithm to the attractive Hubbard model in three dimensions, for different fillings and interaction values. Specifically, we focus on the parameters region in the proximity of the second-order transition to the superconducting and charge-density wave phases, respectively, and calculate (i) their phase-diagrams, (ii) their critical behavior, as well as (iii) the effects of the strong nonlocal correlations on the single-particle properties.
We analyze the impact of the lattice geometry on the thermodynamic transition to magnetically ordered phases in strongly interacting electron systems for various Bravais lattices in three, four, and five dimensions, including both local and nonlocal correlation effects. In a first step we use the dynamical mean-field theory (DMFT), which takes into account purely local correlations, to calculate the magnetic susceptibilities of the Hubbard model on three- (3d-sc), four- (4d-sc), and five- dimensional (5d-sc) simple cubic/hypercubic, as well as on three- dimensional body- (bcc) and face-centered (fcc) cubic lattices, and determine the transition temperature to the corresponding magneticallyordered state. In a second step, we exploit the dynamical vertex approximation ( D Γ A ) , a diagrammatic extension of DMFT, to include the effect of nonlocal correlations which are particularly important in the vicinity of the corresponding phase transition. For the bipartite 3d-sc, 4d-sc, 5d-sc, and bcc lattices nonlocal fluctuations lead to a substantial reduction of the DMFT transition temperature consistent to the overall tendency of mean-field approaches to overestimate the stability of ordered phases. As expected, the magnitude of the difference between the DMFT, being exact in the limit of large connectivity/dimensions, and D Γ A transition temperatures decreases with increasing coordination number. On a more practical perspective, these results also provide a reasonable guidance to evaluate the expected overestimation of the DMFT ordering temperature for different material geometries. For the fcc lattice, on the other hand, the ordered phase observed in DMFT vanishes completely within D Γ A which is consistent with the existence of strong geometric frustration in this lattice.