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.
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 apply a recently introduced hybridization flow functional renormalization group scheme for Anderson-like impurity models as an impurity solver in a dynamical mean-field theory (DMFT) approach to lattice Hubbard models. We present how this scheme is capable of reproducing metallic and insulating solutions of the lattice model. Our setup also offers a numerically rather inexpensive method to calculate two-particle correlation functions. For the paramagnetic Hubbard model on the Bethe lattice in infinite dimensions we calculate the local two-particle vertex for the metallic and the insulating phase. Then we go to a two-site cluster DMFT scheme for the two-dimensional Hubbard model that includes short-range antiferromagnetic fluctuations and obtain the local and nonlocal two-particle vertex functions. We discuss the rich frequency structures of these vertices and compare with the vertex in the single-site solution.
We apply a recently introduced hybridization-flow functional renormalization group scheme for Anderson-like impurity models as impurity solver in a dynamical mean field theory (DMFT) approach to lattice Hubbard models. We present how this scheme is capable of reproducing metallic and insulating solutions of the lattice model. Our setup also offers a numerically rather inexpensive method to calculate two-particle correlation functions. For the paramagnetic Hubbard-model on the Bethe lattice in infinite dimensions we calculate the local two-particle-vertex for the metallic and the insulating phase. Then we go to a two-site cluster-DMFT-scheme for the two-dimensional Hubbard-model that includes short-range antiferromagnetic fluctuations and obtain the local and non-local two-particle-vertex-functions. We discuss the rich frequency structures of these vertices and compare with the vertex in the single-site solution.
Strongly correlated electron systems host a plethora of fascinating physical phenomena and pose formidable challenges in their theoretical analysis. The challenges originate from the inherent complexity of the quantum many-body problem - no classical computer will ever be able to fully simulate these systems - and the lack of an effective single-particle picture, as the strong mutual interactions make it impossible to regard the electrons as independent from each other. As a consequence, most systems of correlated electrons can only be tackled approximately and numerically. In this thesis, we develop a set of numerical methods for strongly correlated electrons, which are inspired by the renormalization group (RG) idea of including degrees of freedom successively from high to low energies. This enables an efficient organization of the diverse fluctuations and is key for an accurate treatment of interacting quantum systems, where collective behavior and composite objects emerge at energy scales far below those of the microscopic constituents. In a first part, we consider the functional renormalization group (fRG), a versatile framework to study the flow of correlation functions upon modulating the underlying action. Though widely used, it has often acted more as a qualitative rather than quantitative method, due a nontransparent approximation induced by truncating the hierarchy of flow equations. We develop an iterative multiloop fRG (mfRG) scheme, which ameliorates this approximation and eliminates many of the drawbacks of fRG experienced hitherto. In particular, it restores the independence of results on the choice of RG regulator and establishes a rigorous relation to the parquet formalism. Furthermore, we show how to derive the flow equations directly from self-consistent many-body relations. This establishes a form of diagrammatic resummations at the two-particle level which circumvents ill-behaved two-particle-irreducible vertices. An application to the prototypical two-dimensional Hubbard model illustrates how our multiloop scheme elevates the fRG approach to correlated electron systems to a quantitative level. Secondly, we employ the numerical renormalization group (NRG), based on the iterative diagonalization of impurity Hamiltonians, in conjunction with the dynamical mean-field theory (DMFT) to describe local correlations in multiorbital systems. Having access to arbitrarily low temperatures and energies, NRG is a unique, real-frequency impurity solver for DMFT. It has been pivotal to the understanding of Hund metals, where strong correlations arise from Hund's rules even at moderate Coulomb repulsion. Building on recent methodological advances, we extend the range of application of DMFT+NRG from orbital-degenerate models to more realistic setups: We first study orbital differentiation in a three-orbital Hund-metal model and unravel key effects of the orbital-selective Mott transition. In a real-materials setting, we then incorporate the bandstructure from density functional theory (DFT) and analyze the archetypal Hund-metal material Sr2RuO4. We particularly follow its RG flow to the Fermi-liquid regime at previously inaccessible low temperatures and generally present DFT+DMFT+NRG as a new computational paradigm for strongly correlated materials. As a side project of our fRG work, we develop an algorithm to count Feynman diagrams from closed many-body relations, which reveals the surprising outcome that totally irreducible contributions are responsible for the factorial growth in the number of diagrams. Additionally, we use NRG to study transport through three-level quantum dots and provide benchmark data for other RG methods, which aim at further describing these systems in nonequilibrium.
We investigate nonequilibrium properties of the single-impurity Anderson model by means of the functional renormalization group (fRG) within Keldysh formalism. We present how the level broadening $\ensuremath{\Gamma}/2$ can be used as flow parameter for the fRG. This choice preserves important aspects of the Fermi-liquid behavior that the model exhibits in case of particle-hole symmetry. An approximation scheme for the Keldysh fRG is developed which accounts for the frequency dependence of the two-particle vertex in a way similar but not equivalent to a recently published approximation to the equilibrium Matsubara fRG. Our method turns out to be a flexible tool for the study of weak to intermediate on-site interactions $U\ensuremath{\lesssim}3\ensuremath{\Gamma}$. In equilibrium we find excellent agreement with numerical RG results for the linear conductance at finite gate voltage, magnetic field, and temperature. In nonequilibrium, our results for the current agree well with time-dependent density-matrix RG. For the nonlinear conductance as function of the bias voltage, we propose reliable results at finite magnetic field and finite temperature. Furthermore, we demonstrate the exponentially small scale of the Kondo temperature to appear in the second-order derivative of the self-energy. We show that the approximation is, however, not able to reproduce the scaling of the effective mass at large interactions.
Technological advances in the fabrication of mesoscopic structures allow for an experimental investigation of low-dimensional interacting many-body systems in nonequilibrium situations. From a theoretical perspective, the accurate description of such systems constitutes a major challenge which has in general not been solved. In light of possible applications, such as quantum computation, but also for a better understanding of correlation effects themselves, it is thus important to put a research focus on methodological advances.The so called functional renormalization group (fRG) has turned out to be a successful and versatile tool in the application to the aforementioned systems. Being an RG method, it has the associated advantages, e.g. a successive integrating out of energy scales and avoiding divergences that perturbation theory would have. In this thesis, we present two paradigmatic extensions of the fRG approach.Typically, the fRG is used to investigate the flow of the generating functional of the one-particle-irreducible (1PI) vertex functions. Typical approximations in this framework bring about certain shortcomings, such as the violation of conservation laws. In order to overcome some of these shortcomings, we formulate the fRG for the generating functional of the two-particle-irreducible (2PI) vertex functions. This step constitutes our first extension of the fRG approach. While typically only nonconserving schemes are generated in 1PI fRG, we can define nonconserving and conserving approximation schemes within 2PI fRG. Partially, these are found to be equivalent to known diagrammatic approximations. We apply the various 2PI fRG schemes to the quantum anharmonic oscillator, as a “toy model”, and to the Anderson impurity model, as a full fermionic quantum many-body problem. Overall, we find that the devised 2PI fRG schemes (even the conserving ones) offer no great advantage over established 1PI methods. As an exception, the nonconserving lowest-order 2PI fRG scheme in which the flow parameter is introduced in the interaction produces better results than 1PI fRG schemes which are comparable in computational effort.Our second extension of the fRG approach is to apply it to models which cannot (directly) be formulated in second quantization. Having a localized spin, the Kondo model is a simple example of such models which nevertheless shows nontrivial (but well-understood, e.g. by Anderson’s Poor Man’s Scaling) correlation effects. It has hardly been treated with fRG in the literature. In this thesis, we present how to overcome obstacles arising in the formulation of the Kondo model in functional integral formalism. Instead of the 1PI fRG, we use the Polchinski fRG which is based on connected amputated Green functions. This makes it possible to make a close connection to Poor Man’s Scaling which employs the concept of the T-matrix. We reproduce the results of Poor Man’s Scaling and also devise a systematic extension to it by keeping the frequency dependence of the flowing amputated Green functions in the RG equations. With this scheme, we do not find a significant qualitative improvement over Poor Man’s Scaling. However, we have layed the groundwork for more successful future extensions of the fRG applied to the Kondo model.
The progress in micro-fabrication technologies makes it possible to address individual, tailored, mesoscopic structures with sizes on a nanometer to micrometer scale. The transport properties of such quasi zero- or one-dimensional systems are considerably influenced by the electron-electron interaction. In theoretical studies this often becomes apparent in the emergence of divergences in a perturbative treatment of the interaction. A successful approach may then be given by renormalization group (RG) methods which are based on the idea of treating energy scales successively from high to low. In the functional RG (fRG) an exact infinite hierarchy of flow equations for multi-particle Green or vertex functions is derived, that allows for systematic approximations on different levels. In numerous studies the fRG has been found to be a powerful tool for the investigation of low dimensional Fermi systems. Usually it is formulated within Matsubara formalism which applies to a system in thermal equilibrium. A central topic of this thesis is the implementation of the fRG within Keldysh formalism, providing the possibility to treat equilibrium and non-equilibrium situations on equal footing. The resulting Keldysh fRG is used in this thesis to explore the transport properties of quantum wires with several barriers and of the single impurity Anderson model. The thesis starts with an introduction to Keldysh formalism and to related concepts. Then the structure of the Green and vertex functions of Keldysh formalism is investigated. This includes studying the consequences of causality, characterizing thermal equilibrium, and formulating corresponding relations for the multi-particle functions. After that, the implementation of the fRG within Keldysh formalism is presented in detail. The flow equations for the vertex functions are derived via their diagrammatic expansion. Two especially useful flow parameters are determined: one of them cuts off the imaginary frequencies that appear as poles of the Fermi functions of the reservoirs. The second flow parameter is an artificially enhanced reservoir hybridization; on one hand it has the advantage to preserve the fluctuation dissipation theorem in equilibrium even in frequency dependent approximations; on the other hand, its use can be more cumbersome. One of the models studied in this thesis is a tight-binding chain with nearest-neighbour interaction. At weak and intermediate interactions it belongs to the class of Luttinger liquids (LLs) whose low energy properties are characterized by power laws. The transport properties of LLs have been studied experimentally for example on carbon nanotubes. However, the presence of transport barriers at the contacts make the interpretation of the experimental data difficult. A theoretical treatment with the fRG permits a rather large flexibility to model the barrier configuration. In this thesis the fRG is first used in its equilibrium Matsubara implementation in order to analyse the temperature dependence of the linear conductance in the presence of two, three or four barriers. It is shown that temperature induces an averaging over quantum phases of propagation, leading to combinations of the power law exponents of the individual barriers. An experimental set-up that allows to verify this behaviour is proposed. In a second step the Keldysh version of the fRG is applied to the problem of transport through a quantum wire at finite bias voltage. For a wire with two high contact barriers it is found that the non-equilibrium particle distribution in the wire is responsible for the appearance of a two-fold power law in the density of states. An additional weakly attached probe lead can be used to make the two different power laws visible in the conductance. The second model which is analysed is the one of the single Anderson impurity, often considered as a paradigm for the investigation of strong local electron correlations. The latter become noticeable in the appearance of the Kondo scale, an energy exponentially small in the interaction strength, that characterizes the low energy physics of the model in the regime of spin fluctuations. On the basis of this model an advanced frequency dependent approximation for the Keldysh fRG is developed in this thesis. Such an approximation can account for subleading effects of the interaction, namely the fact that two-particle interactions can generate decay rates for single-particle states. The derivation of a consistent approximation scheme is presented, that turns out to yield a good description of the model at moderate interaction strength.