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:2410.13646 · Physical Review Research

Certifying steady-state properties of open quantum systems

Luke Mortimer, Donato Farina, G. Di Bello, David Jansen, Andreas Leitherer, Pere Mujal, Antonio Acín

Published Aug 16, 2025• 9 citations
arXivHTMLPDF

Abstract

Estimating the steady-state properties of open many-body quantum systems is a fundamental challenge in quantum science and technologies. In this work, we present a scalable approach based on semidefinite programming to derive certified bounds on the expectation value of an arbitrary observable in the steady state of Lindbladian dynamics. We illustrate our method on a series of many-body systems, including paradigmatic spin-1/2 chains and two-dimensional ladders, considering both equilibrium and nonequilibrium steady states. We benchmark our method with state-of-the-art tensor-network approaches that, unlike our method, are only able to provide estimates, with no guarantee, on steady-state quantities. For the tested models, only modest computational effort is needed to obtain certified nontrivial bounds for system sizes intractable by exact methods. Our method introduces the first general numerical tool for bounding steady-state properties of open quantum dynamics, opening an avenue in the understanding of stable configurations in many-body systems.

DOI
10.1103/hbrt-cn8q

Related papers

Bounding many-body properties under partial information and finite measurement statistics

Luke Mortimer, Leonardo Zambrano, Antonio Acín +1 more

Calculating bounds of properties of many-body quantum systems is of paramount importance, since they guide our understanding of emergent quantum phenomena and complement the insights obtained from estimation methods. Recent semidefinite programming approaches enable probabilistic bounds from finite-shot measurements of easily accessible, yet informationally incomplete, observables. Here we render these methods scalable in the number of qubits by instead utilizing moment-matrix relaxations. After introducing the general formalism, we show how the approach can be adapted with specific knowledge of the system, such as it being the ground state of a given Hamiltonian, possessing specific symmetries or being the steady state of a given Lindbladian. Our approach defines a scalable real-world certification scheme leveraging semidefinite programming relaxations and experimental estimations which, unavoidably, contain shot noise.

quant-pharXiv:2601.10408•Jan 16, 2026•0 citations

Certifying Ground-State Properties of Many-Body Systems

Jie Wang, Jacopo Surace, Irénée Frérot +4 more

A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin- 1 / 2 Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way. Published by the American Physical Society 2024

arXiv:2310.05844•Physical Review X•Jul 11, 2024•15 citations

Local Certification of Many-Body Steady States

Miguel Frías Pérez, Antonio Acín

We present a relaxation-based method to bound expectation values on the steady state of dissipative many-body quantum systems described by master equations of the Lindblad form. Instead of targeting to represent the entire state, we promote the reduced density matrices to our variables and enforce the constraints that are imposed on them by consistency with a global steady state. The resulting constraints have the form of a semidefinite program, which allows us to efficiently bound the values a given expectation value can take in the steady state. Our results show fast convergence of the bounds with the size of the reduced density matrices, giving very competitive predictions for the steady state of several one- and two-dimensional models for an arbitrary number of particles.

quant-pharXiv:2602.06487•Feb 9, 2026•1 citation

Convex Optimization for Nonequilibrium Steady States on a Hybrid Quantum Processor

Jonathan Wei Zhong Lau, Kian Hwee Lim, Kishor Bharti +2 more

Finding the transient and steady state properties of open quantum systems is a central problem in various fields of quantum technologies. Here, we present a quantum-assisted algorithm to determine the steady states of open system dynamics. By reformulating the problem of finding the fixed point of Lindblad dynamics as a feasibility semidefinite program, we bypass several well-known issues with variational quantum approaches to solving for steady states. We demonstrate that our hybrid approach allows us to estimate the steady states of higher dimensional open quantum systems and discuss how our method can find multiple steady states for systems with symmetries.

arXiv:2204.03203•Physical Review Letters•Jun 12, 2023•11 citations

Tensor-Network Framework for Lindbladian Spectra and Steady States

Philipp Westhoff, Mattia Moroder, Ulrich Schollwöck +1 more

