arXiv:1803.03036 · Physical review. B./Physical review. B

Efficient Bethe-Salpeter equation treatment in dynamical mean-field theory

Agnese Tagliavini, Stefan Hummel, Nils Wentzell, Sabine Andergassen, A. Toschi, G. Rohringer

Published Jun 22, 2018• 40 citations

Abstract

We present here two alternative schemes designed to correct the high-frequency truncation errors in the numerical treatment of the Bethe-Salpeter equations. The schemes are applicable to all Bethe-Salpeter calculations with a local two-particle irreducible vertex, which is relevant, e.g., for the dynamical mean-field theory (DMFT) and its diagrammatic extensions. In particular, within a purely diagrammatic framework, we could extend existing algorithms for treating the static case in the particle-hole sector to more general procedures applicable to all bosonic frequencies and all channels. After illustrating the derivation and the theoretical interrelation of the two proposed schemes, these have been applied to the Bethe-Salpeter equations for the auxiliary Anderson impurity models of selected DMFT calculations, where results can be compared against a numerically ``exact'' solution. The successful performance of the proposed schemes suggests that their implementation can significantly improve the accuracy of DMFT calculations at the two-particle level, in particular for more realistic multiorbital calculations where the large number of degrees of freedom substantially restricts the actual frequency range for numerical calculations, as well as---on a broader perspective---of the diagrammatic extensions of DMFT.

DOI: 10.1103/physrevb.97.235140

Related papers

Improved treatment of fermion-boson vertices and Bethe-Salpeter equations in nonlocal extensions of dynamical mean field theory

A. A. Katanin

We reconsider the procedure of calculating fermion-boson vertices and the numerical solution of Bethe-Salpeter equations, used in nonlocal extensions of dynamical mean field theory. Because of the frequency dependence of the vertices, a finite-frequency box for matrix inversions is typically used, which often requires some treatment of the asymptotic behavior of vertices. Recently [J. Kune\ifmmode \check{s}\else \v{s}\fi{}, Phys. Rev. B 83, 085102 (2011); A. Tagliavini et al., Phys. Rev. B 97, 235140 (2018)] it was proposed to split the considered frequency box into smaller and larger ones; in the smaller-frequency box the numerically exact vertices are used, while beyond this box asymptotics of vertices are applied. Yet, this method requires numerical treatment of vertex asymptotics (including corresponding matrix manipulations) in the larger-frequency box and/or knowledge of fermion-boson vertices, which may not be convenient for numerical calculations. In the present paper we derive the formulas which treat analytically the contribution of vertices beyond the chosen frequency box, such that only numerical operations with vertices in the chosen small-frequency box are required. The method is tested on the Hubbard model and can be used in a broad range of applications of nonlocal extensions of dynamical mean field theory.

arXiv:1908.03198Physical review. B./Physical review. BJan 8, 202016 citations

Dual Bethe-Salpeter equation for the multiorbital lattice susceptibility within dynamical mean-field theory

Erik G. C. P. van Loon, Hugo U. R. Strand

Dynamical mean-field theory describes the impact of strong local correlation effects in many-electron systems. While the single-particle spectral function is directly obtained within the formalism, two-particle susceptibilities can also be obtained by solving the Bethe-Salpeter equation. The solution requires handling infinite matrices in Matsubara frequency space. This is commonly treated using a finite frequency cutoff, resulting in slow linear convergence. A decomposition of the two-particle response in local and nonlocal contributions enables a reformulation of the Bethe-Salpeter equation inspired by the dual boson formalism. The reformulation has a drastically improved cubic convergence with respect to the frequency cutoff, considerably facilitating the calculation of susceptibilities in multi-orbital systems. This improved convergence arises from the fact that local contributions can be measured in the impurity solver. The dual Bethe-Salpeter equation uses the fully reducible vertex which is free from vertex divergences. We benchmark the approach on several systems including the spin susceptibility of strontium ruthenate Sr2RuO4 , a strongly correlated Hund's metal with three active orbitals. Published by the American Physical Society 2024

arXiv:2306.05157Physical review. B./Physical review. BApr 23, 20247 citations

High-frequency asymptotics of the vertex function: Diagrammatic parametrization and algorithmic implementation

Nils Wentzell, Gang Li, Agnese Tagliavini +5 more

Vertex functions are a crucial ingredient of several forefront many-body algorithms in condensed matter physics. However, the full treatment of their frequency and momentum dependence severely restricts numerical calculations. A significant advancement requires an efficient treatment of the high-frequency asymptotic behavior of the vertex functions. In this work, we first provide a detailed diagrammatic analysis of the high-frequency structures and their physical interpretation. Based on these insights, we propose a parametrization scheme, which captures the whole high-frequency domain for arbitrary values of the Coulomb interaction and electronic density, and we discuss the details of its algorithmic implementation in many-body solvers based on parquet equations as well as functional renormalization group schemes. Finally, we assess its validity by comparing our results for a single impurity Anderson model with exact diagonalization calculations. The proposed parametrization is pivotal for the algorithmic development of all quantum many-body methods based on vertex functions arising from both local and nonlocal static microscopic interactions as well as effective dynamic interactions which uniformly approach a static value for large frequencies. In this way, our present technique can substantially improve vertex-based diagrammatic approaches including spatial correlations beyond dynamical mean-field theory.

