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

The Journal of Physical Chemistry Letters

Chebyshev Matrix Product States with Canonical Orthogonalization for Spectral Functions of Many-Body Systems

Tong Jiang, Jiajun Ren, Zhigang Shuai

Published Sep 22, 2021• 11 citations
Publisher

Abstract

We propose a method to calculate the spectral functions of many-body systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigenenergies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. For the first time we demonstrate that Chebyshev MPS can be used to deal with ab initio electronic Hamiltonian effectively. We emphasize the strength of coCheMPS to calculate the low excited states of systems with sparse discrete spectrum. We also caution the application for electron-phonon systems with dense density of states.

DOI
10.1021/acs.jpclett.1c02688

Related papers

Chebyshev matrix product states with canonical orthogonalization for spectral functions of strongly correlated systems

Tong Jiang, Jiajun Ren, Zhigang Shuai

We propose a method to calculate the spectral functions of strongly correlated systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigen energies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. We demonstrate for the first time that Chebyshev MPS can be used in quantum chemistry. We also caution the application for electron-phonon system with densed density of states.

•Sep 6, 2021•0 citations

Chebyshev matrix product state approach for spectral functions

Andreas Holzner, Andreas Weichselbaum, Ian P. McCulloch +2 more

We show that recursively generated Chebyshev expansions offer numerically efficient representations for calculating zero-temperature spectral functions of one-dimensional lattice models using matrix product state (MPS) methods. The main features of this Chebyshev matrix product state (CheMPS) approach are as follows: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it can exploit the fact that the latter can be much smaller than the model's many-body bandwidth; (iii) it offers a well-controlled broadening scheme that allows finite-size effects to be either resolved or smeared out, as desired; (iv) it is based on using MPS tools to recursively calculate a succession of Chebyshev vectors $|{t}_{n}\ensuremath{\rangle}$, (v) the entanglement entropies of which were found to remain bounded with increasing recursion order $n$ for all cases analyzed here; and (vi) it distributes the total entanglement entropy that accumulates with increasing $n$ over the set of Chebyshev vectors $|{t}_{n}\ensuremath{\rangle}$, which need not be combined into a single vector. In this way, the growth in entanglement entropy that usually limits density matrix renormalization group (DMRG) approaches is packaged into conveniently manageable units. We present zero-temperature CheMPS results for the structure factor of spin-$\frac{1}{2}$ antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to three benchmark methods, we find that CheMPS (a) yields results comparable in quality to those of correction-vector DMRG, at dramatically reduced numerical cost; (b) agrees well with Bethe ansatz results for an infinite system, within the limitations expected for numerics on finite systems; and (c) can also be applied in the time domain, where it has potential to serve as a viable alternative to time-dependent DMRG (in particular, at finite temperatures). Finally, we present a detailed error analysis of CheMPS for the case of the noninteracting resonant level model.

arXiv:1101.5895•Physical Review B•May 10, 2011•132 citations

Reorthonormalization of Chebyshev matrix product states for dynamical correlation functions

Haidong Xie, Rui-Zhen Huang, Xiaoyang Han +5 more

The Chebyshev expansion offers a numerically efficient and easy-implement algorithm for evaluating dynamic correlation functions using matrix product states (MPS). In this approach, each recursively generated Chebyshev vector is approximately represented by an MPS. However, the recurrence relations of Chebyshev polynomials are broken by the approximation, leading to an error which is accumulated with the increase of the order of expansion. Here we propose a reorthonormalization approach to remove this error introduced in the loss of orthogonality of the Chebyshev polynomials. Our approach, as illustrated by comparison with the exact results for the one-dimensional XY and Heisenberg models, improves significantly the accuracy in the calculation of dynamical correlation functions.

arXiv:1711.09207•Physical review. B./Physical review. B•Feb 7, 2018•27 citations

Chebyshev expansion for impurity models using matrix product states

Martin Ganahl, Patrik Thunström, Frank Verstraete +2 more

We improve a recently developed expansion technique for calculating real-frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost, and in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory. Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group, at lower cost. We also propose an alternative expansion, of $\mathbb{1}\ensuremath{-}exp(\ensuremath{-}\ensuremath{\tau}H)$ instead of the usual $H$, which opens the possibility of using established methods for the time evolution of matrix product states to calculate the spectral functions directly.

arXiv:1403.1209•Physical Review B•Jul 30, 2014•80 citations