Quantum systems coupled to (non-) Markovian environments attract increasing attention due to their peculiar physical properties. Exciting prospects such as unconventional nonequilibrium phases beyond the Mermin-Wagner limit or dissipative state preparation demand a systematic analysis of quantum many-body phases out of equilibrium. Akin to the equilibrium case, this requires the computation of the low-lying eigenstates of Lindbladians, a problem challenging conventional approaches for simulating quantum many-body systems. Here, we undertake a first step to overcome this limitation and introduce a tensor-network-based framework to systematically compute not only steady states, but also low-lying excited states for large, driven quantum many-body systems. Our framework is based on recent advances utilizing complex-time Krylov spaces, and we leverage these ideas to create a toolbox tailored to solve the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems. With an example of the interacting Bose-Hubbard model driven by dissipation-assisted hopping, we demonstrate the high efficiency and accuracy. From a reliable finite-size-scaling analysis of the spectral gap, we find strong evidence for nonlinear hydrodynamic behavior consistent with Kardar-Parisi-Zhang-type superdiffusive relaxation and establish the existence of exponentially accelerated, anomalous relaxation. This method unlocks the capability of spectral analysis of generic open quantum many-body systems suitable also for non-Markovian environments.

arXiv:2509.07709•Physical Review X•Aug 7, 2026•0 citations

Spin- 1/2 XXZ chain coupled to two Lindblad baths: Constructing nonequilibrium steady states from equilibrium correlation functions

Tjark Heitmann, Jonas Richter, Fengping Jin +7 more

State-of-the-art approaches to extract transport coefficients of many-body quantum systems broadly fall into two categories: (i) they target the linear-response regime in terms of equilibrium correlation functions of the closed system; or (ii) they consider an open-system situation typically modeled by a Lindblad equation, where a nonequilibrium steady state emerges from driving the system at its boundaries. While quantitative agreement between (i) and (ii) has been found for selected model and parameter choices, also disagreement has been pointed out in the literature. Studying magnetization transport in the spin-$1/2$ XXZ chain, we here demonstrate that at weak driving, the nonequilibrium steady state in an open system, including its buildup in time, can remarkably be constructed just on the basis of correlation functions in the closed system. We numerically illustrate this direct correspondence of closed-system and open-system dynamics, and show that it allows the treatment of comparatively large open systems, usually only accessible to matrix product state simulations. We also point out potential pitfalls when extracting transport coefficients from nonequilibrium steady states in finite systems.

