Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2404.09444

Generalized Spectral Decomposition for Quantum Impurity Problems

Jun-Bin Wang, Dongchen Huang, Yi‐feng Yang

Published Apr 15, 2024• 0 citations
arXivHTMLPDF

Abstract

Solving quantum impurity problems may advance our understanding of strongly correlated electron physics, but its development in multi-impurity systems has been greatly hindered due to the presence of shared bath. Here, we propose a general operation strategy to disentangle the shared bath into multiple auxiliary baths and relate the problem to a spectral decomposition problem of function matrix for applying the numerical renormalization group (NRG). We prove exactly that such decomposition is possible for models satisfying (block) circulant symmetry, and show how to construct the auxiliary baths for arbitrary impurity configuration by mapping its graph structure to the subgraph of a regular impurity configuration. We further propose an approximate decomposition algorithm to reduce the number of auxiliary baths and save the computational workload. Our work reveals a deep connection between quantum impurity problems and the graph theory, and provides a general scheme to extend the NRG applications for realistic multi-impurity systems.

DOI
10.48550/arxiv.2404.09444

Related papers

Auxiliary-Bath Numerical Renormalization Group Method and Successive Collective Screening in Multi-Impurity Kondo Systems

Danqing Hu, Jiangfan Wang, Yi‐feng Yang

We propose an auxiliary-bath algorithm for the numerical renormalization group (NRG) method to solve multi-impurity models with shared electron baths. The method allows us to disentangle the electron baths into independent Wilson chains to perform standard NRG procedures beyond the widely adopted independent bath approximation. Its application to the 2-impurity model immediately reproduces the well-known even- and odd-parity channels. For 3-impurity Kondo models, we find collective screening of cluster degrees of freedom depending on impurity configurations and clarify the false prediction of a non-Fermi liquid ground state for the $C_3$ symmetric case in previous literature due to improper treatment of disentanglement. Our work highlights the importance of nonlocal spatial correlations due to shared baths and reveals a generic picture of successive collective screening for entropy depletion that is crucial in real correlated systems. Our method greatly expands the applicability of the NRG and opens an avenue for its further development.

arXiv:2305.10297•May 17, 2023•0 citations

Manual for the Flexible DM-NRG code

Örs Legeza, Cătălin Paşcu Moca, Anna I. Tóth +2 more

Quantum impurity models describe interactions between some local degrees of freedom and a continuum of non-interacting fermionic or bosonic states. The investigation of quantum impurity models is a starting point towards the understanding of more complex strongly correlated systems, but quantum impurity models also provide the description of various correlated mesoscopic structures, biological and chemical processes, atomic physics and describe phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson impurity model, or the single- and multi-channel Kondo models. The numerical renormalization group method (NRG) proposed by Wilson in mid 70's has been used in its original form for a longtime as one of the most accurate and powerful methods to deal with quatum impurity problems. Recently, a number of new developments took place: First, a spectral sum-conserving density matrix NRG approach (DM-NRG) has been developed, which has also been generalized for non-Abelian symmetries. In this manual we introduce some of the basic concepts of the NRG method and present recently developed Flexible DM-NRG code. This code uses user-defined non-Abelian symmetries dynamically, computes spectral functions, expectation values of local operators for user-defined impurity models. The code can also use a uniform density of states as well as a user-defined density of states. The current version of the code assumes fermionic bath's and it uses any number of U(1), SU(2) charge SU(2) or Z(2) symmetries. The Flexible DM-NRG code can be downloaded from http://www.phy.bme.hu/~dmnrg

arXiv:0809.3143•Sep 18, 2008•23 citations

Quantum impurities in channel mixing baths

Jin-Guo Liu, Da Wang, Qiang-Hua Wang

We propose a versatile strategy for numerical renormalization group (NRG) solution of general channel-mixing Kondo and Anderson impurity models beyond previous reach. We illustrate the strategy by investigating the quantum phase transitions in models of Anderson impurities coupled to $s$- and $d$-wave superconducting baths. We discuss the effects of nontrivial channel-mixing in such models. Our strategy opens the door toward broad applications of NRG as impurity solver in cluster dynamical mean field theory for strongly correlated electron systems.

arXiv:1509.01461•Physical review. B./Physical review. B•Jan 4, 2016•13 citations

Spectral functions of strongly correlated extended systems via an exact quantum embedding

George H. Booth, Garnet Kin‐Lic Chan

Density matrix embedding theory (DMET) [Phys. Rev. Lett. 109, 186404 (2012)], introduced an approach to quantum cluster embedding methods whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost but was inherently limited by the construction of a bath designed to reproduce ground-state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency-dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the one- (1D) and two-dimensional (2D) Hubbard model, where we obtain zero temperature and thermodynamic limit spectral functions, and show the trivial extension to two-particle Green's functions. This advance therefore extends the scope and applicability of DMET in condensed-matter problems as a computationally tractable route to correlated spectral functions of extended systems and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.

arXiv:1309.2320•Physical Review B•Apr 7, 2015•60 citations

