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:1608.06028 · Physical Review Letters

Steady States of Infinite-Size Dissipative Quantum Chains via Imaginary Time Evolution

Adil A. Gangat, Te I, Ying-Jer Kao

Published Jul 5, 2017• 44 citations
arXivHTMLPDF

Abstract

Directly in the thermodynamic limit, we show how to combine local imaginary and real-time evolution of tensor networks to efficiently and accurately find the nonequilibrium steady states (NESSs) of one-dimensional dissipative quantum lattices governed by a local Lindblad master equation. The imaginary time evolution first bypasses any highly correlated portions of the real-time evolution trajectory by directly converging to the weakly correlated subspace of the NESS, after which, real-time evolution completes the convergence to the NESS with high accuracy. We demonstrate the power of the method with the dissipative transverse field quantum Ising chain. We show that a crossover of an order parameter shown to be smooth in previous finite-size studies remains smooth in the thermodynamic limit.

DOI
10.1103/physrevlett.119.010501

Related papers

Efficient Simulation of Nontrivial Dissipative Spin Chains via Stochastic Unraveling

Andrew Pocklington, Aashish A. Clerk

We present a new technique for efficiently simulating (in polynomial time) a class of one-dimensional (1D) dissipative spin chains that, when mapped to fermions, have quadratic Hamiltonians, with the only nonlinearity coming from Jordan-Wigner strings appearing in the jump operators, despite the fact that these models cannot be mapped to quadratic fermionic master equations. We show that many such Lindblad master equations admit an exact stochastic unraveling, with individual trajectories evolving as Gaussian fermionic states, even though the full master equation describes a system inequivalent to free fermions. This allows one to calculate arbitrary observables efficiently without sign problems, and with bounded sampling complexity. We utilize this new technique to study three paradigmatic dissipative effects: the melting of antiferromagnetic order in the presence of local loss, many-body subradiant phenomenon in systems with correlated loss, and nonequilibrium steady states of a 1D dissipative transverse-field Ising model. Beyond simply providing a powerful numerical technique, our method can also be used to gain both qualitative and quantitative insights into the role of interactions in these models.

arXiv:2503.23469•PRX Quantum•Aug 6, 2025•1 citation

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.12233•Physical Review X•May 14, 2021•29 citations

Dissipative dynamical phase transition as a complex Ising model

Stephen W. Yan, Diego Barberena, Matthew P. A. Fisher +1 more

We investigate a quantum dynamical phase transition induced by the competition between local unitary evolution and dissipation in a qubit chain with a strong, on-site $\mathbb{Z}_2$ symmetry. While the steady-state of this evolution is always maximally-mixed, we show that the dynamical behavior of certain non-local observables on the approach to this steady-state is dictated by a quantum Ising model with a complex transverse-field (cTFIM). A crucial feature of our mapping to the cTFIM is an extensive number of local, weak, $\mathbb{Z}_2$ symmetries possessed by the dynamics, allowing us to investigate these observables analytically and uncover a dynamical phase transition as the relative rate of unitary evolution and dissipation is tuned. We show that the weak-dissipation limit corresponds to a cTFIM with a large magnitude of the imaginary transverse-field, for which the many-body "ground-state" (with smallest real eigenvalue) is gapless, exhibiting quasi-long-range correlations of the local magnetization with a continuously-varying exponent. Correspondingly, the dynamics of the non-local observables show oscillatory behavior with an amplitude decaying exponentially in time. The strong-dissipation limit corresponds to a gapped ferromagnetic phase of the cTFIM, and non-local observables show exponential decay on the approach to equilibrium. This transition in (1+1)-dimensions has a peculiar, "two-sided" nature appearing as either first- or second-order depending on the phase from which the transition is approached, an analytic result which is corroborated by numerical studies. In higher dimensions, we present a field-theoretic understanding of the first-order nature of this transition, when approaching from the ferromagnetic phase of the cTFIM, though the nature of the phase with large imaginary transverse-field remains to be understood.

arXiv:2412.09591•Physical review. B./Physical review. B•Sep 28, 2026•0 citations

Quantum Imaginary Time Evolution on an Infinite 1D Chain

Hao-Ti Hung, Tung Tsao, Ying-Jer Kao

