We present a very efficient solver for the general Anderson impurity problem. It is based on the perturbation around a solution obtained from exact diagonalization using a small number of bath sites. We formulate a perturbation theory which is valid for both weak and strong coupling and interpolates between these limits. Good agreement with numerically exact quantum Monte-Carlo results is found for a single bath site over a wide range of parameters. In particular, the Kondo resonance in the intermediate coupling regime is well reproduced for a single bath site and the lowest order correction. The method is particularly suited for low temperatures and alleviates analytical continuation of imaginary time data due to the absence of statistical noise compared to quantum Monte-Carlo impurity solvers.
Quantum impurity models play an important role in many areas of physics from condensed matter to AMO and quantum information. They are important models for many physical systems but also provide key insights to understanding much more complicated scenarios. In this paper we introduce a simplified method to describe the thermodynamic properties of integrable quantum impurity models. We show this method explicitly using the anisotropic Kondo and the interacting resonant level models. We derive a simplified expression for the free energy of both models in terms of a single physically transparent integral equation which is valid at all temperatures and values of the coupling constants.
We propose a distinct numerical approach to effectively solve the problem of partial diagonalization of the super-large-scale quantum electronic Hamiltonian matrices. The key ingredients of our scheme are the new method for arranging the basis vectors in the computers RAM and the algorithm allowing not to store a matrix in RAM, but to regenerate it on-the-fly during diagonalization procedure. This scheme was implemented in the program, solving the Anderson impurity model in the framework of dynamical mean-field theory (DMFT). The DMFT equations for electronic Hamiltonian with 18 effective orbitals that corresponds to the matrix with the dimension of 2.4 * 10^9 were solved on the distributed memory computational cluster.
We derive an exact mapping from the action of nonequilibrium dynamical mean-field theory (DMFT) to a single-impurity Anderson model (SIAM) with time-dependent parameters, which can be solved numerically by exact diagonalization. The representability of the nonequilibrium DMFT action by a SIAM is established as a rather general property of nonequilibrium Green functions. We also obtain the nonequilibrium DMFT equations using the cavity method alone. We show how to numerically obtain the SIAM parameters using Cholesky or eigenvector matrix decompositions. As an application, we use a Krylov-based time propagation method to investigate the Hubbard model in which the hopping is switched on, starting from the atomic limit. Possible future developments are discussed.
We extend a previously proposed rotation and truncation scheme to optimize quantum Anderson impurity calculations with exact diagonalization [PRB 90, 085102 (2014)] to density-matrix renormalization group (DMRG) calculations. The method reduces the solution of a full impurity problem with virtually unlimited bath sites to that of a small subsystem based on a natural impurity orbital basis set. The later is solved by DMRG in combination with a restricted-active-space truncation scheme. The method allows one to compute Greens functions directly on the real frequency or time axis. We critically test the convergence of the truncation scheme using a one-band Hubbard model solved in the dynamical mean-field theory. The projection is exact in the limit of both infinitely large and small Coulomb interactions. For all parameter ranges the accuracy of the projected solution converges exponentially to the exact solution with increasing subsystem size.
We present an implementation of the hybridization expansion impurity solver which employs sparse matrix exact-diagonalization techniques to compute the time evolution of the local Hamiltonian. This method avoids computationally expensive matrix-matrix multiplications and becomes advantageous over the conventional implementation for models with 5 or more orbitals. In particular, this method will allow the systematic investigation of 7-orbital systems (lanthanide and actinide compounds) within single-site dynamical mean field theory. We illustrate the power and usefulness of our approach with dynamical mean field results for a 5-orbital model which captures some aspects of the physics of the iron based superconductors.