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 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.
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.
In the present Thesis we study the behavior of multi-time correlation functions and of thermodynamical quantities such as heat in open quantum systems undergoing an evolution generally affected by the presence of memory effects, i.e. non-Markovian. In the last decade, a large part of the scientific community in this field has dedicated its efforts to the understanding, precise definition and quantification of non-Markovianity in the quantum realm and now we have at our disposal several benchmark results and a plethora of different estimators that allow to determine the degree of non-Markovianity of a given dynamics. It comes therefore natural to investigate how other different dynamical quantities relate to such estimators also in order to understand the physical implications of memory effects on the statistics of observable quantities. In the first part of this work, a quantitative test of the violation of the so-called quantum regression theorem in presence of a non-Markovian dynamical regime is investigated. The quantum regression theorem represents a procedure that, whenever valid, allows to reconstruct two-time correlation functions of system's operators from the sole knowledge of the dynamics of mean values. It is worth stressing that two-time correlation functions are necessary in order to fully characterize the statistical properties of a quantum system, since they are able to catch aspects of the dynamics, such as fluorescence spectrum, in general not accessible looking at mean values. Despite their relevance however, obtaining two-time correlation functions often represents a formidable task, since the knowledge of the full "system+environment" dynamics is required, a generally too demanding request in the context of open quantum systems theory. The quantum regression theorem represents in this regard the easiest route to determine two-time correlation functions, this highlighting its importance. In this work we show that, in a pure-dephasing spin-boson model, the quantum regression theorem represents a stronger condition than non-Markovianity, in the sense that any presence of memory effects in the reduced dynamics inevitably results in violations to the former. These results have been published in [G.Guarnieri, A. Smirne, B. Vacchini, Phys. Rev. A 90, 022110 (2014)]. The second part of the Thesis is devoted to the characterization of heat ow at the microscopic level in open quantum systems, both finite and infinite dimensional. In particular we begin by studying the time behavior of its mean value in a non-Markovian dynamical regime, showing that, at variance with what happens in the Born-Markov semigroup limiting case, heat can backflow from the environment to the system. After providing a condition for the occurrence of such phenomenon and a measure for its amount for a given dynamics, the relationship with suitable non-Markovianity estimators is sought in two paradigmatic models, namely the spin-boson and the quantum brownian motion. The results, collected in [ G. Guarnieri, C. Uchiyama, B. Vacchini, Phys. Rev. A 93, 012118 (2016); G. Guarnieri, J. Nokkala, R. Schmidt, S. Maniscalco, B. Vacchini, Phys. Rev. A 94, 062101 (2016)], on the one hand allow for the identification of parameter-regions where the heat backflow is absent or maximum. On the other hand they show that the occurrence of heat backflow represents a stricter condition than non-Markovianity, in the sense that non-Markovianity allows for the observation of heat flowing back from the environment to the system and, vice versa, a Markovian dynamics prevents its occurrence. This Thesis concludes with the formulation of a new family of lower bounds to the mean dissipated heat in an environmental-assisted erasure-protocol scenario where Landauer's principle applies. As originally conceived for classical systems, this principle states that every irreversible erasure of information stored in a system inevitably carries along an amount of heat dissipated into the environment which is expended to perform the action. Within the framework recently put forward in [D. Reeb, M. M. Wolf, New J. Phys. 16, 103011 (2014)], which guarantees the validity of Landauer's principle in an open quantum systems scenario, we provide an asymptotically tight family of lower bounds to the dissipated heat which are also valid in the non-equilibrium setting. This construction is applied to an open system consisting of a three-level V-system, in which one transition is externally pumped by a laser field while the other is coupled through an XX-interaction to an environment consisting of a spin chain. Beside calculating all these quantities, an exact solution for the dynamics of such system is also provided. These results are collected in [G. Guarnieri, S. Campbell, J. Goold, S. Pigeon, M. Paternostro, B. Vacchini, in preparation].
The central challenge for describing the dynamics in open quantum systems is that the Hilbert space of typical environments is too large to be treated exactly. In some cases, such as when the environment has a short memory time or only interacts weakly with the system, approximate descriptions of the system are possible. Beyond these, numerically exact methods exist, but these are typically restricted to baths with Gaussian correlations, such as non-interacting bosons. Here we present a numerically-exact method for simulating open quantum systems with arbitrary environments. Our approach automatically reduces the large number of environmental degrees of freedom to those which are most relevant. Specifically, we show how the process tensor -- which describes the effect of the environment -- can be iteratively constructed and compressed using matrix product state techniques. We demonstrate the power of this method by applying it to problems with bosonic, fermionic, and spin environments: electron transport, phonon effects and radiative decay in quantum dots, central spin dynamics, and superradiance. The versatility and efficiency of our automated compression of environments (ACE) method provides a practical general-purpose tool for open quantum systems.
This thesis investigates two thematic lines of research, both underpinned by non-Markovian system-reservoir interactions in quantum optics. The overarching focus is on modelling the open system dynamics in a non-perturbative fashion, broadly on - though not restricted to - instances when the environment is structured. A theory is developed by means of enlarging the open system over environmental degrees of freedom to include memory effects in its dynamics. This is achieved using an established technique that involves mapping a bosonic environment onto a 1D chain of harmonic oscillators. Within this setting, we apply a Heisenberg equation-of-motion approach to derive an exact set coupled differential equations for the open system and a single auxiliary oscillator of the chain. The combined equations are shown to have their interpretation rooted in a quantum Markov stochastic process. Including the auxiliary chain oscillator as part of the original system then enables us to obtain an exact master equation for the enlarged system, avoiding any need for the Born-Markov approximations. Our method is valid for a dissipative two-state system, with cases of multiple excitations and added driving discussed. Separately, we apply the framework of quantum Darwinism to an atom-cavity system, and, subsequently, to a more general multiple-environment model. In both cases, the time-dependent spread of correlations between the open system and fractions of the environment is analysed during the course of the decoherence process. The degree to which information is redundant across different fractions is checked to infer the emergence of classicality. In the second case, we go further and present a decomposition of information in terms of its quantum and classical correlations. A quantitative measure of redundancy is also studied with regard to its ability to witness non-Markovian behaviour. Besides fundamental interest, our results have application to quantum information processing and quantum technologies, keeping in mind the potential beneficial use of non-Markovian effects in reservoir engineering.
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.
Recent advances in quantum technologies and related experiments have created a need for highly accurate, versatile, and computationally efficient simulation techniques for the dynamics of open quantum systems. Long-lived correlation effects (non-Markovianity), system-environment hybridization, and the necessity for accuracy beyond the Born-Markov approximation form particular challenges. Approaches to meet these challenges have been introduced, originating from different fields, such as hierarchical equations of motion, Lindblad-pseudomode formulas, chain-mapping approaches, quantum Brownian motion master equations, stochastic unravelings, and refined quantum master equations. This diversity, while indicative of the field's relevance, has inadvertently led to a fragmentation that hinders cohesive advances and their effective cross-community application to current problems for complex systems. How are different approaches related to each other? What are their strengths and limitations? Here we give a systematic overview and concise discussion addressing these questions. We make use of a unified framework which very conveniently allows to link different schemes and, this way, may also catalyze further progress. In line with the state of the art, this framework is formulated not in a fully reduced space of the system but in an extended state space which in a minimal fashion includes effective reservoir modes. This in turn offers a comprehensive understanding of existing methods, elucidating their physical interpretations, interconnections, and applicability.