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

Nanoscale devices, either biological or artificial, operate in a regime where the usual assumptions of a structureless Markovian bath do not hold. Being able to predict and study the dynamics of such systems is crucial and is usually done by tracing out the bath degrees of freedom, which implies losing information about the environment. To go beyond these approaches we use a numerically exact method relying on a matrix product state representation of the quantum state of a system and its environment to keep track of the bath explicitly. This method is applied to a specific example of interaction that depends on the spatial structure of a system made of two sites. The result is that we predict a non-Markovian dynamics where long-range couplings induce correlations into the environment. The environment dynamics can be naturally extracted from our method and shine a light on long-time feedback effects that are responsible for the observed non-Markovian recurrences in the eigenpopulations of the system.

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.

Understanding the non-equilibrium behavior of quantum systems is a major goal of contemporary physics. Much research is currently focused on the dynamics of many-body systems in lowdimensional lattices following a quench, i.e., a sudden change of parameters. Already such a simple setting poses substantial theoretical challenges for the investigation of the real-time post-quench quantum dynamics. In classical many-body systems the Kolmogorov-Mehl-Johnson-Avrami model describes the phase transformation kinetics of a system that is quenched across a first order phase transition. Here we show that a similar approach can be applied for shedding light on the quench dynamics of an interacting gas of Rydberg atoms, which has become an important experimental platform for the investigation of quantum non-equilibrium effects. We are able to gain an analytic understanding of the time-evolution following a sudden quench from an initial state devoid of Rydberg atoms and identify strikingly different behaviors of the excitation growth in the classical and quantum regimes. Our approach allows us to describe quenches near a non-equilibrium phase transition and provides an approximate analytic solution deep in the quantum domain. arXiv:1709.10383v2 [cond-mat.stat-mech]