Chebyshev pseudosite matrix product state approach for the spectral functions of electron-phonon coupling systems

Pei-Yuan Zhao, Ke Ding, Shuo Yang

The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ coupling system often has a large number of phonon degrees of freedom, whose spectral functions are numerically difficult to compute using matrix product state (MPS) formalisms. To solve this problem, we propose a practical method that combines the Chebyshev MPS and the pseudosite density matrix renormalization group (DMRG) algorithm. The Chebyshev vector is represented by a pseudosite MPS with global $U(1)$ fermion symmetry, which maps ${2}^{{N}_{p}}$ bosonic degrees of freedom onto ${N}_{p}$ pseudosites, each with two states. This approach can handle arbitrary $e$-ph coupling Hamiltonians where pseudosite DMRG performs efficiently. We use this method to study the spectral functions of the doped extended Hubbard-Holstein model in a regime of strong Coulomb repulsion. Key features of the excitation spectra are captured at a modest computational cost. Our results show that weak extended $e$-ph couplings can increase the spectral weight of the holon-folding branch at low phonon frequencies, in agreement with angle-resolved photoemission observations on one-dimensional cuprates.

arXiv:2210.09208•Physical Review Research•Apr 14, 2023•7 citations

Chebyshev matrix product state approach for time evolution

Jad C. Halimeh, Fabian Kolley, Ian P. McCulloch

We present and test a new algorithm for time-evolving quantum many-body systems initially proposed by Holzner et al. [Phys. Rev. B 83, 195115 (2011)]. The approach is based on merging the matrix product state (MPS) formalism with the method of expanding the time-evolution operator in Chebyshev polynomials. We calculate time-dependent observables of a system of hardcore bosons quenched under the Bose-Hubbard Hamiltonian on a one-dimensional lattice. We compare the new algorithm to more standard methods using the MPS architecture. We find that the Chebyshev method gives numerically exact results for small times. However, the reachable times are smaller than the ones obtained with the other state-of-the-art methods. We further extend the new method using a spectral-decomposition-based projective scheme that utilizes an effective bandwidth significantly smaller than the full bandwidth, leading to longer evolution times than the nonprojective method and more efficient information storage, data compression, and less computational effort.

