Solving non-Markovian open quantum systems with multi-channel reservoir coupling

An open quantum system with multiple levels coupled to a bosonic environment at zero temperature is investigated systematically using the non-Markovian quantum-state-diffusion (QSD) method [W. T. Strunz, L. Di\'osi, and N. Gisin, Phys. Rev. Lett. 82, 1801 (1999)]. We have established exact time-local QSD equations for a set of interesting multilevel open systems, including high-spin systems, multiple-transition atomic models, and multilevel atomic models driven by time-dependent external fields. These exact QSD equations have paved the way to evaluate the dynamics of open multilevel atomic systems in the general non-Markovian regimes without any approximation.

We find dynamical invariants for open quantum systems described by the non-Markovian quantum state diffusion (QSD) equation. In stark contrast to closed systems where the dynamical invariant can be identical to the system density operator, these dynamical invariants no longer share the equation of motion for the density operator. Moreover, the invariants obtained with from bi-orthonormal basis can be used to render an exact solution to the QSD equation and the corresponding non-Markovian dynamics without using master equations or numerical simulations. Significantly we show that we can apply these dynamic invariants to reverse-engineering a Hamiltonian that is capable of driving the system to the target state, providing a novel way to design control strategy for open quantum systems.

Open quantum systems and study of decoherence are important for our fundamental understanding of quantum physical phenomena. For practical purposes, a large number of quantum protocols exist that exploit quantum resources, e.g., entanglement, which allows us to go beyond what is possible to achieve by classical means. We combine concepts from open quantum systems and quantum information science and give a proof-of-principle experimental demonstration---with teleportation---that it is possible to implement efficiently a quantum protocol via a non-Markovian open system. The results show that, at the time of implementation of the protocol, it is not necessary to have the quantum resource in the degree of freedom used for the basic protocol---as long as there exists some other degree of freedom or the environment of an open system, which contains useful resources. The experiment is based on a pair of photons, where their polarizations act as open system qubits and frequencies as their environments, while the path degree of freedom of one of the photons represents the state of Alice's qubit to be teleported to Bob's polarization qubit.

Non-Markovian open quantum systems represent the most general dynamics when the quantum system is coupled with a bath environment. The quantum dynamics arising from many important applications are non-Markovian. Although for special cases, such as Hamiltonian evolution and Lindblad evolution, quantum simulation algorithms have been extensively studied, efficient quantum simulations for the dynamics of non-Markovian open quantum systems remain underexplored. The most immediate obstacle for studying such systems is the lack of a universal succinct description of their dynamics. In this work, we fulfill the gap of studying such dynamics by (1) providing a succinct representation of the dynamics of non-Markovian open quantum systems with quantifiable error, and (2) developing an efficient quantum algorithm for simulating such dynamics with cost $${\mathcal {O}}(t\, \textrm{polylog}(t/\epsilon ))$$ for evolution time t and precision $$\epsilon $$ . Our derivation of the succinct representation is based on stochastic Schrödinger equations, which could lead to new alternatives to deal with open quantum systems as well.

We investigate the feasibility of using an effective pulse control method to realize an almost exact state transmission through a uniform chain where the communication channel is immersed in a quantum environment. This pulse control is treated in a nonperturbative fashion, allowing for a more precise view. Using the non-Markovian quantum state diffusion equation, we calculate the dynamics of the transmission fidelity. The results show that the memory effects of the environment can support an effective transmission control that can be boosted under a non-Markovian environment. Our scheme provides another way to realize almost exact state transmission in a spin chain system.

We introduce simple diagrams and rules to calculate the exact time-local master equation for an arbitrary open quantum system interacting with a bath. In it, the time-convolutionless generator is expanded in the system-bath coupling to all orders, and the $n\mathrm{th}$-order contribution is expressed by the superoperators generated from the system operator in the system-bath coupling and the associated $n$-point nonequilibrium correlation functions of the environment operators. This theory generalizes previous attempts to obtain the time-convolutionless generator and provides a general framework to study the non-Markovian, nonequilibrium dynamics of open quantum systems.

We describe a general-purpose computational toolkit for simulating open quantum systems, which provides numerically exact solutions for composites of zero-dimensional quantum systems that may be strongly coupled to multiple, quite general non-Markovian environments. It is based on process tensor matrix product operators (PT-MPOs), which efficiently encapsulate environment influences. The code features implementations of several PT-MPO algorithms, in particular Automated Compression of Environments for general environments comprised of independent modes as well as schemes for generalized spin boson models. The latter includes a divide-and-conquer scheme for periodic PT-MPOs, which enable million time step simulations for realistic models. PT-MPOs can be precalculated and reused for efficiently probing different time-dependent system Hamiltonians. They can also be stacked together and combined to provide numerically complete solutions of small networks of open quantum systems. The code is written in C++ and is fully controllable by configuration files, for which we have developed a versatile and compact human-readable format.

Non-Markovian dynamics arising from the strong coupling of a system to a structured environment is essential in many applications of quantum mechanics and emerging technologies. Deriving an accurate description of general quantum dynamics including memory effects is, however, a demanding task, prohibitive to standard analytical or direct numerical approaches. We present a major release of our open source software package, OQuPy (Open Quantum System in Python), which provides several recently developed numerical methods that address this challenging task. It utilizes the process tensor approach to open quantum systems (OQS) in which a single map, the process tensor, captures all possible effects of an environment on the system. The representation of the process tensor in a tensor network form allows for an exact yet highly efficient description of non-Markovian OQS (NM-OQS). The OQuPy package provides methods to (1) compute the dynamics and multi-time correlations of quantum systems coupled to single and multiple environments, (2) optimize control protocols for NM-OQS, (3) simulate interacting chains of NM-OQS, and (4) compute the mean-field dynamics of an ensemble of NM-OQS coupled to a common central system. Our aim is to provide an easily accessible and extensible tool for researchers of OQS in fields such as quantum chemistry, quantum sensing, and quantum information.

We critically evaluate the most widespread assumption in the theoretical description of coherent control strategies for open quantum systems. We show that, for non-Markovian open systems dynamics, this fixed-dissipator assumption leads to a serious pitfall causing insurmountable difficulties in the effective modelling of the controlled system. We show that the only way to avoid these problems is a full microscopic description of the controlled system in presence of noise. We illustrate our findings with a paradigmatic example.

Using Liouville space and superoperator formalism we consider pure stationary states of open and dissipative quantum systems. We discuss stationary states of open quantum systems, which coincide with stationary states of closed quantum systems. Open quantum systems with pure stationary states of linear oscillator are suggested. We consider stationary states for the Lindblad equation. We discuss bifurcations of pure stationary states for open quantum systems which are quantum analogs of classical dynamical bifurcations.

The fully quantized model of double qubits coupled to a common bath is solved using the quantum state diffusion (QSD) approach in the non-Markovian regime. We have established the explicit time-local non-Markovian QSD equations for the two-qubit dissipative and dephasing models. Diffusive quantum trajectories are applied to the entanglement estimation of two-qubit systems in a non-Markovian regime. In both cases, non-Markovian features of entanglement evolution are revealed through quantum diffusive unravellings in the system state space.

Open quantum systems (OQS) cannot always be described with the Markov approximation, which requires a large separation of system and environment time scales. Here, we give an overview of some of the most important techniques available to tackle the dynamics of an OQS beyond the Markov approximation. Some of these techniques, such as master equations, Heisenberg equations and stochastic methods, are based on solving the reduced OQS dynamics, while others, such as path integral Monte Carlo or chain mapping approaches, are based on solving the dynamics of the full system. We emphasize the physical interpretation and derivation of the various approaches, explore how they are connected and examine how different methods may be suitable for solving different problems.

Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. The non-Markovian environment is often approximated by discrete bosonic modes, thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been previously proposed, the resulting approximation lacks rigorous and general convergence guarantees. In this Letter, we show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that this approximation error typically falls off polynomially with the number of modes. Our results lend rigor to classical and quantum algorithms for approximating non-Markovian dynamics.

The weak-coupled two-level open quantum system described by non-Markovian Time-convolution-less master equation is investigated in this paper. The cut-off frequency, coupling constant and transition frequency, which impact on the system's decay rate, coherence factor and purity, are investigated. The appropriate parameters used in system simulation experiments are determined by comparing analysis results of different values of parameters for the effects of system performance. The control laws used to transfer the system states are designed on the basis of the Lyapunov stability theorem. Numerical simulation experiments are implemented under the MATLAB environment. The features of the free evolution trajectory of the non-Markovian systems and the states transfer from a pure state to a desired pure state under the action of the proposed control laws are studied, respectively. By comparing the experimental results, the effectiveness of the proposed quantum Lyapunov control method applied to the state transfer in non-Markovian open quantum systems is verified. Meanwhile, the influences of different control parameters and cut-off frequencies on the system performance are analyzed.

We investigate the preparation of a target initial state for a two-level (qubit) system from a system-environment equilibrium or correlated state by an external field. The system-environment equilibrium or correlated state results from the inevitable interaction of the system with its environment. An efficient method in an extended auxiliary Liouville space is introduced to describe the dynamics of the non-Markovian open quantum system in the presence of a strong field and an initial system-environment correlation. By using the time evolutions of the population difference, the state trajectory in the Bloch sphere representation and the trace distance between two reduced system states of the open quantum system, the effect of initial system-environment correlations on the preparation of a system state is studied. We introduce an upper bound and a lower bound for the trace distance within our perturbation formalism to describe the diverse behaviors of the dynamics of the trace distance between various correlated states after the system state preparation. These bounds that are much more computable than similar bounds in the literature give a sufficient condition and a necessary condition for the increase of the trace distance and are related to the witnesses of non-Markovianity and initial system-bath correlation.