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.
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.
Stochastic unravelings provide a useful way to represent open quantum system dynamics in terms of pure state realizations, and have been widely studied both from a fundamental and from a computational point of view. They were initially formulated for Markovian dynamics described by the Gorini-Kossakowski-Sudarshan-Lindblad master equation. However, due to recent technological and experimental development, most physical relevant dynamics present temporal correlations beyond the Markov approximation. Such correlations cause decay rates to turn temporarily negative, thus requiring the generalization of stochastic unravelings from Markovian to non-Markovian scenarios. Indeed, many unraveling techniques have been introduced in this regime, and a comprehensive review of the different jump methods is currently missing. In this work, we provide an overview of widely used quantum jump unraveling techniques for non-Markovian systems and also discuss them in terms of their numerical efficiency, divisibility requirements, Hilbert space extension, and measurement interpretation.
Embedding non-Markovian open quantum dynamics into an enlarged Markovian space offers a powerful route to nonperturbative simulations, where the dynamics of the extended space can be governed by multiple distinct Markovian equations. We show that these distinct embeddings arise from different unravelings of Gaussian bath self-energies, generating a family of deterministic, time-local equations for the extended system. Using the Brownian-oscillator spectral density as an illustrative example, we clarify the relationships among existing approaches, including the Hierarchical Equations of Motion (HEOM) and the Lindblad--pseudomode formalism, and demonstrate how this framework enables numerically stable and efficient simulations. This work provides both a transparent theoretical foundation for embedding techniques and a flexible platform for developing new methods to simulate non-Markovian quantum dynamics.
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.
The pseudomode framework provides an exact description of the dynamics of an open quantum system coupled to a non-Markovian environment. Using this framework, the influence of the environment on the system is studied in an equivalent model, where the open system is coupled to a finite number of unphysical pseudomodes that follow a time-local master equation. Building on the insight that this master equation does not need to conserve the hermiticity of the pseudomode state, we here ask for the most general conditions on the master equation that guarantee the correct reproduction of the system's original dynamics. We demonstrate that our generalized approach decreases the number of pseudomodes that are required to model, for example, underdamped environments at finite temperature. We also provide an unraveling of the master equation into quantum jump trajectories of non-Hermitian states, which further facilitates the utilization of the pseudomode technique for numerical calculations by enabling the use of easily parallelizable Monte Carlo simulations. Finally, we show that pseudomodes, despite their unphysical nature, provide a natural picture in which physical processes, such as the creation of system-bath correlations or the exchange of heat, can be studied. Hence, our results pave the way for future investigations of the system-environment interaction leading to a better understanding of open quantum systems far from the Markovian weak-coupling limit. Published by the American Physical Society 2024
The paradigm of considering open quantum systems—i.e., focusing only on the system of interest, and treating the rest of the world as an effective environment—has proven to be a highly effective way to understand a range of quantum systems, across areas of study such as quantum optics, cold atoms, superconducting qubits, and impurities in solid-state systems. A common approach in many of these contexts has been to consider simplified approaches based on the Born and Markov approximations. While these approximations are indeed often appropriate in contexts such as quantum optics, the widespread application of these approximations has been driven more by simplicity than by accuracy. In particular, these Markovian treatments will fail in many cases, such as when coupling to the environment is not weak, when the environment is structured and has resonances, when the system couples to low-frequency modes of the environment, or when the questions of interest involve the propagation of information through the environment. Despite the fact that many real problems are non-Markovian, the Markov approximation is still widely used, as it is often assumed that a fully non-Markovian treatment is too complex to be practical. In this Perspective we discuss a recently developed set of techniques that address this challenge. Centering our discussion around the notion of the process tensor, we demonstrate that the generality of the process tensor concept, coupled with efficient tensor-network methods, opens the door to the description of a wide range of observable non-Markovian processes in a wide range of open quantum systems.
In this thesis, we discuss aspects of the simulation, characterisation, and control of the dynamics of open quantum systems. The latter are systems governed by quantum mechanical laws and interacting with an environment, which is typically much larger than the system itself. According to the standard approach, one is usually able to perform a suitable average over the environmental degrees of freedom, resulting in an effective description of the main system dynamics, which accounts for the effects of the interaction with the surroundings. In standard scenarios, usually studied within the so-called Born-Markov approximation, the system-environment coupling is weak and such that we can perform a neat separation of timescales: the environmental dynamics is assumed to be fast compared to the typical evolution timescale of the system of interest. However, nowadays we are able to inspect physical scenarios where the usual Born-Markov approximation breaks down: the interaction between system and environment can be strong, likely leading to non-negligible memory (non-Markovian) effects. In this work, we discuss some instances in which we cannot work in the standard Born-Markov regime. We introduce and analyse some numerical and analytical techniques to characterise and simulate non-standard scenarios. We show that, even in the Markovian regime, the traditional formulation of thermodynamic irreversibility shows flaws and inconsistencies. We thus use phase-space methods to assess the role of initial correlations shared by the two parties of a bipartite harmonic system in the entropy production rate, including non-Markovian effects. Furthermore, we show that certain interactions enable to redraw the boundaries between the system and the environment in an effective manner, resulting in a new picture where the system degrees of freedom are augmented, while the residual environment is rearranged in such a way that the Born-Markov approximation is recovered. This analytical technique, known as reaction coordinate mapping approach, is employed in this work to show that, upon an accurate choice of the parameter regime, a spin-boson model can serve as a quantum analogue simulator of non-Markovian multiphoton Jaynes-Cummings models. These systematic studies shed light on the thermodynamic characterisation of open quantum systems, as well as on their numerical simulation.