We study the single-impurity Anderson model by means of cluster perturbation theory and the variational cluster approach (VCA). An expression for the VCA grand potential for a system in a noninteracting bath is presented. Results for the single-particle dynamics in different parameter regimes are shown to be in good agreement with established renormalization group results. We aim at a broad and comprehensive overview of the capabilities and shortcomings of the methods. We address the question as to what extent the elusive low-energy properties of the model are reproducible within the framework of VCA. These are furthermore benchmarked against continuous-time quantum Monte Carlo calculations. We also discuss results obtained by an alternative, i.e., self-consistent formulation of VCA, which was introduced recently in the context of nonequilibrium systems.
We study the single-impurity Anderson model by means of cluster perturbation theory and the variational cluster approach (VCA). An expression for the VCA grand potential for a system in a noninteracting bath is presented. Results for the single-particle dynamics in different parameter regimes are shown to be in good agreement with established renormalization group results. We aim at a broad and comprehensive overview of the capabilities and shortcomings of the methods. We address the question as to what extent the elusive low-energy properties of the model are reproducible within the framework of VCA. These are furthermore benchmarked against continuous-time quantum Monte Carlo calculations. We also discuss results obtained by an alternative, i.e., self-consistent formulation of VCA, which was introduced recently in the context of nonequilibrium systems.
We study the single impurity Anderson model by means of the variational cluster approach (VCA). An expression for the VCA grand potential for a system in a noninteracting bath is presented. Results for dynamic correlation functions in different parameter regimes are shown to be in good agreement with renormalization group results. We address the question whether the elusive low energy properties of the model are properly reproduced within the framework of VCA. These are furthermore compared to continuous time quantum Monte Carlo calculations. We also discuss results obtained by an alternative, i. e. self-consistent formulation of VCA, which was introduced recently in the context of non-equilibrium systems.
The variational cluster approach (VCA) is applied to the one-dimensional Hubbard model at zero temperature using clusters (chains) of up to ten sites with full diagonalization and the Lanczos method as cluster solver. Within the framework of the self-energy-functional theory (SFT), different cluster reference systems with and without bath degrees of freedom, in different topologies and with different sets of variational parameters, are considered. Static and one-particle dynamical quantities are calculated for half-filling as a function of U as well as for fixed U as a function of the chemical potential to study the interaction- and filling-dependent metal-insulator (Mott) transition. The recently developed Q-matrix technique is used to compute the SFT grand potential. For benchmarking purposes we compare the VCA results with exact results available from the Bethe ansatz, with essentially exact dynamical density-matrix renormalization-group data, with (cellular) dynamical mean-field theory and full diagonalization of isolated Hubbard chains. Several issues are discussed including convergence of the results with cluster size, the ability of cluster approaches to access the critical regime of the Mott transition, efficiency in the optimization of correlated-site vs bath-site parameters, and of multidimensional parameter optimization. We also study the role of bath sites for the description of excitation properties and as charge reservoirs for the description of filling dependencies. The VCA turns out to be a computationally cheap method which is competitive with established cluster approaches.
The variational cluster approximation (VCA) proposed by M. Potthoff et al. [Phys. Rev. Lett. 91, 206402 (2003)] is extended to electron or spin systems with nonlocal interactions. By introducing more than one source field in the action and employing the Legendre transformation, we derive a generalized self-energy functional with stationary properties. Applying this functional to a proper reference system, we construct the extended VCA (EVCA). In the limit of continuous degrees of freedom for the reference system, EVCA can recover the cluster extension of the extended dynamical mean-field theory (EDMFT). For a system with correlated hopping, the EVCA recovers the cluster extension of the dynamical mean-field theory for correlated hopping. Using a discrete reference system composed of decoupled three-site single impurities, we test the theory for the extended Hubbard model. Quantitatively good results as compared with EDMFT are obtained. We also propose VCA (EVCA) based on clusters with periodic boundary conditions. It has the (extended) dynamical cluster approximation as the continuous limit. A number of related issues are discussed.
The variational cluster approach (VCA) is applied to study spontaneous ferromagnetism in the Hubbard model at zero temperature. We discuss several technical improvements of the numerical implementation of the VCA which become necessary for studies of a ferromagnetically ordered phase, e.g., more accurate techniques to evaluate the variational ground-state energy, improved local as well as global algorithms to find stationary points, and different methods to locate the magnetic phase transition. Using the single-site VCA, i.e., the dynamical impurity approximation (DIA), the ferromagnetic phase diagram of the model in infinite dimensions is worked out. The results are compared with previous dynamical mean-field studies for benchmarking purposes. The DIA results provide a unified picture of ferromagnetism in the infinite-dimensional model by interlinking different parameter regimes that are governed by different mechanisms for ferromagnetic order. Using the DIA and the VCA, we then study ferromagnetism in one-dimensional Hubbard chains with nearest- and next-nearest-neighbor hopping t2. In comparison with previous results from the density-matrix renormalization group, the phase diagram is mapped out as a function of the Hubbard-U, the electron filling, and t2. The stability of the ferromagnetic ground state against local and short-range nonlocal quantum fluctuations is discussed.
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 Anderson Impurity Model is a generic model to describe systems in which one or more energy levels are coupled to an electronic host. On the energy levels the electrons are subject to a Coulomb interaction. Possible applications of this model are magnetic atoms in a metal or quantum dot systems. Already for one impurity it displays very interesting physics including the well known Kondo effect. For two and more impurities the properties of the system in the local moment regime are determined by an interplay of the Kondo effect and the coupling between the impurity spins. Due to the interplay of multiple energy scales, this model poses a serious task for any computational tool. For a single impurity the model can be very efficiently treated using Wilson"s numerical renormalization group (NRG). Since the numerical effort of NRG calculations rises exponentially with the number of impurities, generically it cannot be applied to systems with more than two impurities. The functional renormalization group (fRG) has been very successfully applied to quantum impurity problems in the past. In this work it will be used to calculate the single-particle dynamics for Anderson-like quantum impurity systems in equilibrium. A comparison to NRG in the case of a single impurity will serve as a benchmark before systems with two impurities and quantum dot systems with two and three dots will be studied. The interplay between the Kondo effect and magnetic and orbital ordering will be investigated in these systems. Although the truncation in the presented fRG scheme is motivated by second order perturbation theory, it will be shown that fRG works also in those cases where perturbation theory fails. Generally it will turn out that fRG is capable of treating a broad variety of Anderson impurity systems within the same framework, including asymmetric impurity levels and interaction between multiple impurities.