A full understanding of open quantum systems requires the characterization of both system and environmental properties. However, the complexity of the environmental statistics in the presence of strong system-bath hybridization and long memory effects usually prevents effective non-perturbative methods from going beyond the analysis of the reduced system dynamics. Here we present a model consisting of purified auxiliary bosonic modes to describe, alongside properties of the system, the dynamics of environmental observables for bosonic baths prepared in non-Gaussian initial states. We numerically exemplify this method by simulating non-Markovian multi-photon transfer processes on a coupled cavity waveguide system in the large time delay regime.
A full understanding of open quantum systems requires the characterization of both system and environmental properties. However, the complexity of the environmental statistics in the presence of strong system-bath hybridization and long memory effects usually prevents effective non-perturbative methods from going beyond the analysis of the reduced system dynamics. Here we present a model consisting of purified auxiliary bosonic modes to describe, alongside properties of the system, the dynamics of environmental observables for bosonic baths prepared in non-Gaussian initial states. We numerically exemplify this method by simulating non-Markovian multi-photon transfer processes on a coupled cavity waveguide system in the large time delay regime.
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.
We identify the conditions that guarantee equivalence of the reduced dynamics of an open quantum system (OQS) for two different types of environments-one a continuous bosonic environment leading to a unitary system-environment evolution and the other a discrete-mode bosonic environment resulting in a system-mode (nonunitary) Lindbladian evolution. Assuming initial Gaussian states for the environments, we prove that the two OQS dynamics are equivalent if both the expectation values and two-time correlation functions of the environmental interaction operators are the same at all times for the two configurations. Since the numerical and analytical description of a discrete-mode environment undergoing a Lindbladian evolution is significantly more efficient than that of a continuous bosonic environment in a unitary evolution, our result represents a powerful, nonperturbative tool to describe complex and possibly highly non-Markovian dynamics. As a special application, we recover and generalize the well-known pseudomodes approach to open-system dynamics.
We develop an exact framework for describing the non-Markovian dynamics of an open quantum system interacting with an environment modeled by a generalized spectral density function.The approach relies on mapping the initial system onto an auxiliary configuration, comprising the original open system coupled to a small number of discrete modes, which in turn are each coupled to an independent Markovian reservoir.Based on the connection between the discrete modes and the poles of the spectral density function, we show how expanding the system using the discrete modes allows for the full inclusion non-Markovian effects within an enlarged open system whose dynamics is governed by an exact Lindblad master equation.Initially we apply this result to obtain a generalization of the pseudomode method [B.M. Garraway, Phys. Rev. A. 55, 2290 (1997)] in cases where the spectral density function has a Lorentzian structure.For many other types of spectral density function, we extend our proof to show that an open system dynamics may be modeled physically using discrete modes which admit a non-Hermitian coupling to the system and for such cases determine the equivalent master equation to no longer be of Lindblad form.For applications involving two discrete modes, we demonstrate how to convert between pathological and Lindblad forms of the master equation using the techniques of the pseudomode method.
We outline a non-perturbative approach for simulating the behavior of open quantum systems interacting with a bosonic environment defined by a generalized spectral density function. The method is based on replacing the environment by a set of damped harmonic oscillators - the pseudomodes - thereby forming an enlarged open system whose dynamics is governed by a Markovian master equation. Each pseudomode is connected to one of the poles of the spectral density when analytically continued to the lower-half complex frequency plane. Here, we extend previous results to a completely generic class of open system models, and discuss how our framework can be used as a powerful and versatile tool for analyzing non-Markovian open system dynamics. The effectiveness of the method is demonstrated on the spin-boson model by accurately benchmarking its predictions against numerically exact results.
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
Abstract The treatment of quantum thermodynamic systems beyond weak coupling is of increasing relevance, yet extremely challenging. The evaluation of thermodynamic quantities in strong-coupling regimes requires a nonperturbative knowledge of the bath dynamics, which in turn relies on heavy numerical simulations. To tame these difficulties, considering thermal bosonic baths linearly coupled to the open system, we derive expressions for heat, work, and average system-bath interaction energy that only involve the autocorrelation function of the bath and two-time expectation values of system operators. We then exploit the pseudomode approach, which replaces the physical continuous bosonic bath with a small finite number of damped, possibly interacting, modes, to numerically evaluate these relevant thermodynamic quantities. We show in particular that this method allows for an efficient numerical evaluation of thermodynamic quantities in terms of one-time expectation values of the open system and the pseudomodes. We apply this framework to the investigation of two paradigmatic situations. In the first instance, we study the entropy production for a two-level system (TLS) coupled to an ohmic bath, simulated via interacting pseudomodes, allowing for the presence of time-dependent driving. Secondly, we consider a quantum thermal machine composed of a TLS interacting with two thermal baths at different temperatures, showing that an appropriate sinusoidal modulation of the coupling with the cold bath only is enough to obtain work extraction.