© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1504.06127 · Physical Review A

Matrix-product-operator approach to the nonequilibrium steady state of driven-dissipative quantum arrays

Eduardo Mascarenhas, Hugo Flayac, Vincenzo Savona

Published Aug 19, 2015• 129 citations

Abstract

We develop a numerical procedure to efficiently model the nonequilibrium steady state of one-dimensional arrays of open quantum systems based on a matrix-product operator ansatz for the density matrix. The procedure searches for the null eigenvalue of the Liouvillian superoperator by sweeping along the system while carrying out a partial diagonalization of the single-site stationary problem. It bears full analogy to the density-matrix renormalization-group approach to the ground state of isolated systems, and its numerical complexity scales as a power law with the bond dimension. The method brings considerable advantage when compared to the integration of the time-dependent problem via Trotter decomposition, as it can address arbitrarily long-ranged couplings. Additionally, it ensures numerical stability in the case of weakly dissipative systems thanks to a slow tuning of the dissipation rates along the sweeps. We have tested the method on a driven-dissipative spin chain, under various assumptions for the Hamiltonian, drive, and dissipation parameters, and compared the results to those obtained both by Trotter dynamics and Monte Carlo wave function methods. Accurate and numerically stable convergence was always achieved when applying the method to systems with a gapped Liouvillian and a nondegenerate steady state.

DOI: 10.1103/physreva.92.022116

Related papers

Variational Matrix Product Operators for the Steady State of Dissipative Quantum Systems

Jian Cui, J. I. Cirac, Mari Carmen Bañuls

We present a new variational method based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus, allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations for several dissipative spin models over a wide range of parameters illustrate the performance of the method and show that, indeed, the stationary state is often well described by a MPO of very moderate dimensions.

arXiv:1501.06786•Physical Review Letters•Jun 4, 2015•197 citations

Matrix product operator approach to nonequilibrium Floquet steady states

Zihan Cheng, Andrew C. Potter

We present a numerical method to simulate nonequilibrium Floquet steady states of one-dimensional periodically driven many-body systems coupled to a dissipative bath, based on a matrix product operator ansatz for the Floquet density matrix in frequency space. This method enables access to large systems beyond the reach of exact simulations, while retaining the periodic micromotion information. An excited-state extension of this technique allows computation of the dynamical approach to the steady state. We benchmark our method with a driven-dissipative Ising model and apply it to study the possibility of stabilizing prethermal discrete time-crystalline order by coupling to a cold bath.

arXiv:2206.07740•Physical review. B./Physical review. B•Dec 23, 2022•3 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.12044•Quantum•Sep 17, 2024•6 citations

Superoperator structures and no-go theorems for dissipative quantum phase transitions

Thomas Barthel, Yikang Zhang

In the thermodynamic limit, the steady states of open quantum many-body systems can undergo nonequilibrium phase transitions due to a competition between coherent and driven-dissipative dynamics. Here, we consider Markovian systems and elucidate structures of the Liouville superoperator that generates the time evolution. In many cases of interest, an operator-basis transformation can bring the Liouvillian into a block-triangular form, making it possible to assess its spectrum. The spectral gap sets the asymptotic decay rate. The superoperator structure can be used to bound gaps from below, showing that, in a large class of systems, dissipative phase transitions are actually impossible and that the convergence to steady states follows an exponential temporal decay. Furthermore, when the blocks on the diagonal are Hermitian, the Liouvillian spectra obey Weyl ordering relations. The results apply, for example, to Davies generators and quadratic systems and are also demonstrated for various spin models.

arXiv:2012.05505•Physical Review A•May 26, 2022•10 citations

Driven-dissipative quantum Monte Carlo method for open quantum systems

Alexandra Nagy, Vincenzo Savona

We develop a real-time full configuration-interaction quantum Monte Carlo approach to model driven-dissipative open quantum systems with Markovian system-bath coupling. The method enables stochastic sampling of the Liouville--von Neumann time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the nonequilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce the statistical error. Finally, we demonstrate the efficiency of our approach by applying it to the driven-dissipative two-dimensional $XYZ$ spin-$\frac{1}{2}$ model on a lattice.

arXiv:1802.05931•Physical Review A•May 30, 2018•42 citations

Matrix product simulations of non-equilibrium steady states of quantum spin chains

Tomaž Prosen, Marko Žnidarič

