© 2026 Aqora Quantum S.A.S.

Terms Privacy Legal Notice

arXiv:1904.04822 · Physical Review A

Optimized auxiliary oscillators for the simulation of general open quantum systems

F. Mascherpa, Andrea Smirne, Alejandro D. Somoza, P. Fernández-Acebal, Sandro Donadi, Dario Tamascelli, Susana F. Huelga, Martin Bodo Plenio

Published May 14, 2020• 98 citations

Abstract

A method for the systematic construction of few-body damped harmonic oscillator networks accurately reproducing the effect of general bosonic environments in open quantum systems is presented. Under the sole assumptions of a Gaussian environment and regardless of the system coupled to it, an algorithm to determine the parameters of an equivalent set of interacting damped oscillators obeying a Markovian quantum master equation is introduced. By choosing a suitable coupling to the system and minimizing an appropriate distance between the two-time correlation function of this effective bath and that of the target environment, the error induced in the reduced dynamics of the system is brought under rigorous control. The interactions among the effective modes provide remarkable flexibility in replicating non-Markovian effects on the system even with a small number of oscillators, and the resulting Lindblad equation for the system and the modes may therefore be integrated at a very reasonable computational cost using standard methods for Markovian problems, even in strongly nonperturbative coupling regimes and at arbitrary temperatures including zero. We apply the method to an exactly solvable problem in order to demonstrate its accuracy, and present two studies based on current research in the context of coherent transport in biological aggregates and organic photovoltaics as more realistic examples of its use and potential; performance and versatility are highlighted, and theoretical and numerical advantages over existing methods, as well as possible future improvements, are discussed.

DOI: 10.1103/physreva.101.052108

Related papers

Pseudomode description of general open quantum system dynamics: non-perturbative master equation for the spin-boson model

Graeme Pleasance, Francesco Petruccione

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.

arXiv:2108.05755•Aug 12, 2021•2 citations

Ancilla-train quantum algorithm for simulating non-Markovian open quantum systems

Hans Michael Christensen, Johannes Agerskov, Frederik Nathan

We present a quantum algorithm for simulating open quantum systems coupled to Gaussian environments valid for any configuration and coupling strength. The algorithm is applicable to problems with strongly coupled, or non-Markovian, environments; problems with multiple environments out of mutual equilibrium; and problems with time-dependent Hamiltonians. We show that the algorithm can reproduce the true dynamics of such problems at arbitrary accuracy and, for a broad range of problems, only adds a minor resource cost relative to Trotterized time evolution; the cost is low-degree polynomial in the inverse target accuracy. The algorithm is based on the insight that any Gaussian environment can be represented as a train of ancillary qubits that sequentially interact with the system through a time-local coupling, given by the convolution square root of the bath correlation function; this is a secondary result of our work. In numerical simulations of a spin-boson model, we find accurate simulations can be achieved with moderate ( 10 1 – 10 3 ) ancilla counts in the intermediate to strong-coupling regimes. Our results open up applications of quantum computers for efficient simulation of nonequilibrium and open quantum systems.

arXiv:2509.12717•Physical Review A•Aug 19, 2026•0 citations

Dynamics and spectra of open quantum systems coupled to anharmonic environments

Kian Damezin, Brendon W. Lovett, Erik M. Gauger +1 more

In open quantum systems, the environment is typically modelled as a collection of harmonic oscillators, whereas realistic environments often exhibit unique non-Markovian effects due to the anharmonicity of the environment. Here we present a framework for modelling anharmonic environments in open quantum systems numerically exactly, in which, via a matching criterion, we map the anharmonic environment onto an effective harmonic one. This results in an effective spectral density that is a collection of shifted bare spectra arising from the non-equidistant energy gaps of the anharmonic environment, along with a unique zero-frequency term that effectively acts as static disorder on the system. This is due to the non-zero variance of the diagonal of the bath coupling operator. In tandem, we develop a framework for investigating the effects of environments comprising a few damped anharmonic modes where the matching criterion fails to work due to the non-Gaussianity of the environment. For both continuous and discrete damped anharmonic environments, this work investigates their unique influence on the system's dynamics, including enhanced non-Markovianity compared to a harmonic approximation, novel effects in the absorption spectrum, and the potential for anharmonic environments to enhance energy transport.

quant-pharXiv:2608.27163•Aug 28, 2026•0 citations

Markovian evolution from a novel scheme

Fardin Kheirandish, Zahra Iranshahi, Fatemeh Bakhtiari +1 more

