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:2609.16108

Large-scale quantum simulations of dissipative spin-1/2 Heisenberg chains

João C. Getelina, Andrew Cox, Muhammad Asaduzzaman, Omar Alsheikh, Ryan S. Bennink, James K. Freericks, Alexander F. Kemper

quant-phcond-mat.str-elSubmitted Sep 14, 2026• Announced Sep 16, 2026• 0 citations
arXivHTMLPDF

Abstract

A quantum many-body system coupled to an environment relaxes to a nonequilibrium steady state that can sustain order with no equilibrium counterpart. Computing such steady states is harder than closed-system dynamics as the density matrix problem squares the Hilbert-space dimension, and no free energy selects the steady state. The dissipative spin-1/2 Heisenberg chain is a benchmark example for nonequilibrium steady state physics; various methods have each calculated its phase diagram but do not agree, and a controlled determination at large system size has remained out of reach. Here we simulate the Lindblad dynamics of chains of up to 50 sites on the superconducting processor ibm_kingston -- 100 simultaneously active qubits at up to 1700 entangling-gate depths -- realizing the dissipation via Stinespring dilation. The system's dissipative evolution is a self-correcting mechanism that effectively erases errors, so hardware noise enters only as a weak competing dissipator. We measure static structure factors and resolve ferromagnetic, antiferromagnetic, spin-density-wave, and paramagnetic steady states, mapping the phase diagram with 117 quantum hardware data points across the $J_x$--$J_y$ plane. We uncover a rich non-equilibrium phase diagram of ordered phases with only remnants of the mean-field order, and where sharp transitions give way to the crossovers expected in one dimension. We also find the existence of an incipient (Trotter-induced) spin density wave phase, highlighting the potential of controlled Trotterization as a tool to engineer various magnetic phases in dissipative spin systems. Our quantum simulations largely settle the lingering uncertainty regarding the correct phase diagram of this benchmark system. Moreover, they show that quantum computers are now a feasible tool for addressing scientific questions involving dissipative quantum systems.

Related papers

Enhancing quantum utility: Simulating large-scale quantum spin chains on superconducting quantum computers

Talal Ahmed Chowdhury, Kwangmin Yu, Mahmud Ashraf Shamim +2 more

We present the quantum simulation of the frustrated quantum spin- 1 2 antiferromagnetic Heisenberg spin chain with competing nearest-neighbor ( J 1 ) and next-nearest-neighbor ( J 2 ) exchange interactions in the real superconducting quantum computer with qubits ranging up to 100. In particular, we implement the Hamiltonian with the next-nearest neighbor exchange interaction in conjunction with the nearest-neighbor interaction on IBM's superconducting quantum computer and carry out the time evolution of the spin chain by employing the first-order Trotterization. Furthermore, our implementation of the second-order Trotterization for the isotropic Heisenberg spin chain, involving only nearest-neighbor exchange interaction, enables precise measurement of the expectation values of staggered magnetization observable across a range of up to 100 qubits. Notably, in both cases, our approach results in a constant circuit depth in each Trotter step, independent of the number of qubits. Our demonstration of the accurate measurement of expectation values for the large-scale quantum system using superconducting quantum computers designates the quantum utility of these devices for investigating various properties of many-body quantum systems. This will be a stepping stone to achieving the quantum advantage over classical ones in simulating quantum systems before the fault tolerance quantum era. Published by the American Physical Society 2024

arXiv:2312.12427•Physical Review Research•Jul 25, 2024•29 citations

Quantum Monte Carlo approach to the non-equilibrium steady state of open quantum systems

Alexandra Nagy