We introduce a quantum-circuit algorithm for performing imaginary-time evolution on infinite one-dimensional lattice systems. The method uses a parameterized quantum circuit to represent a uniform matrix product state ansatz. We derive the ITE algorithm using the time-dependent variational principle and employ the quantum Lanczos algorithm to improve the ground-state energy estimate. As a benchmark, we simulate the transverse-field Ising model using both classical simulators and IBM Quantum devices. Our analysis includes a statistical study of the distributions of the cost function and energy density obtained from quantum measurements, illustrating the effects of finite-sampling noise on convergence.

quant-pharXiv:2608.30363•Sep 1, 2026•0 citations

Effective quantum Zeno dynamics in dissipative quantum systems

Vladislav Popkov, Simon Essink, Carlo Presilla A +1 more

We investigate the time evolution of an open quantum system described by a Lindblad master equation with dissipation acting only on a part of the degrees of freedom ${\mathcal{H}}_{0}$ of the system, and targeting a unique dark state in ${\mathcal{H}}_{0}$. We show that, in the Zeno limit of large dissipation, the density matrix of the system traced over the dissipative subspace ${\mathcal{H}}_{0}$, evolves according to another Lindblad dynamics, with renormalized effective Hamiltonian and weak effective dissipation. This behavior is explicitly checked in the case of Heisenberg spin chains with one or both boundary spins strongly coupled to a magnetic reservoir. Moreover, the populations of the eigenstates of the renormalized effective Hamiltonian evolve in time according to a classical Markov dynamics. As a direct application of this result, we propose a computationally efficient exact method to evaluate the nonequilibrium steady state of a general system in the limit of strong dissipation.

arXiv:1806.10422•Physical Review A•Nov 7, 2018•39 citations

Multistability of Driven-Dissipative Quantum Spins

Haggai Landa, Marco Schirò, Grégoire Misguich

We study the dynamics of lattice models of quantum spins one-half, driven by a coherent drive and subject to dissipation. Generically the mean-field limit of these models manifests multistable parameter regions of coexisting steady states with different magnetizations. We introduce an efficient scheme accounting for the corrections to mean field by correlations at leading order, and benchmark this scheme using high-precision numerics based on matrix-product operators in one- and two-dimensional lattices. Correlations are shown to wash the mean-field bistability in dimension one, leading to a unique steady state. In dimension two and higher, we find that multistability is again possible, provided the thermodynamic limit of an infinitely large lattice is taken first with respect to the longtime limit. Variation of the system parameters results in jumps between the different steady states, each showing a critical slowing down in the convergence of perturbations towards the steady state. Experiments with trapped ions can realize the model and possibly answer open questions in the nonequilibrium many-body dynamics of these quantum systems, beyond the system sizes accessible to present numerics.

