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:2209.04956 · Journal of Chemical Theory and Computation

Simulating Open Quantum System Dynamics on NISQ Computers with Generalized Quantum Master Equations

Yuchen Wang, Ellen Mulvihill, Zixuan Hu, Ningyi Lyu, Saurabh Shivpuje, Yudan Liu, Micheline B. Soley, Eitan Geva, Víctor S. Batista, Sabre Kais

Published May 26, 2023• 45 citations
arXivHTMLPDF

Abstract

We present a quantum algorithm based on the generalized quantum master equation (GQME) approach to simulate open quantum system dynamics on noisy intermediate-scale quantum (NISQ) computers. This approach overcomes the limitations of the Lindblad equation, which assumes weak system-bath coupling and Markovity, by providing a rigorous derivation of the equations of motion for any subset of elements of the reduced density matrix. The memory kernel resulting from the effect of the remaining degrees of freedom is used as input to calculate the corresponding non-unitary propagator. We demonstrate how the Sz.-Nagy dilation theorem can be employed to transform the non-unitary propagator into a unitary one in a higher-dimensional Hilbert space, which can then be implemented on quantum circuits of NISQ computers. We validate our quantum algorithm as applied to the spin-boson benchmark model by analyzing the impact of the quantum circuit depth on the accuracy of the results when the subset is limited to the diagonal elements of the reduced density matrix. Our findings demonstrate that our approach yields reliable results on NISQ IBM computers.

DOI
10.1021/acs.jctc.3c00316

Related papers

Simulating Non-Markovian Quantum Dynamics on NISQ Computers Using the Hierarchical Equations of Motion

Xiaohan Dan, Eitan Geva, Víctor S. Batista

Quantum computing offers promising new avenues for tackling the long-standing challenge of simulating the quantum dynamics of complex chemical systems, particularly open quantum systems coupled to external baths. However, simulating such nonunitary dynamics on quantum computers is challenging since quantum circuits are specifically designed to carry out unitary transformations. Furthermore, chemical systems are often strongly coupled to the surrounding environment, rendering the dynamics non-Markovian and beyond the scope of Markovian quantum master equations like Lindblad or Redfield. In this work, we introduce a quantum algorithm designed to simulate non-Markovian dynamics of open quantum systems. Our approach enables the implementation of arbitrary quantum master equations on noisy intermediate-scale quantum (NISQ) computers. We illustrate the method as applied in conjunction with the numerically exact hierarchical equations of motion (HEOM) method. The effectiveness of the resulting quantum HEOM algorithm is demonstrated as applied to simulations of the non-Lindbladian electronic energy and charge transfer dynamics in models of the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran and the Fenna–Matthews–Olson complex.

arXiv:2411.12049•Journal of Chemical Theory and Computation•Feb 14, 2025•22 citations

QFlux: An Open-Source Toolkit for Quantum Dynamics Simulations on Quantum Computers. Part VI - The Generalized Quantum Master Equation

Xiaohan Dan, Pouya Khazaei, Brandon C Allen +8 more

Simulating quantum dynamics at finite temperature in complex chemical systems remains a central challenge in quantum chemistry and materials science. In many cases, it is advantageous to focus on the dynamics of a subsystem of interest, represented by its reduced density matrix, whose interaction with the surrounding environment gives rise to open-system behavior. The Nakajima-Zwanzig generalized quantum master equation (GQME) provides a formally exact framework for capturing such non-Markovian dynamics, where memory effects play a defining role. Here, in Part VI of the QFlux tutorial series, we explore methods for simulating non-Markovian open quantum systems described by the GQME on quantum computers. The approach leverages the Sz.-Nagy dilation theorem to embed memory-dependent, non-unitary evolution into an enlarged unitary space suitable for quantum simulation. Using the spin-boson model as a prototypical example, we demonstrate how to formulate, implement, and analyze GQME dynamics within the QFlux open-source platform. This tutorial provides readers with both the theoretical foundation and practical tools required to study non-Markovian quantum dynamics, bridging fundamental concepts with executable workflows for quantum simulation.

•Jan 30, 2026•1 citation

Variational quantum simulation of time-local quantum master equations via quantum jump

Zhihao Lan, Jie Liu, Zhenyu Li +1 more

