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

Chinese Control Conference (CCC)

Accelerated Simulation of Master Equation for Open Quantum Systems

Bo Yue, Yanfang Wang, Shibei Xue

Published Jul 24, 2023• 0 citations
Publisher

Abstract

Simulation of open quantum systems is undoubtedly critical to design of quantum devices and applications of quantum technologies since dynamics of open quantum systems cannot be solved analytically in general. However, the dimension of open quantum systems grows exponentially as the number of their subsystems increases or as involving complicated environments, which reduces the computational efficiency of the simulation. In this paper, we boost numerical simulation of master equations of open quantum systems from several aspects with which the total time and the occupied memory in the simulation are significantly reduced without losing precision. Our method can be applied to the simulation of non-Markovian quantum systems described by an augmented system model as well. We conduct numerical experiments on simulation of a non-Markovian spin chain. Compared with the algorithm in a well-known toolbox, QuTiP, our method can take 26.3% time and less memory under fair conditions.

DOI
10.23919/ccc58697.2023.10240232

Related papers

Accelerated Simulation of Non-Markovian Quantum Systems Using Quantum Trajectory

Yanfang Wang, Bo Yue, Shibei Xue

Quantum trajectory can be viewed as a numerical simulation method for wavefunctions which consume smaller memories compared to density matrices in master equations. In this paper, we adopt quantum trajectory to simulate the evolution of a non-Markovian quantum system modeled by an augmented system and further use parallel computing to improve efficiency. Compared to the Lindblad master equation solver in QuTiP, it can run more efficiently using less memory. In the example of a 8-spin-chain non-Markovian system, it only takes about 0.57% of the memory footprint and 4.12% of the runtime of the solver for the Lindblad master equation.

Chinese Control Conference (CCC)•Jul 24, 2023•0 citations

Quantum algorithm for simulating the dynamics of an open quantum system

Hefeng Wang, Sahel Ashhab, Franco Nori

In the study of open quantum systems, one typically obtains the decoherence dynamics by solving a master equation. The master equation is derived using knowledge of some basic properties of the system, the environment, and their interaction: One basically needs to know the operators through which the system couples to the environment and the spectral density of the environment. For a large system, it could become prohibitively difficult to even write down the appropriate master equation, let alone solve it on a classical computer. In this paper, we present a quantum algorithm for simulating the dynamics of an open quantum system. On a quantum computer, the environment can be simulated using ancilla qubits with properly chosen single-qubit frequencies and with properly designed coupling to the system qubits. The parameters used in the simulation are easily derived from the parameters of the system $+$ environment Hamiltonian. The algorithm is designed to simulate Markovian dynamics, but it can also be used to simulate non-Markovian dynamics provided that this dynamics can be obtained by embedding the system of interest into a larger system that obeys Markovian dynamics. We estimate the resource requirements for the algorithm. In particular, we show that for sufficiently slow decoherence a single ancilla qubit could be sufficient to represent the entire environment, in principle.

arXiv:1103.3377•Physical Review A•Jun 14, 2011•106 citations

Simulating Non-Markovian Dynamics in Open Quantum Systems

Meng Xu, Vasilii Vadimov, J. T. Stockburger +1 more

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.

quant-pharXiv:2601.02160•Reviews of Modern Physics•Jan 6, 2026•11 citations

Efficient classical simulation of open bosonic quantum systems

Akseli Mäkinen, Joni Ikonen, Takaaki Aoki +3 more

We propose a computationally efficient method to solve the dynamics of operators of bosonic quantum systems coupled to their environments. The method maps the operator under interest to a set of complex-valued functions, and its adjoint master equation to a set of partial differential equations for these functions, which are subsequently solved numerically. In the limit of weak coupling to the environment, the mapping of the operator enables storing the operator efficiently during the simulation, leading to approximately quadratic improvement in the memory consumption compared with the direct approach of solving the adjoint master equation in the number basis, while retaining the computation time comparable. Moreover, the method enables efficient parallelization which allows to optimize for the actual computational time to reach an approximately quadratic speed up, while retaining the memory consumption comparable to the direct approach. We foresee the method to prove useful, e.g., for the verification of the operation of superconducting quantum processors.

arXiv:2203.05526•Mar 10, 2022•0 citations

Simulating non-Markovian open quantum dynamics with neural quantum states

Long Cao, Liwei Ge, Daochi Zhang +6 more

Reducing computational scaling for simulating non-Markovian dissipative dynamics using artificial neural networks is both a major focus and formidable challenge in open quantum systems. To enable neural quantum states (NQSs), we encode environmental memory in dissipatons (quasiparticles with characteristic lifetimes), yielding the dissipaton-embedded quantum master equation (DQME). The resulting NQS-DQME framework achieves compact representation of many-body correlations and non-Markovian memory. Benchmarking against numerically exact hierarchical equations of motion confirms NQS-DQME maintains comparable accuracy while enhancing scalability and interpretability. This methodology opens new paths to explore non-Markovian open quantum dynamics in previously intractable systems.