The Markovian dynamics of open quantum many-body systems are typically governed by the Lindblad master equation, yet obtaining the reduced density matrix for bosonic and fermionic systems remains a formidable numerical challenge due to the exponential growth of the Hilbert space. Here, we introduce a computationally efficient framework for constructing the reduced density matrix by modeling the environment as a copy of the primary system with a monotonically decaying coupling that enforces unidirectional energy flow. Our method directly yields exact solutions that rigorously satisfy the Lindblad master equation for both bosonic and fermionic cases, bypassing the need for costly Liouvillian diagonalization. We validate our approach through applications to paradigmatic models, demonstrating accurate reproduction of dissipative dynamics across a broad parameter regime. This work provides a straightforward and powerful tool for simulating Markovian open-system evolution, with immediate applicability to quantum transport and control problems in many-body physics.

quant-pharXiv:2607.13809•Jul 16, 2026•0 citations

Quantum Dynamics of Interacting oscillators in a thermal medium: A novel scheme

Ronak Moradi, Fardin Kheirandish

We develop an exact analytical framework for studying the nonequilibrium dynamics of interacting quantum harmonic oscillators coupled to a modelled thermal reservoir and driven by external classical fields. The total Hamiltonian is diagonalized via successive Bogoliubov transformations, enabling a nonperturbative treatment of dissipation and driving. For two coupled oscillators with an external source applied to the first, we derive explicit energy expressions and analyze their dependence on coupling and bath parameters. The phase-space structure is characterized through Husimi $Q$-functions for initial separable coherent and number states. We obtain the reduced density matrix elements in the number-state basis and demonstrate that, in the absence of driving and at zero temperature, the reduced density matrix of the main system satisfies the Lindblad master equation. A generating function for the Husimi function is introduced to facilitate computation of higher-order correlations. The analysis is then generalized to $n$ interacting oscillators, with analytical energy expressions provided and the case $n=3$ examined in detail. Our results establish a versatile and exact framework for driven-dissipative quantum systems, with potential applications in quantum thermodynamics, open quantum systems, and many-body physics.

quant-pharXiv:2605.24478•May 26, 2026•0 citations

Convergence of numerical approximations to non-Markovian bosonic gaussian environments

Rahul Trivedi, Daniel Malz, J. I. Cirac

Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. A common approach to the analysis of such systems is to approximate the non-Markovian environment by discrete bosonic modes thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been proposed in previous works [D. Tamascelli et al, PRL (2012), A. W. Chin et al, J. of Math. Phys (2010)], the resulting approximation lacks rigorous convergence guarantees. In this paper, we initiate a rigorous study of the convergence properties of these methods. We show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that, for most physically interesting models of non-Markovian environments, the approximation error falls off polynomially with the number of modes. Our results lend rigor to numerical methods used for approximating non-Markovian quantum dynamics and allow for a quantitative assessment of classical as well as quantum algorithms in simulating non-Markovian quantum systems.

•Jul 15, 2021•1 citation

A method for the systematic construction of few-body damped harmonic oscillator networks accurately reproducing the effect of general bosonic environments in open quantum systems is presented. Under the sole assumptions of a Gaussian environment and regardless of the system coupled to it, an algorithm to determine the parameters of an equivalent set of interacting damped oscillators obeying a Markovian quantum master equation is introduced. By choosing a suitable coupling to the system and minimizing an appropriate distance between the two-time correlation function of this effective bath and that of the target environment, the error induced in the reduced dynamics of the system is brought under rigorous control. The interactions among the effective modes provide remarkable flexibility in replicating non-Markovian effects on the system even with a small number of oscillators, and the resulting Lindblad equation for the system and the modes may therefore be integrated at a very reasonable computational cost using standard methods for Markovian problems, even in strongly nonperturbative coupling regimes and at arbitrary temperatures including zero. We apply the method to an exactly solvable problem in order to demonstrate its accuracy, and present two studies based on current research in the context of coherent transport in biological aggregates and organic photovoltaics as more realistic examples of its use and potential; performance and versatility are highlighted, and theoretical and numerical advantages over existing methods, as well as possible future improvements, are discussed.

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.

