In this contribution to the memorial issue of Göran Lindblad, we investigate the periodically driven Lindblad equation for a two-level system. We analyze the system using both adiabatic diagonalization and numerical simulations of the time-evolution, as well as Floquet theory. Adiabatic diagonalization reveals the presence of exceptional points in the system, which depend on the system parameters. We show how the presence of these exceptional points affects the system evolution, leading to a rapid dephasing at these points and a staircase-like loss of coherence. This phenomenon can be experimentally observed by measuring, for example, the population inversion. We also observe that the presence of exceptional points seems to be related to which underlying Lie algebra the system supports. In the Floquet analysis, we map the time-dependent Liouvillian to a non-Hermitian Floquet Hamiltonian and analyze its spectrum. For weak decay rates, we find a Wannier-Stark ladder spectrum accompanied by corresponding Stark-localized eigenstates. For larger decay rates, the ladders begin to dissolve, and new, less localized states emerge. Additionally, their eigenvalues are exponentially sensitive to perturbations, similar to the skin effect found in certain non-Hermitian Hamiltonians.
In this contribution to the memorial issue of Göran Lindblad, we investigate the periodically driven Lindblad equation for a two-level system. We analyze the system using both adiabatic diagonalization and numerical simulations of the time-evolution, as well as Floquet theory. Adiabatic diagonalization reveals the presence of exceptional points in the system, which depend on the system parameters. We show how the presence of these exceptional points affects the system evolution, leading to a rapid dephasing at these points and a staircase-like loss of coherence. This phenomenon can be experimentally observed by measuring, for example, the population inversion. We also observe that the presence of exceptional points seems to be related to which underlying Lie algebra the system supports. In the Floquet analysis, we map the time-dependent Liouvillian to a non-Hermitian Floquet Hamiltonian and analyze its spectrum. For weak decay rates, we find a Wannier-Stark ladder spectrum accompanied by corresponding Stark-localized eigenstates. For larger decay rates, the ladders begin to dissolve, and new, less localized states emerge. Additionally, their eigenvalues are exponentially sensitive to perturbations, similar to the skin effect found in certain non-Hermitian Hamiltonians.
The dynamics of an isolated quantum system is coherent and unitary. Weak coupling to the environment leads to decoherence, which is traditionally modeled with a Lindblad equation for the system's density matrix. Starting from a pure state, such a system approaches a steady state (mixed or otherwise) in an underdamped or overdamped manner. This transition occurs at an eigenvalue degeneracy of a Lindblad superoperator, called an exceptional point (EP), where corresponding eigenvectors coalesce. Recent years have seen an explosion of interest in creating exceptional points in a truly quantum domain, driven by the enhanced sensitivity and topological features EPs have shown in their classical realizations. Here, we present Floquet analysis of a prototypical qubit whose drive or dissipator strengths are varied periodically. We consider models with a single dissipator that generate global loss (phase damping) or mode-selective loss (spontaneous emission). In all cases, we find that periodic modulations lead to EP lines at small dissipator strengths and a rich EP structure in the parameter space. Our analytical and numerical results show that extending Lindblad Liouvillians to the Floquet domain is a potentially preferred route to accessing exceptional points in the transient dynamics towards the Lindblad steady state.
Abstract We study the approach to the adiabatic limit in periodically driven systems. Specifically focusing on a spin-1/2 in a magnetic field we find that, when the parameters of the Hamiltonian lead to a quasi-degeneracy in the Floquet spectrum, the evolution is not adiabatic even if the frequency of the field is much smaller than the spectral gap of the Hamiltonian. We argue that this is a general phenomenon of periodically driven systems. Although an explanation based on a perturbation theory in ω 0 cannot be given, because of the singularity of the zero frequency limit, we are able to describe this phenomenon by means of a mapping to an extended Hilbert space, in terms of resonances of an effective two-band Wannier–Stark ladder. Remarkably, the phenomenon survives in presence of dissipation towards an environment and can be therefore easily experimentally observed.
Recently, defect production was investigated during nonunitary dynamics due to non-Hermitian Hamiltonian. By ramping up the non-Hermitian coupling linearly in time through an exceptional point, defects are produced in much the same way as approaching a Hermitian critical point. A generalized Kibble-Zurek scaling accounted for the ensuing scaling of the defect density in terms of the speed of the drive and the corresponding critical exponents. Here we extend this setting by adding the recycling term and considering the full Lindbladian time evolution of the problem with quantum jumps. We find that by linearly ramping up the environmental coupling in time, and going beyond the steady-state solution of the Liouvillian, the defect density scales linearly with the speed of the drive for all cases. This scaling is unaffected by the presence of exceptional points of the Liouvillian, which can show up in the transient states. By using a variant of the adiabatic perturbation theory, the scaling of the defect density is determined exactly from a set of algebraic equations. Our study indicates the distinct sensitivity of the Lindbladian time evolution to exceptional points corresponding to steady states and transient states.
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.
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.
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.