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:1111.5931

Ground State and Finite Temperature Lanczos Methods

Peter Prelovšek, J. Bonča

Published Nov 25, 2011• 17 citations
arXivHTMLPDF

Abstract

The present review will focus on recent development of exact-diagonali- zation (ED) methods that use Lanczos algorithm to transform large sparse matrices onto the tridiagonal form. We begin with a review of basic principles of the Lanczos method for computing ground-state static as well as dynamical properties. Next, generalization to finite-temperatures in the form of well established finite-temperature Lanczos method is described. The latter allows for the evaluation of temperatures T>0 static and dynamic quantities within various correlated models. Several extensions and modification of the latter method introduced more recently are analysed. In particular, the low-temperature Lanczos method and the microcanonical Lanczos method, especially applicable within the high-T regime. In order to overcome the problems of exponentially growing Hilbert spaces that prevent ED calculations on larger lattices, different approaches based on Lanczos diagonalization within the reduced basis have been developed. In this context, recently developed method based on ED within a limited functional space is reviewed. Finally, we briefly discuss the real-time evolution of correlated systems far from equilibrium, which can be simulated using the ED and Lanczos-based methods, as well as approaches based on the diagonalization in a reduced basis.

DOI
10.48550/arxiv.1111.5931

Related papers

Microcanonical Lanczos method

X. Zotos

In this article, the development of a novel numerical simulation method based on the Lanczos technique and the microcanonical statistical ensemble is discussed. It allows the evaluation of finite temperature dynamic correlations for quantum many-body systems on lattices larger than hitherto possible. To demonstrate the applicability of the method, a series of recently obtained results for one-dimensional spin and electronic models are reviewed.

The Philosophical Magazine A Journal of Theoretical Experimental and Applied Physics•Mar 10, 2006•3 citations

Exact Diagonalization Techniques for Quantum Spin Systems

Jürgen Schnack

The essence of quantum magnetism lies in correlated many-body states. They determine the physical properties beyond paramagnetism. Even constructions as simple as singlets or triplets cannot be understood without the mathematical concept of superpositions. This is even more evident if time evolution of quantum states and phenomena such as quantum coherenceQuantum coherence are considered. The major approach to assess correlated many-body states in quantum magnetism of molecules is exact diagonalizationDiagonalization of the many-body Hamiltonian describing the physical system. Symmetries are of great help when tackling the numerical task of matrix diagonalization since they allow to decompose a problem of large Hilbert spaceHilbert space dimension into a number of smaller problems. This approach will be presented first. It rests on the identification of the available symmetries—e.g., U(1), SU(2), as well as point groupsPoint group—and their application. However, the applicability of exact diagonalizationDiagonalization is limited by the exponential growth of the Hilbert spaceHilbert space of the spin system with the number of spins: $$\dim \left( {\mathcal{H}} \right) = \mathop \prod \nolimits_{i} \left( {2 s_{i} + 1} \right)$$ . Nowadays, complex Hermitian matricesHermitian matrices of a linear dimension of about 100,000 can be completely diagonalized on supercomputers which limits the number of spins depending on their size $$s_{i}$$ to just a few. I will therefore present a very accurate approximation to exact diagonalization that rests on Lanczos procedures and the concept of typicality. The most popular version, the finite-temperature Lanczos methodFinite-Temperature Lanczos Method (FTLM), allows to evaluate low-energy spectroscopic data, thermal averages as well as time evolutions for systems with Hilbert space dimensions of up to about 1011. Hands-on experience and examples will be provided.

Challenges and advances in computational chemistry and physics•Jan 1, 2023•2 citations

A SIMPLIFIED FORMULATION OF THE LANCZOS TECHNIQUE FOR STRONGLY CORRELATED SYSTEMS

S. Ehika

The basic idea of the Lanczos method is to construct a special basis where the Hamiltonian has a tridiagonal representation. Once in this form, the matrix can be diagonalized easily using standard library subroutines. This approach is known as the standard Lanczos technique (SLT) which still has the problem of diagonalizing large matrix sizes emanating from increase of the size of the Hilbert space with the size of the system. In this current presentation, we develop a simplified formulation of the SLT and then use it to study the rapid convergence to the ground state energy and wavefunction of some finite systems.

•Jan 1, 2009•0 citations

Low-temperature Lanczos method for strongly correlated systems

Markus Aichhorn, Maria Daghofer, Hans Gerd Evertz +1 more

