In this thesis we deal with the description of quantum impurity systems. In order to achieve such an aim, we exploit in the first part of this work the Reduced Density Matrix (RDM) formalism, in the second one the Density Matrix-Numerical Renormalization Group (DM-NRG) method. The first half of the thesis deals with the derivation and the application of out-of-equilibrium and non-perturbative approximations for the dynamics of the Single Impurity Anderson Model (SIAM) in the perturbative, intermediate coupling and Kondo cross-over regimes. In this part of the thesis, the main aim is to develop new technical tools useful to extend the RDM formalism from the well established perturbative regime to the intermediate and, in the future, to the strong coupling regimes. In this part, we analyse the advantages and drawbacks of the so-called Dressed Second Order (DSO), Dressed DSO (DDSO) and the Resonant Tunnelling Approximation (RTA). In the DSO we include processes capable to describe the broadening of the Coulomb blockade peaks but the DSO itself, as the RTA, is not capable to describe the Kondo resonance width properly. In the DDSO we propose a prescription to solve this issue and to describe properly the exponential dependence of the Kondo resonance width. However, in the DDSO the vertex corrections are not included and this yields a spurious conductance peak in the empty orbital regime (Coulomb valley with zero electrons). For this reason we turn our attention to the inclusion of the vertex corrections within the RTA framework. In the second part, we describe transport and thermodynamic properties of carbon-nanotubes (CNTs) based quantum dots using the DM-NRG method. We present results on the spectral function and the linear conductance in cooperation with Catalin Pascu Moca and Gergely Zarand from the Universities of Oradea and Budapest, respectively. In order to model the CNTs we consider an extended Anderson Hamiltonian. Exploiting such an Hamiltonian we study the crossover from the SU(4) to the SU(2) Kondo effect. We demonstrate numerically and analytically that the Kondo temperature scales inversely with respect to the SU(4)-symmetry breaking parameter $\Delta$, tending to the Kondo temperature of the SU(2) Kondo effect in the large $\Delta$ limit. Furthermore we study the CNT quantum dot in the presence of a parallel and perpendicular (with respect to the CNT axis) magnetic field. We find necessary conditions to observe the so-called ``Kondo revivals'' for specific values of the magnetic field. In the final part of this thesis, we model a real world device realised by Jean-Pierre Cleuziou, Ngoc-Viet Nguyen and Wolfgang Wernsdorfer in Grenoble by means of the previously studied Hamiltonian. For this system we search for the best set of parameters to reproduce the experimentally measured linear conductance and we study the underlying Kondo state evaluating the CNT added entropy and specific heat.
In this thesis we deal with the description of quantum impurity systems. In order to achieve such an aim, we exploit in the first part of this work the Reduced Density Matrix (RDM) formalism, in the second one the Density Matrix-Numerical Renormalization Group (DM-NRG) method. The first half of the thesis deals with the derivation and the application of out-of-equilibrium and non-perturbative approximations for the dynamics of the Single Impurity Anderson Model (SIAM) in the perturbative, intermediate coupling and Kondo cross-over regimes. In this part of the thesis, the main aim is to develop new technical tools useful to extend the RDM formalism from the well established perturbative regime to the intermediate and, in the future, to the strong coupling regimes. In this part, we analyse the advantages and drawbacks of the so-called Dressed Second Order (DSO), Dressed DSO (DDSO) and the Resonant Tunnelling Approximation (RTA). In the DSO we include processes capable to describe the broadening of the Coulomb blockade peaks but the DSO itself, as the RTA, is not capable to describe the Kondo resonance width properly. In the DDSO we propose a prescription to solve this issue and to describe properly the exponential dependence of the Kondo resonance width. However, in the DDSO the vertex corrections are not included and this yields a spurious conductance peak in the empty orbital regime (Coulomb valley with zero electrons). For this reason we turn our attention to the inclusion of the vertex corrections within the RTA framework. In the second part, we describe transport and thermodynamic properties of carbon-nanotubes (CNTs) based quantum dots using the DM-NRG method. We present results on the spectral function and the linear conductance in cooperation with Catalin Pascu Moca and Gergely Zarand from the Universities of Oradea and Budapest, respectively. In order to model the CNTs we consider an extended Anderson Hamiltonian. Exploiting such an Hamiltonian we study the crossover from the SU(4) to the SU(2) Kondo effect. We demonstrate numerically and analytically that the Kondo temperature scales inversely with respect to the SU(4)-symmetry breaking parameter Δ, tending to the Kondo temperature of the SU(2) Kondo effect in the large Δ limit. Furthermore we study the CNT quantum dot in the presence of a parallel and perpendicular (with respect to the CNT axis) magnetic field. We find necessary conditions to observe the so-called ``Kondo revivals'' for specific values of the magnetic field. In the final part of this thesis, we model a real world device realised by Jean-Pierre Cleuziou, Ngoc-Viet Nguyen and Wolfgang Wernsdorfer in Grenoble by means of the previously studied Hamiltonian. For this system we search for the best set of parameters to reproduce the experimentally measured linear conductance and we study the underlying Kondo state evaluating the CNT added entropy and specific heat.
The focus of this thesis is set on the calculation of spectral functions for low-dimensional quantum systems. At first the quantum many-body problem is introduced and the related spectral functions of interest are defined. In a next step the applied numerical algorithms are presented: Exact diagonalization & Lanczos algorithm, the numerical renormalization group (NRG) and the density matrix renormalization group (DMRG) in the formulation of matrix product states (MPS). The main focus is set on the calculation of spectral functions with these algorithms. In particular for the DMRG a new algorithm is presented which combines the Lanczos algorithm with matrix product states. The presented algorithms are applied to different problems. The first system denotes a copper crystal with magnetic (cobalt/iron) impurities. These materials are known Kondo-systems, i.e. they show an untypical increase in the resistivity with decreasing temperature. At low temperatures local moments will develop, which will be screened from the surrounding conduction band electrons. The ground state of the system is a long-range strongly correlated many-body state between conduction band electrons and local moments. In the singleparticle impurity spectral function a so-called Kondo resonance at the Fermi energy is formed below the Kondo temperature. In scanning tunneling microscope measurements the spectral function at the surface - the so-called local density of states - of the copper crystals with buried impurities in varying depth was measured. The goal of this study is to verify the existence of long-range Kondo signatures via a comparison of the spectral functions for impurities with different depth below the surface. Therefore in this work numerical simulations are used to calculate the local density of states at the surface. The single impurity Anderson model (SIAM) is applied which shows the Kondo resonance in the impurity spectral function. In order to calculate the local density of states at the surface of the copper crystal the equations of motion of the SIAM and the Dyson equation are used and then solved by a combination of numerical renormalization group and band structure calculations of the copper crystal (linear combination of atomic orbitals). The simulation for impurities in different depth below the surface and the experimental data are in very good agreement and hence support the measurement of a long-range Kondo signature. Furthermore the simulations allow to determine the Kondo temperature of a single impurity. The focus of the second part is set on the one-dimensional Heisenberg model. The spectral function of interest is here the dynamic spin structure factor. The newly developed MPS-Lanczos algorithm within the DMRG is used to calculate the dynamic spin structure factor. The results are compared to analytic results from the Bethe ansatz. In the end the new algorithm is discussed in respect to existing methods within the DMRG to calculate spectral functions.
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 first topic of this thesis is the study of many-body effects in an one-dimensional strongly correlated electronic system - the Kondo lattice model. This system is tackled numerically by means of the density matrix renormalization group, since analytic method, i.e., perturbation the- ory fail due to competing coupling constants. The Kondo lattice model consists of a conduction band of electrons which couple via a spin exchange coupling to a localized spin lattice. We study the spectral properties of the one-dimensional Kondo lattice model as a function of the exchange coupling, the band filling, and the quasimomentum in the ferromagnetic and paramagnetic phases. We compute the dispersion relation of the quasiparticles, their lifetimes, and the Z factor. The exact ground state and the quasiparticle-dispersion relation of the Kondo lattice model with one conduction electron are well known. The quasiparticle could be identified as the spin polaron. Our calculations of the dispersion relation for partial band fillings give a result similar to the one-electron case, which suggests that the quasiparticle in both cases is the spin polaron. We find that the quasiparticle lifetime differs by orders of magnitude between the ferromagnetic and paramagnetic phases and depends strongly on the quasimomentum. Further- more, we study the effects of the Coulomb interaction on the phase diagram, the static magnetic susceptibility and electron spin relaxation. We show that onsite Coulomb interaction supports ferromagnetic order and nearest neighbor Coulomb interaction drives, depending on the elec- tron filling, either a paramagnetic or ferromagnetic order. Furthermore, we calculate electron quasiparticle life times, which can be related to electron spin relaxation and decoherence times, and explain their dependence on the strength of interactions and the electron filling in order to find the sweet spot of parameters where the relaxation time is maximized. We find that effective exchange processes between the electrons dominate the spin relaxation and decoherence rate. In the second topic of this thesis, we numerically calculate the electron transport through carbon nanotube based quantum dot devices. We use a master equation’s approach in first order of the tunneling rate to the leads and an extended constant interaction model to model the carbon nanotube system. This work has been done in collaboration with two experimental groups and we compare their respective experimentally obtained data to our numerical calculations. In both collaborations striking similarity between the numerical data and the experimental data is found. In the first collaboration transport through a carbon nanotube peapod, i.e, a carbon nanotube filled with fullerenes, has been measured. We identify a small hybridization between a fullerene molecule and the surrounding carbon nanotube to be of crucial importance for the understanding of the transport data. In the second collaboration, electron transport through a carbon nanotube rope, i.e., a bundle of carbon nanotubes has been measured. Also here, hybridization between the different nanotubes plays a crucial role. Furthermore, an external magnetic field is applied, which enables the identification of specific spin states of the compound quantum dot system. This might be important for future applications of such devices in spin-dependent electronics.
This thesis contributes to the field of transport through quantum dots. These devices allow for a controlled study of quantum transport and fundamental physical effects, like the Kondo effect. In this thesis we will focus on dots that are well described by generalized Anderson impurity models, where the discrete levels of the quantum dot are tunnel-coupled to fermionic reservoirs. The model parameters, like level energy and width, can be tuned in experiments. Therefore these systems constitute a valuable arena for testing experiment against theory and vice versa. In order to describe these strongly correlated systems, we employ the numerical renormalization group method. This allows us to address both longstanding questions concerning experimental results and new physical phenomena in these fundamental models. This thesis consists of three major projects. The first and most extensive one is concerned with the phase of the transmission amplitude through a quantum dot. Measurements of many-electron quantum dots with small level spacing reveal universal phase behaviour, a result not fully understood for almost 10 years. Recent experiments have seen that, contrarily, for dots with only a few electrons, i.e. large level spacing, the phase depends on the mesoscopic dot parameters. Analyzing a multi-level Anderson model, we show that the generic feature of the two regimes can be reproduced in the regime of overlapping levels or well separated levels, respectively. Thereby the universal character follows from Fano-type antiresonances of the renormalized single-particle levels. Moderate temperature supports the universal character. In the mesoscopic regime, we also investigate the effect of Kondo correlations on the transmission phase. In a second project we analyze a quantum dot coupled to a superconducting reservoir. In contrast to previous belief, the energy resolution of our method is not restricted by the energy scale of the superconducting gap, leading to new insights into the method. The high resolution allows us to resolve sharp peaks in the spectral function that emerge for a certain regime of parameters. A third project deals with a quantum dot coupled to two independent channels, a system known to exhibit non-Fermi liquid behaviour. We investigate the existence of the non-Fermi liquid regime when driving the system out of the Kondo regime by emptying the dot. We find that the extent of the non-Fermi liquid regime strongly depends on the mechanisms that couple impurity and reservoirs but prevent mixing of the latter.
This thesis contributes to the understanding of impurity models. It is divided into two main parts, with a general introduction given in Part I and the research related to it presented in Part II, with the second part being subdivided into two main projects. In the first project, the influence of two many-body effects, the Kondo effect and the Fermi edge singularity, on the absorption and emission spectra of self-assembled quantum dots (QDs) is examined. Whereas the Kondo effect so far was always examined with transport experiments, we show that it has been observed with optical methods for the first time, by comparing experimental data for the absorption line shapes of QDs to calculations with the numerical renormalization group. We continue by examining a QD with strong optical coupling of the energy levels. The resulting interplay of Rabi-oscillations and Kondo effect leads to a new many-body state, a secondary, outer Kondo effect, with Kondo-like correlations between the spin-Kondo and the trion state. The last work regarding optics at QDs addresses the Fermi edge singularity. We show that for QDs this phenomenon can be described numerically on a quantitative level. The second project concerns transport properties of impurity models. First, we present a comprehensive study of the Kondo effect for an InAs-nanowire QD, a system for which the Kondo effect was observed only a few years ago. The second study regarding transport concerns the Kondo effect in bulk metals with magnetic impurities. Although nowadays the Kondo effect is often studied with QDs, it was discovered for iron impurities in noble metals like gold and silver. However, it was unknown for a long time which exact realization of Kondo model describes these systems. We identify the model by comparing numerical calculations for the magnetoresistivity and the dephasing rate for different models to experimental results. The third work about transport concerns the phenomenon that for a fixed type of Kondo model quantities like the magnetoresistivity or the conductivity, respectively, can be scaled onto a universal curve for different parameters, when energies are rescaled with the the Kondo temperature TK, since it is the only relevant low energy scale of the problem. For finite bandwidth, however, different definitions of TK (which coincide in the limit of infinite bandwidth) lead to different TK-values. We show that with a very common definition of TK, finite bandwidth, which is always present at numerical calculations, can deteriorate the universality of rescaled curves, and we offer an alternative definition of TK which ensures proper scaling. In the last study presented in this thesis we calculate the Fermi-liquid coefficients for fully screened multi-channel Kondo models. For temperatures below TK, these models show Fermi-liquid behavior, and the impurity density of states and certain quantities which depend on it, like resistivity, show quadratic dependencies on parameters like temperature or magnetic field, described by the Fermi-liquid coefficients. We calculate these coefficients both analytically and numerically.
In this thesis, we describe methods for solving quantum impurity models out of equilibrium, and consider several applications of these methods. Quantum impurity models describe systems in which a small subsystem (the impurity) is coupled to a large external environment (the bath). These models are useful in describing a wide variety of physical systems, including quantum dots, absorbed adatoms, and molecular junctions. They are also important for the study of lattice models, such as the Hubbard model, where they appear as auxiliary problems within the framework of dynamical mean field theory (DMFT). Many interesting phenomena in these systems, such as transport through a quantum dot, or optical pumping of a Mott insulator, occur away from equilibrium. Theoretical description of these phenomena requires methods that are capable of handling both real time evolution and strong correlation effects presenting a formidable challenge. In this thesis, we describe several methods for the solution of non-equilibrium quantum impurity models which meet this challenge. These methods are based on an expansion in the coupling between the impurity and the bath known as the hybridization expansion. Chapter 1 gives an introduction to the problems addressed in this thesis. Chapter 2 provides a brief overview of Keldysh formalism for nonequilibrium problems. Chapter 3 describes the bare hybridization expansion for the Anderson impurity model as well as its bold extensions. In particular we derive the non-crossing approximation (NCA), the one-crossing approximation (OCA), the inchworm expansion, and a bold expansion for the impurity Green's function. Chapter 4 describes the implementation of a quantum Monte Carlo algorithm based on the inchworm expansion. The rest of the thesis describes applications of these algorithms to problems of physical interest. Chapters 5 and 6 apply the methods previously described to understanding the Kondo effect in quantum dot systems. The first chapter examines the behavior of the Kondo effect under an applied voltage bias, a classic non-equilibrium strongly-correlated problem. The second chapter examines different ways of measuring the cloud of electrons which screen the impurity spin in the Kondo regime. Chapters 7 and 8 apply the same methods to the Hubbard model using the Dynamical mean-field theory (DMFT) mapping. The first chapter studies the Hubbard model in equilibrium using the numerically exact inchworm quantum Monte Carlo method. The second chapter studies the possibility of dynamic switching between the metal and Mott-insulator states of the Hubbard model using the one-crossing approximation (OCA). Finally chapter 9 gives some overall conclusions and suggests future directions for the development of hybridization expansion methods.