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

Using density matrix renormalization group to study open quantum systems

Cheng Guo

Published Jan 1, 2012• 3 citations
Publisher

Abstract

This thesis is concerned with open quantum systems, and more specifically, quantum impurity models. By this we mean a small local quantum system in contact with a macroscopic non-interacting environment. This can be used to model individual impurities in metals and quantum information systems where the influence from the surrounding environment is not negligible. The numerical renormalization group (NRG) is the traditional method to study quantum impurity models. However its application is limited when dealing with real-time dynamics and bosonic systems. In recent years some of NRG's techniques have been introduced to the density matrix renormalization group method (DMRG), which itself is the most powerful numerical method to study one-dimensional quantum systems. The resulting method shows great potential, and this thesis explores and extends the power of the NRG+DMRG combination in treating open quantum systems. We focus mainly on two types of problem. The first is an open quantum system with a time-dependent Hamiltonian, which for example could be the theoretical description of various problems encountered in qubit manipulation. We combine NRG discretization and adaptive time-dependent DMRG (t-DMRG) to study the dissipative Landau-Zener problem. We also use this method to study the quantum decoherence process and the dynamical properties of the telegraph noise model. The results show that the NRG and t-DMRG combination is a fast, accurate and versatile method for such problems. The second type of problem we study involves the quantum critical properties of one- and two-bath spin-boson models. Unlike fermion and spin models, models with bosons are difficult to treat numerically as each boson basis has an infinite number of dimensions. By introducing the optimal boson basis into the variational matrix product state method (VMPS), which is a variant of DMRG, we can deal with an effective local boson basis as large as $10^{10}$. This is the crucial improvement over NRG which can only deal with at most a few dozen local basis states. With this powerful tool we have settled a controversy about the nature of the quantum phase transition in a sub-Ohmic spin-boson model regarding quantum to classical mapping. There, NRG fails to yield the right result due to the highly truncated local boson basis. We also explore the phase diagram of the two-bath spin-boson model and find a new critical phase. We demonstrate that NRG+VMPS with optimal boson basis represents a powerful setting to study quantum impurity models with a bosonic environment.

DOI
10.5282/edoc.15169

Related papers

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

Real-time dynamics of quantum impurity systems in fermionic and bosonic environments

Kleine, Christian

Using the non-perturbative numerical renormalisation group (NRG) we analyse purely fermionic as well as more complex Bose-Fermi quantum impurity systems on their equilibrium properties and beyond: We present an extensive exploration of real-time dynamics for various quenches using the time-dependent NRG. For the purely fermionic system, we use the Single-impurity Anderson model (SIAM) in which the impurity is a localised orbital hybridising with a surrounding conduction band. The model enables charge scattering with a pseudo-gap scattering rate and features a quantum critical point separating the local moment (LM) phase with an unscreened spin moment from a symmetric strong coupling (SSC) phase with a fractionally screened spin moment. The real-time dynamics of the local double occupancy nicely thermalises for quenches regarding the SSC phase. Whereas for quenches within or into the LM phase the double occupancy systematically deviates from thermal equilibrium for long times. A comparison with results of a recently published time-dependent Gutwiller approach reveals clear discrepancies since this approach restricts the formation of the local moment only on the impurity site. We extend the system by an additional bosonic bath leading to the Bose-Fermi Anderson model (BFAM) and study the influence of the additional bosonic degrees of freedom on the real-time dynamics. In equilibrium the bosonic bath features an additional quantum critical point leading to a localised phase once the coupling to a sub-ohmic bosonic bath exceeds a critical value. In general the coupled bosonic bath leads to a renormalisation of the Coulomb interaction. We show that the dynamics is influenced by the specific choice of the bath exponents. For a F-type combination the non-equilibrium dynamics can be completely reproduced by an effective purely fermionic SIAM while for a B-type combination the nature of the bosonic bath and the coupling strength has a profound impact on the electron dynamics which can only be partially understood by an effective SIAM. Furthermore, we investigate the spin dynamics of an initially polarised impurity spin which is coupled to a fermionic and a bosonic bath leading to the Bose-Fermi Kondo model (BFKM). Through the coupling to the fermionic bath the polarisation vanishes with time. We show that in analogy to the spin-boson model (SBM) an initially coupling to the bosonic bath affects the real-time dynamics fundamentally leading to a remaining polarisation for long times.

