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:2507.19323 · The Journal of Chemical Physics

On the discretization error of the discrete generalized quantum master equation

Ruojing Peng, Lachlan P. Lindoy, Joonho Lee

Published Oct 1, 2025• 5 citations
arXivHTMLPDF

Abstract

The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112, 110401 (2014)] can be considered a discrete-time formulation of the Nakajima-Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21, 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. This work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel KN and the continuous-time NZ-QME kernel K(NΔt). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δt → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri's work. The relative performance of various schemes for either computing accurate K(NΔt) from exact dynamics or obtaining accurate dynamics from exact K(NΔt) warrants further investigation.

DOI
10.1063/5.0293127

Related papers

Discrete Generalized Quantum Master Equations

Nancy Makri

Several derivative and integral approximations are explored for discretizing the Nakajima-Zwanzig generalized quantum master equation (NZ-QME or GQME) to obtain discrete quantum master equation (DQME) hierarchies and relations between discrete memory kernel and reduced density matrix (RDM) elements. It is shown that the simplest forward-difference approximation does not allow the reliable determination of the discrete kernel elements, even in the infinitesimal time-step limit, and that discrete kernels obtained in earlier work are flawed, although the procedure can be remedied. The various approximations give rise to DQMEs that differ in structure and in the RDM-kernel relationships. It is shown that the use of a more accurate discretization based on the midpoint derivative and midpoint integral approximations leads to a DQME that exhibits endpoint effects, which reflect the weaker impact of the bath on the RDM during the first time step and which parallel those encountered in the small matrix decomposition of the path integral (SMatPI) with a symmetric factorization of the short-time propagator. The features of the DQME hierarchies and RDM-kernel relations are illustrated through analytical examples involving a simple integrodifferential equation and a scalar GQME model, as well as numerical results for a two-level system (TLS) coupled to a harmonic bath.

Journal of Chemical Theory and Computation•May 6, 2025•10 citations

Comment on "Unified framework for open quantum dynamics with memory"

Nancy Makri, Sohang Kundu, Zhenning Cai +1 more

A recent article by Ivander, Lindoy and Lee [Nature Communications 15, 8087 (2024)] claims to discover the relationship between the generalized quantum master equation (GQME) and the path integral for a system coupled to a harmonic bath. However, this relationship was already established in 2020 by Makri in the context of the small matrix decomposition of the path integral (SMatPI) [J. Chem. Theory and Comput. 16, 4038 (2020)]. The procedure that this article uses in its Supplementary Information (SI) to obtain the various matrices follows the SMatPI decomposition steps for the alternative Trotter ordering. The absence of endpoint effects in the kernel matrices of the discretized GQME expression for the reduced density matrix (RDM) is the consequence of a crude GQME discretization and is not consistent with the SMatPI decomposition of an auxiliary matrix presented in the SI. This form is identical to the transfer tensor method (TTM) of Cerrillo and Cao [Phys. Rev. Lett. 112, 110401 (2014)]. Further, the Dyck path section of this article follows precisely the diagrammatic analysis developed by Wang and Cai in a recent paper [Communications in Computational Physics 36, 389 (2024)]. We elaborate on these three critiques in this Comment.

arXiv:2410.08239•Oct 9, 2024•2 citations

The Transfer Tensor Method: an Analytical Study Case

Marcel Morillas-Rozas, Alberto López-García, Gonzalo Reina Rivero +2 more

The transfer tensor method is a versatile tool for analyzing and propagating general open quantum systems. It captures in a compact manner all memory effects in a non-Markovian system through a straightforward transformation of a set of dynamical maps. Transfer tensors provide the exact convolutional propagator associated with a given time discretization over the past evolution of an open quantum system. Here we show that, for any finite time discretization, the memory kernel of the Nakajima Zwanzig equation deviates from the exact transfer tensors, although both converge in the continuous-time limit, as expected. We examine this behaviour in the context of an analytically solvable model: a two level atom resonant with a lossy cavity in the Jaynes Cummings limit. The atomic dynamics separate into two decoupled degrees of freedom -- the coherence and the population inversion. We derive exact expressions for the dynamical map, the transfer tensors and the memory kernel governing the coherence, and we relate them to their counterparts for the population inversion. As a function of the ratio between the cavity loss rate and the atom-cavity coupling strength, we identify regions of enhanced non-Markovianity in which the system can be described as fully Markovian for certain time-step choices.

quant-pharXiv:2603.08458•The Journal of Chemical Physics•Mar 10, 2026•1 citation

Applicability of transfer tensor method for open quantum system dynamics