We present a modified finite-temperature Lanczos method for the evaluation of dynamical and static quantities of strongly correlated electron systems that complements the finite-temperature method introduced by Jakli\ifmmode \check{c}\else \v{c}\fi{} and Prelov\ifmmode \check{s}\else \v{s}\fi{}ek for low temperatures. Together they allow accurate calculations at any temperature with moderate effort. As an example we calculate the static spin-correlation function and the regular part of the optical conductivity ${\ensuremath{\sigma}}^{\mathrm{reg}}(\ensuremath{\omega})$ of the one-dimensional Hubbard model at half filling and show in detail the connection between the ground-state and finite-temperature method. By using cluster perturbation theory, the finite-temperature spectral function is extended to the infinite system, clearly exhibiting the effects of spin-charge separation.

arXiv:cond-mat/0211165•Physical review. B, Condensed matter•Apr 24, 2003•81 citations

Finite-Temperature Variational Monte Carlo Method for Strongly Correlated Electron Systems

Kensaku Takai, Kota Ido, Takahiro Misawa +2 more

A new computational method for finite-temperature properties of strongly correlated electrons is proposed by extending the variational Monte Carlo method originally developed for the ground state. The method is based on the path integral in the imaginary-time formulation, starting from the infinite-temperature state that is well approximated by a small number of certain random initial states. Lower temperatures are progressively reached by the imaginary-time evolution. The algorithm follows the framework of the quantum transfer matrix and finite-temperature Lanczos methods, but we extends them to treat much larger system sizes without the negative sign problem by optimizing the truncated Hilbert space on the basis of the time-dependent variational principle (TDVP). This optimization algorithm is equivalent to the stochastic reconfiguration (SR) method that has been frequently used for the ground state to optimally truncate the Hilbert space. The obtained finite-temperature states allow an interpretation based on the thermal pure quantum (TPQ) state instead of the conventional canonical-ensemble average. Our method is tested for the one- and two-dimensional Hubbard models and its accuracy and efficiency are demonstrated.

arXiv:1510.05352•Journal of the Physical Society of Japan•Feb 18, 2016•21 citations

Diagonalization- and Numerical Renormalization-Group-Based Methods for Interacting Quantum Systems

R. M. Noack

In these lecture notes, we present a pedagogical review of a number of related numerically exact approaches to quantum many‐body problems. In particular, we focus on methods based on the exact diagonalization of the Hamiltonian matrix and on methods extending exact diagonalization using renormalization group ideas, i.e., Wilson’s Numerical Renormalization Group (NRG) and White’s Density Matrix Renormalization Group (DMRG). These methods are standard tools for the investigation of a variety of interacting quantum systems, especially low‐dimensional quantum lattice models. We also survey extensions to the methods to calculate properties such as dynamical quantities and behavior at finite temperature, and discuss generalizations of the DMRG method to a wider variety of systems, such as classical models and quantum chemical problems. Finally, we briefly review some recent developments for obtaining a more general formulation of the DMRG in the context of matrix product states as well as recent progress in calculating the time evolution of quantum systems using the DMRG and the relationship of the foundations of the method with quantum information theory.

