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:2408.05302 · SciPost Physics

Lindbladian reverse engineering for general non-equilibrium steady states: A scalable null-space approach

Leonardo da Silva Souza, Fernando Iemini

Published Jan 16, 2026• 1 citation
arXivHTMLPDF

Abstract

The study of open system dynamics is of paramount importance both from its fundamental aspects as well as from its potential applications in quantum technologies. In the simpler and most commonly studied case, the dynamics of the system can be described by a Lindblad master equation. However, identifying the Lindbladian that leads to general non-equilibrium steady states (NESS) is usually a non-trivial and challenging task. Here we introduce a method for reconstructing the corresponding Lindbladian master equation given any target NESS, i.e., a Lindbladian Reverse Engineering ( \mathcal{L} ℒ RE) approach. The method maps the reconstruction task to a simple linear problem. Specifically, to the diagonalization of a correlation matrix whose elements are NESS observables and whose size scales linearly (at most quadratically) with the number of terms in the Hamiltonian (Lindblad jump operator) ansatz. The kernel (null-space) of the correlation matrix corresponds to Lindbladian solutions. Moreover, the map defines an iff condition for \mathcal{L} ℒ RE, which works as both a necessary and a sufficient condition; thus, it not only defines, if possible, Lindbladian evolutions leading to the target NESS, but also determines the feasibility of such evolutions in a proposed setup. We illustrate the method in different systems, ranging from bosonic Gaussian systems, dissipative-driven collective spins and random local spin models.

DOI
10.21468/scipostphys.20.1.009

Related papers

Third quantization: a general method to solve master equations for quadratic open Fermi systems

Tomaž Prosen

The Lindblad master equation for an arbitrary quadratic system of n fermions is solved explicitly in terms of diagonalization of a 4 n × 4 n matrix, provided that all Lindblad bath operators are linear in the fermionic variables. The method is applied to the explicit construction of non-equilibrium steady states (NESS) and the calculation of asymptotic relaxation rates in the far from equilibrium problem of heat and spin transport in a nearest neighbour Heisenberg XY spin-1/2 chain in a transverse magnetic field.

arXiv:0801.1257•New Journal of Physics•Apr 15, 2008•502 citations

Non-diagonal Lindblad master equations in quantum reservoir engineering

Diego N. Bernal-García, Lujun Huang, Andrey E. Miroshnichenko +1 more

Reservoir engineering has proven to be a practical approach to control open quantum systems, preserving quantum coherence by appropriately manipulating the reservoir and system-reservoir interactions. In this context, for systems comprised of different parts, it is common to describe the dynamics of a subsystem of interest by performing an adiabatic elimination of the remaining components of the system. This procedure often leads to an effective master equation for the subsystem that is not in the diagonal form of the Gorini-Kossakowski-Lindblad-Sudarshan master equation (here called diagonal Lindblad form). Instead, it has a more general structure (here called non-diagonal Lindblad form), which explicitly reveals the dissipative coupling between the various components of the subsystem. In this work, we present a set of dynamical equations for the first and second moments of the canonical variables for linear Gaussian systems, bosonic and fermionic, described by non-diagonal Lindblad master equations. Our method is efficient and allows one to obtain analytical solutions for the steady state. We supplement our findings with a review of covariance matrix methods, focusing on those related to the measurement of entanglement. Notably, our exploration yields a surprising byproduct: the Duan criterion, commonly applied to bosonic systems for verification of entanglement, is found to be equally valid for fermionic systems. We conclude with a practical example, where we revisit two-mode mechanical entanglement in an optomechanical setup. Our approach, which employs adiabatic elimination for systems governed by time-dependent Hamiltonians, opens the door to examine physical regimes that have not been explored before.

arXiv:2111.04041•Nov 7, 2021•0 citations

Magnetization in a nonequilibrium quantum spin system

X. Z. Zhang