Andrius Gelžinis, Edvardas Rybakovas, Leonas Valkūnas

Accurate simulations of open quantum system dynamics is a long standing issue in the field of chemical physics. Exact methods exist, but are costly, while perturbative methods are limited in their applicability. Recently a new black-box type method, called transfer tensor method (TTM), was proposed [J. Cerrillo and J. Cao, Phys. Rev. Lett. 112, 110401 (2014)]. It allows one to accurately simulate long time dynamics with a numerical cost of solving a time-convolution master equation, provided many initial system evolution trajectories are obtained from some exact method beforehand. The possible time-savings thus strongly depend on the ratio of total versus initial evolution lengths. In this work, we investigate the parameter regimes where an application of TTM would be most beneficial in terms of computational time. We identify several promising parameter regimes. Although some of them correspond to cases when perturbative theories could be expected to perform well, we find that the accuracy of such approaches depends on system parameters in a more complex way than it is commonly thought. We propose that the TTM should be applied whenever system evolution is expected to be long and accuracy of perturbative methods cannot be ensured or in cases when the system under consideration does not correspond to any single perturbative regime.

The Journal of Chemical Physics•Dec 21, 2017•25 citations

Non-Markovian Dynamical Maps: Numerical Processing of Open Quantum Trajectories

Javier Cerrillo, Jianshu Cao

The initial stages of the evolution of an open quantum system encode the key information of its underlying dynamical correlations, which in turn can predict the trajectory at later stages. We propose a general approach based on non-Markovian dynamical maps to extract this information from the initial trajectories and compress it into non-Markovian transfer tensors. Assuming time-translational invariance, the tensors can be used to accurately and efficiently propagate the state of the system to arbitrarily long time scales. The non-Markovian transfer tensor method (TTM) demonstrates the coherent-to-incoherent transition as a function of the strength of quantum dissipation and predicts the noncanonical equilibrium distribution due to the system-bath entanglement. TTM is equivalent to solving the Nakajima-Zwanzig equation and, therefore, can be used to reconstruct the dynamical operators (the system Hamiltonian and memory kernel) from quantum trajectories obtained in simulations or experiments. The concept underlying the approach can be generalized to physical observables with the goal of learning and manipulating the trajectories of an open quantum system.

arXiv:1307.7743•Physical Review Letters•Mar 17, 2014•208 citations

Addendum: Unified framework for open quantum dynamics with memory

Ru‐Wen Peng, Felix Ivander, Lachlan P. Lindoy +1 more

This Addendum presents a detailed analysis of the discretization error in time-integration and time-derivative that appear in the Nakajima-Zwanzig equation. This was brought to our attention by Makri et al. [arXiv:2410.08239]. Our analysis in the Addendum shows that the relationship derived in our earlier work [Nat. Commun. 15, 8087 (2024)] is valid within the choice of discretization and is not contaminated by the discretization error.

