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

Adaptive Lanczos-vector method for dynamic properties within the density matrix renormalization group

Piet E. Dargel, A. Honecker, Robert Peters, R. M. Noack, Thomas Pruschke

Published Apr 22, 2011• 37 citations
arXivHTMLPDF

Abstract

Current widely used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.

DOI
10.1103/physrevb.83.161104

Related papers

Adaptive Lanczos-vector method for dynamic properties within the DMRG

Piet E. Dargel, A. Honecker, Robert Peters +2 more

Current widely-used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.

•Jan 1, 2011•0 citations

Dynamical correlation functions using the density matrix renormalization group

Till D. Kühner, Steven R. White

The density matrix renormalization group (DMRG) method allows for very precise calculations of ground state properties in low-dimensional strongly correlated systems. We investigate two methods to expand the DMRG to calculations of dynamical properties. In the Lanczos vector method the DMRG basis is optimized to represent Lanczos vectors, which are then used to calculate the spectra. This method is fast and relatively easy to implement, but the accuracy at higher frequencies is limited. Alternatively, one can optimize the basis to represent a correction vector for a particular frequency. The correction vectors can be used to calculate the dynamical correlation functions at these frequencies with high accuracy. By separately calculating correction vectors at different frequencies, the dynamical correlation functions can be interpolated and pieced together from these results. For systems with open boundaries we discuss how to construct operators for specific wave vectors using filter functions.

arXiv:cond-mat/9812372•Physical review. B, Condensed matter•Jul 1, 1999•288 citations

Lanczos algorithm with matrix product states for dynamical correlation functions

Piet E. Dargel, Anton Wöllert, A. Honecker +3 more

The density-matrix renormalization group (DMRG) algorithm can be adapted to the calculation of dynamical correlation functions in various ways which all represent compromises between computational efficiency and physical accuracy. In this paper we reconsider the oldest approach based on a suitable Lanczos-generated approximate basis and implement it using matrix product states (MPS) for the representation of the basis states. The direct use of matrix product states combined with an ex post reorthogonalization method allows us to avoid several shortcomings of the original approach, namely the multitargeting and the approximate representation of the Hamiltonian inherent in earlier Lanczos-method implementations in the DMRG framework, and to deal with the ghost problem of Lanczos methods, leading to a much better convergence of the spectral weights and poles. We present results for the dynamic spin structure factor of the spin-1/2 antiferromagnetic Heisenberg chain. A comparison to Bethe ansatz results in the thermodynamic limit reveals that the MPS-based Lanczos approach is much more accurate than earlier approaches at minor additional numerical cost.

arXiv:1203.2523•Physical Review B•May 11, 2012•72 citations

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

Root- N Krylov-space correction vectors for spectral functions with the density matrix renormalization group

Alberto Nocera, G. Alvarez

We propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov-space decomposition to compute the correction vectors. Our approach entails the calculation of the root-$N$ ($N=2$ is the standard square root) of the Hamiltonian propagator using Krylov-space decomposition and repeating this procedure $N$ times to obtain the actual correction vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, whereas achieving better computational performance at large $N$. We apply our method to spin and charge spectral functions of $t\text{\ensuremath{-}}J$ and Hubbard models in the challenging two-leg ladder geometry and provide evidence that the root-$N$ approach reaches a much improved spectral resolution compared to the conventional correction vector.

arXiv:2204.03165•Physical review. B./Physical review. B•Nov 3, 2022•15 citations

Density-matrix algorithm for the calculation of dynamical properties of low-dimensional systems

K. Hallberg

I extend the scope of the density-matrix renormalization-group technique developed by White to the calculation of dynamical correlation functions using the continuous fraction expansion of the Green's function. As an application and performance evaluation I calculate the spin dynamics of the one-dimensional S=1/2 Heisenberg chain.

arXiv:cond-mat/9503094•Physical review. B, Condensed matter•Oct 1, 1995•191 citations
Current widely used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.
Current widely-used approaches to calculate spectral functions using the density-matrix renormalization group in frequency space either necessarily include an artificial broadening (correction-vector method), have limited resolution (time-domain density-matrix renormalization group with Fourier transform method), or are limited to low-energy properties or single dominant modes (original continued fraction method). Here we propose an adaptive Lanczos-vector method to calculate the coefficients of a continued fraction expansion of the spectral function iteratively. We show that one can obtain a very accurate representation of the spectral function very efficiently, and that one can also directly extract the spectral weights and poles for the discrete system.
The density matrix renormalization group (DMRG) method allows for very precise calculations of ground state properties in low-dimensional strongly correlated systems. We investigate two methods to expand the DMRG to calculations of dynamical properties. In the Lanczos vector method the DMRG basis is optimized to represent Lanczos vectors, which are then used to calculate the spectra. This method is fast and relatively easy to implement, but the accuracy at higher frequencies is limited. Alternatively, one can optimize the basis to represent a correction vector for a particular frequency. The correction vectors can be used to calculate the dynamical correlation functions at these frequencies with high accuracy. By separately calculating correction vectors at different frequencies, the dynamical correlation functions can be interpolated and pieced together from these results. For systems with open boundaries we discuss how to construct operators for specific wave vectors using filter functions.
The density-matrix renormalization group (DMRG) algorithm can be adapted to the calculation of dynamical correlation functions in various ways which all represent compromises between computational efficiency and physical accuracy. In this paper we reconsider the oldest approach based on a suitable Lanczos-generated approximate basis and implement it using matrix product states (MPS) for the representation of the basis states. The direct use of matrix product states combined with an ex post reorthogonalization method allows us to avoid several shortcomings of the original approach, namely the multitargeting and the approximate representation of the Hamiltonian inherent in earlier Lanczos-method implementations in the DMRG framework, and to deal with the ghost problem of Lanczos methods, leading to a much better convergence of the spectral weights and poles. We present results for the dynamic spin structure factor of the spin-1/2 antiferromagnetic Heisenberg chain. A comparison to Bethe ansatz results in the thermodynamic limit reveals that the MPS-based Lanczos approach is much more accurate than earlier approaches at minor additional numerical cost.
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 propose a method to compute spectral functions of generic Hamiltonians using the density matrix renormalization group (DMRG) algorithm directly in the frequency domain, based on a modified Krylov-space decomposition to compute the correction vectors. Our approach entails the calculation of the root-NNN (N=2N=2N=2 is the standard square root) of the Hamiltonian propagator using Krylov-space decomposition and repeating this procedure NNN times to obtain the actual correction vector. We show that our method greatly alleviates the burden of keeping a large bond dimension at large target frequencies, a problem found with conventional correction-vector DMRG, whereas achieving better computational performance at large NNN. We apply our method to spin and charge spectral functions of $t\text{\ensuremath{-}}J$ and Hubbard models in the challenging two-leg ladder geometry and provide evidence that the root-NNN approach reaches a much improved spectral resolution compared to the conventional correction vector.