The dynamics described by the non-Hermitian Hamiltonian typically capture the short-term behavior of open quantum systems before quantum jumps occur. In contrast, the long-term dynamics, characterized by the Lindblad master equation, drive the system towards a nonequilibrium steady state (NESS), which is an eigenstate with zero energy of the Liouvillian superoperator, denoted as $\mathcal{L}$. Conventionally, these two types of evolutions exhibit distinct dynamical behaviors. However, in this study we challenge this common belief and demonstrate that the effective non-Hermitian Hamiltonian can accurately represent the long-term dynamics of a critical two-level open quantum system. The criticality of the system arises from the exceptional point of the effective non-Hermitian Hamiltonian. Additionally, the NESS is identical to the coalescent state of the effective non-Hermitian Hamiltonian. We apply this finding to a series of critical open quantum systems and show that a local dissipation channel can induce collective alignment of all spins in the same direction. This direction can be well controlled by modulating the quantum jump operator. The corresponding NESS is a product state and maintains long-time coherence, facilitating quantum control in open many-body systems. This discovery paves the way for a better understanding of the long-term dynamics of critical open quantum systems.

arXiv:2406.00268•Physical review. B./Physical review. B•May 31, 2024•4 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

Arnoldi-Lindblad time evolution: Faster-than-the-clock algorithm for the spectrum of time-independent and Floquet open quantum systems

Fabrizio Minganti, Dolf Huybrechts

The characterization of open quantum systems is a central and recurring problem for the development of quantum technologies. For time-independent systems, an (often unique) steady state describes the average physics once all the transient processes have faded out, but interesting quantum properties can emerge at intermediate timescales. Given a Lindblad master equation, these properties are encoded in the spectrum of the Liouvillian whose diagonalization, however, is a challenge even for small-size quantum systems. Here, we propose a new method to efficiently provide the Liouvillian spectral decomposition. We call this method an Arnoldi-Lindblad time evolution, because it exploits the algebraic properties of the Liouvillian superoperator to efficiently construct a basis for the Arnoldi iteration problem. The advantage of our method is double: (i) It provides a faster-than-the-clock method to efficiently obtain the steady state, meaning that it produces the steady state through time evolution shorter than needed for the system to reach stationarity. (ii) It retrieves the low-lying spectral properties of the Liouvillian with a minimal overhead, allowing to determine both which quantum properties emerge and for how long they can be observed in a system. This method is general and model-independent , and lends itself to the study of large systems where the determination of the Liouvillian spectrum can be numerically demanding but the time evolution of the density matrix is still doable. Our results can be extended to time evolution with a time-dependent Liouvillian. In particular, our method works for Floquet (i.e., periodically driven) systems, where it allows not only to construct the Floquet map for the slow-decaying processes, but also to retrieve the stroboscopic steady state and the eigenspectrum of the Floquet map. Although the method can be applied to any Lindbladian evolution (spin, fermions, bosons, …), for the sake of simplicity we demonstrate the efficiency of our method on several examples of coupled bosonic resonators (as a particular example). Our method outperforms other diagonalization techniques and retrieves the Liouvillian low-lying spectrum even for system sizes for which it would be impossible to perform exact diagonalization.

arXiv:2109.01648•Quantum•Feb 10, 2022•15 citations

Exact solutions of open integrable quantum spin chains

Enej Ilievski

In the thesis we present an analytic approach towards exact description for steady state density operators of nonequilibrium quantum dynamics in the framework of open systems. We employ the so-called quantum Markovian semi-group evolution, i.e. a general form of time-autonomous positivity and trace-preserving dynamical equation for reduced density operators, by only allowing unitarity-breaking dissipative terms acting at the boundaries of a system. Such setup enables to simulate macroscopic reservoirs for different values of effective thermodynamic potentials, causing incoherent transitions between quantum states which are modeled with aid of the Lindblad operators. This serves as a simple minimalistic model for studying quantum transport properties, either in the linear response domain or in more general regimes far from canonical equilibrium. We are mainly exploring possibilities of identifying nonequilibrium situations which are amenable to exact description within matrix product state representation, by exclusively focusing on steady states, i.e. fixed points of the Lindblad equation, of certain prototypic interacting integrable spin chains driven by incoherent polarizing processes. Finally, we define a concept of pseudo-local extensive almost-conserved quantities by allowing a violation of time-invariance up to boundary-localized terms. We elucidate the role of such quantities on non-ergodic behaviour of temporal correlation functions rendering anomalous transport properties. It turns out that such conservation laws can be generated by means of boundary universal quantum transfer operators of the fundamental integrable models.

