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:2609.18806

Fast Evaluation of the Sixth-Order Time-Convolutionless Master-Equation Generator and Beyond

Jiahao Chen, Sirui Chen, Dragomir Davidovic

quant-phphysics.comp-phSubmitted Sep 16, 2026• Announced Sep 17, 2026• 0 citations
v1 · Sep 16, 2026v2 · Oct 20, 2026
arXivHTMLPDF

Abstract

Direct quadrature of the Hadamard-reduced sixth-order time-convolutionless (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 fixed operator coefficients from scalar bath kernels, 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 histories. With the algorithm, evaluating the complete TCL6 time series requires $O(N_t\log^2 N_t)$ operations at fixed system dimension. In addition, we show that TCL$2n$ can be evaluated with $O(N_t\log^{n-1} N_t)$ complexity. A fixed finite Matsubara expansion of the bath correlation function permits $O(N_t\log N_t)$ evaluation at any fixed TCL order. These reductions enable fast long-time simulations of non-Markovian open quantum systems within the regime of validity of the TCL expansion.

Related papers

Sixth-order time-convolutionless master equation and beyond: Late-time resummations, two types of divergences, and the limits of validity

L. Lampert, Srikar Gadamsetty, Shantanu Chaudhary +4 more

Perturbative master equations are essential for modeling open quantum systems but often exhibit late-time divergences when environmental correlations decay algebraically. In this work, we analyze the time-convolutionless (TCL) master equation expanded to order ${\mathrm{TCL}}_{2n}$ and demonstrate that, while van Kampen's cumulants suppress early-time secular growth, they ultimately diverge at long times. To overcome this, we introduce a resummation technique based on the Hadamard trick, which incorporates time integrals directly into the bath's spectral density via elementwise multiplication. This approach establishes a maximum expansion order, ${n}_{\text{max}}=\ensuremath{\lceil}s+1\ensuremath{\rceil}$, and defines a precision limit of the asymptotic states of $O({\ensuremath{\lambda}}^{2\ensuremath{\lceil}s\ensuremath{\rceil}})$, where $s$ is the power-law exponent and $\ensuremath{\lambda}$ is the weak-coupling constant. The resummed master equation features renormalized Bohr frequencies that capture decoherence and spectral overlap effects. In the unbiased spin-boson model, this results in an inflation of the generator at a temperature-independent rate equal to the decoherence rate ${\ensuremath{\nu}}_{2}$ and a finite validity time ${t}_{L}\ensuremath{\approx}(s+1)/{\ensuremath{\nu}}_{2}$. For exponentially decaying correlations, the method recovers a proper Markovian limit below a critical coupling threshold.

arXiv:2406.11088•Physical Review A•Apr 22, 2025•6 citations

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

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

Jiahao Chen, Elyana Crowder, Lian Xiang +1 more

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.

arXiv:2501.04192•APL Quantum•Apr 14, 2025•4 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

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

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
Direct quadrature of the Hadamard-reduced sixth-order time-convolutionless (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 fixed operator coefficients from scalar bath kernels, 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 histories. With the algorithm, evaluating the complete TCL6 time series requires O(Ntlog⁡2Nt)O(N_t\log^2 N_t)O(Nt​log2Nt​) operations at fixed system dimension. In addition, we show that TCL2n2n2n can be evaluated with O(Ntlog⁡n−1Nt)O(N_t\log^{n-1} N_t)O(Nt​logn−1Nt​) complexity. A fixed finite Matsubara expansion of the bath correlation function permits O(Ntlog⁡Nt)O(N_t\log N_t)O(Nt​logNt​) evaluation at any fixed TCL order. These reductions enable fast long-time simulations of non-Markovian open quantum systems within the regime of validity of the TCL expansion.
Perturbative master equations are essential for modeling open quantum systems but often exhibit late-time divergences when environmental correlations decay algebraically. In this work, we analyze the time-convolutionless (TCL) master equation expanded to order TCL2n{\mathrm{TCL}}_{2n}TCL2n​ and demonstrate that, while van Kampen's cumulants suppress early-time secular growth, they ultimately diverge at long times. To overcome this, we introduce a resummation technique based on the Hadamard trick, which incorporates time integrals directly into the bath's spectral density via elementwise multiplication. This approach establishes a maximum expansion order, ${n}_{\text{max}}=\ensuremath{\lceil}s+1\ensuremath{\rceil}$, and defines a precision limit of the asymptotic states of $O({\ensuremath{\lambda}}^{2\ensuremath{\lceil}s\ensuremath{\rceil}})$, where sss is the power-law exponent and $\ensuremath{\lambda}$ is the weak-coupling constant. The resummed master equation features renormalized Bohr frequencies that capture decoherence and spectral overlap effects. In the unbiased spin-boson model, this results in an inflation of the generator at a temperature-independent rate equal to the decoherence rate ${\ensuremath{\nu}}_{2}$ and a finite validity time ${t}_{L}\ensuremath{\approx}(s+1)/{\ensuremath{\nu}}_{2}$. For exponentially decaying correlations, the method recovers a proper Markovian limit below a critical coupling threshold.
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 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.
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.
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.
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.