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:1103.5837 · Physical Review B

Spectral functions for single- and multi-impurity models using density matrix renormalization group

Robert Peters

Published Aug 11, 2011• 29 citations
arXivHTMLPDF

Abstract

This article focuses on the calculation of spectral functions for single- and multi-impurity models using the density matrix renormalization group (DMRG). To calculate spectral functions in DMRG, the correction vector method is presently the most widely used approach. One, however, always obtains Lorentzian convoluted spectral functions, which in applications like the dynamical mean-field theory can lead to wrong results. In order to overcome this restriction, we use chain decompositions in which the resulting effective Hamiltonian can be diagonalized completely to calculate a discrete ``peak'' spectrum. We show that this peak spectrum is a very good approximation to a deconvolution of the correction vector spectral function. Calculating this deconvoluted spectrum directly from the DMRG basis set and operators is the most natural approach, because it uses only information from the system itself. Having calculated this excitation spectrum, one can use an arbitrary broadening to obtain a smooth spectral function or directly analyze the excitations. As a nontrivial test, we apply this method to obtain spectral functions for a model of three coupled Anderson impurities. Although we are focusing in this article on impurity models, the proposed method for calculating the peak spectrum can be easily adapted to usual lattice models.

DOI
10.1103/physrevb.84.075139

Related papers

Dynamic Density-Matrix Renormalization for the Symmetric Single Impurity Anderson Model

Carsten Raas

In this thesis we investigated, implemented, gauged, optimized, and applied a numerical approach given by the dynamic density-matrix renormalization (D-DMRG). The dynamic density-matrix renormalization is the extension of standard DMRG to the calculation of dynamic quantities. D-DMRG provides valuable numerical information on dynamic correlations by computing convolutions of the corresponding spectral densities. The dynamics at zero temperature is determined by computing the expectation values in the local propagator. This can be realized by targeting not only at the ground state and the excited state, but also at the resolvent applied to the excited state. This additional targeted state is called the correction vector. The calculation of the frequency-dependent correction vector is numerically extremely demanding due to the inversion of an almost singular non-hermitian matrix. We compared the performance of several iterative solvers for linear equation systems and managed to stabilize the inversion problem of the D-DMRG by using optimized algorithms. The main limitation of the correction vector D-DMRG is that one cannot obtain data for purely real frequencies but only for frequencies with a certain imaginary part. Hence the extraction of the behavior at purely real frequencies is one of the main problems to be solved in using the D-DMRG. We illustrated how and to which extent such data can be deconvolved to retrieve the wanted spectral densities. We discussed and compared various algorithms to achieve this extraction and presented a non-linear approach from the family of maximum entropy methods. This approach provides a continuous, positive ansatz for the wanted spectral density with the least bias (LB). Even relative abrupt changes of the spectral density can be reproduced satisfactorily. In the vicinity of singularities spurious oscillations occur. The least-bias ansatz can be made more robust towards small numerical inaccuracies and finite-size effects by including besides the entropy functional a xi-functional in the functional to be minimized. We investigated the dynamic propagator of the symmetric single impurity Anderson model by D-DMRG. The Abrikosov-Suhl resonance (ASR) was calculated to gauge the D-DMRG and to demonstrate that features at low energies can be resolved. The low energy scale of the SIAM was reproduced over two orders of magnitude. The Hubbard satellites were investigated. From the broadened D-DMRG raw data we extracted the width and the shifts from the atomic positions. By means of the LB extraction we were able to address the line shape of the satellites. It was found that the Hubbard satellites are strongly asymmetric, the peak is very pronounced and the maximum value is very high. The Hubbard model on the Bethe lattice with infinite coordination number was investigated by means of the dynamic mean-field theory (DMFT) using the D-DMRG as impurity solver. A high-resolution investigation of the electron spectra close to the metal-to-insulator transition in dynamic mean-field theory was presented. The all-numerical, consistent confirmation of a smooth transition at zero temperature is provided.

•Jan 1, 2005•0 citations

Solving the multi-site and multi-orbital Dynamical Mean Field Theory using Density Matrix Renormalization

Y. Núñez-Fernández, K. Hallberg

We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.

arXiv:1711.08745•Nov 23, 2017•1 citation

Spectral functions with the density matrix renormalization group: Krylov-space approach for correction vectors

Alberto Nocera, Gonzalo Álvarez