arXiv:1410.1446•Oct 6, 2014•16 citations
The study of open system dynamics is of paramount importance both from its fundamental aspects as well as from its potential applications in quantum technologies. In the simpler and most commonly studied case, the dynamics of the system can be described by a Lindblad master equation. However, identifying the Lindbladian that leads to general non-equilibrium steady states (NESS) is usually a non-trivial and challenging task. Here we introduce a method for reconstructing the corresponding Lindbladian master equation given any target NESS, i.e., a Lindbladian Reverse Engineering ( \mathcal{L} ℒ RE) approach. The method maps the reconstruction task to a simple linear problem. Specifically, to the diagonalization of a correlation matrix whose elements are NESS observables and whose size scales linearly (at most quadratically) with the number of terms in the Hamiltonian (Lindblad jump operator) ansatz. The kernel (null-space) of the correlation matrix corresponds to Lindbladian solutions. Moreover, the map defines an iff condition for \mathcal{L} ℒ RE, which works as both a necessary and a sufficient condition; thus, it not only defines, if possible, Lindbladian evolutions leading to the target NESS, but also determines the feasibility of such evolutions in a proposed setup. We illustrate the method in different systems, ranging from bosonic Gaussian systems, dissipative-driven collective spins and random local spin models.
The Lindblad master equation for an arbitrary quadratic system of n fermions is solved explicitly in terms of diagonalization of a 4 n × 4 n matrix, provided that all Lindblad bath operators are linear in the fermionic variables. The method is applied to the explicit construction of non-equilibrium steady states (NESS) and the calculation of asymptotic relaxation rates in the far from equilibrium problem of heat and spin transport in a nearest neighbour Heisenberg XY spin-1/2 chain in a transverse magnetic field.
Reservoir engineering has proven to be a practical approach to control open quantum systems, preserving quantum coherence by appropriately manipulating the reservoir and system-reservoir interactions. In this context, for systems comprised of different parts, it is common to describe the dynamics of a subsystem of interest by performing an adiabatic elimination of the remaining components of the system. This procedure often leads to an effective master equation for the subsystem that is not in the diagonal form of the Gorini-Kossakowski-Lindblad-Sudarshan master equation (here called diagonal Lindblad form). Instead, it has a more general structure (here called non-diagonal Lindblad form), which explicitly reveals the dissipative coupling between the various components of the subsystem. In this work, we present a set of dynamical equations for the first and second moments of the canonical variables for linear Gaussian systems, bosonic and fermionic, described by non-diagonal Lindblad master equations. Our method is efficient and allows one to obtain analytical solutions for the steady state. We supplement our findings with a review of covariance matrix methods, focusing on those related to the measurement of entanglement. Notably, our exploration yields a surprising byproduct: the Duan criterion, commonly applied to bosonic systems for verification of entanglement, is found to be equally valid for fermionic systems. We conclude with a practical example, where we revisit two-mode mechanical entanglement in an optomechanical setup. Our approach, which employs adiabatic elimination for systems governed by time-dependent Hamiltonians, opens the door to examine physical regimes that have not been explored before.
The dynamics described by the non-Hermitian Hamiltonian typically capture the short-term behavior of open quantum systems before quantum jumps occur. In contrast, the long-term dynamics, characterized by the Lindblad master equation, drive the system towards a nonequilibrium steady state (NESS), which is an eigenstate with zero energy of the Liouvillian superoperator, denoted as L\mathcal{L}L. Conventionally, these two types of evolutions exhibit distinct dynamical behaviors. However, in this study we challenge this common belief and demonstrate that the effective non-Hermitian Hamiltonian can accurately represent the long-term dynamics of a critical two-level open quantum system. The criticality of the system arises from the exceptional point of the effective non-Hermitian Hamiltonian. Additionally, the NESS is identical to the coalescent state of the effective non-Hermitian Hamiltonian. We apply this finding to a series of critical open quantum systems and show that a local dissipation channel can induce collective alignment of all spins in the same direction. This direction can be well controlled by modulating the quantum jump operator. The corresponding NESS is a product state and maintains long-time coherence, facilitating quantum control in open many-body systems. This discovery paves the way for a better understanding of the long-term dynamics of critical open quantum systems.
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.
The characterization of open quantum systems is a central and recurring problem for the development of quantum technologies. For time-independent systems, an (often unique) steady state describes the average physics once all the transient processes have faded out, but interesting quantum properties can emerge at intermediate timescales. Given a Lindblad master equation, these properties are encoded in the spectrum of the Liouvillian whose diagonalization, however, is a challenge even for small-size quantum systems. Here, we propose a new method to efficiently provide the Liouvillian spectral decomposition. We call this method an Arnoldi-Lindblad time evolution, because it exploits the algebraic properties of the Liouvillian superoperator to efficiently construct a basis for the Arnoldi iteration problem. The advantage of our method is double: (i) It provides a faster-than-the-clock method to efficiently obtain the steady state, meaning that it produces the steady state through time evolution shorter than needed for the system to reach stationarity. (ii) It retrieves the low-lying spectral properties of the Liouvillian with a minimal overhead, allowing to determine both which quantum properties emerge and for how long they can be observed in a system. This method is general and model-independent , and lends itself to the study of large systems where the determination of the Liouvillian spectrum can be numerically demanding but the time evolution of the density matrix is still doable. Our results can be extended to time evolution with a time-dependent Liouvillian. In particular, our method works for Floquet (i.e., periodically driven) systems, where it allows not only to construct the Floquet map for the slow-decaying processes, but also to retrieve the stroboscopic steady state and the eigenspectrum of the Floquet map. Although the method can be applied to any Lindbladian evolution (spin, fermions, bosons, …), for the sake of simplicity we demonstrate the efficiency of our method on several examples of coupled bosonic resonators (as a particular example). Our method outperforms other diagonalization techniques and retrieves the Liouvillian low-lying spectrum even for system sizes for which it would be impossible to perform exact diagonalization.
In the thesis we present an analytic approach towards exact description for steady state density operators of nonequilibrium quantum dynamics in the framework of open systems. We employ the so-called quantum Markovian semi-group evolution, i.e. a general form of time-autonomous positivity and trace-preserving dynamical equation for reduced density operators, by only allowing unitarity-breaking dissipative terms acting at the boundaries of a system. Such setup enables to simulate macroscopic reservoirs for different values of effective thermodynamic potentials, causing incoherent transitions between quantum states which are modeled with aid of the Lindblad operators. This serves as a simple minimalistic model for studying quantum transport properties, either in the linear response domain or in more general regimes far from canonical equilibrium. We are mainly exploring possibilities of identifying nonequilibrium situations which are amenable to exact description within matrix product state representation, by exclusively focusing on steady states, i.e. fixed points of the Lindblad equation, of certain prototypic interacting integrable spin chains driven by incoherent polarizing processes. Finally, we define a concept of pseudo-local extensive almost-conserved quantities by allowing a violation of time-invariance up to boundary-localized terms. We elucidate the role of such quantities on non-ergodic behaviour of temporal correlation functions rendering anomalous transport properties. It turns out that such conservation laws can be generated by means of boundary universal quantum transfer operators of the fundamental integrable models.