arXiv:2404.11093•Physical review. B./Physical review. B•Mar 13, 2026•3 citations

Tensor network methods for non-perturbative dynamics of open quantum systems

Thibaut Lacroix, Adam Burgess, Nicola Lorenzoni +11 more

The description of open quantum system dynamics beyond the perturbative treatment (usually associated with Markovian master equations) is a computationally challenging task due to the unfavorable exponential scaling of memory kernels. Developed over recent decades in the context of quantum information and condensed matter, tensor networks provide both a new formalism and a toolbox for overcoming previous computational bottlenecks. This framework enables the formulation of non-perturbative, numerically exact methods for describing the dynamics of open quantum systems to controllable numerical accuracy. In this review, we present these methods and discuss their commonalities and differences to paint a comprehensive view of the field.

quant-pharXiv:2608.09850•Aug 11, 2026•0 citations
Simulation of open quantum systems is undoubtedly critical to design of quantum devices and applications of quantum technologies since dynamics of open quantum systems cannot be solved analytically in general. However, the dimension of open quantum systems grows exponentially as the number of their subsystems increases or as involving complicated environments, which reduces the computational efficiency of the simulation. In this paper, we boost numerical simulation of master equations of open quantum systems from several aspects with which the total time and the occupied memory in the simulation are significantly reduced without losing precision. Our method can be applied to the simulation of non-Markovian quantum systems described by an augmented system model as well. We conduct numerical experiments on simulation of a non-Markovian spin chain. Compared with the algorithm in a well-known toolbox, QuTiP, our method can take 26.3% time and less memory under fair conditions.
Quantum trajectory can be viewed as a numerical simulation method for wavefunctions which consume smaller memories compared to density matrices in master equations. In this paper, we adopt quantum trajectory to simulate the evolution of a non-Markovian quantum system modeled by an augmented system and further use parallel computing to improve efficiency. Compared to the Lindblad master equation solver in QuTiP, it can run more efficiently using less memory. In the example of a 8-spin-chain non-Markovian system, it only takes about 0.57% of the memory footprint and 4.12% of the runtime of the solver for the Lindblad master equation.
In the study of open quantum systems, one typically obtains the decoherence dynamics by solving a master equation. The master equation is derived using knowledge of some basic properties of the system, the environment, and their interaction: One basically needs to know the operators through which the system couples to the environment and the spectral density of the environment. For a large system, it could become prohibitively difficult to even write down the appropriate master equation, let alone solve it on a classical computer. In this paper, we present a quantum algorithm for simulating the dynamics of an open quantum system. On a quantum computer, the environment can be simulated using ancilla qubits with properly chosen single-qubit frequencies and with properly designed coupling to the system qubits. The parameters used in the simulation are easily derived from the parameters of the system +++ environment Hamiltonian. The algorithm is designed to simulate Markovian dynamics, but it can also be used to simulate non-Markovian dynamics provided that this dynamics can be obtained by embedding the system of interest into a larger system that obeys Markovian dynamics. We estimate the resource requirements for the algorithm. In particular, we show that for sufficiently slow decoherence a single ancilla qubit could be sufficient to represent the entire environment, in principle.
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.
We propose a computationally efficient method to solve the dynamics of operators of bosonic quantum systems coupled to their environments. The method maps the operator under interest to a set of complex-valued functions, and its adjoint master equation to a set of partial differential equations for these functions, which are subsequently solved numerically. In the limit of weak coupling to the environment, the mapping of the operator enables storing the operator efficiently during the simulation, leading to approximately quadratic improvement in the memory consumption compared with the direct approach of solving the adjoint master equation in the number basis, while retaining the computation time comparable. Moreover, the method enables efficient parallelization which allows to optimize for the actual computational time to reach an approximately quadratic speed up, while retaining the memory consumption comparable to the direct approach. We foresee the method to prove useful, e.g., for the verification of the operation of superconducting quantum processors.
Reducing computational scaling for simulating non-Markovian dissipative dynamics using artificial neural networks is both a major focus and formidable challenge in open quantum systems. To enable neural quantum states (NQSs), we encode environmental memory in dissipatons (quasiparticles with characteristic lifetimes), yielding the dissipaton-embedded quantum master equation (DQME). The resulting NQS-DQME framework achieves compact representation of many-body correlations and non-Markovian memory. Benchmarking against numerically exact hierarchical equations of motion confirms NQS-DQME maintains comparable accuracy while enhancing scalability and interpretability. This methodology opens new paths to explore non-Markovian open quantum dynamics in previously intractable systems.
The description of open quantum system dynamics beyond the perturbative treatment (usually associated with Markovian master equations) is a computationally challenging task due to the unfavorable exponential scaling of memory kernels. Developed over recent decades in the context of quantum information and condensed matter, tensor networks provide both a new formalism and a toolbox for overcoming previous computational bottlenecks. This framework enables the formulation of non-perturbative, numerically exact methods for describing the dynamics of open quantum systems to controllable numerical accuracy. In this review, we present these methods and discuss their commonalities and differences to paint a comprehensive view of the field.