arXiv:1905.10349•Physical Review Letters•Jan 29, 2020•97 citations
Directly in the thermodynamic limit, we show how to combine local imaginary and real-time evolution of tensor networks to efficiently and accurately find the nonequilibrium steady states (NESSs) of one-dimensional dissipative quantum lattices governed by a local Lindblad master equation. The imaginary time evolution first bypasses any highly correlated portions of the real-time evolution trajectory by directly converging to the weakly correlated subspace of the NESS, after which, real-time evolution completes the convergence to the NESS with high accuracy. We demonstrate the power of the method with the dissipative transverse field quantum Ising chain. We show that a crossover of an order parameter shown to be smooth in previous finite-size studies remains smooth in the thermodynamic limit.
We present a new technique for efficiently simulating (in polynomial time) a class of one-dimensional (1D) dissipative spin chains that, when mapped to fermions, have quadratic Hamiltonians, with the only nonlinearity coming from Jordan-Wigner strings appearing in the jump operators, despite the fact that these models cannot be mapped to quadratic fermionic master equations. We show that many such Lindblad master equations admit an exact stochastic unraveling, with individual trajectories evolving as Gaussian fermionic states, even though the full master equation describes a system inequivalent to free fermions. This allows one to calculate arbitrary observables efficiently without sign problems, and with bounded sampling complexity. We utilize this new technique to study three paradigmatic dissipative effects: the melting of antiferromagnetic order in the presence of local loss, many-body subradiant phenomenon in systems with correlated loss, and nonequilibrium steady states of a 1D dissipative transverse-field Ising model. Beyond simply providing a powerful numerical technique, our method can also be used to gain both qualitative and quantitative insights into the role of interactions in these models.
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.
We investigate a quantum dynamical phase transition induced by the competition between local unitary evolution and dissipation in a qubit chain with a strong, on-site Z2\mathbb{Z}_2Z2​ symmetry. While the steady-state of this evolution is always maximally-mixed, we show that the dynamical behavior of certain non-local observables on the approach to this steady-state is dictated by a quantum Ising model with a complex transverse-field (cTFIM). A crucial feature of our mapping to the cTFIM is an extensive number of local, weak, Z2\mathbb{Z}_2Z2​ symmetries possessed by the dynamics, allowing us to investigate these observables analytically and uncover a dynamical phase transition as the relative rate of unitary evolution and dissipation is tuned. We show that the weak-dissipation limit corresponds to a cTFIM with a large magnitude of the imaginary transverse-field, for which the many-body "ground-state" (with smallest real eigenvalue) is gapless, exhibiting quasi-long-range correlations of the local magnetization with a continuously-varying exponent. Correspondingly, the dynamics of the non-local observables show oscillatory behavior with an amplitude decaying exponentially in time. The strong-dissipation limit corresponds to a gapped ferromagnetic phase of the cTFIM, and non-local observables show exponential decay on the approach to equilibrium. This transition in (1+1)-dimensions has a peculiar, "two-sided" nature appearing as either first- or second-order depending on the phase from which the transition is approached, an analytic result which is corroborated by numerical studies. In higher dimensions, we present a field-theoretic understanding of the first-order nature of this transition, when approaching from the ferromagnetic phase of the cTFIM, though the nature of the phase with large imaginary transverse-field remains to be understood.
We introduce a quantum-circuit algorithm for performing imaginary-time evolution on infinite one-dimensional lattice systems. The method uses a parameterized quantum circuit to represent a uniform matrix product state ansatz. We derive the ITE algorithm using the time-dependent variational principle and employ the quantum Lanczos algorithm to improve the ground-state energy estimate. As a benchmark, we simulate the transverse-field Ising model using both classical simulators and IBM Quantum devices. Our analysis includes a statistical study of the distributions of the cost function and energy density obtained from quantum measurements, illustrating the effects of finite-sampling noise on convergence.
We investigate the time evolution of an open quantum system described by a Lindblad master equation with dissipation acting only on a part of the degrees of freedom H0{\mathcal{H}}_{0}H0​ of the system, and targeting a unique dark state in H0{\mathcal{H}}_{0}H0​. We show that, in the Zeno limit of large dissipation, the density matrix of the system traced over the dissipative subspace H0{\mathcal{H}}_{0}H0​, evolves according to another Lindblad dynamics, with renormalized effective Hamiltonian and weak effective dissipation. This behavior is explicitly checked in the case of Heisenberg spin chains with one or both boundary spins strongly coupled to a magnetic reservoir. Moreover, the populations of the eigenstates of the renormalized effective Hamiltonian evolve in time according to a classical Markov dynamics. As a direct application of this result, we propose a computationally efficient exact method to evaluate the nonequilibrium steady state of a general system in the limit of strong dissipation.
We study the dynamics of lattice models of quantum spins one-half, driven by a coherent drive and subject to dissipation. Generically the mean-field limit of these models manifests multistable parameter regions of coexisting steady states with different magnetizations. We introduce an efficient scheme accounting for the corrections to mean field by correlations at leading order, and benchmark this scheme using high-precision numerics based on matrix-product operators in one- and two-dimensional lattices. Correlations are shown to wash the mean-field bistability in dimension one, leading to a unique steady state. In dimension two and higher, we find that multistability is again possible, provided the thermodynamic limit of an infinitely large lattice is taken first with respect to the longtime limit. Variation of the system parameters results in jumps between the different steady states, each showing a critical slowing down in the convergence of perturbations towards the steady state. Experiments with trapped ions can realize the model and possibly answer open questions in the nonequilibrium many-body dynamics of these quantum systems, beyond the system sizes accessible to present numerics.