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

TL;DR

This paper reviews and compares three tensor network methods for simulating the entangled dynamics of open quantum systems governed by the Lindblad equation, demonstrating their applications and benchmarking their performance.

Contribution

It provides a comprehensive comparison of quantum trajectories, matrix product density operators, and locally purified tensor networks for open quantum systems simulation.

Findings

Matrix product density operators are best for exciton problems.

All three methods are implemented in the Open Source Matrix Product States package.

A complete set of methods enhances simulation capabilities across scenarios.

Abstract

We present a collection of methods to simulate entangled dynamics of open quantum systems governed by the Lindblad 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 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. 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, push- ing 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. 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 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.