Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Log inSign up
Challenges
Datasets
Workspaces
Papers
Discussions
Leaderboard
Blog
Job Board
Q3AS

© 2026 Aqora Quantum S.A.S.

TermsPrivacyLegal Notice

arXiv:2509.12717 · Physical Review A

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

Hans Michael Christensen, Johannes Agerskov, Frederik Nathan

Published Aug 19, 2026• 0 citations
arXivHTMLPDF

Abstract

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.

DOI
10.1103/dc52-844x

Related papers

Quantum Simulation of Non-Markovian Open Quantum Systems via Low-Rank Kernel Compression and Coupling Dilution

Hengzhun Chen, Zhenning Cai

We present an error-controlled framework for the quantum simulation of non-Markovian open-system dynamics in Gaussian environments. For a system coupled linearly to a Gaussian environment with zero mean, the bath's influence on the reduced dynamics is completely determined by its two-time correlation kernel. This allows us to replace the microscopic environment with an equivalent bath that has the same kernel. On a finite simulation interval, a rank-$N$ approximation of the kernel yields an effective bath of $N$ bosonic modes with time-dependent couplings to the system. To encode these modes in qubits, we first truncate each mode to its lowest $d$ occupation levels. When stronger couplings lead to large occupations, we introduce coupling dilution, which replaces each mode with several more weakly coupled copies while preserving the bath correlation kernel. This distributes excitations across more modes and trades additional modes for smaller local dimensions. We derive explicit reduced-state error bounds and establish polynomial circuit width and depth in simulation time and inverse accuracy under stated spectral-decay and coupling-regularity assumptions, provided elementary unitaries admit polynomial-cost implementations.

quant-pharXiv:2610.29372•Oct 22, 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

Optimized auxiliary oscillators for the simulation of general open quantum systems

F. Mascherpa, Andrea Smirne, Alejandro D. Somoza +5 more

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.

arXiv:1904.04822•Physical Review A•May 14, 2020•98 citations

Efficient quantum algorithm to simulate open systems through a single environmental qubit

Giovanni Di Bartolomeo, Michele Vischi, Tommaso Feri +2 more

We present an efficient algorithm for simulating open quantum systems dynamics described by the Lindblad master equation on quantum computers, addressing key challenges in the field. In contrast to existing approaches, our method achieves two significant advancements. First, we employ a repetition of unitary gates on a set of n system qubits and, remarkably, only a single ancillary bath qubit representing the environment. It follows that, for the typical case of m locality of the Lindblad operators, we reach an exponential improvement of the number of ancilla in terms of m and up to a polynomial improvement in ancilla overhead for large n with respect to other approaches. Although stochasticity is introduced, requiring multiple circuit realizations, the sampling overhead is independent of the system size. Second, we show that, under fixed accuracy conditions, our algorithm enables a reduction in the number of Trotter steps compared to other approaches, substantially decreasing circuit depth. These advancements hold particular significance for near-term quantum computers, where minimizing both width and depth is critical due to inherent noise in their dynamics. Published by the American Physical Society 2024

arXiv:2311.10009•Physical Review Research•Dec 26, 2024•15 citations

Path integral quantum algorithm for simulating non-Markovian quantum dynamics in open quantum systems

Peter Walters, Fei Wang

Research in open quantum system dynamics has received growing interest in recent years, and its ongoing investigation has uncovered many intriguing physics departing from closed systems. A particularly useful application of the open quantum system theory is to simulate quantum dynamical processes in condensed phase materials. As the quantum degree of freedom is constantly under the influence of its thermal environment, the accurate description of the dynamics often requires non-Markovian time evolution at finite temperature. Such calculation is usually quite challenging on classical computers as extensive memory storage is required. In this work, we focus on a quantum system linearly coupled to its harmonic bath and present a path integral based quantum algorithm that time-evolves the reduced density matrix under finite temperature non-Markovian dynamics. To treat the nonunitary time evolution, the Sz.-Nagy dilation scheme is used for the conversion to unitary dynamics that can be implementable on gate-based quantum computers. The modified Hadamard test is then used to retrieve all the information of the reduced density matrix. Complexity analysis shows that the memory requirement has exponential reduction on the quantum machine whereas the runtime complexity stays roughly the same as for the classical computer. It points to the possibility of using this algorithm to simulate multilevel and multisite non-Markovian quantum dynamics that are beyond the reach of classical computers. The algorithm makes no ad hoc assumptions and extends naturally beyond Markovian and weak coupling regimes. We validated the algorithm on the quantum simulator with the spin-boson model and demonstrated its excellent agreement with the classical computer benchmark. Published by the American Physical Society 2024