Hybrid Hamiltonian-diagrammatic quantum impurity solver

Yang Yu, Gaurav Harsha, Lei Zhang +4 more

Quantum impurity models, which describe the coupling between interacting orbitals and a non-interacting bath, play a central role in the physics of strongly correlated electron systems. Solving a quantum impurity model in general requires the use of non-perturbative numerical methods. Hamiltonian-based approaches, which rely on an explicit bath discretization, are typically limited to a small number of bath sites or small entanglement, and diagrammatic methods suffer from sign problems, slow convergence, or diagram truncation approximations. Here we show that these two classes of methods can be combined: augmenting diagrammatic methods with a small auxiliary bath can reduce the residual problem to a regime where low-order perturbation theory is highly accurate and rapidly converging. In a simple benchmark, the precision of the hybrid approach surpasses bold-line calculations by several orders of magnitude; for a strongly interacting two-orbital model with a severe sign problem, convergence is achieved at three orders of magnitude lower computational cost than competing methods; and convergence to the unknown exact result is rapidly accelerated in a difficult realistic problem. Our results establish a practical route to high-precision quantum impurity solutions in correlated quantum systems.

arXiv:2606.11095•Jun 9, 2026•0 citations

Efficient pseudomode representation and complexity of quantum impurity models

Julian Thoenniss, Ilya Vilkoviskiy, Dmitry A. Abanin

Out-of-equilibrium fermionic quantum impurity models (QIM), describing a small interacting system coupled to a continuous fermionic bath, play an important role in condensed matter physics. Solving such models is a computationally demanding task, and a variety of computational approaches are based on finding approximate representations of the bath by a finite number of modes. In this paper, we formulate the problem of finding efficient bath representations as that of approximating a kernel of the bath's Feynman-Vernon influence functional by a sum of complex exponentials, with each term defining a fermionic pseudomode. Under mild assumptions on the analytic properties of the bath spectral density, we provide an analytic construction of pseudomodes, and prove that their number scales polylogarithmically with the maximum evolution time $T$ and the approximation error $\ensuremath{\varepsilon}$. We then demonstrate that the number of pseudomodes can be significantly reduced by an interpolative matrix decomposition (ID). Furthermore, we present a complementary approach, based on constructing rational approximations of the bath's spectral density using the ``AAA'' algorithm, followed by compression with ID. The combination of two approaches yields a pseudomode count scaling as ${N}_{\text{ID}}\ensuremath{\sim}ln(T)ln(1/\ensuremath{\varepsilon})$, and the agreement between the two approaches suggests that the result is close to optimal. Finally, to relate our findings to QIM, we derive an explicit Liouvillian that describes the time evolution of the combined impurity-pseudomodes system. These results establish bounds on the computational resources required for solving out-of-equilibrium QIMs, providing an efficient starting point for tensor-network methods for QIMs.