Strong coupling and environmental memory render many open quantum systems intractable to classical computation. To overcome this barrier, we present a variational quantum algorithm capable of solving generalized form time-local quantum master equations directly on noisy intermediate-scale quantum (NISQ) processors. Our approach combines four key innovations: (1) a pair-vector stochastic Schrödinger equation to unravel time-local quantum master equations, (2) deterministic evolution projected onto parameterized circuits via McLachlan's principle and stochastic jump implementation via singular-value decomposition, (3) using Hadamard tests circuits to measure the jump rates and the element of reduced density matrices, and (4) using no-evolution sampling to sample trajectories. The resulting protocol remains resilient under realistic hardware noise. We validate it on both classical simulators and a 66-qubit superconducting processor for solving the Redfield master equation and fourth-order time-local non-Markovian master equation and find that it can successfully capture key features of complex open quantum dynamics-including non-Markovian oscillations and strong-coupling effects. Our approach establishes a practical pathway toward scalable simulations of open quantum systems in the NISQ era.

The Journal of Chemical Physics•Oct 9, 2025•2 citations

A modified approach for simulating electronically nonadiabatic dynamics via the generalized quantum master equation

Ellen Mulvihill, Alexander Schubert, Xiang Sun +2 more

We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.

The Journal of Chemical Physics•Jan 16, 2019•58 citations

Simulating the dynamics of electronic observables via reduced-dimensionality generalized quantum master equations

Ellen Mulvihill, Eitan Geva

We describe a general-purpose framework for formulating the dynamics of any subset of electronic reduced density matrix elements in terms of a formally exact generalized quantum master equation (GQME). Within this framework, the effect of coupling to the nuclear degrees of freedom, as well as to any projected-out electronic reduced density matrix elements, is captured by a memory kernel and an inhomogeneous term, whose dimensionalities are dictated by the number of electronic reduced density matrix elements included in the subset of interest. We show that the memory kernel and inhomogeneous term within such GQMEs can be calculated from projection-free inputs of the same dimensionality, which can be cast in terms of the corresponding subsets of overall system two-time correlation functions. The applicability and feasibility of such reduced-dimensionality GQMEs is demonstrated on the two-state spin-boson benchmark model. To this end, we compare and contrast the following four types of GQMEs: (1) a full density matrix GQME, (2) a single-population scalar GQME, (3) a populations-only GQME, and (4) a subset GQME for any combination of populations and coherences. Using a method based on the mapping Hamiltonian approach and linearized semiclassical approximation to calculate the projection-free inputs, we find that while single-population GQMEs and subset GQMEs containing only one population are less accurate, they can still produce reasonable results and that the accuracy of the results obtained via the populations-only GQME and a subset GQME containing both populations is comparable to that obtained via the full density matrix GQMEs.