Many-body open quantum systems are exposed to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. As the exact treatment of such models is generally unfeasible, it is favourable to formulate an approximate description by means of the dynamics of the reduced density operator of the system. When studying models with a weak coupling to a memoryless environment, the existence of such solution is granted and the reduced dynamics is governed by the Lindblad quantum master equation. In recent years, open quantum systems have evolved into a major field of studies. Focus of these studies are the characterization of emergent phenomena and dissipative phase transitions, as well as the ongoing debate about whether quantum computing schemes are still hard to simulate classically --- and thus achieve quantum supremacy --- when in presence of some degree of noise-induced decoherence. However, the challenge posed by this task lies in the complexity of the density matrix that increases exponentially with the system size, and the quest for efficient numerical methods is a research field that is still in its infancy. Here, we first develop a real-time full configuration interaction quantum Monte Carlo technique that stems from a class of methods generally known as projector Monte Carlo. The approach enables the stochastic sampling of the Lindblad time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the non-equilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce statistical error. We demonstrate the efficiency of the approach by applying it to the dissipative two-dimensional XYZ spin-1/2 model on a lattice. As the importance sampling of projector approaches is often combined with variational results, we then introduce a novel method that is based on the variational Monte Carlo methods and on a neural network representation of the density matrix. Thanks to the stochastic reconfiguration scheme, the application of the variational principle is translated into the actual integration of the Lindblad quantum master equation. We test the effectiveness of the method by modeling the steady state of the dissipative two-dimensional XYZ spin-1/2 model through a dissipative phase transition, and also the real-time dynamics of the dissipative Ising model. In addition, we discuss the application of the developed methods and the open questions of the field.

•Jan 1, 2020•0 citations

Dissipation Mechanisms and Dissipative Phase Transitions of two coupled Fully Connected Quantum Ising models

Bidyut Dey, Andrea Kazemi Nava, Domenico Giuliano

We study dissipative phase transitions in a system of two coupled fully-connected quantum Ising models interacting with an environment. The dynamics is governed by a Lindblad master equation combining coherent unitary evolution and incoherent dissipative processes, where the unitary part is described within a self-consistent mean-field framework effectively acting on the local Hilbert space of two coupled spins at each site. We analyze two fundamentally different classes of dissipators. In the first case, the jump operators are defined in the instantaneous eigenbasis of the mean-field Hamiltonian and satisfy a detailed-balance condition. In this setting, the relaxation dynamics depends strongly on the quench protocol: a parametric quench of the Hamiltonian leads to conventional relaxation, whereas a temperature quench gives rise to a dynamical phase transition characterized by nonanalytic behavior in time. Yet, in both cases, the system relaxes toward a steady state determined solely by the post-quench parameters and the bath temperature, which closely resembles a thermal Gibbs state of the mean-field Hamiltonian. As a result, the dissipative phase transition occurs at a critical point consistent with the corresponding equilibrium transition. In contrast, when the dissipators are realized via local spin raising and lowering operators, the steady state is genuinely nonequilibrium, leading to a significantly richer phase diagram. In particular, for sufficiently strong system-bath coupling, we observe a reentrant phase featuring a symmetry-broken region bounded by two continuous dissipative phase transitions. Our results evidence how the structure of dissipative processes controls the emergence of equilibrium-like versus genuinely nonequilibrium critical behavior in open quantum systems.

arXiv:2604.27072•Apr 29, 2026•0 citations

Simulating quantum transport via collisional models on a digital quantum computer

Rebecca Erbanni, Xiansong Xu, Tommaso F. Demarie +1 more

Digital quantum computers have the potential to study the dynamics of many-body quantum systems. Nonequilibrium open quantum systems are, however, less straightforward to be implemented. Here we explore the feasibility of studying steady-state transport in strongly interacting many-body quantum systems on a digital quantum computer. To do so, we consider a collisional model representation of the nonequilibrium open dynamics for a boundary-driven $XXZ$ spin chain. More specifically, we investigate how the depth of the quantum circuit is affected by how close we want the steady state to be to the one expected from the underlying master equation. We study the simulation of a boundary-driven spin chain in regimes of weak and strong interactions, which would lead in large systems to diffusive and ballistic dynamics, considering also possible errors in the implementation of the protocol. Last, we analyze the effectiveness of digital simulation via the collisional model of current rectification when the $XXZ$ spin chains are subject to nonuniform magnetic fields and show that, although the circuit depths required to reach steady states are still prohibitive for today's hardware, few collisions are enough to suggest a strong rectifying power.

arXiv:2307.13576•Physical Review A•Sep 28, 2023•9 citations

Probing non-equilibrium dissipative phase transitions with trapped-ion quantum simulators

Casey Haack, Naushad Ahmad Kamar, Daniel A. Paz +2 more

