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

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

Christophe Chatelain

cond-mat.stat-mechquant-phSubmitted Sep 16, 2026• Announced Sep 17, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

Related papers

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

Dynamical simulations of classical stochastic systems using matrix product states

T. H. Johnson, Stephen R. L. Clark, Dieter Jaksch

We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.

arXiv:1006.2639•Physical Review E•Sep 16, 2010•44 citations

One-dimensional many-body entangled open quantum systems with tensor network methods

Daniel Jaschke, Simone Montangero, Lincoln D. Carr

Abstract We present a collection of methods to simulate entangled dynamics of open quantum systems governed by the Lindblad master equation with tensor network methods. Tensor network methods using matrix product states (MPSs) have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the MPS design is tailored for closed one-dimensional systems governed by the Schrödinger equation, the next step for many-body quantum dynamics is the simulation of one-dimensional open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators (MPDOs), and locally purified tensor networks. Selected examples guide possible applications of the methods and serve moreover as a benchmark between the techniques. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose–Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, pushing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that MPDOs are favorable for the exciton problem in our study. All three methods access the same library, i.e., the software package Open Source Matrix Product States , allowing us to have a meaningful comparison between the different approaches based on the selected examples. For example, tensor operations are accessed from the same subroutines and with the same optimization eliminating one possible bias in a comparison of such numerical methods.

arXiv:1804.09796•Quantum Science and Technology•Oct 10, 2018•79 citations

Simulation of Matrix Product States For Dissipation and Thermalization Dynamics of Open Quantum Systems

Souvik Agasti

Abstract We transform the system/reservoir coupling model into a one-dimensional semi-infinite discrete chain through unitary transformation to simulate the open quantum system numerically with the help of time evolving block decimation (TEBD) algorithm. We apply the method to study the dynamics of dissipative systems. We also generate the thermal state of a multimode bath using minimally entangled typical thermal state (METTS) algorithm, and investigate the impact of the thermal bath on an empty system. For both cases, we give an extensive analysis of the impact of the modeling and simulation parameters, and compare the numerics with the analytics.

Journal of Physics Communications•Dec 12, 2019•8 citations

OSMPS: Many-body entangled open quantum systems

Daniel Jaschke, Simone Montangero, Lincoln D. Carr

We present an open source code to simulate entangled dynamics of open quantum systems governed by the Lindblad master equation with tensor network methods. Tensor network methods using matrix product states have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the matrix product state design is tailored for closed one-dimensional systems governed by the Schr\odinger equation, the next step for many-body quantum dynamics is the simulation of one-dimensional open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators, and locally purified tensor networks. The package Open Source Matrix Product States combines these three major methods to simulate one-dimensional open quantum systems with the Lindblad master equation. We access the same underlying tensor network algorithms for all techniques, e.g., the tensor contractions optimized in the same way, allowing us to have a meaningful comparison between the different approaches based on selected examples. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose-Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, pushing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that matrix product density operators are favorable for the exciton problem in our study.

•Apr 25, 2018•0 citations

Superoperators vs. Trajectories for Matrix Product State Simulations of Open Quantum System: A Case Study

Lars Bonnes, Andreas Martin Läuchli

Quantum trajectories and superoperator algorithms implemented within the matrix product state (MPS) framework are powerful tools to simulate the real-time dynamics of open dissipative quantum systems. As for the unitary case, the reachable time-scales as well as system sizes are limited by the (possible) build-up of entanglement entropy. The aforementioned methods constitute complementary approaches how Lindblad master equations can be integrated relying either on a quasi-exact representation of the full density matrix or a stochastic unraveling of the density matrix in terms of pure states. In this work, we systematically benchmark both methods by studying the dynamics of a Bose-Hubbard chain in the presence of local as well as global dephasing. The build-up as well as system-size scaling of entanglement entropy strongly depends on the method and the parameter regime and we discuss the applicability of the methods for these cases as well as study the distribution of observables and time discretization errors that can become a limiting factor for global dissipation.

