A Numerical Renormalization Group approach to Green's Functions for Quantum Impurity Models

A new numerical method for the solution of the Dynamical Mean Field Theory's self-consistent equations is introduced. The method uses the Density Matrix Renormalization Group technique to solve the associated impurity problem. The new algorithm makes no a priori approximations and is only limited by the number of sites that can be considered. We obtain accurate estimates of the critical values of the metal-insulator transitions and provide evidence of substructure in the Hubb...

The many-body problem is usually approached from one of two perspectives: the first originates from an action and is based on Feynman diagrams, the second is centered around a Hamiltonian and deals with quantum states and operators. The connection between results obtained in either way is made through spectral (or Lehmann) representations, well known for two-point correlation functions. Here, we complete this picture by deriving generalized spectral representations for multip...

We show that the Renormalization Group formalism allows to compute with accuracy the zero temperature correlation functions and particle densities of quantum systems.

We propose a versatile strategy for numerical renormalization group solution of general channel-mixing Kondo and Anderson models beyond previous reach, opening the door toward broad applications in protocol non-perturbative machineries, such as dynamical cluster approximation and cluster dynamical mean field theory, for strongly correlated electron systems. We illustrate the strategy by investigating the quantum phase transitions in two quantum impurity models with cases unto...

Quantum impurity problems can be solved using the numerical renormalization group (NRG), which involves discretizing the free conduction electron system and mapping to a `Wilson chain'. It was shown recently that Wilson chains for different electronic species can be interleaved by use of a modified discretization, dramatically increasing the numerical efficiency of the RG scheme [Phys. Rev. B 89, 121105(R) (2014)]. Here we systematically examine the accuracy and efficiency of...

We present a new method for extracting numerically exact imaginary-time Green functions from standard Hirsch-Fye quantum Monte Carlo (HF-QMC) simulations within dynamical mean-field theory. By analytic continuation, angular resolved spectra are obtained without the discretization bias previously associated with HF-QMC results. The method is shown to be accurate even at very low temperatures (T=W/800 for bandwidth W) in the strongly correlated regime.

We propose an auxiliary-bath algorithm for the numerical renormalization group (NRG) method to solve multi-impurity models with shared electron baths. The method allows us to disentangle the electron baths into independent Wilson chains to perform standard NRG procedures beyond the widely adopted independent bath approximation. Its application to the 2-impurity model immediately reproduces the well-known even- and odd-parity channels. For 3-impurity Kondo models, we find succ...

Numerically exact continuous-time Quantum Monte Carlo algorithm for finite fermionic systems with non-local interactions is proposed. The scheme is particularly applicable for general multi-band time-dependent correlations since it does not invoke Hubbard-Stratonovich transformation. The present determinantal grand-canonical method is based on a stochastic series expansion for the partition function in the interaction representation. The results for the Green function and for...

We present a quantum Monte-Carlo algorithm for computing the perturbative expansion in power of the coupling constant $U$ of the out-of-equilibrium Green's functions of interacting Hamiltonians of fermions. The algorithm extends the one presented in Phys. Rev. B 91 245154 (2015), and inherits its main property: it can reach the infinite time (steady state) limit since the computational cost to compute order $U^n$ is uniform versus time; the computing time increases as $2^n$. ...

We present a numerically exact steady-state inchworm Monte Carlo method for nonequilibrium quantum impurity models. Rather than propagating an initial state to long times, the method is directly formulated in the steady-state. This eliminates any need to traverse the transient dynamics and grants access to a much larger range of parameter regimes at vastly reduced computational costs. We benchmark the method on equilibrium Green's functions of quantum dots in the noninteracti...