Colloquium: Non-Markovian dynamics in open quantum systems

Relaxation and decoherence processes are key features of the dynamics of open quantum systems [1]. In the standard approach one tries to develop appropriate master equations for the open system's reduced density matrix ρS which is given by the partial trace taken over the environmental variables coupled to the open system. Invoking the weak-coupling assumption one can formulate in many cases of physical interest a Markovian quantum master equation for the reduced density matrix, expressing the dynamical laws for the irreversible motion of the open system.

A new approach to theory and simulation of the non-Markovian dynamics of open quantum systems is presented. It is based on identification of a parameter which is uniformly bounded on wide time intervals: the occupation of the virtual cloud of quanta. By ‘virtual’ we denote those bath excitations which were emitted by the open system, but eventually will be reabsorbed before any measurement of the bath state. A useful property of the virtual cloud is that the number of its quanta is expected to saturate on long times, since physically this cloud is a (retarded) polarization of the bath around the system. Therefore, the joint state of open system and virtual cloud (we call it dressed state) can be accurately represented in a truncated basis of Fock states, on a wide time scale. At the same time, there can be an arbitrarily large number of the observable quanta (which survive up to measurement), especially if the open system is under driving. However, it turns out that the statistics of the bath-measurement outcomes is classical (in a suitable measurement basis): one can employ a Monte Carlo sampling of these outcomes. Therefore, it is possible to efficiently simulate the dynamics of the observable quantum field. In this work we consider the bath measurement with respect to the coherent states, which yields the Husimi function as the positive (quasi)probability distribution of the outcomes. The joint evolution of the dressed state and the corresponding outcome is called the dressed quantum trajectory. The Monte Carlo sampling of these trajectories yields a stochastic simulation method with promising convergence properties on wide time scales.

A systematic approach to the non-Markovian quantum dynamics of open systems is given by the projection operator techniques of nonequilibrium statistical mechanics. Combining these methods with concepts from quantum information theory and from the theory of positive maps, we derive a class of correlated projection superoperators that take into account in an efficient way statistical correlations between the open system and its environment. The result is used to develop a generalization of the Lindblad theory to the regime of highly non-Markovian quantum processes in structured environments.

We define quantum chaos and integrability in open quantum many-body systems as a dynamical property of single stochastic realizations, referred to as quantum trajectories. This definition relies on the predictions of random matrix theory applied to the subset of the Liouvillian eigenspectrum involved in each quantum trajectory. Our approach, which we name (SSQT), enables a natural distinction between transient and steady-state quantum chaos as general phenomena in open setups. We test the generality and reliability of the SSQT criterion on several dissipative systems, further showing that an open system with a chaotic structure can evolve towards either a chaotic or integrable steady state. We apply our theoretical framework to two driven-dissipative bosonic systems. First, we study the driven-dissipative Bose-Hubbard model, a paradigmatic example of a quantum simulator, clarifying the interplay of integrability, transient, and steady-state chaos across its phase diagram. Our analysis shows the existence of an emergent dissipative quantum chaotic phase, whereas the classical and semiclassical limits display an integrable behavior. In this regime, chaos arises from the quantum and classical fluctuations associated with the dissipation mechanisms. Second, we investigate dissipative quantum chaos in the dispersive readout of a transmon qubit: a measurement technique ubiquitous in superconducting-based quantum hardware. Through the SSQT, we distinguish several regimes where the performance of the measurement instrument can be connected to the integrable or chaotic nature of the underlying driven-dissipative bosonic system. Our work offers a general understanding of the integrable and chaotic dynamics of open quantum systems and paves the way for the investigation of dissipative quantum chaos and its consequences on state-of-the-art noisy intermediate-scale quantum devices.

We discuss two methods of an exact stochastic representation of the non-Markovian quantum dynamics of open systems. The first method employs a pair of stochastic product vectors in the total system's state space, while the second method uses a pair of state vectors in the open system's state space and a random operator acting on the state space of the environment. Both techniques lead to an exact solution of the von Neumann equation for the density matrix of the total system. Employing a spin star model describing a central spin coupled to the bath of surrounding spins, we perform Monte Carlo simulations for both variants of the stochastic dynamics. In addition, we derive an analytical expression for the expectation values of the stochastic dynamics to obtain the exact solution for the density matrix of the central spin.

Information scrambling refers to the phenomenon in which local quantum information in a many-body system becomes dispersed throughout the entire system under unitary evolution. It has been extensively studied in closed quantum systems, where it is quantified by operator size growth, revealing deep connections between condensed matter physics, high-energy physics, and quantum information. However, when extending the study of operator size dynamics to open quantum systems, two different definitions of operator size distributions emerge. These definitions are based on different treatments of the bath. In this work, we aim to establish a unified picture for operator size dynamics in open quantum systems, using the solvable Brownian SYK models at generic system size. In particular, we provide the conditions under which the signature of scrambling transition — discovered using one particular definition — appears in operator size dynamics under the other definition. Additionally, we extend previous studies by exploring finite-size effects that are not captured by the scramblon theory. Our results provide a refined understanding of operator size dynamics in open quantum systems.