arXiv:1411.4831•Nov 18, 2014•14 citations
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.
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.
We adapt the time-evolving block decimation (TEBD) algorithm, originally devised to simulate the dynamics of one-dimensional quantum systems, to simulate the time evolution of nonequilibrium stochastic systems. We describe this method in detail; a system's probability distribution is represented by a matrix product state (MPS) of finite dimension and then its time evolution is efficiently simulated by repeatedly updating and approximately refactorizing this representation. We examine the use of MPS as an approximation method, looking at parallels between the interpretations of applying it to quantum state vectors and probability distributions. In the context of stochastic systems we consider two types of factorization for use in the TEBD algorithm: non-negative matrix factorization (NMF), which ensures that the approximate probability distribution is manifestly non-negative, and the singular value decomposition (SVD). Comparing these factorizations, we find the accuracy of the SVD to be substantially greater than current NMF algorithms. We then apply TEBD to simulate the totally asymmetric simple exclusion process (TASEP) for systems of up to hundreds of lattice sites in size. Using exact analytic results for the TASEP steady state, we find that TEBD reproduces this state such that the error in calculating expectation values can be made negligible even when severely compressing the description of the system by restricting the dimension of the MPS to be very small. Out of the steady state we show for specific observables that expectation values converge as the dimension of the MPS is increased to a moderate size.
Abstract We present a collection of methods to simulate entangled dynamics of open quantum systems governed by the Lindblad master equation with tensor network methods. Tensor network methods using matrix product states (MPSs) have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the MPS design is tailored for closed one-dimensional systems governed by the Schrödinger equation, the next step for many-body quantum dynamics is the simulation of one-dimensional open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators (MPDOs), and locally purified tensor networks. Selected examples guide possible applications of the methods and serve moreover as a benchmark between the techniques. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose–Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, pushing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that MPDOs are favorable for the exciton problem in our study. All three methods access the same library, i.e., the software package Open Source Matrix Product States , allowing us to have a meaningful comparison between the different approaches based on the selected examples. For example, tensor operations are accessed from the same subroutines and with the same optimization eliminating one possible bias in a comparison of such numerical methods.
Abstract We transform the system/reservoir coupling model into a one-dimensional semi-infinite discrete chain through unitary transformation to simulate the open quantum system numerically with the help of time evolving block decimation (TEBD) algorithm. We apply the method to study the dynamics of dissipative systems. We also generate the thermal state of a multimode bath using minimally entangled typical thermal state (METTS) algorithm, and investigate the impact of the thermal bath on an empty system. For both cases, we give an extensive analysis of the impact of the modeling and simulation parameters, and compare the numerics with the analytics.
We present an open source code to simulate entangled dynamics of open quantum systems governed by the Lindblad master equation with tensor network methods. Tensor network methods using matrix product states have been proven very useful to simulate many-body quantum systems and have driven many innovations in research. Since the matrix product state design is tailored for closed one-dimensional systems governed by the Schr\odinger equation, the next step for many-body quantum dynamics is the simulation of one-dimensional open quantum systems. We review the three dominant approaches to the simulation of open quantum systems via the Lindblad master equation: quantum trajectories, matrix product density operators, and locally purified tensor networks. The package Open Source Matrix Product States combines these three major methods to simulate one-dimensional open quantum systems with the Lindblad master equation. We access the same underlying tensor network algorithms for all techniques, e.g., the tensor contractions optimized in the same way, allowing us to have a meaningful comparison between the different approaches based on selected examples. These examples include the finite temperature states of the transverse quantum Ising model, the dynamics of an exciton traveling under the influence of spontaneous emission and dephasing, and a double-well potential simulated with the Bose-Hubbard model including dephasing. We analyze which approach is favorable leading to the conclusion that a complete set of all three methods is most beneficial, pushing the limits of different scenarios. The convergence studies using analytical results for macroscopic variables and exact diagonalization methods as comparison, show, for example, that matrix product density operators are favorable for the exciton problem in our study.
Quantum trajectories and superoperator algorithms implemented within the matrix product state (MPS) framework are powerful tools to simulate the real-time dynamics of open dissipative quantum systems. As for the unitary case, the reachable time-scales as well as system sizes are limited by the (possible) build-up of entanglement entropy. The aforementioned methods constitute complementary approaches how Lindblad master equations can be integrated relying either on a quasi-exact representation of the full density matrix or a stochastic unraveling of the density matrix in terms of pure states. In this work, we systematically benchmark both methods by studying the dynamics of a Bose-Hubbard chain in the presence of local as well as global dephasing. The build-up as well as system-size scaling of entanglement entropy strongly depends on the method and the parameter regime and we discuss the applicability of the methods for these cases as well as study the distribution of observables and time discretization errors that can become a limiting factor for global dissipation.