Frequency-dependent correlations, such as the spectral function and the dynamical structure factor, help illustrate condensed matter experiments. Within the density matrix renormalization group (DMRG) framework, an accurate method for calculating spectral functions directly in frequency is the correction-vector method. The correction vector can be computed by solving a linear equation or by minimizing a functional. This paper proposes an alternative to calculate the correction vector: to use the Krylov-space approach. This paper then studies the accuracy and performance of the Krylov-space approach, when applied to the Heisenberg, the t-J, and the Hubbard models. The cases studied indicate that the Krylov-space approach can be more accurate and efficient than the conjugate gradient, and that the error of the former integrates best when a Krylov-space decomposition is also used for ground state DMRG.

arXiv:1607.03538•Physical review. E•Nov 21, 2016•91 citations

Solving the Multi-site and Multi-orbital Dynamical Mean Field Theory Using Density Matrix Renormalization

Y. Núñez-Fernández, K. Hallberg

We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.

Frontiers in Physics•Feb 26, 2018•9 citations

Dynamical density-matrix renormalization-group method

Eric Jeckelmann

A density-matrix renormalization-group (DMRG) method for calculating dynamical properties and excited states in low-dimensional lattice quantum many-body systems is presented. The method is based on an exact variational principle for dynamical correlation functions and the excited states contributing to them. This dynamical DMRG is an alternate formulation of the correction vector DMRG but is both simpler and more accurate. The finite-size scaling of spectral functions is discussed and a method for analyzing the scaling of dense spectra is described. The key idea of the method is a size-dependent broadening of the spectrum. The dynamical DMRG and the finite-size scaling analysis are demonstrated on the optical conductivity of the one-dimensional Peierls-Hubbard model. Comparisons with analytical results show that the spectral functions of infinite systems can be reproduced almost exactly with these techniques. The optical conductivity of the Mott-Peierls insulator is investigated and it is shown that its spectrum is qualitatively different from the simple spectra observed in Peierls (band) insulators and one-dimensional Mott-Hubbard insulators.

arXiv:cond-mat/0203500•Physical review. B, Condensed matter•Jul 26, 2002•374 citations

DMRG studies of Chebyshev-expanded spectral functions and quantum impurity models

Andreas Holzner

This thesis is concerned with two main topics: first, the advancement of the density matrix renormalization group (DMRG) and, second, its applications. In the first project of this thesis we exploit the common mathematical structure of the numerical renormalization group and the DMRG, namely, matrix product states (MPS), to implement an efficient numerical treatment of a two-lead, multi-level Anderson impurity model. By adopting a star-like geometry, where each species (spin and lead) of conduction electrons is described by its own so-called Wilson chain, instead of a single Wilson chain we achieve a very significant reduction in the numerical resources required to obtain reliable results. Moreover, we show that it is possible to find an "optimal" chain basis, in which chain degrees of freedom of different Wilson chains become effectively decoupled from each other further out on the Wilson chains. This basis turns out to also diagonalize the model's chain-to-chain scattering matrix. In the second project we show that Chebychev expansions offer numerically efficient representations for calculating spectral functions of one-dimensional lattice models using MPS methods. The main features of this Chebychev matrix product state (CheMPS) approach are: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it offers a well-controlled broadening scheme; (iii) it is based on using MPS tools to recursively calculate a succession of Chebychev vectors, (iv) whose entanglement entropies were found to remain bounded with increasing recursion order for all cases analyzed here. We present CheMPS results for the structure factor of spin-1/2 antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to benchmark methods, we find that CheMPS yields results comparable in quality to those of correction vector DMRG, at dramatically reduced numerical cost and agrees well with Bethe Ansatz results for an infinite system, within the limitations expected for numerics on finite systems. Following these technologically focused projects we study the so-called Kondo cloud by means of the DMRG in the third project. The Kondo cloud describes the effect of spatially extended spin-spin correlations of a magnetic moment and the conduction electrons which screen the magnetic moment through the Kondo effect at low temperatures. We focus on the question whether the Kondo screening length, typically assumed to be proportional to the inverse Kondo temperature, can be extracted from the spin-spin correlations. We investigate how perturbations which destroy the Kondo effect, like an applied gate potential or a magnetic field, affect the formation of the screening cloud. In a forth project we address the impact of Quantum (anti-)Zeno physics resulting from repeated single-site resolved observations on the many-body dynamics. We use time-dependent DMRG to obtain the time evolution of the full many-body wave function that is then periodically projected in order to simulate realizations of stroboscopic measurements. For the example of a 1-D lattice of spin-polarized fermions with nearest-neighbor interactions, we find regimes for which many-particle configurations are stabilized and destabilized depending on the interaction strength and the time between observations.