The Journal of Chemical Physics•Jan 11, 2022•20 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
We present a quantum algorithm based on the generalized quantum master equation (GQME) approach to simulate open quantum system dynamics on noisy intermediate-scale quantum (NISQ) computers. This approach overcomes the limitations of the Lindblad equation, which assumes weak system-bath coupling and Markovity, by providing a rigorous derivation of the equations of motion for any subset of elements of the reduced density matrix. The memory kernel resulting from the effect of the remaining degrees of freedom is used as input to calculate the corresponding non-unitary propagator. We demonstrate how the Sz.-Nagy dilation theorem can be employed to transform the non-unitary propagator into a unitary one in a higher-dimensional Hilbert space, which can then be implemented on quantum circuits of NISQ computers. We validate our quantum algorithm as applied to the spin-boson benchmark model by analyzing the impact of the quantum circuit depth on the accuracy of the results when the subset is limited to the diagonal elements of the reduced density matrix. Our findings demonstrate that our approach yields reliable results on NISQ IBM computers.
Quantum computing offers promising new avenues for tackling the long-standing challenge of simulating the quantum dynamics of complex chemical systems, particularly open quantum systems coupled to external baths. However, simulating such nonunitary dynamics on quantum computers is challenging since quantum circuits are specifically designed to carry out unitary transformations. Furthermore, chemical systems are often strongly coupled to the surrounding environment, rendering the dynamics non-Markovian and beyond the scope of Markovian quantum master equations like Lindblad or Redfield. In this work, we introduce a quantum algorithm designed to simulate non-Markovian dynamics of open quantum systems. Our approach enables the implementation of arbitrary quantum master equations on noisy intermediate-scale quantum (NISQ) computers. We illustrate the method as applied in conjunction with the numerically exact hierarchical equations of motion (HEOM) method. The effectiveness of the resulting quantum HEOM algorithm is demonstrated as applied to simulations of the non-Lindbladian electronic energy and charge transfer dynamics in models of the carotenoid-porphyrin-C 60 molecular triad dissolved in tetrahydrofuran and the Fenna–Matthews–Olson complex.
Simulating quantum dynamics at finite temperature in complex chemical systems remains a central challenge in quantum chemistry and materials science. In many cases, it is advantageous to focus on the dynamics of a subsystem of interest, represented by its reduced density matrix, whose interaction with the surrounding environment gives rise to open-system behavior. The Nakajima-Zwanzig generalized quantum master equation (GQME) provides a formally exact framework for capturing such non-Markovian dynamics, where memory effects play a defining role. Here, in Part VI of the QFlux tutorial series, we explore methods for simulating non-Markovian open quantum systems described by the GQME on quantum computers. The approach leverages the Sz.-Nagy dilation theorem to embed memory-dependent, non-unitary evolution into an enlarged unitary space suitable for quantum simulation. Using the spin-boson model as a prototypical example, we demonstrate how to formulate, implement, and analyze GQME dynamics within the QFlux open-source platform. This tutorial provides readers with both the theoretical foundation and practical tools required to study non-Markovian quantum dynamics, bridging fundamental concepts with executable workflows for quantum simulation.
Strong coupling and environmental memory render many open quantum systems intractable to classical computation. To overcome this barrier, we present a variational quantum algorithm capable of solving generalized form time-local quantum master equations directly on noisy intermediate-scale quantum (NISQ) processors. Our approach combines four key innovations: (1) a pair-vector stochastic Schrödinger equation to unravel time-local quantum master equations, (2) deterministic evolution projected onto parameterized circuits via McLachlan's principle and stochastic jump implementation via singular-value decomposition, (3) using Hadamard tests circuits to measure the jump rates and the element of reduced density matrices, and (4) using no-evolution sampling to sample trajectories. The resulting protocol remains resilient under realistic hardware noise. We validate it on both classical simulators and a 66-qubit superconducting processor for solving the Redfield master equation and fourth-order time-local non-Markovian master equation and find that it can successfully capture key features of complex open quantum dynamics-including non-Markovian oscillations and strong-coupling effects. Our approach establishes a practical pathway toward scalable simulations of open quantum systems in the NISQ era.
We present a modified approach for simulating electronically nonadiabatic dynamics based on the Nakajima-Zwanzig generalized quantum master equation (GQME). The modified approach utilizes the fact that the Nakajima-Zwanzig formalism does not require casting the overall Hamiltonian in system-bath form, which is arguably neither natural nor convenient in the case of the Hamiltonian that governs nonadiabatic dynamics. Within the modified approach, the effect of the nuclear degrees of freedom on the time evolution of the electronic reduced density operator is fully captured by a memory kernel super-operator. A methodology for calculating the memory kernel from projection-free inputs is developed. Simulating the electronic dynamics via the modified approach, with a memory kernel obtained using exact or approximate methods, can be more cost effective and/or lead to more accurate results than direct application of those methods. The modified approach is compared to previously proposed GQME-based approaches, and its robustness and accuracy are demonstrated on a benchmark spin-boson model with a memory kernel which is calculated within the Ehrenfest method.
We describe a general-purpose framework for formulating the dynamics of any subset of electronic reduced density matrix elements in terms of a formally exact generalized quantum master equation (GQME). Within this framework, the effect of coupling to the nuclear degrees of freedom, as well as to any projected-out electronic reduced density matrix elements, is captured by a memory kernel and an inhomogeneous term, whose dimensionalities are dictated by the number of electronic reduced density matrix elements included in the subset of interest. We show that the memory kernel and inhomogeneous term within such GQMEs can be calculated from projection-free inputs of the same dimensionality, which can be cast in terms of the corresponding subsets of overall system two-time correlation functions. The applicability and feasibility of such reduced-dimensionality GQMEs is demonstrated on the two-state spin-boson benchmark model. To this end, we compare and contrast the following four types of GQMEs: (1) a full density matrix GQME, (2) a single-population scalar GQME, (3) a populations-only GQME, and (4) a subset GQME for any combination of populations and coherences. Using a method based on the mapping Hamiltonian approach and linearized semiclassical approximation to calculate the projection-free inputs, we find that while single-population GQMEs and subset GQMEs containing only one population are less accurate, they can still produce reasonable results and that the accuracy of the results obtained via the populations-only GQME and a subset GQME containing both populations is comparable to that obtained via the full density matrix GQMEs.
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