We compute the spin susceptibility of the two-dimensional Hubbard model away from half-filling, and analyze the impact of frequency-dependent vertex corrections as obtained from the dynamical mean field theory (DMFT). We find that the local dynamics captured by the DMFT vertex strongly affects nonlocal spin correlations, and thus the momentum dependence of the spin susceptibility. While the widely used random-phase approximation yields commensurate N\'eel-type antiferromagnetism as the dominant instability over a wide doping range, the vertex corrections favor incommensurate ordering wave vectors away from $(\ensuremath{\pi},\ensuremath{\pi})$. Our results indicate that the connection between the magnetic ordering wave vector and the Fermi surface geometry, familiar for weakly interacting systems, can hold in a strongly correlated metal, too.
We compute the spin susceptibility of the two-dimensional Hubbard model away from half-filling, and analyze the impact of frequency-dependent vertex corrections as obtained from the dynamical mean field theory (DMFT). We find that the local dynamics captured by the DMFT vertex strongly affects nonlocal spin correlations, and thus the momentum dependence of the spin susceptibility. While the widely used random-phase approximation yields commensurate N'eel-type antiferromagnetism as the dominant instability over a wide doping range, the vertex corrections favor incommensurate ordering wave vectors away from (\ensuremathπ,\ensuremathπ). Our results indicate that the connection between the magnetic ordering wave vector and the Fermi surface geometry, familiar for weakly interacting systems, can hold in a strongly correlated metal, too.
The treatment of strongly interacting two-dimensional Fermi systems constitutes one of the most challenging problems in the field of condensed matter physics. Many theoretical works focused on the strongly coupled Hubbard model, since it is expected to capture the most important physics of cuprate superconductors. Due to our methodological improvements and understanding of the frequency dependence of the two-particle vertex function, the functional renormalization group combined with the dynamical mean field theory can now be used to study competing correlations in the strongly interacting regime. While limited to moderate interaction strengths, the functional renormalization group describes efficiently systems with a hierarchy of different energy scales and competing correlations. For instance, it provides definite evidence for d-wave superconductivity in the two-dimensional Hubbard model at moderate coupling. In a first project, we study the frequency dependence of the one-particle irreducible vertex function generated by the functional renormalization group flow. The frequency dependence, which becomes singular for strong interactions, appears to be important already for moderate couplings, and it cannot be represented by separate channels each depending on a single linear combination of frequencies. For strongly interacting systems, the dynamical mean field theory captures strong local correlation effects nonperturbatively. With this approximation, we study the impact of local correlations on the magnetic susceptibility. The local dynamics strongly affects the spin response function and its momentum dependence. In contrast to the widely used random-phase approximation with self-energy corrections that predicts Néel antiferromagnetic order, the local vertex corrections favor an incommensurate order similar to the ordering instability predicted by the Fermi surface geometry, as for weakly interacting systems. The dynamical mean field theory is also used as starting point for the functional renormalization group flow. We demonstrate that, due to our improvements for the parametrization of the vertex function, this approach is actually able to access the strong coupling physics. Moreover, we derive a flow scheme that conserves the local contributions and reduces the truncation error of the flow equations. In the strongly interacting regime, we capture strong d-wave pairing correlations driven by magnetic fluctuations with a mechanism similar to the one observed in the weakly interacting system.
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.
Using the dynamical mean-field theory (DMFT) as a ``booster rocket'', the functional renormalization group (fRG) can be upgraded from a weak-coupling method to a powerful computation tool for strongly interacting fermion systems. The strong local correlations are treated nonperturbatively by the DMFT, while the fRG flow can be formulated such that it is driven exclusively by nonlocal correlations, which are more amenable to approximations. We show that the full frequency dependence of the two-particle vertex needs to be taken into account in this approach, and demonstrate that this is actually possible---in spite of the singular frequency dependence of the vertex at strong coupling. We are thus able to present results obtained from the DMFT-boosted fRG for the two-dimensional Hubbard model in the strongly interacting regime. We find strong antiferromagnetic correlations from half filling to 18% hole doping and at the lowest temperature we can access, a sizable d-wave pairing interaction driven by magnetic correlations at the edge of the antiferromagnetic regime.
Away from half filling and above its Neel temperature, the three-dimensional Hubbard model is a paramagnetic metal with antiferromagnetic spin fluctuations, where controlled momentum-resolved results have been lacking. Using connected-determinant diagrammatic Monte Carlo, we compute the momentum-resolved self-energy for n = 0.4-0.95, T/t = 0.05-0.33 and U/t = 4-8, and the full-zone spin susceptibility at twelve of these points, focusing on U/t = 4. At U/t = 4 and toward half filling the spin fluctuations become strong: for n = 0.7-0.9 at T/t = 0.14 they enhance the susceptibility peak up to sevenfold and raise the momentum-dependent part of the scattering rate to four to five times its second-order value. They are driven by an effective spin interaction screened to about 2U/3 and nearly uniform over the zone. Yet the peak height is set by the slightly smaller value of this interaction at the peak wave vector, and near magnetic order this small difference changes the height by a large factor. Two-particle self-consistent theory reproduces the peak because its sum rule fixes its single constant interaction near this value, but for n = 0.7-0.9 at T/t = 0.14 its self-energy overestimates the local scattering rate and underestimates the momentum-dependent part. Though strong, the fluctuations are short ranged, with a correlation length of at most about two lattice spacings, and for n = 0.7-0.9 and T/t >= 0.14 the scattering rate near the Fermi surface varies by at most about 1% at fixed band energy, without hot spots. At U/t = 4 dynamical mean-field theory (DMFT) reproduces the local scattering rate and quasiparticle weight within a few percent, and the nonlocal part it misses is approximately a function of the band energy. Strong spin fluctuations thus coexist with nearly isotropic quasiparticles, and the results provide a benchmark for extensions of DMFT and for spin-fluctuation theories.
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.
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.