We generalize the theory of flow equations to open quantum systems focusing on Lindblad master equations. We introduce and discuss three different generators of the flow that transform a linear non-Hermitian operator into a diagonal one. We first test our dissipative flow equations on a generic matrix and on a physical problem with a driven-dissipative single fermionic mode. We then move to problems with many fermionic modes and discuss the interplay between coherent (disordered) dynamics and localized losses. Our method can also be applied to non-Hermitian Hamiltonians.
We generalize the theory of flow equations to open quantum systems focusing on Lindblad master equations. We introduce and discuss three different generators of the flow that transform a linear non-Hermitian operator into a diagonal one. We first test our dissipative flow equations on a generic matrix and on a physical problem with a driven-dissipative single fermionic mode. We then move to problems with many fermionic modes and discuss the interplay between coherent (disordered) dynamics and localized losses. Our method can also be applied to non-Hermitian Hamiltonians.
We compare two approaches to open quantum systems, namely, the non-Hermitian dynamics and the Lindblad master equation. In order to deal with more general dissipative phenomena, we propose the unified master equation that combines the characteristics of both of these approaches. This allows us to assess the differences between them as well as to clarify which observed features come from the Lindblad or the non-Hermitian part, when it comes to experiment. Using a generic two-mode single-atom laser system as a practical example, we analytically solve the dynamics of the normalized density matrix operator. We study the two-level model in a number of cases (depending on parameters and types of dynamics), compute different observables and study their physical properties. It turns out that one can not only able to describe the different types of damping in dissipative quantum optical systems but also mimic the undamped anharmonic oscillatory phenomena which happen in quantum systems with more than two levels (while staying within the framework of the analytically simple two-mode approximation).
The study of many-body quantum physics has been recently challenged by the appearance of an increasing number of experimental platforms where genuine quantum phenomena take place notwithstanding the presence of an environment and of dissipation. Exciton-polaritons, lossy atomic and molecular gases, cavity QED arrays, arrays of trapped ions and Rydberg atoms, are only few prominent examples of a long list. The present work is ideally divided into two parts, the former is devoted to methodological results; the latter to the theoretical characterization of experimental data. Whereas the notion of equilibrium has been a fruitful guide to the development of standard many-body physics, these setups are inherently out of equilibrium and their modelisation requires the introduction of new theoretical tools. Several methods for addressing this dissipative out-of-equilibrium dynamics have been proposed, based for instance on quantum trajectories, truncated Wigner expansions, tensor networks, machine learning, but the solution of many-body physics for open quantum systems remains a formidable task. Clearly, techniques developed in the framework of Hamiltonian closed systems are a continuous source of inspiration for novel developments, and in this work we present the generalization of one such technique, the so-called flow equations, to the dissipative framework where the generator of the dynamics is a linear non-Hermitian operator; notable examples include non-Hermitian Hamiltonians and Lindblad master equations. In the second part of this work we focus on two-body lossy dynamics in strongly-correlated cold atomic gases. In particular, we derive an effective description for the long-time dynamics of a spin ½ fermionic system by extending a pre-existing theory for the bosonic case.
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.
We show how to solve a large class of Lindblad master equations for noninteracting particles on L sites. In the first portion we concentrate on bosonic particles, while in the second we will address fermionic particles. In both cases we show how to reduce the problem to diagonalizing an $L\ifmmode\times\else\texttimes\fi{}L$ non-Hermitian matrix. In particular, for boundary dissipative driving of a uniform chain, the matrix is a tridiagonal bordered Toeplitz matrix which can be solved analytically for the normal master modes and their relaxation rates (rapidities). In the regimes in which an analytical solution cannot be found, our approach can still provide a speedup in the numerical evaluation. For bosonic particles, we use this numerical method to study the relaxation gap at nonequilibrium phase transitions in a boundary driven bosonic ladder with synthetic gauge fields. We conclude by showing how to construct the nonequilibrium steady state. The analysis for fermionic particles closely follows that of bosons, but with important differences due to the different commutation rules.
Open quantum systems provide an essential theoretical basis for the development of novel quantum technologies, since any real quantum system inevitably interacts with its environment. Lindblad master equations capture the effect of Markovian environments. Closed quantum systems can be treated using flow equations with the particle conserving generator. We generalize this generator to non-Hermitian matrices and open quantum systems governed by Lindbladians, comparing our results with recently proposed generators by Rosso et al. (arXiv:2007.12044). In comparison, we find that our advocated generator provides an efficient flow with good accuracy in spite of truncations.
We examine a standard scheme to obtain the non-Hermitian Hamiltonian (NHH) from the Lindblad master equation by neglecting its jump term, and propose an alternative approach to address its limitations. The NHH obtained by the conventional scheme fails to provide a good approximation for fermionic many-body systems, even on short timescales. To resolve this issue, we present a framework called the local NHH formalism, which describes the loss process in each individual mode locally. This formalism is applicable to general dissipative fermionic systems and remains consistent with the underlying Lindblad master equation at the level of the equations of motion of local observables. The local NHH formalism also provides a convenient framework for spectral analysis, compared to the Lindblad master equation. As an illustration, we consider a fermionic superfluid subject to one-body loss and find the loss-induced population increase. The conventional NHH fails to capture these unique phenomena.