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.
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. 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. We demonstrate the efficiency of our method on several examples of coupled bosonic systems.
We study driven and open quantum systems described by a time-periodic Lindblad master equation. In closed systems, the stroboscopic dynamics can always be described by an effective time-independent Floquet Hamiltonian; this idea is the basis of Floquet engineering. However, in the presence of dissipation, the existence of an effective time-independent Floquet Lindbladian is not guaranteed due to the non-unitary nature of the evolution. Using Floquet theory, we construct a well-defined time-independent Floquet Lindbladian in an extended Sambe-Liouville space, transforming the initial time-dependent problem to a static and non-Hermitian eigenvalue problem. For harmonic driving, we introduce a matrix continued fraction method to nonperturbatively resum multiphoton processes and construct an effective Floquet Lindbladian acting only on the physical Liouville space. Compared to other high-frequency expansions, this method has the advantage of providing the whole infinite series expansion at once. Using a resolvent formalism, we show how to obtain a spectral Floquet representation of correlation functions of an open quantum system. As an application, we consider a dissipating two-level system in a linearly polarized field and calculate its resonance fluorescence spectrum. Furthermore, we consider a parametrically driven quantum dot with pump and loss for which we calculate its spectral function and current-voltage characteristics.
Markovian open quantum systems are well modeled by the Lindblad Master Equation (ME)
dtdρt=Lρt, where
L is a linear (super-)operator and
ρt is the system state, a positive matrix. When designing or characterizing a quantum system, one is usually interested in the steady state
ρ∞ (such that
Lρ∞=0), the first few excited states, and trajectories
t↦ρt. In finite dimension,
ρt is an
n×n matrix,
L thus typically costs
n4 to store explicitly as a dense matrix, and
O(n6) to diagonalize or invert exactly, making standard linear algebraic techniques expensive for large systems. However,
L usually costs only
O(n3) to apply. This makes iterative methods appealing, but they do not work without a good preconditioner. In this article, our main observation is that a part of the Lindblad equation, corresponding to the so-called no-jump evolution
S, can be inverted efficiently. Using this inverse map, we introduce an auxiliary completely positive trace-preserving (CPTP) map
Φ whose fixed point is directly related to
ρ∞, all the other eigenvalues having smaller magnitude. The map
Φ is thus well suited to iterative methods, and
ρ∞ can be found in a few Arnoldi iterations. Using the same inverse map
S−1 as preconditioner, we compute the low-lying spectrum efficiently via shift-invert Arnoldi, and, as a proof of concept, build an implicit time integrator that is competitive on stiff systems in the low-precision regime. For the steady-state and low excited states problems, our methods scale like
O(n3) per iteration and offer state-of-the-art performance on CPU and GPU.
In this paper, we investigate periodically driven open quantum systems within the framework of Floquet-Lindblad master equations. Specifically, we discuss Lindblad master equations in the presence of a coherent, time-periodic driving and establish their general spectral features. We also clarify the notions of transient and nondecaying solutions from this spectral perspective, and then prove that any physical system described by a Floquet-Lindblad equation must have at least one physical nonequilibrium steady state (NESS), corresponding to an eigenoperator of the Floquet-Lindblad evolution superoperator
UF with unit eigenvalue. Since the Floquet-Lindblad formalism encapsulates the entire information regarding the NESS, it in principle enables us to obtain nonlinear effects to all orders at once. The Floquet-Lindblad formalism thus provides a powerful tool for studying driven-dissipative solid-state systems, which we illustrate by deriving the nonlinear optical response of a simple two-band model of an insulating solid and comparing it with prior results established through Keldysh techniques.
Floquet engineering of isolated systems is often based on the concept of the effective time-independent Floquet Hamiltonian, which describes the stroboscopic evolution of a periodically driven quantum system in steps of the driving period and which is routinely obtained analytically using high-frequency expansions. The generalization of these concepts to open quantum systems described by a Markovian master equation of Lindblad type turns out to be nontrivial: On the one hand, already for a two-level system two different phases can be distinguished, where the effective time-independent Floquet generator (describing the stroboscopic evolution) is either again Markovian and of Lindblad type or not. On the other hand, even though in the high-frequency regime a Lindbladian Floquet generator (Floquet Lindbladian) is numerically found to exist, this behavior is, curiously, not correctly reproduced within analytical high-frequency expansions. Here, we demonstrate that a proper Floquet Lindbladian can still be obtained from a high-frequency expansion, when treating the problem in a suitably chosen rotating frame. Within this approach, we can then also describe the transition to a phase at lower driving frequencies, where no Floquet Lindbladian exists, and show that the emerging non-Markovianity of the Floquet generator can entirely be attributed to the micromotion of the open driven system.
The stroboscopic evolution of a time-periodically driven isolated quantum system can always be described by an effective time-independent Hamiltonian. Whether this concept can be generalized to open Floquet systems, described by a Markovian Lindblad master equation with a time-periodic generator (Lindbladian), remains an open question. By using a two-level system as a model, we explicitly show the existence of two well-defined parameter regions. In one region the stroboscopic evolution can be described by a Markovian master equation with a time-independent Floquet Lindbladian. In the other it cannot; but here the one-cycle evolution operator can be reproduced with an effective non-Markovian master equation that is homogeneous but nonlocal in time. Interestingly, we find that the boundary between the phases depends on when during the period the evolution is stroboscopically monitored. This reveals the nontrivial role played by the micromotion in the dynamics of open Floquet systems.