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.
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.
The nonequilibrium effects of dissipation and drive play a key role in an immense variety of nanoscale and mesoscale quantum systems. To understand the behaviour of open quantum systems, we need accurate methods that capture the influence of the environment on the system, while managing the exponentially large Hilbert space required to describe the system. Tensor network algorithms offer an efficient way to approach this challenge. In this thesis, we develop and apply tensor network techniques to study the dynamics and steady states of various open quantum systems. The first part of the thesis focuses on the driven dissipative many body physics in coupled cavity arrays described by Born-Markov master equations. We extend transfer matrix product operator methods to Liouvillian dynamics, and utilize them to compute dynamical correlation functions and fluorescence spectrum of an infinite coupled cavity array in 1D. We also investigate thermalization, and observe the emergence of a quasi-thermal steady state with a negative effective temperature. In another study, we use infinite projected entangled pair state (iPEPS) methods to compute steady states of coupled cavity lattices in 2D. We find that a straightforward adaptation of iPEPS to Liouvillian dynamics is unstable, contradicting a recent publication in the field. The second part investigates more general systems involving strong couplings and structured environments that induce non-Markovian dynamics. We develop a powerful time-evolving matrix product operator (TEMPO) algorithm that builds on Feynman-Vernon influence functional formalism, and uses matrix product states (MPS) to represent the temporal non-Markovian correlations efficiently. We apply TEMPO to study the localization phase transition of the spin-boson model and the dynamics of two spatially separated two-level systems coupled to a common environment. Finally, we propose the Toblerone TEMPO algorithm, which extends TEMPO to many-body systems interacting with general bosonic environments.
We introduce a general numerical method to compute the dynamics and multitime correlations of chains of quantum systems, where each system may couple strongly to a structured environment. The method combines the process tensor formalism for general (possibly non-Markovian) open quantum systems with time-evolving block decimation for one-dimensional chains. It systematically reduces the numerical complexity originating from system-environment correlations before integrating them into the full many-body problem, making a wide range of applications numerically feasible. We illustrate the power of this method by studying two examples. First, we study the thermalization of individual spins of a short XYZ Heisenberg chain with strongly coupled thermal leads. Our results confirm the complete thermalization of the chain when coupled to a single bath, and they reveal distinct effective temperatures in low-, mid-, and high-frequency regimes when the chain is placed between a hot and a cold bath. Second, we study the dynamics of diffusion in a longer XY chain, when each site couples to its own bath.
Being able to study the dynamics of quantum systems interacting with several environments is important in many settings ranging from quantum chemistry to quantum thermodynamics, through out-of-equilibrium systems. For such problems tensor network-based methods are state-of-the-art approaches for performing numerically exact simulations. However, to be used efficiently in this multi-environment and non-perturbative context, these methods require an optimized choice of the topology of the wave-function Ansätze. This is often done by analysing cross-correlations between different system and environment degrees of freedom. Here, we show for canonical model Hamiltonians that simple orderings of bosonic environmental modes, which enable the joint {System + Environments} state to be written as a matrix product state, considerably reduce the bond dimension required for convergence despite introducing long-ranged interactions. These results suggest that complex correlation analyses for tweaking tensor networks topology (e.g. entanglement renormalization) are usually not necessary and that tree tensor network states are sub-optimal compared to simple matrix product states in many important models of physical open systems.
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.
A better understanding of dissipation is crucial for understanding real-world quantum systems. Indeed, all quantum systems experience interactions with an (often) uncontrollable outside environment that can lead to a decay of excited state populations and a loss of quantum coherences. The study of dissipation is timely as the development of next-generation nanoscale quantum technologies is on its way, and the existence of non-trivial quantum effects in biological systems is being seriously investigated. However, descriptions of dissipation in quantum systems are reduced (most of the time) to time-local approaches and (everywhere) to space-local independent environments. These simplifying assumptions do render analytic and numerical calculations possible, yet they get rid of a breadth of physical processes that can alter radically the quantum systems' dynamics. In this thesis, building on a numerically exact tensor networks method, we developed a technique able to handle spatio-temporal correlations between a quantum system and bosonic (i.e. vibrational, electromagnetic, magnons, etc.) environments. With this method we studied the signalling process - a form of information backflow - in quantum systems, and uncovered how it can induce non-trivial dynamics, and be leveraged to populate otherwise inaccessible excited states. We also evidenced the ability of 'non-local' environment reorganisation, induced by system-environment interactions, to radically change the nature of the thermodynamically favoured system ground state. The new phenomenology of physical processes, resulting from considering quantum systems interacting with a common environment, has important consequences for the design of nanodevices as it gives access to new control, sensing and cross-talk mechanisms. In another vein, these results might also give us a new framework to study and interpret (quantum?) effects in the biological realm.