arXiv:2303.00430•Physical review. B./Physical review. B•Nov 27, 2023•6 citations
Estimating the steady-state properties of open many-body quantum systems is a fundamental challenge in quantum science and technologies. In this work, we present a scalable approach based on semidefinite programming to derive certified bounds on the expectation value of an arbitrary observable in the steady state of Lindbladian dynamics. We illustrate our method on a series of many-body systems, including paradigmatic spin-1/2 chains and two-dimensional ladders, considering both equilibrium and nonequilibrium steady states. We benchmark our method with state-of-the-art tensor-network approaches that, unlike our method, are only able to provide estimates, with no guarantee, on steady-state quantities. For the tested models, only modest computational effort is needed to obtain certified nontrivial bounds for system sizes intractable by exact methods. Our method introduces the first general numerical tool for bounding steady-state properties of open quantum dynamics, opening an avenue in the understanding of stable configurations in many-body systems.
Calculating bounds of properties of many-body quantum systems is of paramount importance, since they guide our understanding of emergent quantum phenomena and complement the insights obtained from estimation methods. Recent semidefinite programming approaches enable probabilistic bounds from finite-shot measurements of easily accessible, yet informationally incomplete, observables. Here we render these methods scalable in the number of qubits by instead utilizing moment-matrix relaxations. After introducing the general formalism, we show how the approach can be adapted with specific knowledge of the system, such as it being the ground state of a given Hamiltonian, possessing specific symmetries or being the steady state of a given Lindbladian. Our approach defines a scalable real-world certification scheme leveraging semidefinite programming relaxations and experimental estimations which, unavoidably, contain shot noise.
A ubiquitous problem in quantum physics is to understand the ground-state properties of many-body systems. Confronted with the fact that exact diagonalization quickly becomes impossible when increasing the system size, variational approaches are typically employed as a scalable alternative: Energy is minimized over a subset of all possible states and then different physical quantities are computed over the solution state. Despite remarkable success, rigorously speaking, all that variational methods offer are upper bounds on the ground-state energy. On the other hand, so-called relaxations of the ground-state problem based on semidefinite programming represent a complementary approach, providing lower bounds to the ground-state energy. However, in their current implementation, neither variational nor relaxation methods offer provable bound on other observables in the ground state beyond the energy. In this work, we show that the combination of the two classes of approaches can be used to derive certifiable bounds on the value of any observable in the ground state, such as correlation functions of arbitrary order, structure factors, or order parameters. We illustrate the power of this approach in paradigmatic examples of 1D and 2D spin- 1 / 2 Heisenberg models. To improve the scalability of the method, we exploit the symmetries and sparsity of the considered systems to reach sizes of hundreds of particles at much higher precision than previous works. Our analysis therefore shows how to obtain certifiable bounds on many-body ground-state properties beyond energy in a scalable way. Published by the American Physical Society 2024
We present a relaxation-based method to bound expectation values on the steady state of dissipative many-body quantum systems described by master equations of the Lindblad form. Instead of targeting to represent the entire state, we promote the reduced density matrices to our variables and enforce the constraints that are imposed on them by consistency with a global steady state. The resulting constraints have the form of a semidefinite program, which allows us to efficiently bound the values a given expectation value can take in the steady state. Our results show fast convergence of the bounds with the size of the reduced density matrices, giving very competitive predictions for the steady state of several one- and two-dimensional models for an arbitrary number of particles.
Finding the transient and steady state properties of open quantum systems is a central problem in various fields of quantum technologies. Here, we present a quantum-assisted algorithm to determine the steady states of open system dynamics. By reformulating the problem of finding the fixed point of Lindblad dynamics as a feasibility semidefinite program, we bypass several well-known issues with variational quantum approaches to solving for steady states. We demonstrate that our hybrid approach allows us to estimate the steady states of higher dimensional open quantum systems and discuss how our method can find multiple steady states for systems with symmetries.
Quantum systems coupled to (non-) Markovian environments attract increasing attention due to their peculiar physical properties. Exciting prospects such as unconventional nonequilibrium phases beyond the Mermin-Wagner limit or dissipative state preparation demand a systematic analysis of quantum many-body phases out of equilibrium. Akin to the equilibrium case, this requires the computation of the low-lying eigenstates of Lindbladians, a problem challenging conventional approaches for simulating quantum many-body systems. Here, we undertake a first step to overcome this limitation and introduce a tensor-network-based framework to systematically compute not only steady states, but also low-lying excited states for large, driven quantum many-body systems. Our framework is based on recent advances utilizing complex-time Krylov spaces, and we leverage these ideas to create a toolbox tailored to solve the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems. With an example of the interacting Bose-Hubbard model driven by dissipation-assisted hopping, we demonstrate the high efficiency and accuracy. From a reliable finite-size-scaling analysis of the spectral gap, we find strong evidence for nonlinear hydrodynamic behavior consistent with Kardar-Parisi-Zhang-type superdiffusive relaxation and establish the existence of exponentially accelerated, anomalous relaxation. This method unlocks the capability of spectral analysis of generic open quantum many-body systems suitable also for non-Markovian environments.
State-of-the-art approaches to extract transport coefficients of many-body quantum systems broadly fall into two categories: (i) they target the linear-response regime in terms of equilibrium correlation functions of the closed system; or (ii) they consider an open-system situation typically modeled by a Lindblad equation, where a nonequilibrium steady state emerges from driving the system at its boundaries. While quantitative agreement between (i) and (ii) has been found for selected model and parameter choices, also disagreement has been pointed out in the literature. Studying magnetization transport in the spin-1/21/21/2 XXZ chain, we here demonstrate that at weak driving, the nonequilibrium steady state in an open system, including its buildup in time, can remarkably be constructed just on the basis of correlation functions in the closed system. We numerically illustrate this direct correspondence of closed-system and open-system dynamics, and show that it allows the treatment of comparatively large open systems, usually only accessible to matrix product state simulations. We also point out potential pitfalls when extracting transport coefficients from nonequilibrium steady states in finite systems.