We show that the numerical renormalization group is a viable multi-band impurity solver for dynamical mean-field theory (DMFT), offering unprecedented real-frequency spectral resolution at arbitrarily low energies and temperatures. We use it to obtain a numerically exact DMFT solution to the Hund metal problem for a three-band model on a Bethe lattice at 1/3 filling. The ground state is a Fermi liquid. The one-particle spectral function undergoes a coherence-incoherence crossover with increasing temperature, with spectral weight being transferred from low to high energies. Further, it exhibits a strong particle-hole asymmetry. In the incoherent regime, the self-energy displays approximate power-law behavior for positive frequencies only. The spin and orbital spectral functions show "spin-orbital separation": spin screening occurs at much lower energies than orbital screening. The renormalization group flows clearly reveal the relevant physics at all energy scales.
We show that the numerical renormalization group is a viable multi-band impurity solver for dynamical mean-field theory (DMFT), offering unprecedented real-frequency spectral resolution at arbitrarily low energies and temperatures. We use it to obtain a numerically exact DMFT solution to the Hund metal problem for a three-band model on a Bethe lattice at 1/3 filling. The ground state is a Fermi liquid. The one-particle spectral function undergoes a coherence-incoherence crossover with increasing temperature, with spectral weight being transferred from low to high energies. Further, it exhibits a strong particle-hole asymmetry. In the incoherent regime, the self-energy displays approximate power-law behavior for positive frequencies only. The spin and orbital spectral functions show "spin-orbital separation": spin screening occurs at much lower energies than orbital screening. The renormalization group flows clearly reveal the relevant physics at all energy scales.
For a long time strong electronic correlations in metals have mainly been associated with Mottness, the proximity to a Mott metal-insulator transition (MIT), where large Coulomb interactions induce the localization of charges. However, triggered by the discovery of the iron-based superconductors about ten years ago, it was realized that multi-orbital materials with only moderate Coulomb but sizeable Hund’s rule interactions – so-called Hund metals – allow for a distinct screening mechanism towards strong correlations: Hundness. Here, Hund’s rule constrains the spin rather than the charge dynamics. This discovery led to a vividly debated fundamental issue in the field of strongly correlated condensed matter systems, which is the main topic of the present thesis: what is the origin of strong correlations in the normal phase of Hund metals, Mottness or Hundness? And what are their decisive fingerprints? The goal of this dissertation is twofold. First, we present and advance our method: the numerical renormalization group (NRG) as viable real-frequency multi-band impurity solver for dynamical mean-field theory (DMFT), a common approach to tackle strongly correlated systems. Second, we apply DMFT+NRG to shed light on the Hund-metal problem raised above. In the first part of this thesis we present our state-of-the-art NRG solver, which offers direct access to data with unprecedented real-frequency spectral resolution at arbitrarily low energies and temperatures in contrast to commonly used Quantum Monte Carlo solvers. It is based on matrix product states and exploits non-abelian symmetries to reduce numerical costs. In the case of orbital symmetry, this allows us to treat multi-band models with more than two bands, and thus to tackle the Hund-metal problem for the first time with NRG. For multi-band models without orbital symmetry, an “interleaved” scheme of NRG (iNRG) was recently developed, dramatically increasing the numerical efficiency. Remarkably, the accuracy of iNRG is comparable to standard NRG, as we reveal in a detailed study. This finding establishes iNRG as a promising DMFT solver for material-specific model simulations. In the second part of this thesis we study a minimal toy model for Hund metals with DMFT+NRG, the orbital-symmetric three-band Hubbard-Hund model (3HHM) close to a lattice filling of 1/3. Our major insight is “spin-orbital separation” (SOS), a Hund’s-ruleinduced two-stage Kondo-type screening process, in which orbital screening occurs at much higher energies than spin screening. In Hund metals, i.e. far from a MIT phase boundary, SOS thus causes large electron masses by strongly reducing the coherence scale below which a Fermi liquid is formed. Further, it opens up a broad incoherent and strongly particle-hole asymmetric intermediate energy regime that reaches up to bare excitation scales. This SOS regime shows fractional power-law behavior and is characterized by resilient “Hund quasiparticles” with itinerant orbital degrees of freedom coupled non-trivially to quasi-free large spins. At zero temperature, the local density of states exhibits a two-tier quasiparticle peak on top of a broad incoherent background. In contrast, in Mott-correlated metals, i.e. close to the MIT phase boundary, the SOS regime becomes negligibly small and the Hubbard bands are well separated. These findings lead to distinct signatures of Hundness and Mottness in the temperature dependence of ARPES spectra, static local susceptibilities, resistivity, thermopower and entropy, many of which were also found in realistic simulations of the archetypal Hund- and Mott-correlated materials, Sr2RuO4 and V2O3. In summary, we provide evidence that and elucidate how Hundness evokes strong correlation effects in Hund metals. This might help to better interpret experimental results and guide superconducting theories.
Hund metals have attracted attention in recent years due to their unconventional superconductivity, which supposedly originates from non-Fermi-liquid (NFL) properties of the normal state.When studying Hund metals using dynamical mean-field theory, one arrives at a self-consistent "Hund impurity problem" involving a multiorbital quantum impurity with nonzero Hund coupling interacting with a metallic bath.If its spin and orbital degrees of freedom are screened at different energy scales, T sp < T orb , the intermediate energy window is governed by a novel NFL fixed point, whose nature had not yet been clarified.We resolve this problem by providing an analytical solution of a paradigmatic example of a Hund impurity problem, involving two spin and three orbital degrees of freedom.To this end, we combine a state-ofthe-art implementation of the numerical renormalization group, capable of exploiting non-Abelian symmetries, with a generalization of Affleck and Ludwig's conformal field theory (CFT) approach for multichannel Kondo models.We characterize the NFL fixed point of Hund metals in detail for a Kondo model with an impurity forming an SUð2Þ × SUð3Þ spin-orbital multiplet, tuned such that the NFL energy window is very wide.The impurity's spin and orbital susceptibilities then exhibit striking power-law behavior, which we explain using CFT arguments.We find excellent agreement between CFT predictions and numerical renormalization group results.Our main physical conclusion is that the regime of spinorbital separation, where orbital degrees of freedom have been screened but spin degrees of freedom have not, features anomalously strong local spin fluctuations: the impurity susceptibility increases as χ imp sp ∼ ω -γ , with γ > 1.
The natural orbitals renormalization group (NORG) has previously been proposed as an efficient numerical method for solving zero-temperature properties of multisite and multiorbital quantum impurity systems. Here, we implement the NORG as an impurity solver for dynamical mean-field theory (DMFT). In comparison with the exact diagonalization method, the NORG method can treat much more bath sites in an impurity model to which the DMFT maps a lattice model and can find accurate zero-temperature Matsubara and low-frequency retarded Green's functions. We demonstrate the effectiveness of this method on a two-orbital Hubbard model on the Bethe lattice and find successfully the orbital selective Mott transition with a Kondo resonance peak in the wide band and two holon-doublon bound state excitation peaks in the narrow band.
Many correlated metallic materials are described by Landau Fermi-liquid theory at low energies, but for Hund metals the Fermi-liquid coherence scale T_{FL} is found to be surprisingly small. In this Letter, we study the simplest impurity model relevant for Hund metals, the three-channel spin-orbital Kondo model, using the numerical renormalization group (NRG) method and compute its global phase diagram. In this framework, T_{FL} becomes arbitrarily small close to two new quantum critical points that we identify by tuning the spin or spin-orbital Kondo couplings into the ferromagnetic regimes. We find quantum phase transitions to a singular Fermi-liquid or a novel non-Fermi-liquid phase. The new non-Fermi-liquid phase shows frustrated behavior involving alternating overscreenings in spin and orbital sectors, with universal power laws in the spin (ω^{-1/5}), orbital (ω^{1/5}) and spin-orbital (ω^{1}) dynamical susceptibilities. These power laws, and the NRG eigenlevel spectra, can be fully understood using conformal field theory arguments, which also clarify the nature of the non-Fermi-liquid phase.
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.
Numerical studies on strongly correlated fermionic systems are very complicated and still provide essential problems. The main reason is the exponential growth of the un- derlying Hilbert state space with the system size and the fermionic sign problem for Monte Carlo studies. Among the most widely employed numerical techniques for study- ing two-dimensional quantum many-body systems are cluster extensions of the dynamical mean-field theory (DMFT), e.g. dynamical cluster approximation (DCA). They map an infinitely large multi-dimensional lattice problem to a one-dimensional impurity problem. In 2015 it was shown that the density matrix renormalisation group (DMRG) used as an impurity solver for DMFT (DMFT+DMRG) on the imaginary-frequency axis allows to solve multi-site and multi-band problems extremely fast compared to other solvers. Within this thesis, we further develop this DMRG+DMFT approach to apply the method on real material settings. The step from artificial, completely degenerate multi-band mod- els with simple dispersion relations on a Bethe lattice, studied in 2015, to systems with realistic band structures and lifted degeneracies involves more challenges than originally suspected. In this thesis, we will first recapitulate relevant methods for our approach like matrix prod- uct states, the density matrix renormalisation group and several time evolution methods. In this context we will present several improvements ranging from optimised time evo- lutions to entanglement based optimisations of tensor networks. Second, we will present a very detailed description of the dynamical mean field theory. We will focus on both methodological aspects and implementation details. This chapter is intended to allow other researcher to implement their own DMFT code using DMRG as an impurity solver. Third, we will discuss three different models to show the extent of problems DMRG+ DMFT is able to solve. We will focus on multi-site DCA calculations in the case of the two-dimensional Hubbard model and show that DMRG allows to tackle systems with intermediate interaction strengths at low temperatures, which are unsolvable with other solvers. In the second case, the real material Sr2VO4, we will show the first two-site DCA results for a realistic three-band model. In contrast to assumptions, partly reintroducing the momentum dependence of the self-energy does not improve agreement between exper- imental observations and theoretical results. Finally, we will move on to another realistic three-band model, which describes Sr2RuO4, to show how to deal with the influence of spin-orbit coupling on DMFT. We will present the first low-temperature results for this material and will confirm previous results of simplified model calculations.