arXiv:2610.12017

Model Order Reduction for Parametric Dissipative Quantum Systems

Bernard Haasdonk, Benjamin Stamm, Zhuoyao Zeng

Published Oct 8, 2026• 0 citations

Abstract

We consider reduced basis approximations of non-equilibrium steady states of parameter-dependent Lindbladians, whose ambient matrix space grows exponentially with system size. We give a unified structural view of the problem through its characterizations as an eigenvalue problem or as a trace-constrained source problem. Building on these characterizations, we propose a residual-driven reduced basis framework for both formulations. Under the assumption of a unique steady state, we derive residual-based state-error estimates, and establish (sub-)exponential convergence of the greedy algorithms under real-analytic parameter dependence of the problem data. Numerical experiments on boundary-driven Fermi--Hubbard and dissipative XYZ chains support the analysis and suggest that the trace-constrained greedy method offers the best overall balance of accuracy and computational cost.

DOI: 10.48550/arxiv.2610.12017

Related papers

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.07709Physical Review XAug 7, 20260 citations

Approximate Reduced Lindblad Dynamics via Algebraic and Adiabatic Methods

Tommaso Grigoletto, Alain Sarlette, Francesco Ticozzi +1 more

We present an algebraic framework for approximate model reduction of Markovian open quantum dynamics that guarantees complete positivity and trace preservation by construction. First, we show that projecting a Lindblad generator on its center manifold -- the space spanned by eigenoperators with purely imaginary eigenvalue -- yields an asymptotically exact reduced quantum dynamical semigroup whose dynamics is unitary, with exponentially decaying transient error controlled by the generator's spectral gap. Second, for analytic perturbations of a Lindblad generator with a tractable center manifold, we propose a perturbative reduction that keeps the reduced space fixed at the unperturbed center manifold. The resulting generator is shown to remain a valid Lindbladian; explicit finite-time error bounds that quantify leakage from the unperturbed center sector are provided. We further clarify the connection to adiabatic elimination methods, by both showing how the algebraic reduction can be directly related to a first-order adiabatic-elimination and by providing sufficient conditions under which the latter method can be applied while preserving complete positivity. We showcase the usefulness of our techniques in dissipative many-body quantum systems exhibiting non-stationary long-time dynamics.

quant-pharXiv:2603.11982Quantum Science and TechnologyMar 13, 20260 citations

Steady-state solution methods for open quantum optical systems

Paul D. Nation

We discuss the numerical solution methods available when solving for the steady-state density matrix of a time-independent open quantum optical system, where the system operators are expressed in a suitable basis representation as sparse matrices. In particular, we focus on the difficulties posed by the non-Hermitian structure of the Lindblad super operator, and the numerical techniques designed to mitigate these pitfalls. In addition, we introduce a doubly iterative inverse-power method that can give reduced memory and runtime requirements in situations where other iterative methods are limited due to poor bandwidth and profile reduction. The relevant methods are demonstrated on several prototypical quantum optical systems where it is found that iterative methods based on iLU factorization using reverse Cuthill-Mckee ordering tend to outperform other solution techniques in terms of both memory consumption and runtime as the size of the underlying Hilbert space increases. For eigenvalue solving, Krylov iterations using the stabilized bi-conjugate gradient method outperform generalized minimal residual methods. In contrast, minimal residual methods work best for solvers based on direct LU decomposition. This work serves as a guide for solving the steady-state density matrix of an arbitrary quantum optical system, and points to several avenues of future research that will extend the applicability of these classical algorithms in absence of a quantum computer.

arXiv:1504.06768Apr 25, 20157 citations

A method to investigate many-body Liouvillian spectra beyond the steady state

Philipp Westhoff

Open quantum systems, where quantum dynamics is influenced by (non-)Markovian environments, have attracted increasing attention due to their rich physical properties. Understanding phenomena such as non-equilibrium phases, anomalous thermalization, and dissipative state preparation requires the knowledge of the low-lying eigenvalues and eigenstates of the Lindbladian. In this thesis, we introduce a framework that systematically computes not only steady states but also low-lying excited states in large, dissipative, interacting quantum many-body systems. The framework is based on tensor-network methods and utilizes recent advances in complex-time Krylov spaces. Extending these ideas to the challenging non-Hermitian eigenvalue problem ubiquitous in open quantum systems enables simulations at unprecedented Hilbert space dimensions. We employ the framework to analyze a dissipative state preparation protocol for Bose-Einstein condensates (BECs) on optical lattices. Analytical and numerical studies allow us to characterize the steady state and perform a finite-size scaling of the dissipative gap, uncovering strong evidence for non-linear hydrodynamics described by the Kardar-Parisi-Zhang universality class. Furthermore, we show analytically that the dissipative state preparation can be substantially accelerated via the quantum Mpemba effect. Our approach exploits weak symmetries to analytically identify a class of simple, experimentally realizable states that converge exponentially faster to the steady state than typical random initializations. Overall, this work establishes a versatile framework for the spectral analysis of generic open quantum many-body systems, enabling the analysis of complex collective phenomena and additionally paving the way toward faster dissipative preparation of highly entangled states in analog quantum simulators.

Jan 1, 20250 citations

Reduced Basis Method for Driven-Dissipative Quantum Systems

Hans Christiansen, V. Baran, Jens Paaske

Reduced basis methods provide an efficient way of mapping out phase diagrams of strongly correlated many-body quantum systems. The method relies on using the exact solutions at select parameter values to construct a low-dimensional basis from which observables can be efficiently and reliably computed throughout the parameter space. Here, we show that this method can be generalized to driven-dissipative Markovian systems, allowing efficient calculations of observables in the transient and steady states. A subsequent distillation of the reduced basis vectors according to their explained variances allows for an unbiased exploration of the most pronounced parameter dependencies indicative of phase boundaries in the thermodynamic limit.

arXiv:2505.05460Physical Review LettersOct 9, 20252 citations

Exact metastability in a class of driven-dissipative quantum many-body systems

David D. Noachtar, Aashish A. Clerk

Metastability in many-body quantum systems and its associated exponentially-long timescales have been the subject of considerable recent interest. Here, we focus on a class of driven-dissipative many-body open quantum systems described by a Lindbladian having hidden time-reversal symmetry (a form of quantum detailed balance). Examples include boundary-driven interacting spin chains, bosonic lattice models and driven-dissipative collective spin models. We suggest that for such systems, slow timescales in the vicinity of a dissipative first-order phase transition can be analytically predicted using a special purification of the non-equilibrium steady state. We show the accuracy of our conjecture through detailed studies of a dissipative transverse-field Ising model with collective and local decay, and a driven-dissipative nonlinear cavity model. Our results allow quantitative insights into metastability and slow dynamics for a range of systems, including cases where semiclassical or path-integral instanton approaches are intractable.

quant-pharXiv:2606.07736Jun 9, 20260 citations