arXiv:1610.06520Physical review. B./Physical review. BAug 3, 2020104 citations

Improved effective vertices in the multiorbital two-particle self-consistent method from dynamical mean-field theory

Karim Zantout, Steffen Backes, Aleksandar Razpopov +2 more

In this paper, we present a multiorbital form of the two-particle self-consistent approach (TPSC) where the effective local and static irreducible interaction vertices are determined by means of the dynamical mean-field theory (DMFT). This approach replaces the approximate ansatz equations for the double-occupations $\ensuremath{\langle}{n}_{\ensuremath{\alpha},\ensuremath{\sigma}}^{}{n}_{\ensuremath{\beta},{\ensuremath{\sigma}}^{\ensuremath{'}}}^{}\ensuremath{\rangle}$ by sampling them directly for the same model using DMFT. Compared to the usual Hartree-Fock-like ansatz, this leads to more accurate local vertices in the weakly correlated regime and provides access to stronger correlated systems that were previously out of reach. This approach is extended by replacing the local component of the TPSC self-energy by the DMFT impurity self-energy, which results in an improved self-energy that incorporates strong local correlations but retains a nontrivial momentum dependence. We find that this combination of TPSC and DMFT provides a significant improvement over the multiorbital formulation of TPSC as it allows to determine the components of the spin vertex without artificial symmetry assumptions and opens the possibility to include the transversal particle-hole channel. The new approach is also able to remove unphysical divergences in the charge vertices in TPSC. We find a general trend that lower temperatures can be accessed in the calculation. Benchmarking single-particle quantities, such as the local spectral function with other many-body methods, we find significant improvement in the more strongly correlated regime.

arXiv:2211.01400Physical review. B./Physical review. BJun 1, 202312 citations

One-particle irreducible functional approach: A route to diagrammatic extensions of the dynamical mean-field theory

G. Rohringer, A. Toschi, Hartmut Hafermann +3 more

We present an approach which is based on the one-particle irreducible (1PI) generating functional formalism and includes electronic correlations on all length scales beyond the local correlations of dynamical mean-field theory (DMFT). This formalism allows us to unify aspects of the dynamical vertex approximation (D$\ensuremath{\Gamma}$A) and the dual fermion (DF) scheme, yielding a consistent formulation of nonlocal correlations at the one- and two-particle level beyond DMFT within the functional integral formalism. In particular, the considered approach includes one-particle reducible contributions from the three- and more-particle vertices in the dual fermion approach, as well as some diagrams not included in the ladder version of D$\ensuremath{\Gamma}$A. To demonstrate the applicability and physical content of the 1PI approach, we compare the diagrammatics of 1PI, DF, and D$\ensuremath{\Gamma}$A, as well as the numerical results of these approaches for the half-filled Hubbard model in two dimensions.

arXiv:1301.7546Physical Review BSep 9, 201388 citations

Dynamical Mean-Field Theory Studies of Correlated Electrons

Yang Yu

Correlated electron systems exhibit a variety of intriguing emergent phenomena, including Mott physics, pseudogap behavior, and unconventional superconductivity. One of the most widely used numerical approaches for investigating these systems is dynamical mean-field theory (DMFT). This dissertation addresses two related problems: how to develop better quantum impurity solvers for DMFT and how to use DMFT results to gain a better understanding of correlated electron systems. The first part of this dissertation develops numerical methods for solving quantum impurity problems, which constitute the central computational task in DMFT calculations. In Chap. 4, we develop a method that combines a tensor-train numerical integration scheme with the inchworm hybridization expansion for quantum impurity problems. Benchmarks on spinless, single-orbital, and multi-orbital impurity models demonstrate the convergence of the Green’s functions to exact or reference results and identify the regimes in which tensor-train compression of Feynman diagrams and the inchworm expansion itself are effective. In Chap. 5, we introduce the concept of counterterms within the hybridization expansion and apply it to a bare hybridization expansion solver that relies on a sum-of-exponentials integration technique. Tests on spinless, multi-orbital, and realistic embedding impurity problems show that counterterms can significantly reduce the required perturbative order and improve the efficiency of hybridization-expansion calculations. The second part of this dissertation applies the cluster extension of DMFT to the single-band Hubbard model. In Chap. 6, we introduce a fluctuation decomposition for the self-energy that eliminates the Fierz ambiguity encountered in conventional fluctuation diagnostics of the Hubbard model. The resulting analysis shows that spin fluctuations make an important contribution to the pseudogap self-energy, while multi-boson contributions are also necessary for a complete description. In Chap. 7, we analyze the spin-fermion coupling and its role in d-wave pairing. Our calculations show that, for antinodal fermions in the pseudogap regime, their coupling to antiferromagnetic spin fluctuations is strongly enhanced relative to the non-interacting limit. A modified Bethe–Salpeter equation analysis further shows that this enhanced spin–fermion coupling boosts the spin-fluctuation-induced contribution to the d-wave pairing susceptibility.

Jan 1, 20260 citations