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:2501.04192 · APL Quantum

Benchmarking TCL4: Assessing the usability and reliability of fourth-order approximations

Jiahao Chen, Elyana Crowder, Lian Xiang, Dragomir M Davidović

Published Apr 14, 2025• 4 citations
arXivHTMLPDF

Abstract

The non-Markovian dynamics of an open quantum system can be rigorously derived using the Feynman–Vernon influence functional approach. Although this formalism is exact, practical numerical implementations often require compromises. The time-convolutionless (TCL) master equation offers an exact framework, yet its application typically relies on a perturbative expansion of both the time forward and time backward state propagators. Due to the significant computational effort involved—and the scarcity of analytical solutions for most open quantum systems—the fourth-order perturbative TCL generator (TCL4) has only been benchmarked on a limited range of systems and parameter spaces. Recent advancements, however, have made the computation of TCL4 faster and more accessible. In this paper, we benchmark the TCL4 master equation against numerically exact methods for the biased spin-boson model. We focus on the regime near critical bath coupling where perturbative master equations are expected to become inaccurate. Our findings reveal that the TCL4 approach is most reliable at low temperatures and more efficient than the numerical exact methods. This study aims to delineate the conditions under which the TCL4 perturbative master equation enhances the accuracy of the TCL2.

DOI
10.1063/5.0255350

Related papers

Asymptotic TCL4 Generator for the Spin-Boson Model: Analytical Derivation and Benchmarking

Kumar, Prem, Athulya, K. P., Ghosh, Sibasish

The spin-boson model is a widely used model for understanding the properties of a two-level open quantum system. Accurately describing its dynamics often requires going beyond the weak system-environment coupling approximation. However, calculating the higher-order generators of such a dynamics, with a system-environment coupling that is not too weak, has been known to be challenging, both numerically and analytically. This work presents the analytical derivation of the complete fourth-order time-convolutionless (TCL) generator for a generic spin-boson model, accurate up to 4th order in the system-environment coupling parameter, under the assumption that the environmental spectral density is an odd function of frequency. In the case of a semiconductor double-quantum-dot system, our results reveal corrections to the dynamics that may become physically significant in some parameter regimes. Furthermore, we report that the widely used second-order TCL master equation tends to overestimate the non-Markovianity of a dynamics over a large parameter regime. Within the regime of its applicability, our results provide a computational advantage over numerically exact techniques. The accuracy of the fourth-order TCL generator is rigorously benchmarked against specialized analytical calculations done for the Ohmic spectral density with Drude cutoff and against the numerically exact Hierarchical Equations of Motion technique.

arXiv:2506.17009•Jun 20, 2025•0 citations

Exact generator and its high order expansions in time-convolutionless generalized master equation: Applications to spin-boson model and excitation energy transfer

Yanying Liu, Yaming Yan, Meng Xu +2 more

The time-convolutionless (TCL) quantum master equation provides a powerful tool to simulate reduced dynamics of a quantum system coupled to a bath. The key quantity in the TCL master equation is the so-called kernel or generator, which describes effects of the bath degrees of freedom. Since the exact TCL generators are usually hard to calculate analytically, most applications of the TCL generalized master equation have relied on approximate generators using second and fourth order perturbative expansions. By using the hierarchical equation of motion (HEOM) and extended HEOM methods, we present a new approach to calculating the exact TCL generator and its high order perturbative expansions. The new approach is applied to the spin-boson model with different sets of parameters, to investigate the convergence of the high order expansions of the TCL generator. We also discuss circumstances where the exact TCL generator becomes singular for the spin-boson model, and a model of excitation energy transfer in the Fenna-Matthews-Olson complex.

arXiv:1806.07540•Chinese Journal of Chemical Physics•Aug 1, 2018•15 citations

Expansion of time-convolutionless non-Markovian quantum master equations: A case study using the Fano-Anderson model

Tim Alhäuser, Heinz-Peter Breuer

We explore the performance of the time-convolutionless (TCL) projection operator technique using the Fano-Anderson model as a test case. Comparing the exact TCL master equation with an expansion in powers of the strength of the system-environment coupling, we analyze the transient dynamics as well as the steady-state behavior. For a Lorentzian spectral density we demonstrate that the dimensionless expansion parameter corresponds to the ratio of the environmental correlation time to the relaxation time of the system, and we derive the convergence radius for the TCL expansion, which is seen to depend on the ratio of detuning and width of the spectral density. We further study the quantum non-Markovianity of the model based on the evolution of the Bures distance between quantum states and how it is represented by the second and fourth order of the expansion. Our results highlight both the strengths and the limitations of the TCL formalism in capturing key features of open quantum systems and, in particular, the challenges of accurately describing strongly coupled systems and non-Markovian dynamics.

quant-pharXiv:2604.22010•The Journal of Chemical Physics•Apr 27, 2026•0 citations

Open quantum systems beyond Fermi's golden rule: Diagrammatic expansion of the steady-state time-convolutionless master equations

Michael S. Ferguson, Oded Zilberberg, G. Blatter