Physical Review Research•Feb 1, 2024•13 citations

Efficiently Solving Spin-Boson Dynamics Via Non-Markovian Quantum Trajectories

Zeng-Zhao Li, Cho‐Tung Yip, Hai-Yao Deng +4 more

We develop a systematic and efficient approach for numerically solving the non-Markovian quantum state diffusion equations for open quantum systems coupled to an environment up to arbitrary orders of noises or coupling strengths. As an important application, we consider a real-time simulation of a spin-boson model in a strong coupling regime that is difficult to deal with using conventional methods. We show that the non-Markovian stochastic Schrodinger equation can be efficiently implemented as a real--time simulation for this model, so as to give an accurate description of spin-boson dynamics beyond the rotating-wave approximation.

•Jun 6, 2014•0 citations
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.
We present an error-controlled framework for the quantum simulation of non-Markovian open-system dynamics in Gaussian environments. For a system coupled linearly to a Gaussian environment with zero mean, the bath's influence on the reduced dynamics is completely determined by its two-time correlation kernel. This allows us to replace the microscopic environment with an equivalent bath that has the same kernel. On a finite simulation interval, a rank-NNN approximation of the kernel yields an effective bath of NNN bosonic modes with time-dependent couplings to the system. To encode these modes in qubits, we first truncate each mode to its lowest ddd occupation levels. When stronger couplings lead to large occupations, we introduce coupling dilution, which replaces each mode with several more weakly coupled copies while preserving the bath correlation kernel. This distributes excitations across more modes and trades additional modes for smaller local dimensions. We derive explicit reduced-state error bounds and establish polynomial circuit width and depth in simulation time and inverse accuracy under stated spectral-decay and coupling-regularity assumptions, provided elementary unitaries admit polynomial-cost implementations.
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.
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 present an efficient algorithm for simulating open quantum systems dynamics described by the Lindblad master equation on quantum computers, addressing key challenges in the field. In contrast to existing approaches, our method achieves two significant advancements. First, we employ a repetition of unitary gates on a set of n system qubits and, remarkably, only a single ancillary bath qubit representing the environment. It follows that, for the typical case of m locality of the Lindblad operators, we reach an exponential improvement of the number of ancilla in terms of m and up to a polynomial improvement in ancilla overhead for large n with respect to other approaches. Although stochasticity is introduced, requiring multiple circuit realizations, the sampling overhead is independent of the system size. Second, we show that, under fixed accuracy conditions, our algorithm enables a reduction in the number of Trotter steps compared to other approaches, substantially decreasing circuit depth. These advancements hold particular significance for near-term quantum computers, where minimizing both width and depth is critical due to inherent noise in their dynamics. Published by the American Physical Society 2024
Research in open quantum system dynamics has received growing interest in recent years, and its ongoing investigation has uncovered many intriguing physics departing from closed systems. A particularly useful application of the open quantum system theory is to simulate quantum dynamical processes in condensed phase materials. As the quantum degree of freedom is constantly under the influence of its thermal environment, the accurate description of the dynamics often requires non-Markovian time evolution at finite temperature. Such calculation is usually quite challenging on classical computers as extensive memory storage is required. In this work, we focus on a quantum system linearly coupled to its harmonic bath and present a path integral based quantum algorithm that time-evolves the reduced density matrix under finite temperature non-Markovian dynamics. To treat the nonunitary time evolution, the Sz.-Nagy dilation scheme is used for the conversion to unitary dynamics that can be implementable on gate-based quantum computers. The modified Hadamard test is then used to retrieve all the information of the reduced density matrix. Complexity analysis shows that the memory requirement has exponential reduction on the quantum machine whereas the runtime complexity stays roughly the same as for the classical computer. It points to the possibility of using this algorithm to simulate multilevel and multisite non-Markovian quantum dynamics that are beyond the reach of classical computers. The algorithm makes no ad hoc assumptions and extends naturally beyond Markovian and weak coupling regimes. We validated the algorithm on the quantum simulator with the spin-boson model and demonstrated its excellent agreement with the classical computer benchmark. Published by the American Physical Society 2024
We develop a systematic and efficient approach for numerically solving the non-Markovian quantum state diffusion equations for open quantum systems coupled to an environment up to arbitrary orders of noises or coupling strengths. As an important application, we consider a real-time simulation of a spin-boson model in a strong coupling regime that is difficult to deal with using conventional methods. We show that the non-Markovian stochastic Schrodinger equation can be efficiently implemented as a real--time simulation for this model, so as to give an accurate description of spin-boson dynamics beyond the rotating-wave approximation.