Abstract The Lindblad equation describes the time evolution of a density matrix of a quantum mechanical system. Stationary solutions are obtained by time-averaging the solution, which will in general depend on the initial state. We provide an analytical expression for the steady states of the Lindblad equation using the quantum jump unraveling, a version of an ergodic theorem, and the stationary probabilities of the corresponding discrete-time Markov chains. Our result is valid when the number of states appearing the in quantum trajectory is finite. The classical case of a Markov jump-process is recovered as a special case, and differences between the two are discussed.
Abstract The Lindblad equation describes the time evolution of a density matrix of a quantum mechanical system. Stationary solutions are obtained by time-averaging the solution, which will in general depend on the initial state. We provide an analytical expression for the steady states of the Lindblad equation using the quantum jump unraveling, a version of an ergodic theorem, and the stationary probabilities of the corresponding discrete-time Markov chains. Our result is valid when the number of states appearing the in quantum trajectory is finite. The classical case of a Markov jump-process is recovered as a special case, and differences between the two are discussed.
Abstract The dynamics of Gaussian states for open quantum systems described by Lindblad equations can be solved analytically for systems with quadratic Hamiltonians and linear Lindbladians, showing the familiar phenomena of dissipation and decoherence. It is well known that the Lindblad dynamics can be expressed as an ensemble average over stochastic pure-state dynamics, which can be interpreted as individual experimental implementations, where the form of the stochastic dynamics depends on the measurement setup. Here we consider quantum-jump and stochastic Schrödinger dynamics for initially Gaussian states. While both unravellings converge to the same Lindblad dynamics when averaged, the individual dynamics can differ qualitatively. For the stochastic Schrödinger equation, Gaussian states remain Gaussian during the evolution, with stochastic differential equations governing the evolution of the phase-space centre and a deterministic evolution of the covariance matrix. In contrast to this, individual pure-state dynamics arising from the quantum-jump evolution do not remain Gaussian in general. Applying results developed in the non-Hermitian context for Hagedorn wavepackets, we formulate a method to generate quantum-jump trajectories that is described entirely in terms of the evolution of an underlying Gaussian state. To illustrate the behaviours of the different unravellings in comparison to the Lindblad dynamics, we consider two examples in detail, which can be largely treated analytically, a harmonic oscillator subject to position measurement and a damped harmonic oscillator. In both cases, we highlight the differences as well as the similarities of the stochastic Schrödinger and the quantum-jump dynamics.
This work is concerned with determination of the steady-state structure of time-independent Lindblad master equations, especially those possessing more than one steady state. The approach here is to treat Lindblad systems as generalizations of unitary quantum mechanics, extending the intuition of symmetries and conserved quantities to the dissipative case. We combine and apply various results to obtain an exhaustive characterization of the infinite-time behavior of Lindblad evolution, including both the structure of the infinite-time density matrix and its dependence on initial conditions. The effect of the environment in the infinite-time limit can therefore be tracked exactly for arbitrary state initialization and without knowledge of dynamics at intermediate time. As a consequence, sufficient criteria for determining the steady state of a Lindblad master equation are obtained. These criteria are knowledge of the initial state, a basis for the steady-state subspace, and all conserved quantities. We give examples of two-qubit dissipation and single-mode d-photon absorption where all quantities are determined analytically. Applications of these techniques to quantum information, computation, and feedback control are discussed.
We introduce jump-time unraveling as a distinct description of open quantum systems. As our starting point, we consider quantum jump trajectories, which emerge, physically, from continuous quantum measurements, or, formally, from the unraveling of Markovian quantum master equations. If the stochastically evolving quantum trajectories are ensemble-averaged at specific times, the resulting quantum states are solutions to the associated quantum master equation. We demonstrate that quantum trajectories can also be ensemble-averaged at specific jump counts. The resulting jump-time-averaged quantum states are then solutions to a discrete, deterministic evolution equation, with time replaced by the jump count. This jump-time evolution represents a trace-preserving quantum dynamical map if and only if the associated open system does not exhibit dark states. In the presence of dark states, on the other hand, jump-time-averaged states may decay into the dark states and the jump-time evolution may eventually terminate. Jump-time-averaged quantum states and the associated jump-time evolution are operationally accessible in continuous measurement schemes, when quantum jumps are detected and used to trigger the readout measurements. We illustrate the jump-time evolution with the examples of a two-level system undergoing relaxation or dephasing, a damped harmonic oscillator, and a free particle exposed to collisional decoherence.
We show that all nonrelativistic quantum processes, whether open or closed, are either unitary or probabilistic unitary, i.e., probabilistic combinations of unitary evolutions. This means that, for open quantum systems, their continuous dynamics can always be described by the Lindblad master equation with all jump operators being unitary. We call this formalism the probabilistic unitary formulation of open quantum system dynamics. This formalism is shown to be exact under all cases, and does not rely on any assumptions other than the continuity and differentiability of the density matrix. Moreover, it requires as few as $d\ensuremath{-}1$ jump operators, instead of ${d}^{2}\ensuremath{-}1$, to describe the open dynamics in the most general case, where d is the dimension of Hilbert space of the system. Importantly, different from the conventional Lindblad master equation, this formalism is state dependent, meaning that the Hamiltonian, jump operators, and rates in general all depend on the current state of the density matrix. Hence one needs to know the explicit expression of the density matrix in order to write down the probabilistic unitary master equation explicitly. Experimentally, the formalism provides a scheme to control a quantum state to evolve along designed nonunitary quantum trajectories, and can be potentially useful in quantum computing and quantum control scenes since only unitary resources are needed for implementation.
The Lindblad equation determines the time evolution of the density operator of open quantum systems. While valid for any system size, its use is, in practice, restricted to prototype/surrogate models with the aim of tackling specific aspects of the overall quantum complexity of a multi-atomic system. Path integral molecular dynamics (PIMD) instead provides static and dynamical quantum statistical averages of physical observables for systems in equilibrium composed of up to thousands of atoms over timescales up to nanoseconds, under the condition that short-time quantum coherence is not relevant for the properties of interest. PIMD relies on the well-established technique of molecular dynamics (MD) with its associated classical trajectories. However, it cannot describe a direct time evolution of a system and its convergence to a stationary state in situations out of equilibrium. In this work, we analyze the link between the Lindblad equation and PIMD; specifically, we will discuss how PIMD can actually be used to calculate the time evolution of ensemble-averaged physical observables and their convergence to a stationary state for situations out of equilibrium, bypassing the need of explicitly solving the Lindblad equation. Yet, at the same time, the Lindblad equation and PIMD are linked to one another through a formal relation of equivalence, which provides an argument for the consistency of PIMD results, namely the positivity of the density operator at any time. A numerical study of a prototype system, which is of interest in chemical physics, will be used to showcase the method.
For the efficient simulation of open quantum systems, we often use quantum jump trajectories given by pure states that evolve stochastically to unravel the dynamics of the underlying master equation. In the Markovian regime, when the dynamics is described by a Gorini-Kossakowski-Sudarshan-Lindblad (GKSL) master equation, this procedure is known as Monte Carlo wave function approach. However, beyond ultraweak system-bath coupling, the dynamics of the system is not described by an equation of GKSL type, but rather by the Redfield equation, which can be brought into pseudo-Lindblad form. Here, negative dissipation strengths prohibit the conventional approach. To overcome this problem, we propose a pseudo-Lindblad quantum trajectory (PLQT) unraveling. It does not require an effective extension of the state space, like other approaches, except for the addition of a single classical bit. We test the PLQT for the eternal non-Markovian master equation for a single qubit and an interacting Fermi-Hubbard chain coupled to a thermal bath and discuss its computational effort compared to solving the full master equation.