The iTEBD algorithm beyond unitary evolution

TL;DR

This paper extends the iTEBD algorithm to simulate a wider range of one-dimensional quantum and classical systems, including open dynamics, thermal states, and 2D partition functions, beyond its original unitary evolution scope.

Contribution

The authors generalize the iTEBD algorithm to handle arbitrary tensor network operators, enabling new applications in open quantum systems, thermal states, and 2D classical models.

Findings

Successfully simulates open system dynamics in 1D

Calculates thermal states and partition functions in 1D and 2D

Enhances the iTEBD algorithm's versatility for tensor network applications

Abstract

The infinite time-evolving block decimation (iTEBD) algorithm [Phys. Rev. Lett. 98, 070201 (2007)] allows to simulate unitary evolution and to compute the ground state of one-dimensional quantum lattice systems in the thermodynamic limit. Here we extend the algorithm to tackle a much broader class of problems, namely the simulation of arbitrary one-dimensional evolution operators that can be expressed as a (translationally invariant) tensor network. Relatedly, we also address the problem of finding the dominant eigenvalue and eigenvector of a one-dimensional transfer matrix that can be expressed in the same way. New applications include the simulation, in the thermodynamic limit, of open (i.e. master equation) dynamics and thermal states in 1D quantum systems, as well as calculations with partition functions in 2D classical systems, on which we elaborate. The present extension of the algorithm also plays a prominent role in the infinite projected entangled-pair states (iPEPS) approach to infinite 2D quantum lattice systems.