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

Multistability of Driven-Dissipative Quantum Spins

Haggai Landa, Marco Schirò, Grégoire Misguich

Published Jan 29, 2020• 97 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/physrevlett.124.043601

Related papers

Correlation-induced steady states and limit cycles in driven dissipative quantum systems

Haggai Landa, Marco Schirò, Grégoire Misguich

We study a driven-dissipative model of spins one-half (qubits) on a lattice with nearest-neighbor interactions. Focusing on the role of spatially extended spin-spin correlations in determining the phases of the system, we characterize the spatial structure of the correlations in the steady state, as well as their temporal dynamics. In dimension one we use essentially exact matrix-product-operator simulations on large systems, and pushing these calculations to dimension two, we obtain accurate results on small cylinders. We also employ an approximation scheme based on solving the dynamics of the mean field dressed by the feedback of quantum fluctuations at leading order. This approach allows us to study the effect of correlations in large lattices with over one hundred thousand spins, as the spatial dimension is increased up to five. In dimension two and higher we find two new states that are stabilized by quantum correlations and do not exist in the mean-field limit of the model. One of these is a steady state with mean magnetization values that lie between the two bistable mean-field values and whose correlation functions have properties reminiscent of both. The correlation length of the new phase diverges at a critical point, beyond which we find emerging a new limit cycle state with the magnetization and correlators oscillating periodically in time.

arXiv:2001.05474•Physical review. B./Physical review. B•Aug 3, 2020•16 citations

Beyond mean-field bistability in driven-dissipative lattices: Bunching-antibunching transition and quantum simulation

Juan José Mendoza-Arenas, Stephen R. L. Clark, Simone Felicetti +4 more

In the present work we investigate the existence of multiple nonequilibrium steady states in a coherently driven $XY$ lattice of dissipative two-level systems. A commonly used mean-field ansatz, in which spatial correlations are neglected, predicts a bistable behavior with a sharp shift between low- and high-density states. In contrast one-dimensional matrix product methods reveal these effects to be artifacts of the mean-field approach, with both disappearing once correlations are taken fully into account. Instead, a bunching-antibunching transition emerges. This indicates that alternative approaches should be considered for higher spatial dimensions, where classical simulations are currently infeasible. Thus we propose a circuit QED quantum simulator implementable with current technology to enable an experimental investigation of the model considered.

arXiv:1510.06651•Physical Review A•Feb 16, 2016•89 citations

Exact multistability and dissipative time crystals in interacting fermionic lattices

Hadiseh Alaeian, Berislav Buča

Abstract The existence of multistability in quantum systems beyond the mean-field approximation remains an intensely debated open question. Quantum fluctuations are finite-size corrections to the mean-field as the full exact solution is unobtainable and they usually destroy the multistability present on the mean-field level. Here, by identifying and using exact modulated dynamical symmetries in a driven-dissipative fermionic chain we exactly prove multistability in the presence of quantum fluctuations. Further, unlike common cases in our model, rather than destroying multistability, the quantum fluctuations themselves exhibit multistability, which is absent on the mean-field level for our systems. Moreover, the studied model acquires additional thermodynamic dynamical symmetries that imply persistent periodic oscillations, constituting the first case of a boundary time crystal,to the best of our knowledge, a genuine extended many-body quantum system with the previous cases being only in emergent single- or few-body models. The model can be made into a dissipative time crystal in the limit of large dissipation (i.e. the persistent oscillations are stabilized by the dissipation) making it both a boundary and dissipative time crystal.

Communications Physics•Dec 7, 2022•22 citations

Transitions de phases dissipatives dans les systèmes quantiques ouverts sur reseaux.

Florent Storme

This thesis is about the theory of open quantum many-body physics with a particular focus on driven-dissipative photonic and spin lattices. After a review of the main physical platforms and theoretical concepts, we describe our original results. In Chapter 2, we present a theory of the single-mode Kerr model with a time-dependent pump. This model describes a single-mode optical cavity with a third-order optical nonlinearity. In the regime of parameters where the semiclassical analysis shows bistability, we find that a dynamic hysteresis loop appears in the exact solution. The hysteresis area as a function of the sweep time shows a double power-law decay where the second exponent is independent of the system parameters. We show how these effects are related to the emergence of a dissipative phase transition. We also describe concisely the experimental results which have recently confirmed such theoretical predictions. In Chapter 3, we present the corner-space renormalization method. In order to obtain the steady-state density matrix of the lattice, we solve the Lindblad master equation in a subspace of the Hilbert space (the "corner"). The states spanning the corner space are selected iteratively using eigenvectors of the density matrix of smaller lattice systems, merging in real space two lattices at each iteration and selecting M pairs of states by maximizing their joint probability. The accuracy of the results is then improved by increasing the dimension M of the corner space until convergence is reached. The method has been benchmarked on a two-dimensional Bose-Hubbard model with coherent driving. The strength and limitations of the method are critically discussed.In Chapter 4, we investigate a dissipative phase transition in the two-dimensional anisotropic Heisenberg XYZ model. Using the corner-space renormalization method we present a finite-size analysis of steady-state observables. In particular, we show the critical behaviour of the magnetic susceptibility, the entropy growth and the entanglement witnesses, providing a first evaluation of the critical exponents characterizing the transition. A study of the dynamics of finite-size systems is also consistent with a critical slowing down. For comparison, we present the corresponding analysis for one-dimensional arrays, showing the absence of criticality due to the reduced dimension.In Chapter 5, we explore the non-equilibrium photonic phases of a dissipative Bose-Hubbard model with incoherent pumping of coupled two-level systems. Within a Gutzwiller mean-field approach, we determine the steady-state phase diagram of the system. We predict a second-order phase transition between an incompressible Mott-like phase and a coherent delocalized phase.