arXiv:1507.01226•Physical Review B•Sep 14, 2015•31 citations
We propose a method to calculate the spectral functions of many-body systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigenenergies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. For the first time we demonstrate that Chebyshev MPS can be used to deal with ab initio electronic Hamiltonian effectively. We emphasize the strength of coCheMPS to calculate the low excited states of systems with sparse discrete spectrum. We also caution the application for electron-phonon systems with dense density of states.
We propose a method to calculate the spectral functions of strongly correlated systems by Chebyshev expansion in the framework of matrix product states coupled with canonical orthogonalization (coCheMPS). The canonical orthogonalization can improve the accuracy and efficiency significantly because the orthogonalized Chebyshev vectors can provide an ideal basis for constructing the effective Hamiltonian in which the exact recurrence relation can be retained. In addition, not only the spectral function but also the excited states and eigen energies can be directly calculated, which is usually impossible for other MPS-based methods such as time-dependent formalism or correction vector. The remarkable accuracy and efficiency of coCheMPS over other methods are demonstrated by calculating the spectral functions of spin chain and ab initio hydrogen chain. We demonstrate for the first time that Chebyshev MPS can be used in quantum chemistry. We also caution the application for electron-phonon system with densed density of states.
We show that recursively generated Chebyshev expansions offer numerically efficient representations for calculating zero-temperature spectral functions of one-dimensional lattice models using matrix product state (MPS) methods. The main features of this Chebyshev matrix product state (CheMPS) approach are as follows: (i) it achieves uniform resolution over the spectral function's entire spectral width; (ii) it can exploit the fact that the latter can be much smaller than the model's many-body bandwidth; (iii) it offers a well-controlled broadening scheme that allows finite-size effects to be either resolved or smeared out, as desired; (iv) it is based on using MPS tools to recursively calculate a succession of Chebyshev vectors $|{t}_{n}\ensuremath{\rangle}$, (v) the entanglement entropies of which were found to remain bounded with increasing recursion order nnn for all cases analyzed here; and (vi) it distributes the total entanglement entropy that accumulates with increasing nnn over the set of Chebyshev vectors $|{t}_{n}\ensuremath{\rangle}$, which need not be combined into a single vector. In this way, the growth in entanglement entropy that usually limits density matrix renormalization group (DMRG) approaches is packaged into conveniently manageable units. We present zero-temperature CheMPS results for the structure factor of spin-12\frac{1}{2}21​ antiferromagnetic Heisenberg chains and perform a detailed finite-size analysis. Making comparisons to three benchmark methods, we find that CheMPS (a) yields results comparable in quality to those of correction-vector DMRG, at dramatically reduced numerical cost; (b) agrees well with Bethe ansatz results for an infinite system, within the limitations expected for numerics on finite systems; and (c) can also be applied in the time domain, where it has potential to serve as a viable alternative to time-dependent DMRG (in particular, at finite temperatures). Finally, we present a detailed error analysis of CheMPS for the case of the noninteracting resonant level model.
The Chebyshev expansion offers a numerically efficient and easy-implement algorithm for evaluating dynamic correlation functions using matrix product states (MPS). In this approach, each recursively generated Chebyshev vector is approximately represented by an MPS. However, the recurrence relations of Chebyshev polynomials are broken by the approximation, leading to an error which is accumulated with the increase of the order of expansion. Here we propose a reorthonormalization approach to remove this error introduced in the loss of orthogonality of the Chebyshev polynomials. Our approach, as illustrated by comparison with the exact results for the one-dimensional XY and Heisenberg models, improves significantly the accuracy in the calculation of dynamical correlation functions.
We improve a recently developed expansion technique for calculating real-frequency spectral functions of any one-dimensional model with short-range interactions, by postprocessing computed Chebyshev moments with linear prediction. This can be achieved at virtually no cost, and in sharp contrast to existing methods based on the dampening of the moments, improves the spectral resolution rather than lowering it. We validate the method for the exactly solvable resonating level model and the single impurity Anderson model. It is capable of resolving sharp Kondo resonances, as well as peaks within the Hubbard bands when employed as an impurity solver for dynamical mean-field theory. Our method works at zero temperature and allows for arbitrary discretization of the bath spectrum. It achieves similar precision as the dynamical density matrix renormalization group, at lower cost. We also propose an alternative expansion, of $\mathbb{1}\ensuremath{-}exp(\ensuremath{-}\ensuremath{\tau}H)$ instead of the usual HHH, which opens the possibility of using established methods for the time evolution of matrix product states to calculate the spectral functions directly.
The electron-phonon $(e\text{\ensuremath{-}}\mathrm{ph})$ coupling system often has a large number of phonon degrees of freedom, whose spectral functions are numerically difficult to compute using matrix product state (MPS) formalisms. To solve this problem, we propose a practical method that combines the Chebyshev MPS and the pseudosite density matrix renormalization group (DMRG) algorithm. The Chebyshev vector is represented by a pseudosite MPS with global U(1)U(1)U(1) fermion symmetry, which maps 2Np{2}^{{N}_{p}}2Np​ bosonic degrees of freedom onto Np{N}_{p}Np​ pseudosites, each with two states. This approach can handle arbitrary eee-ph coupling Hamiltonians where pseudosite DMRG performs efficiently. We use this method to study the spectral functions of the doped extended Hubbard-Holstein model in a regime of strong Coulomb repulsion. Key features of the excitation spectra are captured at a modest computational cost. Our results show that weak extended eee-ph couplings can increase the spectral weight of the holon-folding branch at low phonon frequencies, in agreement with angle-resolved photoemission observations on one-dimensional cuprates.
We present and test a new algorithm for time-evolving quantum many-body systems initially proposed by Holzner et al. [Phys. Rev. B 83, 195115 (2011)]. The approach is based on merging the matrix product state (MPS) formalism with the method of expanding the time-evolution operator in Chebyshev polynomials. We calculate time-dependent observables of a system of hardcore bosons quenched under the Bose-Hubbard Hamiltonian on a one-dimensional lattice. We compare the new algorithm to more standard methods using the MPS architecture. We find that the Chebyshev method gives numerically exact results for small times. However, the reachable times are smaller than the ones obtained with the other state-of-the-art methods. We further extend the new method using a spectral-decomposition-based projective scheme that utilizes an effective bandwidth significantly smaller than the full bandwidth, leading to longer evolution times than the nonprojective method and more efficient information storage, data compression, and less computational effort.