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.
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.
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.
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.
In this thesis, we tackle various problems in strongly correlated electron systems, which can be addressed properly via the non-perturbative dynamical mean field theory (DMFT) approach using the continuous time quantum Monte Carlo method as an impurity solver. First, we revisit the old Nagaoka ferromagnetism problem in the U=oo Hubbard model and study the stability of the ferromagnetic state as a function of the temperature, the doping level, and the next-nearest-neighbor hopping t'. We then address the nature of the Mott transition in the two-dimensional Hubbard model at half-filling using cluster DMFT. Cluster DMFT can incorporate the short-range correlations beyond DMFT by extending the spatial range in which correlations are treated exactly to a finite cluster size. The non-local correlations reduce substantially the critical interaction U and modify the shape of the transition lines in the phase diagram. We then concentrate on the calculation of two-particle response functions from the ab initio perspective by means of computing the one-particle excitation spectrum using the combination of the density functional theory (DFT) and DMFT and extracting the two-particle irreducible vertex function from a local two-particle Green's function computed within DMFT. In particular, we derive the equations for calculating the magnetic/charge susceptibility and the pairing susceptibility in superconductivity. This approach is applied to the Hubbard model and the periodic Anderson model and we determine the phase diagram of magnetism and superconductivity in these models. We show that the superconducting phase is indeed stable near the magnetic phase where the pairing interaction mediated by spin fluctuations is dominantly enhanced. The non-local correlation effect to superconductivity is also discussed using the dual fermion approach and the dynamical vertex approximation. We finally apply the vertex function approach within DFT+DMFT to a Fe-based superconductor, BaFe2As2, and compute the dynamical magnetic susceptibility in this material. Our calculation results show a good agreement with the magnetic excitation spectra observed in a neutron scattering experiment. The response function calculation method derived in this thesis can capture both a localized and an itinerant nature of collective excitations in strongly correlated electron systems.
We study the Anderson-Hubbard model to examine the interplay between the disordered potential and the Hubbard interactions. We use the dynamical mean-field theory (DMFT) with an exact diagonalization as an impurity solver to effectively deal with the quantum fluctuations due to the local interactions. In order to characterize the paramagnetic metal-insulator transition in disordered systems, we compute the typical local density of states, as well as the arithmetically-averaged local density of states. Our approach is found to reproduce the characteristic features of the metal-insulator transition in infinite dimensions, such as the existence of the disorder-driven metallic phase for moderate values of disorder strength and interactions, as observed in earlier studies. In two dimensions, the system turns out to exhibit more insulating behavior when we increase the cluster size in the cellular DMFT. The metallic phase, which is induced from weak Mott insulators by the introduction of a disordered potential, is not observed in two dimensions.
We propose an approach describing correlated electronic systems within the two-particle irreducible functional renormalization-group (2PI-fRG) approach, which uses initial self-energy and the two-particle irreducible vertices, obtained from (extended) dynamical mean-field theory [E(DMFT)]. Using the 2PI-fRG approach allows us to treat both local and nonlocal interactions. In the case of purely local interaction the corresponding equations have similar (although not identical) structure to the earlier developed DMFT+fRG (DMF2RG) approach. In a simplest truncation, neglecting scale dependence of the two-particle irreducible vertices, we reproduce the results for the two-particle vertices/susceptibilities in the ladder approximation of the dual-boson (DB) or the dynamic vertex approximation approach; in more sophisticated truncations the method allows us to consider nonlocal corrections to the self-energy as well as the interplay of charge and spin correlations. The proposed scheme is tested on the two-dimensional standard and extended $U\text{\ensuremath{-}}V$ half-filled Hubbard models. For the standard Hubbard model we obtain nonlocal self-energy, which is in agreement with numerical studies; for the extended Hubbard model we find the boundary of charge instability, which agrees well with the results of the DB approach. We also ascertain that the effect of spin correlations on electron interaction in the charge channel, not considered previously in the DB approach, only slightly reduces critical next-nearest-neighbor interaction of charge instability of the extended Hubbard model at the considered finite low temperature, yielding better agreement with dynamic cluster approximation. The considered method is rather general and can be applied to study various phenomena in strongly correlated electronic systems.
Using the dynamical mean-field theory (DMFT) as a ``booster rocket'', the functional renormalization group (fRG) can be upgraded from a weak-coupling method to a powerful computation tool for strongly interacting fermion systems. The strong local correlations are treated nonperturbatively by the DMFT, while the fRG flow can be formulated such that it is driven exclusively by nonlocal correlations, which are more amenable to approximations. We show that the full frequency dependence of the two-particle vertex needs to be taken into account in this approach, and demonstrate that this is actually possible---in spite of the singular frequency dependence of the vertex at strong coupling. We are thus able to present results obtained from the DMFT-boosted fRG for the two-dimensional Hubbard model in the strongly interacting regime. We find strong antiferromagnetic correlations from half filling to 18% hole doping and at the lowest temperature we can access, a sizable d-wave pairing interaction driven by magnetic correlations at the edge of the antiferromagnetic regime.