The density matrix renormalization group method is introduced in energy space to study Anderson impurity models. The method allows for calculations in the thermodynamic limit and is advantageous for studying not only the dynamical properties, but also the quantum entanglement of the ground state at the vicinity of an impurity quantum phase transition. This method is applied to obtain numerically exactly the ground-state phase diagram of the single-impurity Anderson model on the honeycomb lattice at half-filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase, but no Kondo screening phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model using the method introduced here. Although the high-energy excitations are obviously different, we find that the ground-state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single-impurity Anderson model on the honeycomb lattice, thus providing a quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally threefold degenerate at the valence fluctuating point. This should be contrasted with the degeneracy of the energy spectrum because the ground state is found to be always singlet. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum phase transition.
The density matrix renormalization group method is introduced in energy space to study Anderson impurity models. The method allows for calculations in the thermodynamic limit and is advantageous for studying not only the dynamical properties, but also the quantum entanglement of the ground state at the vicinity of an impurity quantum phase transition. This method is applied to obtain numerically exactly the ground-state phase diagram of the single-impurity Anderson model on the honeycomb lattice at half-filling. The calculation of local static quantities shows that the phase diagram contains two distinct phases, the local moment (LM) phase and the asymmetric strong coupling (ASC) phase, but no Kondo screening phase. These results are supported by the local spin and charge excitation spectra, which exhibit qualitatively different behavior in these two phases and also reveal the existence of the valence fluctuating point at the phase boundary. For comparison, we also study the low-energy effective pseudogap Anderson model using the method introduced here. Although the high-energy excitations are obviously different, we find that the ground-state phase diagram and the asymptotically low-energy excitations are in good quantitative agreement with those for the single-impurity Anderson model on the honeycomb lattice, thus providing a quantitative justification for the previous studies based on low-energy approximate approaches. Furthermore, we find that the lowest entanglement level is doubly degenerate for the LM phase, whereas it is singlet for the ASC phase and is accidentally threefold degenerate at the valence fluctuating point. This should be contrasted with the degeneracy of the energy spectrum because the ground state is found to be always singlet. Our results therefore clearly demonstrate that the low-lying entanglement spectrum can be used to determine with high accuracy the phase boundary of the impurity quantum phase transition.
We introduce a block Lanczos (BL) recursive technique to construct quasi-one-dimensional models, suitable for density-matrix renormalization group (DMRG) calculations, from single- as well as multiple-impurity Anderson models in any spatial dimensions. This new scheme, named BL-DMRG method, allows us to calculate not only local but also spatially dependent static and dynamical quantities of the ground state for general Anderson impurity models without losing elaborate geometrical information of the lattice. We show that the BL-DMRG method can be easily extended to treat a multiorbital Anderson impurity model where not only inter- and intraorbital Coulomb interactions but also Hund's coupling and pair hopping interactions are included. We also show that the symmetry adapted BL bases can be utilized, when it is appropriate, to reduce the computational cost. As a demonstration, we apply the BL-DMRG method to three different models for graphene with a structural defect and with a single hydrogen or fluorine absorbed, where a single Anderson impurity is coupled to conduction electrons in the honeycomb lattice. These models include (i) a single adatom on the honeycomb lattice, (ii) a substitutional impurity in the honeycomb lattice, and (iii) an effective model for a single carbon vacancy in graphene. Our analysis of the local dynamical magnetic susceptibility and the local density of states at the impurity site reveals that, for the particle-hole symmetric case at half-filling of electron density, the ground state of model (i) behaves as an isolated magnetic impurity with no Kondo screening, while the ground state of the other two models forms a spin-singlet state where the impurity moment is screened by the conduction electrons. We also calculate the real-space dependence of the spin-spin correlation functions between the impurity site and the conduction sites for these three models. Our results clearly show that, reflecting the presence or absence of unscreened magnetic moment at the impurity site, the spin-spin correlation functions decay as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}3}$, differently from the noninteracting limit ($\ensuremath{\propto}$ ${r}^{\ensuremath{-}2}$), for model (i) and as $\ensuremath{\propto}$ ${r}^{\ensuremath{-}4}$, exactly the same as the noninteracting limit, for models (ii) and (iii) in the asymptotic r, where r is the distance between the impurity site and the conduction site. Finally, based on our results, we shed light on recent experiments on graphene where the formation of local magnetic moments as well as the Kondo-like behavior have been observed.
We use the numerical renormalization group method to investigate the spectral properties of a single-impurity Anderson model with a gap $\ensuremath{\delta}$ across the Fermi level in the conduction-electron spectrum. For any finite $\ensuremath{\delta}>0$, at half filling the ground state of the system is always a doublet. Away from half filling a quantum phase transition (QPT) occurs as function of the gap value $\ensuremath{\delta}$, and the system evolves from the strong-coupling (SC) Kondo-type state, corresponding to $\ensuremath{\delta}<{\ensuremath{\delta}}_{C}$ toward a localized moment (LM) regime for $\ensuremath{\delta}>{\ensuremath{\delta}}_{C}$. The opening of the gap leads to the formation of one (two) bound states when the system is in the SC (LM) regime. The evolution across the QPT of their positions and the corresponding weights together with the dynamic properties of the model are investigated.
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.
We study the ground state properties of a single magnetic impurity in the honeycomb lattice. To understand the non-perturbative properties, we use the density matrix renormalization group method where the honeycomb lattice model is mapped onto an effective one-dimensional model by means of the Lanczos recursion technique. We apply this method to three different models: (i) a magnetic adatom on the honeycomb lattice, (ii) a substitutional magnetic impurity in the honeycomb lattice, and (iii) a defect in the honeycomb lattice. Here, we focus only on the particle-hole symmetric cases and find that the ground state of model (i) is qualitatively different from the others: the magnetic impurity spin behaves like an isolated spin in model (i) while it is screened in models (ii) and (iii) by conduction electrons to form a Kondo singlet like state.
Quantum impurity models are the prototypical examples of quantum many-body dynamics which manifests in their spectral and transport properties. Single channel Anderson(and Kondo model) leads to the Fermi liquid ground state in the strong coupling regime which corresponds to a stable infrared fixed point at which the quantum impurity gets completely screened by the conduction electrons. Quantum impurity models with non-trivial density of states exhibit quantum phase transition and this quantum criticality lies in the universality class of local quantum critical systems. In this paper, we report first study of the flow equation renormalization of gapped Kondo model which has gapped density of states, the gap being at the Fermi level. Flow equation renormalization group method has proved to be one of the very robust renormalization methods to study Kondo physics both in equilibrium as well as in non-equilibrium. Here we have shown that this method can also be employed to study local quantum criticality. We have calculated the flow equations for the Kondo coupling and solved them for various values of the gap parameter and we find that there is suppression of Kondo divergence as gap is increased which signifies that as gap is increased, renormalization flow goes away from the strong coupling fixed point. We have also calculated the spin susceptibility and we find that as gap is increased, susceptibility goes over to the Curie behaviour and hence confirming the renormalization flow towards the local moment fixed point.
We investigate static and dynamical ground-state properties of the two-impurity Anderson model at half filling in the limit of vanishing impurity separation using the dynamical density-matrix renormalization group method. In the weak-coupling regime, we find a quantum phase transition as a function of inter-impurity hopping driven by the charge degrees of freedom. For large values of the local Coulomb repulsion, the transition is driven instead by a competition between local and non-local magnetic correlations. We find evidence that, in contrast to the usual phenomenological picture, it seems to be the bare effective exchange interactions which trigger the observed transition.