Open quantum many-body systems with controllable dissipation can exhibit novel features in their dynamics and steady states. A paradigmatic example is the dissipative transverse field Ising model. It has been shown recently that the steady state of this model with all-to-all interactions is genuinely non-equilibrium near criticality, exhibiting a modified time-reversal symmetry and violating the fluctuation-dissipation theorem. Experimental study of such non-equilibrium steady-state phase transitions is however lacking. Here we propose realistic experimental setups and measurement schemes for current trapped-ion quantum simulators to demonstrate this phase transition, where controllable dissipation is engineered via a continuous weak optical pumping laser. With extensive numerical calculations, we show that strong signatures of this dissipative phase transition and its non-equilibrium properties can be observed with a small system size across a wide range of system parameters. In addition, we show that the same signatures can also be seen if the dissipation is instead achieved via Floquet dynamics with periodic and probabilistic resetting of the spins. Dissipation engineered in this way may allow the simulation of more general types of driven-dissipative systems or facilitate the dissipative preparation of useful many-body entangled states.

arXiv:2311.06199•Nov 10, 2023•0 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
A quantum many-body system coupled to an environment relaxes to a nonequilibrium steady state that can sustain order with no equilibrium counterpart. Computing such steady states is harder than closed-system dynamics as the density matrix problem squares the Hilbert-space dimension, and no free energy selects the steady state. The dissipative spin-1/2 Heisenberg chain is a benchmark example for nonequilibrium steady state physics; various methods have each calculated its phase diagram but do not agree, and a controlled determination at large system size has remained out of reach. Here we simulate the Lindblad dynamics of chains of up to 50 sites on the superconducting processor ibm_kingston -- 100 simultaneously active qubits at up to 1700 entangling-gate depths -- realizing the dissipation via Stinespring dilation. The system's dissipative evolution is a self-correcting mechanism that effectively erases errors, so hardware noise enters only as a weak competing dissipator. We measure static structure factors and resolve ferromagnetic, antiferromagnetic, spin-density-wave, and paramagnetic steady states, mapping the phase diagram with 117 quantum hardware data points across the JxJ_xJx​--JyJ_yJy​ plane. We uncover a rich non-equilibrium phase diagram of ordered phases with only remnants of the mean-field order, and where sharp transitions give way to the crossovers expected in one dimension. We also find the existence of an incipient (Trotter-induced) spin density wave phase, highlighting the potential of controlled Trotterization as a tool to engineer various magnetic phases in dissipative spin systems. Our quantum simulations largely settle the lingering uncertainty regarding the correct phase diagram of this benchmark system. Moreover, they show that quantum computers are now a feasible tool for addressing scientific questions involving dissipative quantum systems.
We present the quantum simulation of the frustrated quantum spin- 1 2 antiferromagnetic Heisenberg spin chain with competing nearest-neighbor ( J 1 ) and next-nearest-neighbor ( J 2 ) exchange interactions in the real superconducting quantum computer with qubits ranging up to 100. In particular, we implement the Hamiltonian with the next-nearest neighbor exchange interaction in conjunction with the nearest-neighbor interaction on IBM's superconducting quantum computer and carry out the time evolution of the spin chain by employing the first-order Trotterization. Furthermore, our implementation of the second-order Trotterization for the isotropic Heisenberg spin chain, involving only nearest-neighbor exchange interaction, enables precise measurement of the expectation values of staggered magnetization observable across a range of up to 100 qubits. Notably, in both cases, our approach results in a constant circuit depth in each Trotter step, independent of the number of qubits. Our demonstration of the accurate measurement of expectation values for the large-scale quantum system using superconducting quantum computers designates the quantum utility of these devices for investigating various properties of many-body quantum systems. This will be a stepping stone to achieving the quantum advantage over classical ones in simulating quantum systems before the fault tolerance quantum era. Published by the American Physical Society 2024
Many-body open quantum systems are exposed to an essentially uncontrollable environment that acts as a source of decoherence and dissipation. As the exact treatment of such models is generally unfeasible, it is favourable to formulate an approximate description by means of the dynamics of the reduced density operator of the system. When studying models with a weak coupling to a memoryless environment, the existence of such solution is granted and the reduced dynamics is governed by the Lindblad quantum master equation. In recent years, open quantum systems have evolved into a major field of studies. Focus of these studies are the characterization of emergent phenomena and dissipative phase transitions, as well as the ongoing debate about whether quantum computing schemes are still hard to simulate classically --- and thus achieve quantum supremacy --- when in presence of some degree of noise-induced decoherence. However, the challenge posed by this task lies in the complexity of the density matrix that increases exponentially with the system size, and the quest for efficient numerical methods is a research field that is still in its infancy. Here, we first develop a real-time full configuration interaction quantum Monte Carlo technique that stems from a class of methods generally known as projector Monte Carlo. The approach enables the stochastic sampling of the Lindblad time evolution of the density matrix thanks to a massively parallel algorithm, thus providing estimates of observables on the non-equilibrium steady state. We present the underlying theory and introduce an initiator technique and importance sampling to reduce statistical error. We demonstrate the efficiency of the approach by applying it to the dissipative two-dimensional XYZ spin-1/2 model on a lattice. As the importance sampling of projector approaches is often combined with variational results, we then introduce a novel method that is based on the variational Monte Carlo methods and on a neural network representation of the density matrix. Thanks to the stochastic reconfiguration scheme, the application of the variational principle is translated into the actual integration of the Lindblad quantum master equation. We test the effectiveness of the method by modeling the steady state of the dissipative two-dimensional XYZ spin-1/2 model through a dissipative phase transition, and also the real-time dynamics of the dissipative Ising model. In addition, we discuss the application of the developed methods and the open questions of the field.
We study dissipative phase transitions in a system of two coupled fully-connected quantum Ising models interacting with an environment. The dynamics is governed by a Lindblad master equation combining coherent unitary evolution and incoherent dissipative processes, where the unitary part is described within a self-consistent mean-field framework effectively acting on the local Hilbert space of two coupled spins at each site. We analyze two fundamentally different classes of dissipators. In the first case, the jump operators are defined in the instantaneous eigenbasis of the mean-field Hamiltonian and satisfy a detailed-balance condition. In this setting, the relaxation dynamics depends strongly on the quench protocol: a parametric quench of the Hamiltonian leads to conventional relaxation, whereas a temperature quench gives rise to a dynamical phase transition characterized by nonanalytic behavior in time. Yet, in both cases, the system relaxes toward a steady state determined solely by the post-quench parameters and the bath temperature, which closely resembles a thermal Gibbs state of the mean-field Hamiltonian. As a result, the dissipative phase transition occurs at a critical point consistent with the corresponding equilibrium transition. In contrast, when the dissipators are realized via local spin raising and lowering operators, the steady state is genuinely nonequilibrium, leading to a significantly richer phase diagram. In particular, for sufficiently strong system-bath coupling, we observe a reentrant phase featuring a symmetry-broken region bounded by two continuous dissipative phase transitions. Our results evidence how the structure of dissipative processes controls the emergence of equilibrium-like versus genuinely nonequilibrium critical behavior in open quantum systems.
Digital quantum computers have the potential to study the dynamics of many-body quantum systems. Nonequilibrium open quantum systems are, however, less straightforward to be implemented. Here we explore the feasibility of studying steady-state transport in strongly interacting many-body quantum systems on a digital quantum computer. To do so, we consider a collisional model representation of the nonequilibrium open dynamics for a boundary-driven XXZXXZXXZ spin chain. More specifically, we investigate how the depth of the quantum circuit is affected by how close we want the steady state to be to the one expected from the underlying master equation. We study the simulation of a boundary-driven spin chain in regimes of weak and strong interactions, which would lead in large systems to diffusive and ballistic dynamics, considering also possible errors in the implementation of the protocol. Last, we analyze the effectiveness of digital simulation via the collisional model of current rectification when the XXZXXZXXZ spin chains are subject to nonuniform magnetic fields and show that, although the circuit depths required to reach steady states are still prohibitive for today's hardware, few collisions are enough to suggest a strong rectifying power.
Open quantum many-body systems with controllable dissipation can exhibit novel features in their dynamics and steady states. A paradigmatic example is the dissipative transverse field Ising model. It has been shown recently that the steady state of this model with all-to-all interactions is genuinely non-equilibrium near criticality, exhibiting a modified time-reversal symmetry and violating the fluctuation-dissipation theorem. Experimental study of such non-equilibrium steady-state phase transitions is however lacking. Here we propose realistic experimental setups and measurement schemes for current trapped-ion quantum simulators to demonstrate this phase transition, where controllable dissipation is engineered via a continuous weak optical pumping laser. With extensive numerical calculations, we show that strong signatures of this dissipative phase transition and its non-equilibrium properties can be observed with a small system size across a wide range of system parameters. In addition, we show that the same signatures can also be seen if the dissipation is instead achieved via Floquet dynamics with periodic and probabilistic resetting of the spins. Dissipation engineered in this way may allow the simulation of more general types of driven-dissipative systems or facilitate the dissipative preparation of useful many-body entangled states.
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.