arXiv:1612.00656 · Nature Communications

A simple tensor network algorithm for two-dimensional steady states

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

Published Oct 30, 2017• 134 citations

Abstract

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.

DOI: 10.1038/s41467-017-01511-6

Related papers

A simple tensor network algorithm for two-dimensional steady states

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

Jan 1, 2017127 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, 20163 citations

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

Conor Mc Keever, Marzena Hanna Szymańska

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.

arXiv:2012.12233Physical Review XMay 14, 202129 citations

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.08258Physical Review LettersFeb 22, 201984 citations

Tensor-network-based variational Monte Carlo approach to the non-equilibrium steady state of open quantum systems

Dawid A. Hryniuk, Marzena Hanna Szymańska

We introduce a novel method of efficiently simulating the non-equilibrium steady state of large many-body open quantum systems with highly non-local interactions, based on a variational Monte Carlo optimization of a matrix product operator ansatz. Our approach outperforms and offers several advantages over comparable algorithms, such as an improved scaling of the computational cost with respect to the bond dimension for periodic systems. We showcase the versatility of our approach by studying the phase diagrams and correlation functions of the dissipative quantum Ising model with collective dephasing and long-ranged power law interactions for spin chains of up to N=100 spins.

arXiv:2405.12044QuantumSep 17, 20246 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, 20200 citations