We present a quantum algorithm for simulating open quantum systems coupled to Gaussian environments valid for any configuration and coupling strength. The algorithm is applicable to problems with strongly coupled, or non-Markovian, environments; problems with multiple environments out of mutual equilibrium; and problems with time-dependent Hamiltonians. We show that the algorithm can reproduce the true dynamics of such problems at arbitrary accuracy and, for a broad range of problems, only adds a minor resource cost relative to Trotterized time evolution; the cost is low-degree polynomial in the inverse target accuracy. The algorithm is based on the insight that any Gaussian environment can be represented as a train of ancillary qubits that sequentially interact with the system through a time-local coupling, given by the convolution square root of the bath correlation function; this is a secondary result of our work. In numerical simulations of a spin-boson model, we find accurate simulations can be achieved with moderate ( 10 1 – 10 3 ) ancilla counts in the intermediate to strong-coupling regimes. Our results open up applications of quantum computers for efficient simulation of nonequilibrium and open quantum systems.

In open quantum systems, the environment is typically modelled as a collection of harmonic oscillators, whereas realistic environments often exhibit unique non-Markovian effects due to the anharmonicity of the environment. Here we present a framework for modelling anharmonic environments in open quantum systems numerically exactly, in which, via a matching criterion, we map the anharmonic environment onto an effective harmonic one. This results in an effective spectral density that is a collection of shifted bare spectra arising from the non-equidistant energy gaps of the anharmonic environment, along with a unique zero-frequency term that effectively acts as static disorder on the system. This is due to the non-zero variance of the diagonal of the bath coupling operator. In tandem, we develop a framework for investigating the effects of environments comprising a few damped anharmonic modes where the matching criterion fails to work due to the non-Gaussianity of the environment. For both continuous and discrete damped anharmonic environments, this work investigates their unique influence on the system's dynamics, including enhanced non-Markovianity compared to a harmonic approximation, novel effects in the absorption spectrum, and the potential for anharmonic environments to enhance energy transport.

The Markovian dynamics of open quantum many-body systems are typically governed by the Lindblad master equation, yet obtaining the reduced density matrix for bosonic and fermionic systems remains a formidable numerical challenge due to the exponential growth of the Hilbert space. Here, we introduce a computationally efficient framework for constructing the reduced density matrix by modeling the environment as a copy of the primary system with a monotonically decaying coupling that enforces unidirectional energy flow. Our method directly yields exact solutions that rigorously satisfy the Lindblad master equation for both bosonic and fermionic cases, bypassing the need for costly Liouvillian diagonalization. We validate our approach through applications to paradigmatic models, demonstrating accurate reproduction of dissipative dynamics across a broad parameter regime. This work provides a straightforward and powerful tool for simulating Markovian open-system evolution, with immediate applicability to quantum transport and control problems in many-body physics.

We develop an exact analytical framework for studying the nonequilibrium dynamics of interacting quantum harmonic oscillators coupled to a modelled thermal reservoir and driven by external classical fields. The total Hamiltonian is diagonalized via successive Bogoliubov transformations, enabling a nonperturbative treatment of dissipation and driving. For two coupled oscillators with an external source applied to the first, we derive explicit energy expressions and analyze their dependence on coupling and bath parameters. The phase-space structure is characterized through Husimi

Q

-functions for initial separable coherent and number states. We obtain the reduced density matrix elements in the number-state basis and demonstrate that, in the absence of driving and at zero temperature, the reduced density matrix of the main system satisfies the Lindblad master equation. A generating function for the Husimi function is introduced to facilitate computation of higher-order correlations. The analysis is then generalized to

n

interacting oscillators, with analytical energy expressions provided and the case

n=3

examined in detail. Our results establish a versatile and exact framework for driven-dissipative quantum systems, with potential applications in quantum thermodynamics, open quantum systems, and many-body physics.

Non-Markovian effects are important in modeling the behavior of open quantum systems arising in solid-state physics, quantum optics as well as in study of biological and chemical systems. A common approach to the analysis of such systems is to approximate the non-Markovian environment by discrete bosonic modes thus mapping it to a Lindbladian or Hamiltonian simulation problem. While systematic constructions of such modes have been proposed in previous works [D. Tamascelli et al, PRL (2012), A. W. Chin et al, J. of Math. Phys (2010)], the resulting approximation lacks rigorous convergence guarantees. In this paper, we initiate a rigorous study of the convergence properties of these methods. We show that under some physically motivated assumptions on the system-environment interaction, the finite-time dynamics of the non-Markovian open quantum system computed with a sufficiently large number of modes is guaranteed to converge to the true result. Furthermore, we show that, for most physically interesting models of non-Markovian environments, the approximation error falls off polynomially with the number of modes. Our results lend rigor to numerical methods used for approximating non-Markovian quantum dynamics and allow for a quantitative assessment of classical as well as quantum algorithms in simulating non-Markovian quantum systems.