Non-Markovian dynamics of time-fractional open quantum systems

We present a general analysis of the role of initial correlations between the open system and an environment on quantum dynamics of the open system.

In the dynamics of open quantum systems, the backflow of information to the reduced system under study has been suggested as the actual physical mechanism inducing memory and thus leading to non-Markovian quantum dynamics. To this aim, the trace-distance revivals between distinct evolved system states have been shown to be subordinated to the establishment of system-environment correlations or changes in the environmental state. We show that this interpretation can be substantiated also for a class of entropic quantifiers. We exploit a suitably regularized version of Umegaki's quantum relative entropy, known as telescopic relative entropy, that is tightly connected to the quantum Jensen-Shannon divergence. In particular, we derive general upper bounds on the telescopic relative entropy revivals conditioned and determined by the formation of correlations and changes in the environment. We illustrate our findings by means of examples, considering in particular the Jaynes-Cummings model and a phase covariant dynamics.

This paper discusses several methods for describing the dynamics of open quantum systems, where the environment of the open system is infinite-dimensional. These are purifications, phase space forms, master equation and liouville equation forms. The main contribution is in using Feynman-Kac formalisms to describe the infinite-demsional components.

Quantum simulation is a powerful tool to study the properties of quantum systems. The dynamics of open quantum systems are often described by completely positive (CP) maps, for which several quantum simulation schemes exist. Such maps, however, represent only a subset of a larger class of maps: the general dynamical maps which are linear, Hermitian preserving, and trace preserving but not necessarily positivity preserving. Here we present a simulation scheme for these general dynamical maps, which occur when the underlying system-reservoir model undergoes entangling (and thus non-Markovian) dynamics. Such maps also arise as the inverse of CP maps, which are commonly used in error mitigation. We illustrate our simulation scheme on an IBM quantum processor, demonstrating its ability to recover the initial state of a Lindblad evolution. This paves the way for a novel form of quantum error mitigation. Our scheme only requires one ancilla qubit as an overhead and a small number of one and two qubit gates. Consequently, we expect it to be of practical use in near-term quantum devices.

We characterize (non)Markovian dynamics of open quantum systems. Two recently proposed measures of non-Markovianity are analyzed: one based on the concept of divisibility of the dynamical map and the other one based on distinguishability of quantum states. The characterization of the correspondingg enerators in the Heisenberg picture is provided as well.

Recently, a "unified" quantum master equationwas derived and shown to be of the Gorini-Kossakowski-Lindblad-Sudarshan form. This equationdescribes the dynamics of open quantum systems in a manner that forgoes the full secular approximation and retains the impact of coherences between eigenstates close in energy. We implement full counting statistics with the unified quantum master equationto investigate the statistics of energy currents through open quantum systems with nearly degenerate levels. We show that, in general, this equationgives rise to dynamics that satisfy fluctuation symmetry, a sufficient condition for the Second Law of Thermodynamics at the level of average fluxes. For systems with nearly degenerate energy levels, such that coherences build up, the unified equation is simultaneously thermodynamically consistent and more accurate than the fully secular master equation. We exemplify our results for a "V" system facilitating energy transport between two thermal baths at different temperatures. We compare the statistics of steady-state heat currents through this system as predicted by the unified equationto those given by the Redfield equation, which is less approximate but, in general, not thermodynamically consistent. We also compare results to the secular equation, where coherences are entirely abandoned. We find that maintaining coherences between nearly degenerate levels is essential to properly capture the current and its cumulants. On the other hand, the relative fluctuations of the heat current, which embody the thermodynamic uncertainty relation, display inconsequential dependence on quantum coherences.

Quantum mechanics must be regarded as open systems. On one hand, this is due to the fact that, like in classical physics, any realistic system is subjected to a coupling to an uncontrollable environment which influences it in a non-negligible way. The theory of open quantum systems thus plays a major role in many applications of quantum physics since perfect isolation of quantum system is not possible and since a complete microscopic description or control of the environment degrees of freedom is not feasible or only partially so [1]. Practical considerations therefore force one to seek for a simpler, effectively probabilistic description in terms of an open system. There is a close physical and mathematical connection between the evolution of an open system, the state changes induced by quantum measurements, and the classical notion of a stochastic process. The paper provides a bibliographic review of this interrelations, it shows the mathematical equivalence between markovian master equation and generalized piecewise deterministic processes [1] and it introduces the open system in an open observed environment model.

We present a renormalization group (RG) method which allows for an analytical study of the transient dynamics of open quantum systems on all time scales. Whereas oscillation frequencies and decay rates of exponential time evolution follow from the fixed point positions, the long-time behavior of pre-exponential functions is related to the scaling behavior around the fixed points. We show that certain terms of the RG flow are only cut off by inverse time, which leads to a difference between infrared and ultraviolet scaling. An evaluation for the ohmic spin boson model at weak damping reveals significant deviations from previous predictions in the long-time regime. We propose that weak coupling problems for stationary quantities can in principle turn into strong coupling ones for the determination of the long-time behavior.