In order to model realistic quantum devices it is necessary to simulate quantum systems strongly coupled to their environment. To date, most understanding of open quantum systems is restricted either to weak system–bath couplings or to special cases where specific numerical techniques become effective. Here we present a general and yet exact numerical approach that efficiently describes the time evolution of a quantum system coupled to a non-Markovian harmonic environment. Our method relies on expressing the system state and its propagator as a matrix product state and operator, respectively, and using a singular value decomposition to compress the description of the state as time evolves. We demonstrate the power and flexibility of our approach by numerically identifying the localisation transition of the Ohmic spin-boson model, and considering a model with widely separated environmental timescales arising for a pair of spins embedded in a common environment.

Revealing memory effects in phasecovariant quantum master equations J Teittinen, H Lyyra, B Sokolov et al. AbstractStrong coupling between a system and its environment leads to the emergence of non-Markovian dynamics, which cannot be described by a time-local master equation. One way to capture such dynamics is to use numerical real-time path integrals, where assuming a finite bath memory time enables manageable simulation scaling. However, by comparing to the exactly soluble independent boson model, we show that the presence of transient negative decay rates in the exact dynamics can result in simulations with unphysical exponential growth of density matrix elements when the finite memory approximation is used. We therefore reformulate this approximation in such a way that the exact dynamics are reproduced identically and then apply our new method to the spin-boson model with superohmic environmental coupling, commonly used to model phonon environments, but which cannot be solved exactly. Our new method allows us to easily access parameter regimes where we find revivals in population dynamics which are due to non-Markovian backflow of information from the bath to the system.

No abstract

The dynamics of a wide range of technologically important quantum systems are dominated by their interaction with just a few environmental modes. Such highly structured environments give rise to long-lived bath correlations that induce complex dynamics which are very difficult to simulate. These difficulties are further aggravated when spatial correlations between different parts of the system are important. By modeling the dynamics of a pair of two-level quantum systems in a common, structured, environment we show that a recently developed numerical approach, the timeevolving matrix product operator, is capable of accurate simulation under exactly these conditions. We find that tuning the separation to match the wavelength of the dominant environmental modes can drastically modify the system dynamics. To further explore this behavior, we show that the full dynamics of the bath can be calculated directly from those of the system, thus allowing us to develop intuition for the complex system dynamics observed.

No abstract

No abstract

Tracing out the environmental degrees of freedom is a necessary procedure when simulating open quantum systems. While being an essential step in deriving a tractable master equation it represents a loss of information. In situations where there is strong interplay between the system and environmental degrees of freedom this loss makes understanding the dynamics challenging. These dynamics, when viewed in isolation, have no time-local description: they are non-Markovian and memory effects induce complex features that are difficult to interpret. To address this problem, we here show how to use system correlations, calculated by any method, to infer any correlation function of a Gaussian environment, so long as the coupling between system and environment is linear. This not only allows reconstruction of the full dynamics of both system and environment, but also opens avenues into studying the effect of a system on its environment. In order to obtain accurate bath dynamics, we exploit a numerically exact approach to simulating the system dynamics, which is based on the construction and contraction of a tensor network that represents the process tensor of this open quantum system. Using this we are able to find any system correlation function exactly. To demonstrate the applicability of our method we show how heat moves between different modes of a bosonic bath when coupled to a two-level system that is subject to an off-resonant drive.

Tracing out the environmental degrees of freedom is a necessary procedure when simulating open quantum systems. While being an essential step in deriving a tractable master equation it represents a loss of information. In situations where there is strong interplay between the system and environmental degrees of freedom this loss makes understanding the system's dynamics challenging. These dynamics, when viewed in isolation, are non-Markovian and memory effects induce complex features that are difficult to interpret. Here we exploit a numerically exact approach to simulating the system dynamics of the spin-boson model, which is based on the construction and contraction of tensor network that represents the process tensor of this open quantum system. We are then able to find any system correlation function exactly. We show that we can use these to infer any correlation function of a Gaussian environment, so long as the coupling between system and envinroment is linear. This not only allows reconstruction of the full dynamics of both system and environment, but also opens avenues into studying the effect of a system on its environment. To demonstrate the applicability of our method we show how heat moves between different modes of a bosonic bath when coupled to a two-level system that is subject to an off-resonant drive.