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:1901.05793 · Physical review. B./Physical review. B

Matrix product states approaches to operator spreading in ergodic quantum systems

Kévin Hémery, Frank Pollmann, David J. Luitz

Published Sep 11, 2019• 27 citations
arXivHTMLPDF

Abstract

We review different matrix-product-state (MPS) approaches to study the spreading of operators in generic nonintegrable quantum systems. As a common ground to all methods, we quantify this spreading by means of the Frobenius norm of the commutator of a spreading operator with a local operator, which is usually referred to as the out-of-time-order correlation (OTOC) function. We compare two approaches based on matrix-product states in the Schr\"odinger picture: the time-dependent block decimation (TEBD) and the time-dependent variational principle (TDVP), as well as TEBD based on matrix-product operators directly in the Heisenberg picture. The results of all methods are compared to numerically exact results using Krylov space exact time evolution. We find that for the Schr\"odinger picture, the TDVP algorithm performs better than the TEBD algorithm.Moreover, the tails of the OTOC are accurately obtained both by TDVP MPS and TEBD MPO. They are in very good agreement with exact results at short times, and appear to be converged in bond dimensions even at longer times. However, the growth and saturation regimes are not well captured by either of the methods.

DOI
10.1103/physrevb.100.104303

Related papers

Tensor network approaches to operator spreading in ergodic quantum systems

Kévin Hémery, Frank Pollmann, David J. Luitz

We review different tensor network approaches to study the spreading of operators in generic nonintegrable quantum systems. As a common ground to all methods, we quantify this spreading by means of the Frobenius norm of the commutator of a spreading operator with a local operator, which is usually referred to as the out of time order correlation (OTOC) function. We compare two approaches based on matrix-product states in the Schrodinger picture: the time dependent block decimation (TEBD) and the time dependent variational principle (TDVP), as well as TEBD based on matrix-product operators directly in the Heisenberg picture. The results of all methods are compared to numerically exact results using Krylov space exact time evolution. We find that for the Schrodinger picture the TDVP algorithm performs better than the TEBD algorithm. Moreover the tails of the OTOC are accurately obtained both by TDVP MPS and TEBD MPO. They are in very good agreement with exact results at short times, and appear to be converged in bond dimension even at longer times. However the growth and saturation regimes are not well captured by both methods.

•Jan 17, 2019•2 citations

Time-evolution methods for matrix-product states

Sebastian Paeckel, Thomas Köhler, Andreas Swoboda +3 more

Matrix-product states have become the de facto standard for the representation of one-dimensional quantum many body states. During the last few years, numerous new methods have been introduced to evaluate the time evolution of a matrix-product state. Here, we will review and summarize the recent work on this topic as applied to finite quantum systems. We will explain and compare the different methods available to construct a time-evolved matrix-product state, namely the time-evolving block decimation, the MPO WI,II method, the global Krylov method, the local Krylov method and the one- and two-site time-dependent variational principle. We will also apply these methods to four different representative examples of current problem settings in condensed matter physics.

arXiv:1901.05824•Annals of Physics•Oct 31, 2019•763 citations

Operator front broadening in chaotic and integrable quantum chains

Javier Lopez-Piqueres, Brayden Ware, Sarang Gopalakrishnan +1 more

Operator spreading under unitary time evolution has attracted a lot of attention recently as a way to probe many-body quantum chaos. While quantities such as out-of-time-ordered correlators (OTOCs) do distinguish interacting from noninteracting systems, it has remained unclear to what extent they can truly diagnose chaotic vs integrable dynamics in many-body quantum systems. Here, we analyze operator spreading in generic one-dimensional many-body quantum systems using a combination of matrix product operator (MPO) and analytical techniques, focusing on the operator right weight. First, we show that while small bond dimension MPOs allow one to capture the exponentially decaying tail of the operator front, in agreement with earlier results, they lead to significant quantitative and qualitative errors for the actual front---defined by the maximum of the right weight. We find that while the operator front broadens diffusively in both integrable and chaotic interacting spin chains, the precise shape and scaling of the height of the front in integrable systems is anomalous for all accessible times. We interpret these results using a quasiparticle picture. This provides a sharp, although rather subtle, signature of many-body quantum chaos in the operator front.

