Strong correlations arising from interactions between particles are responsible for many phenomena in condensed matter physics, from superconductivity to heavy fermions. An accurate theoretical description of such complex systems requires to solve so-called many-body problems, which are inherently difficult to tackle due to the exponential growth of the underlying Fock space with the size of the system. Usually, single-particle descriptions or perturbation theory at first orders fail to provide satisfactory results in correlated regimes, and pertinent Ansätze or cutting edge numerical methods are required to solvethese problems.Among the wide variety of strongly correlated electronic systems, quantum impurity systems are of central interest. These systems, where a few localized degrees of freedom that undergo strong interactions hybridize with a larger system of free particles, have been widely studied for decades. Beyond their relevance to study dilute dynamic impurities in metals, quantum dots, or more complex nanostructures, they are exploited in a broader class of problems through dynamical mean field theory (DMFT), that is designed for fully interacting lattice problems. Despite the effort devoted to studying of these problems,some questions remain open and new approaches are required to address them.In this thesis, quantum impurity problems are tackled with the use of natural orbitals, a representation rooted in the framework of quantum chemistry. The numerical renormalization group, the standard non-perturbative method to treat these systems, is employed to demonstrate the hierarchical structure of correlations unveiled by natural orbitals. An ansatz for the wavefunction of the ground state is then proposed, separating explicitely a small set of correlated orbitals to a large uncorrelated one, leaving the problem to a few-body complexity. This simplification brought new insights for simulating quantumimpurity problems, and a new algorithm is developed to generate an optimized set of correlated orbitals independently of existing methods, going beyond their usual limitations.In the first place, the algorithm is used to compute the screening cloud in large one-dimensional chains and in a realistic 2D square lattice coupled to a localized adatom. This observable is of great interest to understand the screening mechanism of the magnetic moment of the impurity that plays a central role in transport properties, for unidimensional systems or in bulk metals experiencing RKKY interactions. Spatial correlations were found to be highly non-symmetric, breaking the s-wave picture. In one direction, corre-lations vanish way faster than the usual isotropic correlation length, while on the other one they spread further than this same expected length. This anisotropy suggests that impurities in bulk systems may interact differently depending on their relative positions with respect to the lattice.Taking further advantage of the efficiency of this new algorithm allowed us to perform a statistical study on large systems to investigate correlations around a quantum impurity in a disordered host, averaging over a large amount of realizations, which was not possible with previous approaches. The question motivating this study was the fate of the correlation cloud spreading into the disordered electronic bath. Although spatial correlations are expected to vanish when the disorder increases, long range correlationssurvive to some configurations of the potential, inducing a bimodal structure of probability distribution functions of relevant scales.Finally, these disordered problems are investigated in the non-interacting situation with random matrix theory, which reproduced qualitatively the results obtained before. The joint distribution of one-particle orbitals energies and amplitudes is calculated, which allows to compute any disordered averaged local correlation functions, and calculations in the large-N limit are presented.
Strong correlations arising from interactions between particles are responsible for many phenomena in condensed matter physics, from superconductivity to heavy fermions. An accurate theoretical description of such complex systems requires to solve so-called many-body problems, which are inherently difficult to tackle due to the exponential growth of the underlying Fock space with the size of the system. Usually, single-particle descriptions or perturbation theory at first orders fail to provide satisfactory results in correlated regimes, and pertinent Ansätze or cutting edge numerical methods are required to solvethese problems.Among the wide variety of strongly correlated electronic systems, quantum impurity systems are of central interest. These systems, where a few localized degrees of freedom that undergo strong interactions hybridize with a larger system of free particles, have been widely studied for decades. Beyond their relevance to study dilute dynamic impurities in metals, quantum dots, or more complex nanostructures, they are exploited in a broader class of problems through dynamical mean field theory (DMFT), that is designed for fully interacting lattice problems. Despite the effort devoted to studying of these problems,some questions remain open and new approaches are required to address them.In this thesis, quantum impurity problems are tackled with the use of natural orbitals, a representation rooted in the framework of quantum chemistry. The numerical renormalization group, the standard non-perturbative method to treat these systems, is employed to demonstrate the hierarchical structure of correlations unveiled by natural orbitals. An ansatz for the wavefunction of the ground state is then proposed, separating explicitely a small set of correlated orbitals to a large uncorrelated one, leaving the problem to a few-body complexity. This simplification brought new insights for simulating quantumimpurity problems, and a new algorithm is developed to generate an optimized set of correlated orbitals independently of existing methods, going beyond their usual limitations.In the first place, the algorithm is used to compute the screening cloud in large one-dimensional chains and in a realistic 2D square lattice coupled to a localized adatom. This observable is of great interest to understand the screening mechanism of the magnetic moment of the impurity that plays a central role in transport properties, for unidimensional systems or in bulk metals experiencing RKKY interactions. Spatial correlations were found to be highly non-symmetric, breaking the s-wave picture. In one direction, corre-lations vanish way faster than the usual isotropic correlation length, while on the other one they spread further than this same expected length. This anisotropy suggests that impurities in bulk systems may interact differently depending on their relative positions with respect to the lattice.Taking further advantage of the efficiency of this new algorithm allowed us to perform a statistical study on large systems to investigate correlations around a quantum impurity in a disordered host, averaging over a large amount of realizations, which was not possible with previous approaches. The question motivating this study was the fate of the correlation cloud spreading into the disordered electronic bath. Although spatial correlations are expected to vanish when the disorder increases, long range correlationssurvive to some configurations of the potential, inducing a bimodal structure of probability distribution functions of relevant scales.Finally, these disordered problems are investigated in the non-interacting situation with random matrix theory, which reproduced qualitatively the results obtained before. The joint distribution of one-particle orbitals energies and amplitudes is calculated, which allows to compute any disordered averaged local correlation functions, and calculations in the large-N limit are presented.
Dynamical quantum impurities in metals induce electronic correlations in real space that are difficult to simulate due to their multiscale nature, so that only $s\text{\ensuremath{-}}\mathrm{wave}$ scattering in clean metallic hosts has been investigated so far. However, screening clouds should show anisotropy due to lack of full rotational invariance in two- and three-dimensional lattices, while inherent disorder will also induce spatial inhomogeneities. To tackle these challenges, we present an efficient and robust algorithm based on the recursive generation of natural orbitals defined as eigenvectors of the truncated single-particle density matrix. This method provides well-converged many-body wave functions on lattices with up to tens of thousands of sites, bypassing some limitations of other approaches. The algorithm is put to the test by investigating the charge screening cloud around an interacting resonant level, both on clean and disordered lattices, achieving accurate spatial resolution from short to long distances. We thus demonstrate strong anisotropy of spatial correlations around an adatom in the half-filled square lattice. Taking advantage of the efficiency of the algorithm, we further compute the disorder-induced distribution of Kondo temperatures over several thousands of random realizations, at the same time gaining access to the full spatial profile of the screening cloud in each sample. While the charge screening cloud is typically shortened due to the polarization of the impurity by the disorder potential, we surprisingly find that rare disorder configurations preserve the long-range nature of Kondo correlations in the electronic bath.
In recent years, nano-sized systems involving local charge or spin states have received wide interests because of their potential application in emerging fields such as quantum information and quantum computation. For instance, organometallic molecular complexes may serve as building blocks of quantum storage devices, because the spin-unpaired d or f electrons at transition metal centers may be employed to construct spin qubits. Moreover, if the system contains strong electron-electron interaction, the involving local quantum states are subject to prominent electron correlation effects (such as the Kondo effect). Theoretically, spatially confined nano-systems are often described by quantum impurity models. Thus the accurate prediction of the intrinsic properties of quantum impurity systems and the deep understanding on the response and evolution of local quantum states under external fields or in dissipative environment are fundamentally important for the design and fabrication of quantum devices. The accurate characterization of quantum coherence, correlation, and entanglement in quantum impurity systems remains a great challenge. Enormous efforts have been made to achieve this goal. A variety of theoretical methods have been developed, including the numerical renormalization group method, the quantum Monte Carlo method, and many others. However, all the existing methods are subject to certain limitations regarding accuracy or efficiency. Therefore, we choose to view this problem from a new perspective—the perspective of open quantum systems. Over the past decade, we have developed a formally exact quantum dissipation theory, the hierarchical equations of motion (HEOM) theory, for fermionic open systems. The HEOM theory captures the combined effects of system-environment dissipation, many-body interaction, and non-Markovian memory in a nonperturbative manner. It is capable of addressing static and dynamic responses of system observables in both equilibrium and nonequilibrium situations. We have implemented the HEOM method in our self-designed computer program HEOM-QUICK. We have also devised a series of advanced algorithms which substantially enhance the numerical efficiency of HEOM. The HEOM-QUICK program thus provides an accurate, efficient, and versatile theoretical tool for the investigation of strongly correlated quantum impurity systems. We have applied the HEOM method to study a variety of problems associated with quantum impurity systems. These include the intrinsic properties of quantum impurity models and lattice models, the transient current response of quantum dots to ac voltages, the thermopower and local heating effect in nonequilibrium quantum dots, and the tuning of local spin states in adsorbed molecular magnets. In particular, the HEOM method has been combined with density functional theory (DFT) method, and thus allows for first-principles-based simulation on the precise tuning of local spin states in adsorbed magnetic molecules. Numerical simulations have discovered or reproduced many important quantum phenomena that are essential for relevant experiments, including the Kondo effect, Mott metal-insulator transition, quantum memristive effect, etc. This paper gives a brief overview of our development of the HEOM method for fermionic open systems, as well as its applications to various strongly correlated quantum impurity systems.
Correlated electron systems exhibit a variety of intriguing emergent phenomena, including Mott physics, pseudogap behavior, and unconventional superconductivity. One of the most widely used numerical approaches for investigating these systems is dynamical mean-field theory (DMFT). This dissertation addresses two related problems: how to develop better quantum impurity solvers for DMFT and how to use DMFT results to gain a better understanding of correlated electron systems. The first part of this dissertation develops numerical methods for solving quantum impurity problems, which constitute the central computational task in DMFT calculations. In Chap. 4, we develop a method that combines a tensor-train numerical integration scheme with the inchworm hybridization expansion for quantum impurity problems. Benchmarks on spinless, single-orbital, and multi-orbital impurity models demonstrate the convergence of the Green’s functions to exact or reference results and identify the regimes in which tensor-train compression of Feynman diagrams and the inchworm expansion itself are effective. In Chap. 5, we introduce the concept of counterterms within the hybridization expansion and apply it to a bare hybridization expansion solver that relies on a sum-of-exponentials integration technique. Tests on spinless, multi-orbital, and realistic embedding impurity problems show that counterterms can significantly reduce the required perturbative order and improve the efficiency of hybridization-expansion calculations. The second part of this dissertation applies the cluster extension of DMFT to the single-band Hubbard model. In Chap. 6, we introduce a fluctuation decomposition for the self-energy that eliminates the Fierz ambiguity encountered in conventional fluctuation diagnostics of the Hubbard model. The resulting analysis shows that spin fluctuations make an important contribution to the pseudogap self-energy, while multi-boson contributions are also necessary for a complete description. In Chap. 7, we analyze the spin-fermion coupling and its role in d-wave pairing. Our calculations show that, for antinodal fermions in the pseudogap regime, their coupling to antiferromagnetic spin fluctuations is strongly enhanced relative to the non-interacting limit. A modified Bethe–Salpeter equation analysis further shows that this enhanced spin–fermion coupling boosts the spin-fluctuation-induced contribution to the d-wave pairing susceptibility.
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.
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.
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.