Steady-state observables, such as occupation numbers and currents, are crucial experimental signatures in open quantum systems.The time-convolutionless (TCL) master equation, which is both exact and time-local, is an ideal candidate for the perturbative computation of such observables.We develop a diagrammatic approach to evaluate the steady-state TCL generator based on operators rather than superoperators.We obtain the steady-state occupation numbers, extend our formulation to the calculation of currents, and provide a simple physical interpretation of the diagrams.We further benchmark our method on a single noninteracting level coupled to Fermi reservoirs, where we recover the exact expansion to next-to-leading order.The low number of diagrams appearing in our formulation makes the extension to higher orders accessible.Combined, these properties make the steady-state time-convolutionless master equation an effective tool for the calculation of steady-state properties in open quantum systems.

arXiv:2010.09838•Physical Review Research•May 19, 2021•26 citations

Fast Evaluation of the Sixth-Order Time-Convolutionless (TCL6) Master-Equation Generator

Jiahao Chen, Sirui Chen, Dragomir Davidovic

Higher-order time-convolutionless (TCL) master equations provide systematic corrections for non-Markovian dynamics in open quantum systems, but their time-ordered integrals become increasingly costly to evaluate. Direct quadrature of the Hadamard-reduced sixth-order (TCL6) generator at all $N_t$ sampled times requires $O(N_t^3)$ operations at fixed system dimension. We derive an exact reduction for a finite-dimensional system coupled through one Hermitian operator to a stationary centered Gaussian bath. By separating operator coefficients from scalar time-dependent spectral densities, the complete TCL6 time series is reduced to cumulative sums and first-order recurrences, one-dimensional causal convolutions, and an exact dyadic recursion for interlocked sums. With the algorithm, evaluating the complete TCL6 time series requires $O(N_t\log^2N_t)$ operations. These reductions enable fast, practical simulations of non-Markovian dynamics involving multiphonon and multiphoton processes in open quantum systems.

quant-pharXiv:2609.18806•Sep 17, 2026•0 citations

Exact calculation of the time convolutionless master equation generator: Application to the nonequilibrium resonant level model

Lyran Kidon, Eli Y. Wilner, Eran Rabani

The generalized quantum master equation provides a powerful tool to describe the dynamics in quantum impurity models driven away from equilibrium. Two complementary approaches, one based on Nakajima-Zwanzig-Mori time-convolution (TC) and the other on the Tokuyama-Mori time-convolutionless (TCL) formulations provide a starting point to describe the time-evolution of the reduced density matrix. A key in both approaches is to obtain the so called "memory kernel" or "generator," going beyond second or fourth order perturbation techniques. While numerically converged techniques are available for the TC memory kernel, the canonical approach to obtain the TCL generator is based on inverting a super-operator in the full Hilbert space, which is difficult to perform and thus, nearly all applications of the TCL approach rely on a perturbative scheme of some sort. Here, the TCL generator is expressed using a reduced system propagator which can be obtained from system observables alone and requires the calculation of super-operators and their inverse in the reduced Hilbert space rather than the full one. This makes the formulation amenable to quantum impurity solvers or to diagrammatic techniques, such as the nonequilibrium Green's function. We implement the TCL approach for the resonant level model driven away from equilibrium and compare the time scales for the decay of the generator with that of the memory kernel in the TC approach. Furthermore, the effects of temperature, source-drain bias, and gate potential on the TCL/TC generators are discussed.

