We present a perturbation theory for non-Markovian quantum-state diffusion (QSD), the theory of diffusive quantum trajectories for open systems in a bosonic environment [Phys. Rev. A 58, 1699 (1998)]. We establish a systematic expansion in the ratio between the environmental correlation time and the typical system time scale. The leading order recovers the Markov theory, so here we concentrate on the next-order correction corresponding to first-order non-Markovian master equations. These perturbative equations greatly simplify the general non-Markovian QSD approach, and allow for efficient numerical simulations beyond the Markov approximation. Furthermore, we show that each perturbative scheme for QSD naturally gives rise to a perturbative scheme for the master equation which we study in some detail. Analytical and numerical examples are presented, including the quantum Brownian motion model.
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 tutorial, I will review some recent progresses in non-Markovian dynamics of quantum open systems. I will be focused with non-Markovian quantum trajectories and non-Markovian master equation approaches. In particular, I will show how to derive an exact non-Markovian master equation from the corresponding quantum trajectories. I will also show how to use quantum trajectories to derive the well-known Kraus operators for pure dephasing noise. The applications of quantum trajectories to quantum information and decoherence will be briefly reviewed.
The description of open quantum system dynamics beyond the perturbative treatment (usually associated with Markovian master equations) is a computationally challenging task due to the unfavorable exponential scaling of memory kernels. Developed over recent decades in the context of quantum information and condensed matter, tensor networks provide both a new formalism and a toolbox for overcoming previous computational bottlenecks. This framework enables the formulation of non-perturbative, numerically exact methods for describing the dynamics of open quantum systems to controllable numerical accuracy. In this review, we present these methods and discuss their commonalities and differences to paint a comprehensive view of the field.
An open quantum system is a quantum system that interacts with some environment whose degrees of freedom have been coarse grained away. This model describes non-equilibrium processes more general than scattering-matrix formulations. Furthermore, the microscopically-derived environment provides a model of noise, dissipation and decoherence far more general than Markovian (white noise) models. The latter are fully characterized by Lindblad equations and can be motivated phenomenologically. Non-Markovian processes consistently account for backreaction with the environment and can incorporate effects such as finite temperature and spatial correlations. We consider linear systems with bilinear coupling to the environment, or quantum Brownian motion, and nonlinear systems with weak coupling to the environment. For linear systems we provide exact solutions with analytical results for a variety of spectral densities. Furthermore, we point out an important mathematical subtlety which led to incorrect master-equation coefficients in earlier derivations, given nonlocal dissipation. For nonlinear systems we provide perturbative solutions by translating the formalism of canonical perturbation theory into the context of master equations. It is shown that unavoidable degeneracy causes an unfortunate reduction in accuracy between perturbative master equations and their solutions. We also extend the famous theorem of Lindblad, Gorini, Kossakowski and Sudarshan on completely positivity to non-Markovian master equations. Our application is primarily to model atoms interacting via a common electromagnetic field. The electromagnetic field contains correlations in both space and time, which are related to its relativistic (photon-mediated) nature. As such, atoms residing in the same field experience different environmental effects depending upon their relative position and orientation. Our more accurate solutions were necessary to assess sudden death of entanglement at zero temperature. In contrast to previous claims, we found that all initial states of two-level atoms undergo finite-time disentanglement. We were also able to access regimes which cannot be described by Lindblad equations and other simpler methods, such as near resonance. Finally we revisit the infamous Abraham-Lorentz force, wherein a single particle in motion experiences backreaction from the electromagnetic field. This leads to a number of well-known problems including pre-acceleration and runaway solutions. We found a more a more-suitable open-system treatment of the nonrelativistic particle to be perfectly causal and dissipative without any extraneous requirements for finite size of the particle, weak coupling to the field, etc..