arXiv:2103.13414•Physical review. B./Physical review. B•Sep 14, 2021•18 citations

A first introduction to Matrix Product State algorithms for the integration of Lindblad equation

Christophe Chatelain

In this introductory review, we present and compare four algorithms for the numerical integration of the Lindblad equation for one-dimensional quantum lattice systems. All four methods are based on Matrix Product State representations and can be viewed as extensions of the Time-Evolving Block Decimation (TEBD) algorithm to open quantum systems. Two approaches directly integrate the vectorized Lindblad equation, one of them explicitly enforcing the positivity of the density matrix. The other two rely on stochastic unravelings of the Lindblad equation, namely the quantum trajectory and quantum state diffusion approaches. We discuss the principles, numerical implementation, accuracy, and computational efficiency of the different methods, and benchmark them against an exactly solvable free fermion model.

cond-mat.stat-mechCross-listedarXiv:2609.18841•Sep 17, 2026•0 citations

Time-dependent variational principle in matrix-product state manifolds: Pitfalls and potential

Benedikt Kloss, Yevgeny Bar Lev, David R. Reichman

We study the applicability of the time-dependent variational principle in matrix-product state manifolds for the long time description of quantum interacting systems. By studying integrable and nonintegrable systems for which the long time dynamics are known we demonstrate that convergence of long time observables is subtle and needs to be examined carefully. Remarkably, for the disordered nonintegrable system we consider the long time dynamics are in good agreement with the rigorously obtained short time behavior and with previous obtained numerically exact results, suggesting that at least in this case, the apparent convergence of this approach is reliable. Our study indicates that, while great care must be exercised in establishing the convergence of the method, it may still be asymptotically accurate for a class of disordered nonintegrable quantum systems.

arXiv:1710.09378•Physical review. B./Physical review. B•Jan 26, 2018•87 citations

Comparison of MPS based real time evolution algorithms for Anderson Impurity Models

Daniel Bauernfeind, Markus Aichhorn, Hans Gerd Evertz

We perform a detailed comparison of two Matrix Product States (MPS) based time evolution algorithms for Anderson Impurity Models. To describe the bath, we use both the star-geometry as well as the commonly employed Wilson chain geometry. For each bath geometry, we use either the Time Dependent Variational Principle (TDVP) or the Time Evolving Block Decimation (TEBD) to perform the time evolution. To apply TEBD for the star-geometry, we use a specially adapted algorithm that can deal with the long-range coupling terms. Analyzing the major sources of errors, one expects them to be proportional to the system size for all algorithms. Surprisingly, we find errors independent of system size except for TEBD in chain geometry. Additionally, we show that the right combination of bath representation and time evolution algorithm is important. While TDVP in chain geometry is a very precise approach, TEBD in star geometry is much faster, such that for a given accuracy it is superior to TDVP in chain geometry. This makes the adapted version of TEBD in star geometry the most efficient method to solve impurity problems.