A time-dependent density matrix renormalization group method with a matrix product ansatz is employed for explicit computation of non-equilibrium steady state density operators of several integrable and non-integrable quantum spin chains, which are driven far from equilibrium by means of Markovian couplings to external baths at the two ends. It is argued that even though the time evolution cannot be simulated efficiently due to fast entanglement growth, the steady states in and out of equilibrium can be typically accurately approximated, with the result that chains of length of the order of n≈100 spins are accessible. Our results are demonstrated by performing explicit simulations of steady states and calculations of energy/spin densities/currents in several problems of heat and spin transport in quantum spin chains. A previously conjectured relation between quantum chaos and normal transport is re-confirmed with high accuracy for much larger systems.

arXiv:0811.4188•Journal of Statistical Mechanics Theory and Experiment•Feb 1, 2009•239 citations

We develop a numerical procedure to efficiently model the nonequilibrium steady state of one-dimensional arrays of open quantum systems based on a matrix-product operator ansatz for the density matrix. The procedure searches for the null eigenvalue of the Liouvillian superoperator by sweeping along the system while carrying out a partial diagonalization of the single-site stationary problem. It bears full analogy to the density-matrix renormalization-group approach to the ground state of isolated systems, and its numerical complexity scales as a power law with the bond dimension. The method brings considerable advantage when compared to the integration of the time-dependent problem via Trotter decomposition, as it can address arbitrarily long-ranged couplings. Additionally, it ensures numerical stability in the case of weakly dissipative systems thanks to a slow tuning of the dissipation rates along the sweeps. We have tested the method on a driven-dissipative spin chain, under various assumptions for the Hamiltonian, drive, and dissipation parameters, and compared the results to those obtained both by Trotter dynamics and Monte Carlo wave function methods. Accurate and numerically stable convergence was always achieved when applying the method to systems with a gapped Liouvillian and a nondegenerate steady state.

We present a new variational method based on the matrix product operator (MPO) ansatz, for finding the steady state of dissipative quantum chains governed by master equations of the Lindblad form. Instead of requiring an accurate representation of the system evolution until the stationary state is attained, the algorithm directly targets the final state, thus, allowing for a faster convergence when the steady state is a MPO with small bond dimension. Our numerical simulations for several dissipative spin models over a wide range of parameters illustrate the performance of the method and show that, indeed, the stationary state is often well described by a MPO of very moderate dimensions.

We present a numerical method to simulate nonequilibrium Floquet steady states of one-dimensional periodically driven many-body systems coupled to a dissipative bath, based on a matrix product operator ansatz for the Floquet density matrix in frequency space. This method enables access to large systems beyond the reach of exact simulations, while retaining the periodic micromotion information. An excited-state extension of this technique allows computation of the dynamical approach to the steady state. We benchmark our method with a driven-dissipative Ising model and apply it to study the possibility of stabilizing prethermal discrete time-crystalline order by coupling to a cold bath.

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.

In the thermodynamic limit, the steady states of open quantum many-body systems can undergo nonequilibrium phase transitions due to a competition between coherent and driven-dissipative dynamics. Here, we consider Markovian systems and elucidate structures of the Liouville superoperator that generates the time evolution. In many cases of interest, an operator-basis transformation can bring the Liouvillian into a block-triangular form, making it possible to assess its spectrum. The spectral gap sets the asymptotic decay rate. The superoperator structure can be used to bound gaps from below, showing that, in a large class of systems, dissipative phase transitions are actually impossible and that the convergence to steady states follows an exponential temporal decay. Furthermore, when the blocks on the diagonal are Hermitian, the Liouvillian spectra obey Weyl ordering relations. The results apply, for example, to Davies generators and quadratic systems and are also demonstrated for various spin models.

We develop a real-time full configuration-interaction quantum Monte Carlo approach to model driven-dissipative open quantum systems with Markovian system-bath coupling. The method enables stochastic sampling of the Liouville--von Neumann time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the nonequilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce the statistical error. Finally, we demonstrate the efficiency of our approach by applying it to the driven-dissipative two-dimensional

XYZ

spin-

\frac{1}{2}

model on a lattice.

A time-dependent density matrix renormalization group method with a matrix product ansatz is employed for explicit computation of non-equilibrium steady state density operators of several integrable and non-integrable quantum spin chains, which are driven far from equilibrium by means of Markovian couplings to external baths at the two ends. It is argued that even though the time evolution cannot be simulated efficiently due to fast entanglement growth, the steady states in and out of equilibrium can be typically accurately approximated, with the result that chains of length of the order of n≈100 spins are accessible. Our results are demonstrated by performing explicit simulations of steady states and calculations of energy/spin densities/currents in several problems of heat and spin transport in quantum spin chains. A previously conjectured relation between quantum chaos and normal transport is re-confirmed with high accuracy for much larger systems.