We provide a framework for understanding dynamical metastability in open many-body systems of free bosons, whereby the dynamical stability properties of the system in the infinite-size (thermodynamic) limit may sharply differ from those of any finite-size truncation, and anomalous transient dynamics may arise. By leveraging pseudospectral techniques, we trace the discrepancy between asymptotic and transient dynamics to the non-normality of the underlying quadratic bosonic Lindbladian (QBL) generator and show that two distinct flavors of dynamical metastability can arise. QBLs exhibiting type I dynamical metastability, previously discussed in the context of anomalous transient amplification [Phys. Rev. Lett. 127, 245701 (2021)], are dynamically unstable in the infinite-size limit yet stable once open boundaries are imposed. In contrast, type II dynamically metastable QBLs, which we uncover in this work, are dynamically stable for infinite size but become unstable under open boundary conditions for arbitrary finite system size. We exhibit representative models for both types of metastability in the dissipative, as well as the closed-system (Hamiltonian) settings, and analyze distinctive behavior they can engender. We show that dynamical metastability manifests itself in the generation of entanglement entropy by way of a transient which reflects the stability phase of the infinite (rather than the actual finite) system and, as a result, is directly tied to the emergence of supervolume entanglement scaling in type I systems. Finally, we demonstrate how, even in Hermitian, and especially in highly non-normal regimes, the spectral properties of an infinite-size QBL are reflected in the linear response functions of the corresponding finite QBLs by way of resonant pseudospectral modes.We provide a framework for understanding dynamical metastability in open many-body systems of free bosons, whereby the dynamical stability properties of the system in the infinite-size (thermodynamic) limit may sharply differ from those of any finite-size truncation, and anomalous transient dynamics may arise. By leveraging pseudospectral techniques, we trace the discrepancy between asymptotic and transient dynamics to the non-normality of the underlying quadratic bosonic Lindbladian (QBL) generator and show that two distinct flavors of dynamical metastability can arise. QBLs exhibiting type I dynamical metastability, previously discussed in the context of anomalous transient amplification [Phys. Rev. Lett. 127, 245701 (2021)], are dynamically unstable in the infinite-size limit yet stable once open boundaries are imposed. In contrast, type II dynamically metastable QBLs, which we uncover in this work, are dynamically stable for infinite size but become unstable under open boundary conditions for arbitrary finite system size. We exhibit representative models for both types of metastability in the dissipative, as well as the closed-system (Hamiltonian) settings, and analyze distinctive behavior they can engender. We show that dynamical metastability manifests itself in the generation of entanglement entropy by way of a transient which reflects the stability phase of the infinite (rather than the actual finite) system and, as a result, is directly tied to the emergence of supervolume entanglement scaling in type I systems. Finally, we demonstrate how, even in Hermitian, and especially in highly non-normal regimes, the spectral properties of an infinite-size QBL are reflected in the linear response functions of the corresponding finite QBLs by way of resonant pseudospectral modes.

arXiv:2405.08873Physical Review ASep 9, 202411 citations

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

Leonardo da Silva Souza, Fernando Iemini

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 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.

arXiv:2408.05302SciPost PhysicsJan 16, 20261 citation

Further Bosonic Representations

WORLD SCIENTIFIC eBooksApr 27, 20150 citations

Quantum Gibbs sampling through the detectability lemma

Di Fang, Jianfeng Lu, Yu Tong +1 more

Gibbs state preparation is an important subroutine in quantum computing. In this work we use the detectability lemma to improve Gibbs state preparation. Specifically, we design new Gibbs state preparation methods that do not rely on simulating Lindbladian evolution, thus avoiding the overhead from it. For local Lindbladians consisting of $M$ terms, this approach reduces the cost by a factor of $O(M)$ . We also combine the detectability lemma operator and quantum singular value transformation to implement ground state projection operators of frustration-free Hamiltonians, resulting in a quadratic speedup in the spectral gap dependence. Applying this method to Lindbladians for the Gibbs state of local commuting Hamiltonians, we achieve quadratically better dependence on the Lindbladian spectral gap.Gibbs state preparation is an important subroutine in quantum computing. In this work we use the detectability lemma to improve Gibbs state preparation. Specifically, we design new Gibbs state preparation methods that do not rely on simulating Lindbladian evolution, thus avoiding the overhead from it. For local Lindbladians consisting of $M$ terms, this approach reduces the cost by a factor of $O(M)$. We also combine the detectability lemma operator and quantum singular value transformation to implement ground state projection operators of frustration-free Hamiltonians, resulting in a quadratic speedup in the spectral gap dependence. Applying this method to Lindbladians for the Gibbs state of local commuting Hamiltonians, we achieve quadratically better dependence on the Lindbladian spectral gap.

quant-pharXiv:2604.07214Apr 9, 20260 citations

Extreme bosonic linear channels

Alexander Semenovich Holevo

Theoretical and Mathematical PhysicsFeb 1, 201316 citations

M

terms, this approach reduces the cost by a factor of

O(M)

. We also combine the detectability lemma operator and quantum singular value transformation to implement ground state projection operators of frustration-free Hamiltonians, resulting in a quadratic speedup in the spectral gap dependence. Applying this method to Lindbladians for the Gibbs state of local commuting Hamiltonians, we achieve quadratically better dependence on the Lindbladian spectral gap.