Nature Communications•Aug 12, 2025•2 citations
The transfer tensor method (TTM) [Cerrillo and Cao, Phys. Rev. Lett. 112, 110401 (2014)] can be considered a discrete-time formulation of the Nakajima-Zwanzig quantum master equation (NZ-QME) for modeling non-Markovian quantum dynamics. A recent paper [Makri, J. Chem. Theory Comput. 21, 5037 (2025)] raised concerns regarding the consistency of the TTM discretization, particularly a spurious term at the initial time t = 0. This work presents a detailed analysis of the discretization structure of the TTM, clarifying the origin of the initial-time correction and establishing a consistent relationship between the TTM discrete-time memory kernel KN and the continuous-time NZ-QME kernel K(NΔt). This relationship is validated numerically using the spin-boson model, demonstrating convergence of reconstructed memory kernels and accurate dynamical evolution as Δt → 0. While the TTM provides a consistent discretization, we note that alternative schemes are also viable, such as the midpoint derivative/midpoint integral scheme proposed in Makri's work. The relative performance of various schemes for either computing accurate K(NΔt) from exact dynamics or obtaining accurate dynamics from exact K(NΔt) warrants further investigation.
Several derivative and integral approximations are explored for discretizing the Nakajima-Zwanzig generalized quantum master equation (NZ-QME or GQME) to obtain discrete quantum master equation (DQME) hierarchies and relations between discrete memory kernel and reduced density matrix (RDM) elements. It is shown that the simplest forward-difference approximation does not allow the reliable determination of the discrete kernel elements, even in the infinitesimal time-step limit, and that discrete kernels obtained in earlier work are flawed, although the procedure can be remedied. The various approximations give rise to DQMEs that differ in structure and in the RDM-kernel relationships. It is shown that the use of a more accurate discretization based on the midpoint derivative and midpoint integral approximations leads to a DQME that exhibits endpoint effects, which reflect the weaker impact of the bath on the RDM during the first time step and which parallel those encountered in the small matrix decomposition of the path integral (SMatPI) with a symmetric factorization of the short-time propagator. The features of the DQME hierarchies and RDM-kernel relations are illustrated through analytical examples involving a simple integrodifferential equation and a scalar GQME model, as well as numerical results for a two-level system (TLS) coupled to a harmonic bath.
A recent article by Ivander, Lindoy and Lee [Nature Communications 15, 8087 (2024)] claims to discover the relationship between the generalized quantum master equation (GQME) and the path integral for a system coupled to a harmonic bath. However, this relationship was already established in 2020 by Makri in the context of the small matrix decomposition of the path integral (SMatPI) [J. Chem. Theory and Comput. 16, 4038 (2020)]. The procedure that this article uses in its Supplementary Information (SI) to obtain the various matrices follows the SMatPI decomposition steps for the alternative Trotter ordering. The absence of endpoint effects in the kernel matrices of the discretized GQME expression for the reduced density matrix (RDM) is the consequence of a crude GQME discretization and is not consistent with the SMatPI decomposition of an auxiliary matrix presented in the SI. This form is identical to the transfer tensor method (TTM) of Cerrillo and Cao [Phys. Rev. Lett. 112, 110401 (2014)]. Further, the Dyck path section of this article follows precisely the diagrammatic analysis developed by Wang and Cai in a recent paper [Communications in Computational Physics 36, 389 (2024)]. We elaborate on these three critiques in this Comment.
The transfer tensor method is a versatile tool for analyzing and propagating general open quantum systems. It captures in a compact manner all memory effects in a non-Markovian system through a straightforward transformation of a set of dynamical maps. Transfer tensors provide the exact convolutional propagator associated with a given time discretization over the past evolution of an open quantum system. Here we show that, for any finite time discretization, the memory kernel of the Nakajima Zwanzig equation deviates from the exact transfer tensors, although both converge in the continuous-time limit, as expected. We examine this behaviour in the context of an analytically solvable model: a two level atom resonant with a lossy cavity in the Jaynes Cummings limit. The atomic dynamics separate into two decoupled degrees of freedom -- the coherence and the population inversion. We derive exact expressions for the dynamical map, the transfer tensors and the memory kernel governing the coherence, and we relate them to their counterparts for the population inversion. As a function of the ratio between the cavity loss rate and the atom-cavity coupling strength, we identify regions of enhanced non-Markovianity in which the system can be described as fully Markovian for certain time-step choices.
Accurate simulations of open quantum system dynamics is a long standing issue in the field of chemical physics. Exact methods exist, but are costly, while perturbative methods are limited in their applicability. Recently a new black-box type method, called transfer tensor method (TTM), was proposed [J. Cerrillo and J. Cao, Phys. Rev. Lett. 112, 110401 (2014)]. It allows one to accurately simulate long time dynamics with a numerical cost of solving a time-convolution master equation, provided many initial system evolution trajectories are obtained from some exact method beforehand. The possible time-savings thus strongly depend on the ratio of total versus initial evolution lengths. In this work, we investigate the parameter regimes where an application of TTM would be most beneficial in terms of computational time. We identify several promising parameter regimes. Although some of them correspond to cases when perturbative theories could be expected to perform well, we find that the accuracy of such approaches depends on system parameters in a more complex way than it is commonly thought. We propose that the TTM should be applied whenever system evolution is expected to be long and accuracy of perturbative methods cannot be ensured or in cases when the system under consideration does not correspond to any single perturbative regime.
The initial stages of the evolution of an open quantum system encode the key information of its underlying dynamical correlations, which in turn can predict the trajectory at later stages. We propose a general approach based on non-Markovian dynamical maps to extract this information from the initial trajectories and compress it into non-Markovian transfer tensors. Assuming time-translational invariance, the tensors can be used to accurately and efficiently propagate the state of the system to arbitrarily long time scales. The non-Markovian transfer tensor method (TTM) demonstrates the coherent-to-incoherent transition as a function of the strength of quantum dissipation and predicts the noncanonical equilibrium distribution due to the system-bath entanglement. TTM is equivalent to solving the Nakajima-Zwanzig equation and, therefore, can be used to reconstruct the dynamical operators (the system Hamiltonian and memory kernel) from quantum trajectories obtained in simulations or experiments. The concept underlying the approach can be generalized to physical observables with the goal of learning and manipulating the trajectories of an open quantum system.