arXiv:1510.04369•The Journal of Chemical Physics•Dec 21, 2015•52 citations
The non-Markovian dynamics of an open quantum system can be rigorously derived using the Feynman–Vernon influence functional approach. Although this formalism is exact, practical numerical implementations often require compromises. The time-convolutionless (TCL) master equation offers an exact framework, yet its application typically relies on a perturbative expansion of both the time forward and time backward state propagators. Due to the significant computational effort involved—and the scarcity of analytical solutions for most open quantum systems—the fourth-order perturbative TCL generator (TCL4) has only been benchmarked on a limited range of systems and parameter spaces. Recent advancements, however, have made the computation of TCL4 faster and more accessible. In this paper, we benchmark the TCL4 master equation against numerically exact methods for the biased spin-boson model. We focus on the regime near critical bath coupling where perturbative master equations are expected to become inaccurate. Our findings reveal that the TCL4 approach is most reliable at low temperatures and more efficient than the numerical exact methods. This study aims to delineate the conditions under which the TCL4 perturbative master equation enhances the accuracy of the TCL2.
The spin-boson model is a widely used model for understanding the properties of a two-level open quantum system. Accurately describing its dynamics often requires going beyond the weak system-environment coupling approximation. However, calculating the higher-order generators of such a dynamics, with a system-environment coupling that is not too weak, has been known to be challenging, both numerically and analytically. This work presents the analytical derivation of the complete fourth-order time-convolutionless (TCL) generator for a generic spin-boson model, accurate up to 4th order in the system-environment coupling parameter, under the assumption that the environmental spectral density is an odd function of frequency. In the case of a semiconductor double-quantum-dot system, our results reveal corrections to the dynamics that may become physically significant in some parameter regimes. Furthermore, we report that the widely used second-order TCL master equation tends to overestimate the non-Markovianity of a dynamics over a large parameter regime. Within the regime of its applicability, our results provide a computational advantage over numerically exact techniques. The accuracy of the fourth-order TCL generator is rigorously benchmarked against specialized analytical calculations done for the Ohmic spectral density with Drude cutoff and against the numerically exact Hierarchical Equations of Motion technique.
The time-convolutionless (TCL) quantum master equation provides a powerful tool to simulate reduced dynamics of a quantum system coupled to a bath. The key quantity in the TCL master equation is the so-called kernel or generator, which describes effects of the bath degrees of freedom. Since the exact TCL generators are usually hard to calculate analytically, most applications of the TCL generalized master equation have relied on approximate generators using second and fourth order perturbative expansions. By using the hierarchical equation of motion (HEOM) and extended HEOM methods, we present a new approach to calculating the exact TCL generator and its high order perturbative expansions. The new approach is applied to the spin-boson model with different sets of parameters, to investigate the convergence of the high order expansions of the TCL generator. We also discuss circumstances where the exact TCL generator becomes singular for the spin-boson model, and a model of excitation energy transfer in the Fenna-Matthews-Olson complex.
We explore the performance of the time-convolutionless (TCL) projection operator technique using the Fano-Anderson model as a test case. Comparing the exact TCL master equation with an expansion in powers of the strength of the system-environment coupling, we analyze the transient dynamics as well as the steady-state behavior. For a Lorentzian spectral density we demonstrate that the dimensionless expansion parameter corresponds to the ratio of the environmental correlation time to the relaxation time of the system, and we derive the convergence radius for the TCL expansion, which is seen to depend on the ratio of detuning and width of the spectral density. We further study the quantum non-Markovianity of the model based on the evolution of the Bures distance between quantum states and how it is represented by the second and fourth order of the expansion. Our results highlight both the strengths and the limitations of the TCL formalism in capturing key features of open quantum systems and, in particular, the challenges of accurately describing strongly coupled systems and non-Markovian dynamics.
Steady-state observables, such as occupation numbers and currents, are crucial experimental signatures in open quantum systems.The time-convolutionless (TCL) master equation, which is both exact and time-local, is an ideal candidate for the perturbative computation of such observables.We develop a diagrammatic approach to evaluate the steady-state TCL generator based on operators rather than superoperators.We obtain the steady-state occupation numbers, extend our formulation to the calculation of currents, and provide a simple physical interpretation of the diagrams.We further benchmark our method on a single noninteracting level coupled to Fermi reservoirs, where we recover the exact expansion to next-to-leading order.The low number of diagrams appearing in our formulation makes the extension to higher orders accessible.Combined, these properties make the steady-state time-convolutionless master equation an effective tool for the calculation of steady-state properties in open quantum systems.
Higher-order time-convolutionless (TCL) master equations provide systematic corrections for non-Markovian dynamics in open quantum systems, but their time-ordered integrals become increasingly costly to evaluate. Direct quadrature of the Hadamard-reduced sixth-order (TCL6) generator at all NtN_tNt​ sampled times requires O(Nt3)O(N_t^3)O(Nt3​) operations at fixed system dimension. We derive an exact reduction for a finite-dimensional system coupled through one Hermitian operator to a stationary centered Gaussian bath. By separating operator coefficients from scalar time-dependent spectral densities, the complete TCL6 time series is reduced to cumulative sums and first-order recurrences, one-dimensional causal convolutions, and an exact dyadic recursion for interlocked sums. With the algorithm, evaluating the complete TCL6 time series requires O(Ntlog⁡2Nt)O(N_t\log^2N_t)O(Nt​log2Nt​) operations. These reductions enable fast, practical simulations of non-Markovian dynamics involving multiphonon and multiphoton processes in open quantum systems.
The generalized quantum master equation provides a powerful tool to describe the dynamics in quantum impurity models driven away from equilibrium. Two complementary approaches, one based on Nakajima-Zwanzig-Mori time-convolution (TC) and the other on the Tokuyama-Mori time-convolutionless (TCL) formulations provide a starting point to describe the time-evolution of the reduced density matrix. A key in both approaches is to obtain the so called "memory kernel" or "generator," going beyond second or fourth order perturbation techniques. While numerically converged techniques are available for the TC memory kernel, the canonical approach to obtain the TCL generator is based on inverting a super-operator in the full Hilbert space, which is difficult to perform and thus, nearly all applications of the TCL approach rely on a perturbative scheme of some sort. Here, the TCL generator is expressed using a reduced system propagator which can be obtained from system observables alone and requires the calculation of super-operators and their inverse in the reduced Hilbert space rather than the full one. This makes the formulation amenable to quantum impurity solvers or to diagrammatic techniques, such as the nonequilibrium Green's function. We implement the TCL approach for the resonant level model driven away from equilibrium and compare the time scales for the decay of the generator with that of the memory kernel in the TC approach. Furthermore, the effects of temperature, source-drain bias, and gate potential on the TCL/TC generators are discussed.