•Jan 1, 2015•0 citations

Strongly Correlated Models in One Dimension from a Density Matrix Renormalization Group Perspective

Stefan Rommer

The density matrix renormalization group (DMRG) is a numerical method for studying low dimensional strongly correlated models. The study of these models is one of the most active areas of condensed matter physics today and has been spurred not only by the discovery of high temperature superconductivity but also by experiments on new compounds. In this thesis we study fundamentals of the density matrix renormalization group method as well as applications of the method to Heisenberg spin chains. The introductory text is meant as a brief overview of the DMRG as well as Heisenberg spin systems and provides an introduction to the papers. In papers I and II we investigate the density matrix renormalization group with the aim to better understand the underlying structure of the method. We find that if the renormalization converges to a fixed point in the thermodynamic limit, the resulting basis states can be written as a matrix product form. We also show an explicit relation of the DMRG to variational calculations using the fact that the matrix product states can be rederived through a variational ansatz without any reference to the DMRG calculation. It is then shown how the matrix product states can be used to define a set of Bloch waves for describing elementary excitations of the Heisenberg spin-1 chain. The low lying energy spectrum is calculated. In papers III and IV we investigate the thermodynamics of impurities in a spin-1/2 chain. For this purpose we have shown that a recently developed transfer matrix DMRG method can be used to calculate impurity properties. Both impurity corrections to thermodynamic properties and local properties at finite temperature can be obtained. The method gives quicker and more accurate results than quantum Monte Carlo. The numerical results show good agreement with field theory calculations for a number of impurity configurations and the cross-over function for the magnetic susceptibility from high- to low-temperature behavior shows the expected data collapse. Our results also give quantitative predictions for experiments on doped spin-1/2 chains.

•Jan 1, 1998•0 citations

Numerical renormalization group method for quantum impurity systems

R. Bulla, Theo A. Costi, Thomas Pruschke

In the early 1970s, Wilson developed the concept of a fully nonperturbative renormalization group transformation. When applied to the Kondo problem, this numerical renormalization group (NRG) method gave for the first time the full crossover from the high-temperature phase of a free spin to the low-temperature phase of a completely screened spin. The NRG method was later generalized to a variety of quantum impurity problems. The purpose of this review is to give a brief introduction to the NRG method, including some guidelines for calculating physical quantities, and to survey the development of the NRG method and its various applications over the last 30 years. These applications include variants of the original Kondo problem such as the non-Fermi-liquid behavior in the two-channel Kondo model, dissipative quantum systems such as the spin-boson model, and lattice systems in the framework of the dynamical mean-field theory.

arXiv:cond-mat/0701105•Reviews of Modern Physics•Apr 2, 2008•1,475 citations

Density matrix renormalization group study of a quantum impurity model with Landau-Zener time-dependent Hamiltonian

Cheng Guo, Andreas Weichselbaum, Stefan K. Kehrein +2 more

We use the adaptive time-dependent density matrix renormalization group method (t-DMRG) to study the nonequilibrium dynamics of a benchmark quantum impurity system which has a time-dependent Hamiltonian. This model is a resonant-level model, obtained by a mapping from a certain Ohmic spin-boson model describing the dissipative Landau-Zener transition. We map the resonant-level model onto a Wilson chain, then calculate the time-dependent occupation ${n}_{d}(t)$ of the resonant level. We compare t-DMRG results with exact results at zero temperature and find very good agreement. We also give a physical interpretation of the numerical results.

arXiv:0810.2900•Physical Review B•Mar 30, 2009•20 citations

Numerical renormalization group for quantum impurities in a bosonic bath

R. Bulla, Hyunjung Lee, Ning-Hua Tong +1 more