arXiv:cond-mat/0510321•AIP conference proceedings•Jan 1, 2005•118 citations
The present review will focus on recent development of exact-diagonali- zation (ED) methods that use Lanczos algorithm to transform large sparse matrices onto the tridiagonal form. We begin with a review of basic principles of the Lanczos method for computing ground-state static as well as dynamical properties. Next, generalization to finite-temperatures in the form of well established finite-temperature Lanczos method is described. The latter allows for the evaluation of temperatures T>0 static and dynamic quantities within various correlated models. Several extensions and modification of the latter method introduced more recently are analysed. In particular, the low-temperature Lanczos method and the microcanonical Lanczos method, especially applicable within the high-T regime. In order to overcome the problems of exponentially growing Hilbert spaces that prevent ED calculations on larger lattices, different approaches based on Lanczos diagonalization within the reduced basis have been developed. In this context, recently developed method based on ED within a limited functional space is reviewed. Finally, we briefly discuss the real-time evolution of correlated systems far from equilibrium, which can be simulated using the ED and Lanczos-based methods, as well as approaches based on the diagonalization in a reduced basis.
The essence of quantum magnetism lies in correlated many-body states. They determine the physical properties beyond paramagnetism. Even constructions as simple as singlets or triplets cannot be understood without the mathematical concept of superpositions. This is even more evident if time evolution of quantum states and phenomena such as quantum coherenceQuantum coherence are considered. The major approach to assess correlated many-body states in quantum magnetism of molecules is exact diagonalizationDiagonalization of the many-body Hamiltonian describing the physical system. Symmetries are of great help when tackling the numerical task of matrix diagonalization since they allow to decompose a problem of large Hilbert spaceHilbert space dimension into a number of smaller problems. This approach will be presented first. It rests on the identification of the available symmetries—e.g., U(1), SU(2), as well as point groupsPoint group—and their application. However, the applicability of exact diagonalizationDiagonalization is limited by the exponential growth of the Hilbert spaceHilbert space of the spin system with the number of spins: $$\dim \left( {\mathcal{H}} \right) = \mathop \prod \nolimits_{i} \left( {2 s_{i} + 1} \right)$$ . Nowadays, complex Hermitian matricesHermitian matrices of a linear dimension of about 100,000 can be completely diagonalized on supercomputers which limits the number of spins depending on their size sis_{i}si​ to just a few. I will therefore present a very accurate approximation to exact diagonalization that rests on Lanczos procedures and the concept of typicality. The most popular version, the finite-temperature Lanczos methodFinite-Temperature Lanczos Method (FTLM), allows to evaluate low-energy spectroscopic data, thermal averages as well as time evolutions for systems with Hilbert space dimensions of up to about 1011. Hands-on experience and examples will be provided.
The basic idea of the Lanczos method is to construct a special basis where the Hamiltonian has a tridiagonal representation. Once in this form, the matrix can be diagonalized easily using standard library subroutines. This approach is known as the standard Lanczos technique (SLT) which still has the problem of diagonalizing large matrix sizes emanating from increase of the size of the Hilbert space with the size of the system. In this current presentation, we develop a simplified formulation of the SLT and then use it to study the rapid convergence to the ground state energy and wavefunction of some finite systems.
We present a modified finite-temperature Lanczos method for the evaluation of dynamical and static quantities of strongly correlated electron systems that complements the finite-temperature method introduced by Jakli\ifmmode \check{c}\else \v{c}\fi{} and Prelov\ifmmode \check{s}\else \v{s}\fi{}ek for low temperatures. Together they allow accurate calculations at any temperature with moderate effort. As an example we calculate the static spin-correlation function and the regular part of the optical conductivity ${\ensuremath{\sigma}}^{\mathrm{reg}}(\ensuremath{\omega})$ of the one-dimensional Hubbard model at half filling and show in detail the connection between the ground-state and finite-temperature method. By using cluster perturbation theory, the finite-temperature spectral function is extended to the infinite system, clearly exhibiting the effects of spin-charge separation.
A new computational method for finite-temperature properties of strongly correlated electrons is proposed by extending the variational Monte Carlo method originally developed for the ground state. The method is based on the path integral in the imaginary-time formulation, starting from the infinite-temperature state that is well approximated by a small number of certain random initial states. Lower temperatures are progressively reached by the imaginary-time evolution. The algorithm follows the framework of the quantum transfer matrix and finite-temperature Lanczos methods, but we extends them to treat much larger system sizes without the negative sign problem by optimizing the truncated Hilbert space on the basis of the time-dependent variational principle (TDVP). This optimization algorithm is equivalent to the stochastic reconfiguration (SR) method that has been frequently used for the ground state to optimally truncate the Hilbert space. The obtained finite-temperature states allow an interpretation based on the thermal pure quantum (TPQ) state instead of the conventional canonical-ensemble average. Our method is tested for the one- and two-dimensional Hubbard models and its accuracy and efficiency are demonstrated.
In these lecture notes, we present a pedagogical review of a number of related numerically exact approaches to quantum many‐body problems. In particular, we focus on methods based on the exact diagonalization of the Hamiltonian matrix and on methods extending exact diagonalization using renormalization group ideas, i.e., Wilson’s Numerical Renormalization Group (NRG) and White’s Density Matrix Renormalization Group (DMRG). These methods are standard tools for the investigation of a variety of interacting quantum systems, especially low‐dimensional quantum lattice models. We also survey extensions to the methods to calculate properties such as dynamical quantities and behavior at finite temperature, and discuss generalizations of the DMRG method to a wider variety of systems, such as classical models and quantum chemical problems. Finally, we briefly review some recent developments for obtaining a more general formulation of the DMRG in the context of matrix product states as well as recent progress in calculating the time evolution of quantum systems using the DMRG and the relationship of the foundations of the method with quantum information theory.