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:2012.12233 · Physical Review X

Stable iPEPO Tensor-Network Algorithm for Dynamics of Two-Dimensional Open Quantum Lattice Models

Conor Mc Keever, Marzena Hanna Szymańska

Published May 14, 2021• 29 citations
arXivHTMLPDF

Abstract

Being able to accurately describe the dynamics steady states of driven and/or dissipative but quantum correlated lattice models is of fundamental importance in many areas of science: from quantum information to biology.An efficient numerical simulation of large open systems in two spatial dimensions is a challenge.In this work, we develop a tensor network method, based on an infinite projected entangled pair operator ansatz, applicable directly in the thermodynamic limit.We incorporate techniques of finding optimal truncations of enlarged network bonds by optimizing an objective function appropriate for open systems.Comparisons with numerically exact calculations, both for the dynamics and the steady state, demonstrate the power of the method.In particular, we consider dissipative transverse quantum Ising, driven-dissipative hard-core boson, and dissipative anisotropic XY models in non-mean-field limits, proving able to capture substantial entanglement in the presence of dissipation.Our method enables us to study regimes that are accessible to current experiments but lie well beyond the applicability of existing techniques.

DOI
10.1103/physrevx.11.021035

Related papers

Tensor Network Annealing Algorithm for Two-Dimensional Thermal States

Augustine Kshetrimayum, Matteo Rizzi, Jens Eisert +1 more

Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.

arXiv:1809.08258•Physical Review Letters•Feb 22, 2019•84 citations

Efficient Time Evolution of 2D Open-Quantum Lattice Models with Long-Range Interactions using Tensor Networks

Jack Dunham, Marzena H. Szymańska

Simulating many-body open quantum systems is an extremely challenging problem, with methods often restricted to either models with nearest-neighbor interactions or semi-classical approximations. In particular, modeling two-dimensional systems with realistic long-range interactions, in addition to dissipation, is of vital importance to the development of modern quantum computing and simulation platforms. In this paper, we present a construction of the time-evolution operator, as a projected entangled pair operator (denoted tePEPO), that can be used to evolve a tensor network ansatz through time. Interactions beyond nearest-neighbor, including interactions between sites not collinear in the lattice, can be represented efficiently as a tePEPO. Furthermore, we obtain approximations to realistic radial long-range interactions decaying with a power-law, that give accurate results with small tePEPO bond dimension. Finally, we consider a physical example of a Rydberg atom Hamiltonian with long-range dipolar interactions, and show evidence of a dipole-dipole blockading effect in presence of dissipation. This work demonstrates the applicability of tensor networks to two-dimensional systems widely studied in experiments, but previously inaccessible to non-semi-classical methods.

quant-pharXiv:2512.01781•Dec 2, 2025•0 citations

Tensor network simulations of open quantum systems

Dainius Kilda

The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.

•Mar 3, 2020•0 citations

A simple tensor network algorithm for two-dimensional steady states

Augustine Kshetrimayum, Hendrik Weimer, Román Orús

Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.

arXiv:1612.00656•Nature Communications•Oct 30, 2017•134 citations

A simple tensor network algorithm for two-dimensional steady states

Augustine Kshetrimayum, Hendrik Weimer, Román Orús

Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.

•Jan 1, 2017•127 citations

A simple tensor network algorithm for 2d steady states

Augustine Kshetrimayum, Hendrik Weimer, Román Orús

Understanding dissipation in 2D quantum many-body systems is a remarkably difficult open challenge. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady-states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin-1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.

•Dec 2, 2016•3 citations
Being able to accurately describe the dynamics steady states of driven and/or dissipative but quantum correlated lattice models is of fundamental importance in many areas of science: from quantum information to biology.An efficient numerical simulation of large open systems in two spatial dimensions is a challenge.In this work, we develop a tensor network method, based on an infinite projected entangled pair operator ansatz, applicable directly in the thermodynamic limit.We incorporate techniques of finding optimal truncations of enlarged network bonds by optimizing an objective function appropriate for open systems.Comparisons with numerically exact calculations, both for the dynamics and the steady state, demonstrate the power of the method.In particular, we consider dissipative transverse quantum Ising, driven-dissipative hard-core boson, and dissipative anisotropic XY models in non-mean-field limits, proving able to capture substantial entanglement in the presence of dissipation.Our method enables us to study regimes that are accessible to current experiments but lie well beyond the applicability of existing techniques.
Tensor network methods have become a powerful class of tools to capture strongly correlated matter, but methods to capture the experimentally ubiquitous family of models at finite temperature beyond one spatial dimension are largely lacking. We introduce a tensor network algorithm able to simulate thermal states of two-dimensional quantum lattice systems in the thermodynamic limit. The method develops instances of projected entangled pair states and projected entangled pair operators for this purpose. It is the key feature of this algorithm to resemble the cooling down of the system from an infinite temperature state until it reaches the desired finite-temperature regime. As a benchmark, we study the finite-temperature phase transition of the Ising model on an infinite square lattice, for which we obtain remarkable agreement with the exact solution. We then turn to study the finite-temperature Bose-Hubbard model in the limits of two (hard-core) and three bosonic modes per site. Our technique can be used to support the experimental study of actual effectively two-dimensional materials in the laboratory, as well as to benchmark optical lattice quantum simulators with ultracold atoms.
Simulating many-body open quantum systems is an extremely challenging problem, with methods often restricted to either models with nearest-neighbor interactions or semi-classical approximations. In particular, modeling two-dimensional systems with realistic long-range interactions, in addition to dissipation, is of vital importance to the development of modern quantum computing and simulation platforms. In this paper, we present a construction of the time-evolution operator, as a projected entangled pair operator (denoted tePEPO), that can be used to evolve a tensor network ansatz through time. Interactions beyond nearest-neighbor, including interactions between sites not collinear in the lattice, can be represented efficiently as a tePEPO. Furthermore, we obtain approximations to realistic radial long-range interactions decaying with a power-law, that give accurate results with small tePEPO bond dimension. Finally, we consider a physical example of a Rydberg atom Hamiltonian with long-range dipolar interactions, and show evidence of a dipole-dipole blockading effect in presence of dissipation. This work demonstrates the applicability of tensor networks to two-dimensional systems widely studied in experiments, but previously inaccessible to non-semi-classical methods.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
Understanding dissipation in 2D quantum many-body systems is an open challenge which has proven remarkably difficult. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin 1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.
Understanding dissipation in 2D quantum many-body systems is a remarkably difficult open challenge. Here we show how numerical simulations for this problem are possible by means of a tensor network algorithm that approximates steady-states of 2D quantum lattice dissipative systems in the thermodynamic limit. Our method is based on the intuition that strong dissipation kills quantum entanglement before it gets too large to handle. We test its validity by simulating a dissipative quantum Ising model, relevant for dissipative systems of interacting Rydberg atoms, and benchmark our simulations with a variational algorithm based on product and correlated states. Our results support the existence of a first order transition in this model, with no bistable region. We also simulate a dissipative spin-1/2 XYZ model, showing that there is no re-entrance of the ferromagnetic phase. Our method enables the computation of steady states in 2D quantum lattice systems.