•Jan 1, 2012•1 citation
This article focuses on the calculation of spectral functions for single- and multi-impurity models using the density matrix renormalization group (DMRG). To calculate spectral functions in DMRG, the correction vector method is presently the most widely used approach. One, however, always obtains Lorentzian convoluted spectral functions, which in applications like the dynamical mean-field theory can lead to wrong results. In order to overcome this restriction, we use chain decompositions in which the resulting effective Hamiltonian can be diagonalized completely to calculate a discrete ``peak'' spectrum. We show that this peak spectrum is a very good approximation to a deconvolution of the correction vector spectral function. Calculating this deconvoluted spectrum directly from the DMRG basis set and operators is the most natural approach, because it uses only information from the system itself. Having calculated this excitation spectrum, one can use an arbitrary broadening to obtain a smooth spectral function or directly analyze the excitations. As a nontrivial test, we apply this method to obtain spectral functions for a model of three coupled Anderson impurities. Although we are focusing in this article on impurity models, the proposed method for calculating the peak spectrum can be easily adapted to usual lattice models.
In this thesis we investigated, implemented, gauged, optimized, and applied a numerical approach given by the dynamic density-matrix renormalization (D-DMRG). The dynamic density-matrix renormalization is the extension of standard DMRG to the calculation of dynamic quantities. D-DMRG provides valuable numerical information on dynamic correlations by computing convolutions of the corresponding spectral densities. The dynamics at zero temperature is determined by computing the expectation values in the local propagator. This can be realized by targeting not only at the ground state and the excited state, but also at the resolvent applied to the excited state. This additional targeted state is called the correction vector. The calculation of the frequency-dependent correction vector is numerically extremely demanding due to the inversion of an almost singular non-hermitian matrix. We compared the performance of several iterative solvers for linear equation systems and managed to stabilize the inversion problem of the D-DMRG by using optimized algorithms. The main limitation of the correction vector D-DMRG is that one cannot obtain data for purely real frequencies but only for frequencies with a certain imaginary part. Hence the extraction of the behavior at purely real frequencies is one of the main problems to be solved in using the D-DMRG. We illustrated how and to which extent such data can be deconvolved to retrieve the wanted spectral densities. We discussed and compared various algorithms to achieve this extraction and presented a non-linear approach from the family of maximum entropy methods. This approach provides a continuous, positive ansatz for the wanted spectral density with the least bias (LB). Even relative abrupt changes of the spectral density can be reproduced satisfactorily. In the vicinity of singularities spurious oscillations occur. The least-bias ansatz can be made more robust towards small numerical inaccuracies and finite-size effects by including besides the entropy functional a xi-functional in the functional to be minimized. We investigated the dynamic propagator of the symmetric single impurity Anderson model by D-DMRG. The Abrikosov-Suhl resonance (ASR) was calculated to gauge the D-DMRG and to demonstrate that features at low energies can be resolved. The low energy scale of the SIAM was reproduced over two orders of magnitude. The Hubbard satellites were investigated. From the broadened D-DMRG raw data we extracted the width and the shifts from the atomic positions. By means of the LB extraction we were able to address the line shape of the satellites. It was found that the Hubbard satellites are strongly asymmetric, the peak is very pronounced and the maximum value is very high. The Hubbard model on the Bethe lattice with infinite coordination number was investigated by means of the dynamic mean-field theory (DMFT) using the D-DMRG as impurity solver. A high-resolution investigation of the electron spectra close to the metal-to-insulator transition in dynamic mean-field theory was presented. The all-numerical, consistent confirmation of a smooth transition at zero temperature is provided.
We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.
Frequency-dependent correlations, such as the spectral function and the dynamical structure factor, help illustrate condensed matter experiments. Within the density matrix renormalization group (DMRG) framework, an accurate method for calculating spectral functions directly in frequency is the correction-vector method. The correction vector can be computed by solving a linear equation or by minimizing a functional. This paper proposes an alternative to calculate the correction vector: to use the Krylov-space approach. This paper then studies the accuracy and performance of the Krylov-space approach, when applied to the Heisenberg, the t-J, and the Hubbard models. The cases studied indicate that the Krylov-space approach can be more accurate and efficient than the conjugate gradient, and that the error of the former integrates best when a Krylov-space decomposition is also used for ground state DMRG.
We implement an efficient numerical method to calculate response functions of complex impurities based on the Density Matrix Renormalization Group (DMRG) and use it as the impurity-solver of the Dynamical Mean Field Theory (DMFT). This method uses the correction vector to obtain precise Green's functions on the real frequency axis at zero temperature. By using a self-consistent bath configuration with very low entanglement, we take full advantage of the DMRG to calculate dynamical response functions paving the way to treat large effective impurities such as those corresponding to multi-orbital interacting models and multi-site or multi-momenta clusters. This method leads to reliable calculations of non-local self energies at arbitrary dopings and interactions and at any energy scale.
A density-matrix renormalization-group (DMRG) method for calculating dynamical properties and excited states in low-dimensional lattice quantum many-body systems is presented. The method is based on an exact variational principle for dynamical correlation functions and the excited states contributing to them. This dynamical DMRG is an alternate formulation of the correction vector DMRG but is both simpler and more accurate. The finite-size scaling of spectral functions is discussed and a method for analyzing the scaling of dense spectra is described. The key idea of the method is a size-dependent broadening of the spectrum. The dynamical DMRG and the finite-size scaling analysis are demonstrated on the optical conductivity of the one-dimensional Peierls-Hubbard model. Comparisons with analytical results show that the spectral functions of infinite systems can be reproduced almost exactly with these techniques. The optical conductivity of the Mott-Peierls insulator is investigated and it is shown that its spectrum is qualitatively different from the simple spectra observed in Peierls (band) insulators and one-dimensional Mott-Hubbard insulators.
This thesis is concerned with two main topics: first, the advancement of the density matrix renormalization group (DMRG) and, second, its applications. In the first project of this thesis we exploit the common mathematical structure of the numerical renormalization group and the DMRG, namely, matrix product states (MPS), to implement an efficient numerical treatment of a two-lead, multi-level Anderson impurity model. By adopting a star-like geometry, where each species (spin and lead) of conduction electrons is described by its own so-called Wilson chain, instead of a single Wilson chain we achieve a very significant reduction in the numerical resources required to obtain reliable results. Moreover, we show that it is possible to find an "optimal" chain basis, in which chain degrees of freedom of different Wilson chains become effectively decoupled from each other further out on the Wilson chains. This basis turns out to also diagonalize the model's chain-to-chain scattering matrix. In the second project we show that Chebychev expansions offer numerically efficient representations for calculating spectral functions of one-dimensional lattice models using MPS methods. The main features of this Chebychev matrix product state (CheMPS) approach are: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it offers a well-controlled broadening scheme; (iii) it is based on using MPS tools to recursively calculate a succession of Chebychev vectors, (iv) whose entanglement entropies were found to remain bounded with increasing recursion order for all cases analyzed here. We present CheMPS results for the structure factor of spin-1/2 antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to benchmark methods, we find that CheMPS yields results comparable in quality to those of correction vector DMRG, at dramatically reduced numerical cost and agrees well with Bethe Ansatz results for an infinite system, within the limitations expected for numerics on finite systems. Following these technologically focused projects we study the so-called Kondo cloud by means of the DMRG in the third project. The Kondo cloud describes the effect of spatially extended spin-spin correlations of a magnetic moment and the conduction electrons which screen the magnetic moment through the Kondo effect at low temperatures. We focus on the question whether the Kondo screening length, typically assumed to be proportional to the inverse Kondo temperature, can be extracted from the spin-spin correlations. We investigate how perturbations which destroy the Kondo effect, like an applied gate potential or a magnetic field, affect the formation of the screening cloud. In a forth project we address the impact of Quantum (anti-)Zeno physics resulting from repeated single-site resolved observations on the many-body dynamics. We use time-dependent DMRG to obtain the time evolution of the full many-body wave function that is then periodically projected in order to simulate realizations of stroboscopic measurements. For the example of a 1-D lattice of spin-polarized fermions with nearest-neighbor interactions, we find regimes for which many-particle configurations are stabilized and destabilized depending on the interaction strength and the time between observations.