arXiv:2409.08816•Physical review. B./Physical review. B•Aug 13, 2025•5 citations
Solving quantum impurity problems may advance our understanding of strongly correlated electron physics, but its development in multi-impurity systems has been greatly hindered due to the presence of shared bath. Here, we propose a general operation strategy to disentangle the shared bath into multiple auxiliary baths and relate the problem to a spectral decomposition problem of function matrix for applying the numerical renormalization group (NRG). We prove exactly that such decomposition is possible for models satisfying (block) circulant symmetry, and show how to construct the auxiliary baths for arbitrary impurity configuration by mapping its graph structure to the subgraph of a regular impurity configuration. We further propose an approximate decomposition algorithm to reduce the number of auxiliary baths and save the computational workload. Our work reveals a deep connection between quantum impurity problems and the graph theory, and provides a general scheme to extend the NRG applications for realistic multi-impurity systems.
We propose an auxiliary-bath algorithm for the numerical renormalization group (NRG) method to solve multi-impurity models with shared electron baths. The method allows us to disentangle the electron baths into independent Wilson chains to perform standard NRG procedures beyond the widely adopted independent bath approximation. Its application to the 2-impurity model immediately reproduces the well-known even- and odd-parity channels. For 3-impurity Kondo models, we find collective screening of cluster degrees of freedom depending on impurity configurations and clarify the false prediction of a non-Fermi liquid ground state for the C3C_3C3​ symmetric case in previous literature due to improper treatment of disentanglement. Our work highlights the importance of nonlocal spatial correlations due to shared baths and reveals a generic picture of successive collective screening for entropy depletion that is crucial in real correlated systems. Our method greatly expands the applicability of the NRG and opens an avenue for its further development.
Quantum impurity models describe interactions between some local degrees of freedom and a continuum of non-interacting fermionic or bosonic states. The investigation of quantum impurity models is a starting point towards the understanding of more complex strongly correlated systems, but quantum impurity models also provide the description of various correlated mesoscopic structures, biological and chemical processes, atomic physics and describe phenomena such as dissipation or dephasing. Prototypes of these models are the Anderson impurity model, or the single- and multi-channel Kondo models. The numerical renormalization group method (NRG) proposed by Wilson in mid 70's has been used in its original form for a longtime as one of the most accurate and powerful methods to deal with quatum impurity problems. Recently, a number of new developments took place: First, a spectral sum-conserving density matrix NRG approach (DM-NRG) has been developed, which has also been generalized for non-Abelian symmetries. In this manual we introduce some of the basic concepts of the NRG method and present recently developed Flexible DM-NRG code. This code uses user-defined non-Abelian symmetries dynamically, computes spectral functions, expectation values of local operators for user-defined impurity models. The code can also use a uniform density of states as well as a user-defined density of states. The current version of the code assumes fermionic bath's and it uses any number of U(1), SU(2) charge SU(2) or Z(2) symmetries. The Flexible DM-NRG code can be downloaded from http://www.phy.bme.hu/~dmnrg
We propose a versatile strategy for numerical renormalization group (NRG) solution of general channel-mixing Kondo and Anderson impurity models beyond previous reach. We illustrate the strategy by investigating the quantum phase transitions in models of Anderson impurities coupled to sss- and ddd-wave superconducting baths. We discuss the effects of nontrivial channel-mixing in such models. Our strategy opens the door toward broad applications of NRG as impurity solver in cluster dynamical mean field theory for strongly correlated electron systems.
Density matrix embedding theory (DMET) [Phys. Rev. Lett. 109, 186404 (2012)], introduced an approach to quantum cluster embedding methods whereby the mapping of strongly correlated bulk problems to an impurity with finite set of bath states was rigorously formulated to exactly reproduce the entanglement of the ground state. The formalism provided similar physics to dynamical mean-field theory at a tiny fraction of the cost but was inherently limited by the construction of a bath designed to reproduce ground-state, static properties. Here, we generalize the concept of quantum embedding to dynamic properties and demonstrate accurate bulk spectral functions at similarly small computational cost. The proposed spectral DMET utilizes the Schmidt decomposition of a response vector, mapping the bulk dynamic correlation functions to that of a quantum impurity cluster coupled to a set of frequency-dependent bath states. The resultant spectral functions are obtained on the real-frequency axis, without bath discretization error, and allows for the construction of arbitrary dynamic correlation functions. We demonstrate the method on the one- (1D) and two-dimensional (2D) Hubbard model, where we obtain zero temperature and thermodynamic limit spectral functions, and show the trivial extension to two-particle Green's functions. This advance therefore extends the scope and applicability of DMET in condensed-matter problems as a computationally tractable route to correlated spectral functions of extended systems and provides a competitive alternative to dynamical mean-field theory for dynamic quantities.
Quantum impurity models, which describe the coupling between interacting orbitals and a non-interacting bath, play a central role in the physics of strongly correlated electron systems. Solving a quantum impurity model in general requires the use of non-perturbative numerical methods. Hamiltonian-based approaches, which rely on an explicit bath discretization, are typically limited to a small number of bath sites or small entanglement, and diagrammatic methods suffer from sign problems, slow convergence, or diagram truncation approximations. Here we show that these two classes of methods can be combined: augmenting diagrammatic methods with a small auxiliary bath can reduce the residual problem to a regime where low-order perturbation theory is highly accurate and rapidly converging. In a simple benchmark, the precision of the hybrid approach surpasses bold-line calculations by several orders of magnitude; for a strongly interacting two-orbital model with a severe sign problem, convergence is achieved at three orders of magnitude lower computational cost than competing methods; and convergence to the unknown exact result is rapidly accelerated in a difficult realistic problem. Our results establish a practical route to high-precision quantum impurity solutions in correlated quantum systems.
Out-of-equilibrium fermionic quantum impurity models (QIM), describing a small interacting system coupled to a continuous fermionic bath, play an important role in condensed matter physics. Solving such models is a computationally demanding task, and a variety of computational approaches are based on finding approximate representations of the bath by a finite number of modes. In this paper, we formulate the problem of finding efficient bath representations as that of approximating a kernel of the bath's Feynman-Vernon influence functional by a sum of complex exponentials, with each term defining a fermionic pseudomode. Under mild assumptions on the analytic properties of the bath spectral density, we provide an analytic construction of pseudomodes, and prove that their number scales polylogarithmically with the maximum evolution time TTT and the approximation error $\ensuremath{\varepsilon}$. We then demonstrate that the number of pseudomodes can be significantly reduced by an interpolative matrix decomposition (ID). Furthermore, we present a complementary approach, based on constructing rational approximations of the bath's spectral density using the ``AAA'' algorithm, followed by compression with ID. The combination of two approaches yields a pseudomode count scaling as ${N}_{\text{ID}}\ensuremath{\sim}ln(T)ln(1/\ensuremath{\varepsilon})$, and the agreement between the two approaches suggests that the result is close to optimal. Finally, to relate our findings to QIM, we derive an explicit Liouvillian that describes the time evolution of the combined impurity-pseudomodes system. These results establish bounds on the computational resources required for solving out-of-equilibrium QIMs, providing an efficient starting point for tensor-network methods for QIMs.