We present a detailed description of the recently proposed numerical renormalization group method for models of quantum impurities coupled to a bosonic bath. Specifically, the method is applied to the spin-boson model, both in the Ohmic and sub-Ohmic cases. We present various results for static as well as dynamic quantities and discuss details of the numerical implementation, e.g., the discretization of a bosonic bath with arbitrary continuous spectral density, the suitable choice of a finite basis in the bosonic Hilbert space, and questions of convergence with respect to truncation parameters. The method is shown to provide high-accuracy data over the whole range of model parameters and temperatures, which are in agreement with exact results and other numerical data from the literature.

arXiv:cond-mat/0407559•Physical Review B•Jan 25, 2005•174 citations
This thesis is concerned with open quantum systems, and more specifically, quantum impurity models. By this we mean a small local quantum system in contact with a macroscopic non-interacting environment. This can be used to model individual impurities in metals and quantum information systems where the influence from the surrounding environment is not negligible. The numerical renormalization group (NRG) is the traditional method to study quantum impurity models. However its application is limited when dealing with real-time dynamics and bosonic systems. In recent years some of NRG's techniques have been introduced to the density matrix renormalization group method (DMRG), which itself is the most powerful numerical method to study one-dimensional quantum systems. The resulting method shows great potential, and this thesis explores and extends the power of the NRG+DMRG combination in treating open quantum systems. We focus mainly on two types of problem. The first is an open quantum system with a time-dependent Hamiltonian, which for example could be the theoretical description of various problems encountered in qubit manipulation. We combine NRG discretization and adaptive time-dependent DMRG (t-DMRG) to study the dissipative Landau-Zener problem. We also use this method to study the quantum decoherence process and the dynamical properties of the telegraph noise model. The results show that the NRG and t-DMRG combination is a fast, accurate and versatile method for such problems. The second type of problem we study involves the quantum critical properties of one- and two-bath spin-boson models. Unlike fermion and spin models, models with bosons are difficult to treat numerically as each boson basis has an infinite number of dimensions. By introducing the optimal boson basis into the variational matrix product state method (VMPS), which is a variant of DMRG, we can deal with an effective local boson basis as large as 101010^{10}1010. This is the crucial improvement over NRG which can only deal with at most a few dozen local basis states. With this powerful tool we have settled a controversy about the nature of the quantum phase transition in a sub-Ohmic spin-boson model regarding quantum to classical mapping. There, NRG fails to yield the right result due to the highly truncated local boson basis. We also explore the phase diagram of the two-bath spin-boson model and find a new critical phase. We demonstrate that NRG+VMPS with optimal boson basis represents a powerful setting to study quantum impurity models with a bosonic environment.
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
Using the non-perturbative numerical renormalisation group (NRG) we analyse purely fermionic as well as more complex Bose-Fermi quantum impurity systems on their equilibrium properties and beyond: We present an extensive exploration of real-time dynamics for various quenches using the time-dependent NRG. For the purely fermionic system, we use the Single-impurity Anderson model (SIAM) in which the impurity is a localised orbital hybridising with a surrounding conduction band. The model enables charge scattering with a pseudo-gap scattering rate and features a quantum critical point separating the local moment (LM) phase with an unscreened spin moment from a symmetric strong coupling (SSC) phase with a fractionally screened spin moment. The real-time dynamics of the local double occupancy nicely thermalises for quenches regarding the SSC phase. Whereas for quenches within or into the LM phase the double occupancy systematically deviates from thermal equilibrium for long times. A comparison with results of a recently published time-dependent Gutwiller approach reveals clear discrepancies since this approach restricts the formation of the local moment only on the impurity site. We extend the system by an additional bosonic bath leading to the Bose-Fermi Anderson model (BFAM) and study the influence of the additional bosonic degrees of freedom on the real-time dynamics. In equilibrium the bosonic bath features an additional quantum critical point leading to a localised phase once the coupling to a sub-ohmic bosonic bath exceeds a critical value. In general the coupled bosonic bath leads to a renormalisation of the Coulomb interaction. We show that the dynamics is influenced by the specific choice of the bath exponents. For a F-type combination the non-equilibrium dynamics can be completely reproduced by an effective purely fermionic SIAM while for a B-type combination the nature of the bosonic bath and the coupling strength has a profound impact on the electron dynamics which can only be partially understood by an effective SIAM. Furthermore, we investigate the spin dynamics of an initially polarised impurity spin which is coupled to a fermionic and a bosonic bath leading to the Bose-Fermi Kondo model (BFKM). Through the coupling to the fermionic bath the polarisation vanishes with time. We show that in analogy to the spin-boson model (SBM) an initially coupling to the bosonic bath affects the real-time dynamics fundamentally leading to a remaining polarisation for long times.
The density matrix renormalization group (DMRG) is a numerical method for studying low dimensional strongly correlated models. The study of these models is one of the most active areas of condensed matter physics today and has been spurred not only by the discovery of high temperature superconductivity but also by experiments on new compounds. In this thesis we study fundamentals of the density matrix renormalization group method as well as applications of the method to Heisenberg spin chains. The introductory text is meant as a brief overview of the DMRG as well as Heisenberg spin systems and provides an introduction to the papers. In papers I and II we investigate the density matrix renormalization group with the aim to better understand the underlying structure of the method. We find that if the renormalization converges to a fixed point in the thermodynamic limit, the resulting basis states can be written as a matrix product form. We also show an explicit relation of the DMRG to variational calculations using the fact that the matrix product states can be rederived through a variational ansatz without any reference to the DMRG calculation. It is then shown how the matrix product states can be used to define a set of Bloch waves for describing elementary excitations of the Heisenberg spin-1 chain. The low lying energy spectrum is calculated. In papers III and IV we investigate the thermodynamics of impurities in a spin-1/2 chain. For this purpose we have shown that a recently developed transfer matrix DMRG method can be used to calculate impurity properties. Both impurity corrections to thermodynamic properties and local properties at finite temperature can be obtained. The method gives quicker and more accurate results than quantum Monte Carlo. The numerical results show good agreement with field theory calculations for a number of impurity configurations and the cross-over function for the magnetic susceptibility from high- to low-temperature behavior shows the expected data collapse. Our results also give quantitative predictions for experiments on doped spin-1/2 chains.
In the early 1970s, Wilson developed the concept of a fully nonperturbative renormalization group transformation. When applied to the Kondo problem, this numerical renormalization group (NRG) method gave for the first time the full crossover from the high-temperature phase of a free spin to the low-temperature phase of a completely screened spin. The NRG method was later generalized to a variety of quantum impurity problems. The purpose of this review is to give a brief introduction to the NRG method, including some guidelines for calculating physical quantities, and to survey the development of the NRG method and its various applications over the last 30 years. These applications include variants of the original Kondo problem such as the non-Fermi-liquid behavior in the two-channel Kondo model, dissipative quantum systems such as the spin-boson model, and lattice systems in the framework of the dynamical mean-field theory.
We use the adaptive time-dependent density matrix renormalization group method (t-DMRG) to study the nonequilibrium dynamics of a benchmark quantum impurity system which has a time-dependent Hamiltonian. This model is a resonant-level model, obtained by a mapping from a certain Ohmic spin-boson model describing the dissipative Landau-Zener transition. We map the resonant-level model onto a Wilson chain, then calculate the time-dependent occupation nd(t){n}_{d}(t)nd​(t) of the resonant level. We compare t-DMRG results with exact results at zero temperature and find very good agreement. We also give a physical interpretation of the numerical results.
We present a detailed description of the recently proposed numerical renormalization group method for models of quantum impurities coupled to a bosonic bath. Specifically, the method is applied to the spin-boson model, both in the Ohmic and sub-Ohmic cases. We present various results for static as well as dynamic quantities and discuss details of the numerical implementation, e.g., the discretization of a bosonic bath with arbitrary continuous spectral density, the suitable choice of a finite basis in the bosonic Hilbert space, and questions of convergence with respect to truncation parameters. The method is shown to provide high-accuracy data over the whole range of model parameters and temperatures, which are in agreement with exact results and other numerical data from the literature.