The quantitative control of the dynamic correlations of single impurity Anderson models is essential in several very active fields. We analyze the one-particle Green function with a constant energy resolution by dynamic density-matrix renormalization. In contrast to other approaches, sharp dominant resonances at high energies are found. Their origin and importance are discussed.
The quantitative control of the dynamic correlations of single impurity Anderson models is essential in several very active fields. We analyze the one-particle Green function with a constant energy resolution by dynamic density-matrix renormalization. In contrast to other approaches, sharp dominant resonances at high energies are found. Their origin and importance are discussed.
A dynamic density-matrix renormalization group approach to the spectral properties of quantum impurity problems is presented. The method is demonstrated on the spectral density of the flat-band symmetric single-impurity Anderson model. We show that this approach provides the impurity spectral density for all frequencies and coupling strengths. In particular, Hubbard satellites at high energy can be obtained with a good resolution. The main difficulties are the necessary discretization of the host band hybridized with the impurity and the resolution of sharp spectral features such as the Abrikosov–Suhl resonance.
The impurity Green function and dynamical susceptibilities for the two-channel Anderson impurity model are calculated. An exact expression for the self-energy of the impurity Green function as ratio of two correlation functions is given. The imaginary part of the self-energy scales as $\sqrt{\ensuremath{\mid}\ensuremath{\omega}∕{T}_{K}\ensuremath{\mid}}$ for $T\ensuremath{\rightarrow}0$, serving as a hallmark for non-Fermi behavior. The many-body resonance is pinned to a universal value of $1∕(2\ensuremath{\pi}\ensuremath{\Delta})$ at $\ensuremath{\omega}=0$ for arbitrary local occupation in contrary to the single-channel Anderson impurity model. Its shape becomes increasingly more symmetric for the Kondo regimes of the model. The dynamical spin and channel susceptibilities are governed by two energy scales TK and Th, and their ${\ensuremath{\chi}}^{\ensuremath{''}}(\ensuremath{\omega})$ approach a constant value for $\ensuremath{\omega}\ensuremath{\rightarrow}0$.
An essential role for the description of correlated quantum many-particle systems is played by impurity problems. These consist of a small number of localized orbitals, where the electrons are subject to a Coulomb-interaction, coupled to noninteracting bath degrees of freedom. Impurity models are used for the description of magnetic impurities in metals, of quantum points in nanostructures, and in the context of the dynamical mean field theory (DMFT). For their solution a large number of numerical methods exists. We mention here quantum Monte Carlo methods, the numerical renormalization group, or the exact diagonalization method. In this thesis we introduce a new approach to impurity problems, which is based on the functional renormalization group (fRG). In contrast to the methods mentioned above this approach is approximative, but it can be applied with comparatively lower numerical effort. The aim of this thesis is to examine the advantages and limits of this method. An important role in this connection is played by the single impurity Anderson model (SIAM), which is already understood very well with other methods and is seen as a prototype for models with strong local correlations. The spin fluctuations of this model are governed by the so-called Kondo energy scale and a central question of this thesis is whether this scale can be reproduced with our approach. We begin with a short introduction to the Green's function formalism for correlated quantum many particle systems. Then the fRG flow equations for the one-particle irreducible vertex functions are deduced, and several approximation schemes are discussed. In chapter 4 we introduce the fRG scheme used in this thesis. At first the SIAM is mapped to a semi-infinite chain, in which the interacting orbital is given by the first site of this chain. The system is then subdivided into two parts. The first part, which is called 'core', consists of the interacting orbital and the first L bath sites. The remaining bath sites form the second part. At the beginning both parts are decoupled from each other such that the core can be solved exactly. Starting with this exact solution the coupling between the core and the remaining bath sites is switched on slowly in a renormalization group flow. This way it is possible to calculate the one-particle and two-particle correlation functions on the correlated site. We call this flow scheme 'hybridization flow'. The flow equations are formulated in an effective theory on the first bath site outside the core, which turns out to be advantageous compared to other implementations of the hybridization flow. In chapter 5 this hybridization flow method is applied to the SIAM with semi-elliptic bath density of states. We discuss the differences which arise for different core-sizes (L = 0,1,2,3) and for different truncations of the fRG flow equations. In the cases L=0 and L=2 the local Fermi liquid properties of the SIAM and the associated Kondo scale can not be reproduced in both approximation schemes. For the core sizes L=1 and L=3 the ground state of the core is already a spin singlet state in the beginning of the flow and the Fermi liquid properties are reproduced. The Kondo energy scale, which determines the width of the central quasi-particle peak and the size of the local spin-susceptibility, is however not accurately resolved. This is partly due to the finite temperature used in our implementation such that for larger interaction strength the temperature is above the Kondo energy scale. The application of the hybridization flow method to quantum impurity problems in the context of the DMFT is discussed in chapter 6. In this approach lattice models with a local interaction, like the Hubbard model, are mapped to an effective impurity problem for a single lattice site embedded in a dynamical mean field representing the influence of the other electrons. This scheme is exact in infinitely many dimensions and a nonperturbative approximation method for finite-dimensional systems. In particular the Mott metal-insulator transition in the Hubbard model is successfully described by this method. One class of quantities, which is easily accessible from the hybridization flow scheme, is given by the local two-particle correlation functions. These play an essential role in non-local extensions of the DMFT. Therefore we mainly focus on these quantities. At first we apply our flow scheme to the Hubbard model on the Bethe lattice in infinite dimensions and calculate the local one-particle irreducible vertex functions in the insulating and the metallic phase in good agreement with previous calculations, which use the exact diagonalization method as impurity solver. After this we extend our flow scheme to a cluster DMFT method, which includes short-range antiferromagnetic correlations and calculate again the local and non-local vertex functions for the two-dimensional Hubbard model. We compare these to the results of the single-site DMFT.
The self-energy method for quantum impurity models expresses the correlation part of the self-energy in terms of the ratio of two Green's functions and allows for a more accurate calculation of equilibrium spectral functions than is possible directly from the one-particle Green's function [Bulla et al., J. Phys.: Condens. Matter 10, 8365 (1998)], for example, within the numerical renormalization group method. In addition, the self-energy itself is a central quantity required in the dynamical mean field theory of strongly correlated lattice models. Here, we show how to generalize the self-energy method to the time-dependent situation for the prototype model of strong correlations, the Anderson impurity model. We use the equation-of-motion method to obtain closed expressions for the local Green's function in terms of a time-dependent correlation self-energy, with the latter being given as a ratio of a one-particle time-dependent Green's function and a higher-order correlation function. We benchmark this self-energy approach to time-dependent spectral functions against the direct approach within the time-dependent numerical renormalization group method. The self-energy approach improves the accuracy of time-dependent spectral function calculations, and the closed-form expressions for the Green's function allow for a clear picture of the time evolution of spectral features at the different characteristic time scales. The self-energy approach is of potential interest also for other quantum impurity solvers for real-time evolution, including time-dependent density matrix renormalization group and continuous-time quantum Monte Carlo techniques.
The paper addresses the many-body quasiparticle dynamics of the Anderson impurity model at finite temperatures in the framework of the equation-of-motion method. We find a new exact identity relating the one-particle and many-particle Green’s functions. Using this identity we present a consistent and general scheme for a construction of generalized mean fields (elastic scattering corrections) and self-energy (inelastic scattering) in terms of Dyson equation. A new approach for the complex expansion for the single-particle propagator in terms of Coulomb repulsion U and hybridization V is proposed. Using the exact identity, the essentially new many-body dynamical solution of SIAM has been derived. This approach offers a new way for the systematic construction of the approximative interpolating dynamical solutions of the strongly correlated electron systems.
In this thesis we investigated, implemented, gauged, optimized, and applied a numerical approach given by the dynamic density-matrix renormalization (D-DMRG). The dynamic density-matrix renormalization is the extension of standard DMRG to the calculation of dynamic quantities. D-DMRG provides valuable numerical information on dynamic correlations by computing convolutions of the corresponding spectral densities. The dynamics at zero temperature is determined by computing the expectation values in the local propagator. This can be realized by targeting not only at the ground state and the excited state, but also at the resolvent applied to the excited state. This additional targeted state is called the correction vector. The calculation of the frequency-dependent correction vector is numerically extremely demanding due to the inversion of an almost singular non-hermitian matrix. We compared the performance of several iterative solvers for linear equation systems and managed to stabilize the inversion problem of the D-DMRG by using optimized algorithms. The main limitation of the correction vector D-DMRG is that one cannot obtain data for purely real frequencies but only for frequencies with a certain imaginary part. Hence the extraction of the behavior at purely real frequencies is one of the main problems to be solved in using the D-DMRG. We illustrated how and to which extent such data can be deconvolved to retrieve the wanted spectral densities. We discussed and compared various algorithms to achieve this extraction and presented a non-linear approach from the family of maximum entropy methods. This approach provides a continuous, positive ansatz for the wanted spectral density with the least bias (LB). Even relative abrupt changes of the spectral density can be reproduced satisfactorily. In the vicinity of singularities spurious oscillations occur. The least-bias ansatz can be made more robust towards small numerical inaccuracies and finite-size effects by including besides the entropy functional a xi-functional in the functional to be minimized. We investigated the dynamic propagator of the symmetric single impurity Anderson model by D-DMRG. The Abrikosov-Suhl resonance (ASR) was calculated to gauge the D-DMRG and to demonstrate that features at low energies can be resolved. The low energy scale of the SIAM was reproduced over two orders of magnitude. The Hubbard satellites were investigated. From the broadened D-DMRG raw data we extracted the width and the shifts from the atomic positions. By means of the LB extraction we were able to address the line shape of the satellites. It was found that the Hubbard satellites are strongly asymmetric, the peak is very pronounced and the maximum value is very high. The Hubbard model on the Bethe lattice with infinite coordination number was investigated by means of the dynamic mean-field theory (DMFT) using the D-DMRG as impurity solver. A high-resolution investigation of the electron spectra close to the metal-to-insulator transition in dynamic mean-field theory was presented. The all-numerical, consistent confirmation of a smooth transition at zero temperature is provided.