•Nov 24, 2017•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.

•Jan 1, 2017•127 citations

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

João C. Getelina, Andrew Cox, Muhammad Asaduzzaman +4 more

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.

quant-pharXiv:2609.16108•Sep 16, 2026•0 citations
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.
We study a driven-dissipative model of spins one-half (qubits) on a lattice with nearest-neighbor interactions. Focusing on the role of spatially extended spin-spin correlations in determining the phases of the system, we characterize the spatial structure of the correlations in the steady state, as well as their temporal dynamics. In dimension one we use essentially exact matrix-product-operator simulations on large systems, and pushing these calculations to dimension two, we obtain accurate results on small cylinders. We also employ an approximation scheme based on solving the dynamics of the mean field dressed by the feedback of quantum fluctuations at leading order. This approach allows us to study the effect of correlations in large lattices with over one hundred thousand spins, as the spatial dimension is increased up to five. In dimension two and higher we find two new states that are stabilized by quantum correlations and do not exist in the mean-field limit of the model. One of these is a steady state with mean magnetization values that lie between the two bistable mean-field values and whose correlation functions have properties reminiscent of both. The correlation length of the new phase diverges at a critical point, beyond which we find emerging a new limit cycle state with the magnetization and correlators oscillating periodically in time.
In the present work we investigate the existence of multiple nonequilibrium steady states in a coherently driven XYXYXY lattice of dissipative two-level systems. A commonly used mean-field ansatz, in which spatial correlations are neglected, predicts a bistable behavior with a sharp shift between low- and high-density states. In contrast one-dimensional matrix product methods reveal these effects to be artifacts of the mean-field approach, with both disappearing once correlations are taken fully into account. Instead, a bunching-antibunching transition emerges. This indicates that alternative approaches should be considered for higher spatial dimensions, where classical simulations are currently infeasible. Thus we propose a circuit QED quantum simulator implementable with current technology to enable an experimental investigation of the model considered.
Abstract The existence of multistability in quantum systems beyond the mean-field approximation remains an intensely debated open question. Quantum fluctuations are finite-size corrections to the mean-field as the full exact solution is unobtainable and they usually destroy the multistability present on the mean-field level. Here, by identifying and using exact modulated dynamical symmetries in a driven-dissipative fermionic chain we exactly prove multistability in the presence of quantum fluctuations. Further, unlike common cases in our model, rather than destroying multistability, the quantum fluctuations themselves exhibit multistability, which is absent on the mean-field level for our systems. Moreover, the studied model acquires additional thermodynamic dynamical symmetries that imply persistent periodic oscillations, constituting the first case of a boundary time crystal,to the best of our knowledge, a genuine extended many-body quantum system with the previous cases being only in emergent single- or few-body models. The model can be made into a dissipative time crystal in the limit of large dissipation (i.e. the persistent oscillations are stabilized by the dissipation) making it both a boundary and dissipative time crystal.
This thesis is about the theory of open quantum many-body physics with a particular focus on driven-dissipative photonic and spin lattices. After a review of the main physical platforms and theoretical concepts, we describe our original results. In Chapter 2, we present a theory of the single-mode Kerr model with a time-dependent pump. This model describes a single-mode optical cavity with a third-order optical nonlinearity. In the regime of parameters where the semiclassical analysis shows bistability, we find that a dynamic hysteresis loop appears in the exact solution. The hysteresis area as a function of the sweep time shows a double power-law decay where the second exponent is independent of the system parameters. We show how these effects are related to the emergence of a dissipative phase transition. We also describe concisely the experimental results which have recently confirmed such theoretical predictions. In Chapter 3, we present the corner-space renormalization method. In order to obtain the steady-state density matrix of the lattice, we solve the Lindblad master equation in a subspace of the Hilbert space (the "corner"). The states spanning the corner space are selected iteratively using eigenvectors of the density matrix of smaller lattice systems, merging in real space two lattices at each iteration and selecting M pairs of states by maximizing their joint probability. The accuracy of the results is then improved by increasing the dimension M of the corner space until convergence is reached. The method has been benchmarked on a two-dimensional Bose-Hubbard model with coherent driving. The strength and limitations of the method are critically discussed.In Chapter 4, we investigate a dissipative phase transition in the two-dimensional anisotropic Heisenberg XYZ model. Using the corner-space renormalization method we present a finite-size analysis of steady-state observables. In particular, we show the critical behaviour of the magnetic susceptibility, the entropy growth and the entanglement witnesses, providing a first evaluation of the critical exponents characterizing the transition. A study of the dynamics of finite-size systems is also consistent with a critical slowing down. For comparison, we present the corresponding analysis for one-dimensional arrays, showing the absence of criticality due to the reduced dimension.In Chapter 5, we explore the non-equilibrium photonic phases of a dissipative Bose-Hubbard model with incoherent pumping of coupled two-level systems. Within a Gutzwiller mean-field approach, we determine the steady-state phase diagram of the system. We predict a second-order phase transition between an incompressible Mott-like phase and a coherent delocalized phase.
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.
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.