One of the most important and accurate approaches for studying interacting many-electron correlated systems is the quantum renormalization group (RG) procedure, upon which we propose a concept in the framework of natural orbitals so that we can generalize the RG into general orbital space. To demonstrate its effectiveness, we apply this new approach on a quantum cluster-impurity model with four impurities in comparison with the exact diagonalization method. We also find a shortcut to dramatically improving this RG procedure. A profound impact on the cluster dynamical mean-field theory and multi-impurity/orbital Kondo problems is further discussed.
One of the most important and accurate approaches for studying interacting many-electron correlated systems is the quantum renormalization group (RG) procedure, upon which we propose a concept in the framework of natural orbitals so that we can generalize the RG into general orbital space. To demonstrate its effectiveness, we apply this new approach on a quantum cluster-impurity model with four impurities in comparison with the exact diagonalization method. We also find a shortcut to dramatically improving this RG procedure. A profound impact on the cluster dynamical mean-field theory and multi-impurity/orbital Kondo problems is further discussed.
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.
Renormalization group methods become natural for quantum many-body systems where there is an emergence of new energy scales like as Kondo scale in case of Kondo model. It was the Numerical Renormalization Group which gave the full quantitative understanding of the strong coupling fixed point of Kondo model and also of Anderson Impurity model. However, it was Andersons poor man scaling method which gave the analytical insights into the scaling behaviour of quantum impurity models. That makes poor man scaling method very important, and hence it has been used extensively in carrying out renormalization studies of quantum impurity models. This method has been generalized and also applied to other quantum many-body Hamiltonians. In this paper, we give a detailed presentation of poor man scaling method and apply it to study the scaling behaviour of Kondo model and Anderson impurity model as well. Though there are many reviews of perturbative renormalization of Kondo model, in this paper we have also presented a detailed discussion about scaling analysis of Anderson impurity model. We also present another method called flow equation method which is a generalization of poor man scaling method and has recently been applied to study the quantum many-body dynamics both in equilibrium and out-of equilibrium. We also present a numerical procedure to extract Kondo scale from flow equations of Kondo model. The importance of this procedure is that it can be applied to calculate the Kondo scale for situations where poor mans scaling method can not be applied.
Generalized quantum impurity models—which feature a few localized and strongly correlated degrees of freedom coupled to itinerant conduction electrons—describe diverse physical systems, from magnetic moments in metals to nanoelectronics quantum devices such as quantum dots or single-molecule transistors. Correlated materials can also be understood as self-consistent impurity models through dynamical mean-field theory. Accurate simulation of such models is challenging, especially at low temperatures, due to many-body effects from electronic interactions, resulting in strong renormalization. In particular, the interplay between local impurity complexity and Kondo physics is highly nontrivial. A common approach, which we further develop in this work, is to consider instead a simpler effective impurity model that still captures the low-energy physics of interest. The mapping from a bare to an effective model is typically done perturbatively, but even this can be difficult for complex systems, and the resulting effective model parameters can nevertheless be quite inaccurate. Here we develop a nonperturbative, unsupervised machine learning approach to systematically obtain low-energy effective impurity-type models, based on the renormalization-group framework. The method is shown to be general and flexible, as well as accurate and systematically improvable. We benchmark the method against exact results for the Anderson impurity model, and we provide an outlook for more complex models beyond the reach of existing methods. Published by the American Physical Society 2024
iNORG: An open-source quantum impurity solver package based on the natural orbitals renormalization group
In the context of dynamical mean-field theory (DMFT) calculations for strongly correlated electron systems, quantum impurity solvers play a central computational role in treating correlated lattice models and realistic materials. Consequently, developing efficient and robust quantum impurity solvers remains a key challenge. In this paper, we present an open-source quantum impurity solver package based on the natural orbitals renormalization group (NORG) method, dubbed iNORG. This software delivers high accuracy with reduced computational cost by optimizing the bath representation using natural orbitals and incorporating advanced features such as efficient Hilbert space selection and efficient algorithms for computing Green's functions. We first introduce the basic principle of the NORG method and then discuss the implementation details. The software framework, major features, and installation procedure for iNORG are explained as well. Finally, several simple examples are presented to demonstrate the usage of iNORG.
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.
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.