We present an efficient exact diagonalization scheme for the extended dynamical mean-field theory and apply it to the extended Hubbard model on the square lattice with nonlocal charge-charge interactions. Our solver reproduces the phase diagram of this approximation with good accuracy. Details on the numerical treatment of the large Hilbert space of the auxiliary Holstein-Anderson impurity problem are provided. Benchmarks with a numerically exact strong-coupling continuous-time quantum Monte Carlo solver show better convergence behavior of the exact diagonalization in the deep insulator. Special attention is given to possible effects due to the discretization of the bosonic bath. We discuss the quality of real axis spectra and address the question of screening in the Mott insulator within extended dynamical mean-field theory.
We present an efficient exact diagonalization scheme for the extended dynamical mean-field theory and apply it to the extended Hubbard model on the square lattice with nonlocal charge-charge interactions. Our solver reproduces the phase diagram of this approximation with good accuracy. Details on the numerical treatment of the large Hilbert space of the auxiliary Holstein-Anderson impurity problem are provided. Benchmarks with a numerically exact strong-coupling continuous-time quantum Monte Carlo solver show better convergence behavior of the exact diagonalization in the deep insulator. Special attention is given to possible effects due to the discretization of the bosonic bath. We discuss the quality of real axis spectra and address the question of screening in the Mott insulator within extended dynamical mean-field theory.
We here present how a self-consistent solution of the dynamical mean field theory equations can be obtained using exact diagonalization of an Anderson impurity model with accuracies comparable to those found using renormalization group or quantum Monte Carlo methods. We show how one can solve a correlated quantum impurity coupled to several hundred uncorrelated bath sites, using a restricted active basis set. The number of bath sites determines the resolution of the obtained spectral function, which consists of peaks roughly spaced by the band width divided by the number of bath sites. The self-consistency cycles are fully performed on the real frequency axis and expressed as numerical stable matrix operations. The same impurity solver has been used on Ligand Field and finite size cluster calculations and is capable of treating involved Hamiltonians including the full rotational invariant Coulomb interaction, spin-orbit coupling and low symmetry crystal-fields. The proposed method allows for the calculation of a variety of correlation functions at little extra cost.
To explore correlated electrons in the presence of local and nonlocal disorder, the Blackman-Esterling-Berk method for averaging over off-diagonal disorder is implemented into dynamical mean-field theory using tensor notation. The impurity model combining disorder and correlations is solved using the recently developed fork tensor-product state solver, which allows one to calculate the single particle spectral functions on the real-frequency axis. In the absence of off-diagonal hopping, we establish exact bounds of the spectral function of the noninteracting Bethe lattice with coordination number Z. In the presence of interaction, the Mott insulating paramagnetic phase of the one-band Hubbard model is computed at zero temperature in alloys with site- and off-diagonal disorder. When the Hubbard U parameter is increased, transitions from an alloy band insulator through a correlated metal into a Mott insulating phase are found to take place.
The dynamical mean-field theory (DMFT) maps a correlated lattice problem onto an impurity problem of a single correlated site coupled to an uncorrelated bath. Most implementations solve the DMFT equations using quantum Monte-Carlo sampling on the imaginary time and frequency (Matsubara) axis. We will here review alternative methods using exact diagonalization, i.e., representing the many-body ground state of the impurity as a sum over Slater determinants and calculating Green’s functions using iterative Lanczos procedures. The advantage being that these methods have no sign problem, can handle involved multi-orbital Hamiltonians (low crystal symmetry, spin-orbit coupling) and – when working completely on the real axis – do not need a mathematically ill-posed analytical continuation. The disadvantage of traditional implementations of exact diagonalization has been the exponential scaling of the calculation problem as a function of number of bath discretization points. In the last part we will review how recent advances in exact diagonalization can evade the exponential barrier thereby increasing the number of bath discretization points to reach the thermodynamic limit.
In this paper we explore the use of an equation of motion decoupling method as an impurity solver to be used in conjunction with the dynamical mean field self-consistency condition for the solution of lattice models. We benchmark the impurity solver against exact diagonalization, and apply the method to study the infinite U Hubbard model, the periodic Anderson model and the pd model. This simple and numerically efficient approach yields the spectra expected for strongly correlated materials, with a quasiparticle peak and a Hubbard band. It works in a large range of parameters, and therefore can be used for the exploration of real materials using the local density approximation and dynamical mean field theory.
The dynamical mean-field theory (DMFT), which is successful in the study of strongly correlated fermions, was recently extended to boson systems [K. Byczuk and D. Vollhardt, Phys. Rev. B 77, 235106 (2008)]. In this paper, we employ the bosonic DMFT to study the Bose-Hubbard model which describes on-site interacting bosons in a lattice. Using exact diagonalization as the impurity solver, we get the DMFT solutions for the Green's function, the occupation density, as well as the condensate fraction on a Bethe lattice. Various phases are identified: the Mott insulator, the Bose-Einstein condensed (BEC) phase, and the normal phase. At finite temperatures, we obtain the crossover between the Mott-type regime and the normal phase, as well as the BEC-to-normal phase transition. Phase diagrams on the $\ensuremath{\mu}/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane and on the $T/U\text{\ensuremath{-}}\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}/U$ plane are produced ($\stackrel{\ifmmode \tilde{}\else \~{}\fi{}}{t}$ is the scaled hopping amplitude). We compare our results with the previous ones, and discuss the implication of these results to experiments.
The development of polynomial cost solvers for correlated quantum impurity models, with controllable errors, is a central challenge in quantum many-body physics, where these models find applications ranging from nanoscience to the dynamical mean-field theory (DMFT). Here, we describe how configuration interaction (CI) approximations to exact diagonalization (ED) may be used as solvers in DMFT. CI approximations retain the main advantages of ED, such as the ability to treat general interactions and off-diagonal hybridizations and to obtain real spectral information, but are of polynomial cost. Furthermore, their errors can be controlled by monitoring the convergence of physical quantities as a function of the CI hierarchy. Using benchmark DMFT applications, such as single-site DMFT of the one-dimensional Hubbard model and $2\ifmmode\times\else\texttimes\fi{}2$ cluster DMFT of the two-dimensional Hubbard model, we show that CI approximations allow us to obtain near-exact ED results for a tiny fraction of the cost. This is true over the entire range of interaction strengths including ``difficult'' regimes, such as in the pseudogap phase of the two-dimensional Hubbard model. We use the ability of CI approximations to treat large numbers of orbitals to demonstrate convergence of the bath representation in the $2\ifmmode\times\else\texttimes\fi{}2$ cluster DMFT using a 24-bath orbital representation. CI approximations thus form a promising route to extend ED to problems that are currently difficult to study using other solvers such as continuous-time quantum Monte Carlo, including impurity models with large numbers of orbitals and general interactions.