arXiv:1906.09077•Jun 21, 2019•4 citations
We review different matrix-product-state (MPS) approaches to study the spreading of operators in generic nonintegrable quantum systems. As a common ground to all methods, we quantify this spreading by means of the Frobenius norm of the commutator of a spreading operator with a local operator, which is usually referred to as the out-of-time-order correlation (OTOC) function. We compare two approaches based on matrix-product states in the Schr"odinger picture: the time-dependent block decimation (TEBD) and the time-dependent variational principle (TDVP), as well as TEBD based on matrix-product operators directly in the Heisenberg picture. The results of all methods are compared to numerically exact results using Krylov space exact time evolution. We find that for the Schr"odinger picture, the TDVP algorithm performs better than the TEBD algorithm.Moreover, the tails of the OTOC are accurately obtained both by TDVP MPS and TEBD MPO. They are in very good agreement with exact results at short times, and appear to be converged in bond dimensions even at longer times. However, the growth and saturation regimes are not well captured by either of the methods.
We review different tensor network approaches to study the spreading of operators in generic nonintegrable quantum systems. As a common ground to all methods, we quantify this spreading by means of the Frobenius norm of the commutator of a spreading operator with a local operator, which is usually referred to as the out of time order correlation (OTOC) function. We compare two approaches based on matrix-product states in the Schrodinger picture: the time dependent block decimation (TEBD) and the time dependent variational principle (TDVP), as well as TEBD based on matrix-product operators directly in the Heisenberg picture. The results of all methods are compared to numerically exact results using Krylov space exact time evolution. We find that for the Schrodinger picture the TDVP algorithm performs better than the TEBD algorithm. Moreover the tails of the OTOC are accurately obtained both by TDVP MPS and TEBD MPO. They are in very good agreement with exact results at short times, and appear to be converged in bond dimension even at longer times. However the growth and saturation regimes are not well captured by both methods.
Matrix-product states have become the de facto standard for the representation of one-dimensional quantum many body states. During the last few years, numerous new methods have been introduced to evaluate the time evolution of a matrix-product state. Here, we will review and summarize the recent work on this topic as applied to finite quantum systems. We will explain and compare the different methods available to construct a time-evolved matrix-product state, namely the time-evolving block decimation, the MPO WI,II method, the global Krylov method, the local Krylov method and the one- and two-site time-dependent variational principle. We will also apply these methods to four different representative examples of current problem settings in condensed matter physics.
Operator spreading under unitary time evolution has attracted a lot of attention recently as a way to probe many-body quantum chaos. While quantities such as out-of-time-ordered correlators (OTOCs) do distinguish interacting from noninteracting systems, it has remained unclear to what extent they can truly diagnose chaotic vs integrable dynamics in many-body quantum systems. Here, we analyze operator spreading in generic one-dimensional many-body quantum systems using a combination of matrix product operator (MPO) and analytical techniques, focusing on the operator right weight. First, we show that while small bond dimension MPOs allow one to capture the exponentially decaying tail of the operator front, in agreement with earlier results, they lead to significant quantitative and qualitative errors for the actual front---defined by the maximum of the right weight. We find that while the operator front broadens diffusively in both integrable and chaotic interacting spin chains, the precise shape and scaling of the height of the front in integrable systems is anomalous for all accessible times. We interpret these results using a quasiparticle picture. This provides a sharp, although rather subtle, signature of many-body quantum chaos in the operator front.
In this introductory review, we present and compare four algorithms for the numerical integration of the Lindblad equation for one-dimensional quantum lattice systems. All four methods are based on Matrix Product State representations and can be viewed as extensions of the Time-Evolving Block Decimation (TEBD) algorithm to open quantum systems. Two approaches directly integrate the vectorized Lindblad equation, one of them explicitly enforcing the positivity of the density matrix. The other two rely on stochastic unravelings of the Lindblad equation, namely the quantum trajectory and quantum state diffusion approaches. We discuss the principles, numerical implementation, accuracy, and computational efficiency of the different methods, and benchmark them against an exactly solvable free fermion model.
We study the applicability of the time-dependent variational principle in matrix-product state manifolds for the long time description of quantum interacting systems. By studying integrable and nonintegrable systems for which the long time dynamics are known we demonstrate that convergence of long time observables is subtle and needs to be examined carefully. Remarkably, for the disordered nonintegrable system we consider the long time dynamics are in good agreement with the rigorously obtained short time behavior and with previous obtained numerically exact results, suggesting that at least in this case, the apparent convergence of this approach is reliable. Our study indicates that, while great care must be exercised in establishing the convergence of the method, it may still be asymptotically accurate for a class of disordered nonintegrable quantum systems.
We perform a detailed comparison of two Matrix Product States (MPS) based time evolution algorithms for Anderson Impurity Models. To describe the bath, we use both the star-geometry as well as the commonly employed Wilson chain geometry. For each bath geometry, we use either the Time Dependent Variational Principle (TDVP) or the Time Evolving Block Decimation (TEBD) to perform the time evolution. To apply TEBD for the star-geometry, we use a specially adapted algorithm that can deal with the long-range coupling terms. Analyzing the major sources of errors, one expects them to be proportional to the system size for all algorithms. Surprisingly, we find errors independent of system size except for TEBD in chain geometry. Additionally, we show that the right combination of bath representation and time evolution algorithm is important. While TDVP in chain geometry is a very precise approach, TEBD in star geometry is much faster, such that for a given accuracy it is superior to TDVP in chain geometry. This makes the adapted version of TEBD in star geometry the most efficient method to solve impurity problems.