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

Universal structure of computing moments for exact quantum dynamics: Application to arbitrary system–bath couplings

Ruihao Bi, Wei Liu, Wenjie Dou

Published Jun 11, 2025• 6 citations
arXivHTMLPDF

Abstract

We introduce a general procedure for computing higher-order moments of correlation functions in open quantum systems, extending the scope of our recent work on Memory Kernel Coupling Theory (MKCT) [Liu et al., arXiv:2407.01923 (2024)]. This approach is demonstrated for arbitrary system-bath coupling that can be expressed as a polynomial, HSB=V̂(α0+α1q̂+α2q̂2+⋯) , where we show that the recursive commutators of a system operator obey a universal hierarchy. Exploiting this structure, the higher-order moments are obtained by evaluating the expectation values of the system and bath operators separately, with bath expectation values derived from the derivatives of a generating function. We further apply MKCT to compute the dipole autocorrelation function for the spin-boson model with both linear and quadratic coupling, achieving agreement with the hierarchical equations of motion approach. Our findings suggest a promising path toward accurate dynamics for complex open quantum systems.

DOI
10.1063/5.0273707

Related papers

From higher-order moments to time correlation functions in strongly correlated systems: A DMRG-based memory kernel coupling theory

Yunhao Liu, Wenjie Dou

We introduce a hybrid approach for computing dynamical observables in strongly correlated systems using higher-order moments. This method integrates memory kernel coupling theory (MKCT) with the density matrix renormalization group (DMRG), extending our recent work on MKCT [Liu et al., Phys. Rev. Lett. 135, 148001(2025)] to strongly correlated systems. The method establishes that correlation functions can be derived from the moments. Within our framework, operators and wavefunctions are represented as matrix product operators (MPOs) and matrix product states (MPSs), respectively. Crucially, the repeated application of the Liouville operator is achieved through an iterative procedure analogous to the DMRG algorithm itself. We demonstrate the effectiveness and efficiency of MKCT-DMRG by computing the spectral function of the Hubbard model. Furthermore, we successfully apply the method to compute the electronic friction in the Hubbard-Holstein model. In all cases, the results show excellent agreement with time-dependent DMRG (TD-DMRG) benchmarks. The advantage of MKCT-DMRG over TD-DMRG is computational efficiency, which avoids expensive real-time propagation in TD-DMRG. These findings establish MKCT-DMRG as a promising and accurate framework for simulating challenging dynamical properties in strongly correlated quantum systems.

arXiv:2509.13140•The Journal of Chemical Physics•Dec 1, 2025•3 citations

Memory Kernel Coupling Theory: Obtaining Time Correlation Function from Higher-Order Moments

Wei Liu, Yu Su, Yao Wang +1 more

Dynamical observables can often be described by time correlation functions (TCFs). However, efficiently calculating TCFs for complex quantum systems is a significant challenge, which generally requires solving the full dynamics of the systems. This Letter presents the memory kernel coupling theory (MKCT), a general formalism for evaluating TCFs. The MKCT builds upon Mori's memory kernel formalism for TCFs. Our theory further decomposes the memory kernel into auxiliary kernels. Rapid decay of auxiliary kernels allows us to truncate the coupled equations of motion with high accuracy. Notably, only higher-order moments are sufficient as the input for obtaining TCFs. While this formalism is general, we carry out the numerical demonstration for typical open quantum systems-the spin-boson model and the single impurity Anderson model.

Physical Review Letters•Aug 29, 2025•7 citations

Generalized quantum master equation from memory kernel coupling theory

Rui-Hao Bi, Wei Liu, Wenjie Dou

The generalized quantum master equation provides a powerful framework for non-Markovian dynamics of open quantum systems. However, the accurate and efficient evaluation of the memory kernel remains a challenge. In this work, we introduce a comprehensive tensorial extension to the Memory Kernel Coupling Theory (MKCT) to overcome this bottleneck. By elevating the original scalar formalism to a tensorial framework, the extended MKCT enables the calculation of general expectation values and cross-correlation functions. We demonstrate the numerical accuracy and efficiency of this method across multiple benchmark systems: capturing transient populations and coherences in the spin-boson model, resolving the excitonic absorption spectrum of the Fenna-Matthews-Olson complex, and simulating charge mobility in one-dimensional lattice models. These successful applications establish the tensorial MKCT as a highly efficient tool for investigating complex dynamics in open quantum systems.

physics.chem-phCross-listedarXiv:2603.01458•Journal of Chemical Theory and Computation•Mar 3, 2026•1 citation

Memory Kernel Coupling Theory: Obtain Time Correlation Function from Higher-order Moments

Wei Liu, Yu Su, Yao Wang +1 more

Dynamical observables can often be described by time correlation functions (TCFs). However, efficiently calculating TCFs for complex quantum systems is a significant challenge, which generally requires solving the full dynamics of the systems. This Letter presents the memory kernel coupling theory (MKCT), a general formalism for evaluating TCFs. The MKCT builds upon Mori's memory kernel formalism for TCFs. Our theory further decomposes the memory kernel into auxiliary kernels. Rapid decay of auxiliary kernels allows us to truncate the coupled equations of motion with high accuracy. Notably, only higher-order moments are sufficient as the input for obtaining TCFs. While this formalism is general, we carry out the numerical demonstration for a typical open quantum system--the spin-boson model.

arXiv:2407.01923•Jul 2, 2024•1 citation

An Exact Operator Formulation of Statistical Quasiparticle Theory for Nonlinear System-Bath Coupling

Yu Su, Yao Wang

We develop an operator formulation of statistical quasiparticle theory for open quantum systems with nonlinear system-bath coupling, extending beyond the linear coupling setting of the standard Gaussian influence functional formulation. For Gaussian baths, the statistical quasiparticles, termed dissipatons, are explicitly constructed as collective operators in the thermofield bath space. Our construction provides a microscopic operator basis for the dissipaton algebra and generates hierarchical equations of motion directly from the Liouville equation. Nonlinear interactions are incorporated through a recursive generalized Wick's theorem, enabling automatic generation of hierarchical equations for polynomial couplings. The resulting framework connects explicit bath operator representations to nonperturbative, non-Markovian dynamics and unifies the treatment of linear and nonlinear system-bath interactions.

quant-pharXiv:2609.24603•Sep 22, 2026•0 citations

A unified stochastic formulation of dissipative quantum dynamics. I. Generalized hierarchical equations

Chang‐Yu Hsieh, Jianshu Cao

We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.

arXiv:1701.05709•The Journal of Chemical Physics•Jan 2, 2018•85 citations
We introduce a general procedure for computing higher-order moments of correlation functions in open quantum systems, extending the scope of our recent work on Memory Kernel Coupling Theory (MKCT) [Liu et al., arXiv:2407.01923 (2024)]. This approach is demonstrated for arbitrary system-bath coupling that can be expressed as a polynomial, HSB=V̂(α0+α1q̂+α2q̂2+⋯) , where we show that the recursive commutators of a system operator obey a universal hierarchy. Exploiting this structure, the higher-order moments are obtained by evaluating the expectation values of the system and bath operators separately, with bath expectation values derived from the derivatives of a generating function. We further apply MKCT to compute the dipole autocorrelation function for the spin-boson model with both linear and quadratic coupling, achieving agreement with the hierarchical equations of motion approach. Our findings suggest a promising path toward accurate dynamics for complex open quantum systems.
We introduce a hybrid approach for computing dynamical observables in strongly correlated systems using higher-order moments. This method integrates memory kernel coupling theory (MKCT) with the density matrix renormalization group (DMRG), extending our recent work on MKCT [Liu et al., Phys. Rev. Lett. 135, 148001(2025)] to strongly correlated systems. The method establishes that correlation functions can be derived from the moments. Within our framework, operators and wavefunctions are represented as matrix product operators (MPOs) and matrix product states (MPSs), respectively. Crucially, the repeated application of the Liouville operator is achieved through an iterative procedure analogous to the DMRG algorithm itself. We demonstrate the effectiveness and efficiency of MKCT-DMRG by computing the spectral function of the Hubbard model. Furthermore, we successfully apply the method to compute the electronic friction in the Hubbard-Holstein model. In all cases, the results show excellent agreement with time-dependent DMRG (TD-DMRG) benchmarks. The advantage of MKCT-DMRG over TD-DMRG is computational efficiency, which avoids expensive real-time propagation in TD-DMRG. These findings establish MKCT-DMRG as a promising and accurate framework for simulating challenging dynamical properties in strongly correlated quantum systems.
Dynamical observables can often be described by time correlation functions (TCFs). However, efficiently calculating TCFs for complex quantum systems is a significant challenge, which generally requires solving the full dynamics of the systems. This Letter presents the memory kernel coupling theory (MKCT), a general formalism for evaluating TCFs. The MKCT builds upon Mori's memory kernel formalism for TCFs. Our theory further decomposes the memory kernel into auxiliary kernels. Rapid decay of auxiliary kernels allows us to truncate the coupled equations of motion with high accuracy. Notably, only higher-order moments are sufficient as the input for obtaining TCFs. While this formalism is general, we carry out the numerical demonstration for typical open quantum systems-the spin-boson model and the single impurity Anderson model.
The generalized quantum master equation provides a powerful framework for non-Markovian dynamics of open quantum systems. However, the accurate and efficient evaluation of the memory kernel remains a challenge. In this work, we introduce a comprehensive tensorial extension to the Memory Kernel Coupling Theory (MKCT) to overcome this bottleneck. By elevating the original scalar formalism to a tensorial framework, the extended MKCT enables the calculation of general expectation values and cross-correlation functions. We demonstrate the numerical accuracy and efficiency of this method across multiple benchmark systems: capturing transient populations and coherences in the spin-boson model, resolving the excitonic absorption spectrum of the Fenna-Matthews-Olson complex, and simulating charge mobility in one-dimensional lattice models. These successful applications establish the tensorial MKCT as a highly efficient tool for investigating complex dynamics in open quantum systems.
Dynamical observables can often be described by time correlation functions (TCFs). However, efficiently calculating TCFs for complex quantum systems is a significant challenge, which generally requires solving the full dynamics of the systems. This Letter presents the memory kernel coupling theory (MKCT), a general formalism for evaluating TCFs. The MKCT builds upon Mori's memory kernel formalism for TCFs. Our theory further decomposes the memory kernel into auxiliary kernels. Rapid decay of auxiliary kernels allows us to truncate the coupled equations of motion with high accuracy. Notably, only higher-order moments are sufficient as the input for obtaining TCFs. While this formalism is general, we carry out the numerical demonstration for a typical open quantum system--the spin-boson model.
We develop an operator formulation of statistical quasiparticle theory for open quantum systems with nonlinear system-bath coupling, extending beyond the linear coupling setting of the standard Gaussian influence functional formulation. For Gaussian baths, the statistical quasiparticles, termed dissipatons, are explicitly constructed as collective operators in the thermofield bath space. Our construction provides a microscopic operator basis for the dissipaton algebra and generates hierarchical equations of motion directly from the Liouville equation. Nonlinear interactions are incorporated through a recursive generalized Wick's theorem, enabling automatic generation of hierarchical equations for polynomial couplings. The resulting framework connects explicit bath operator representations to nonperturbative, non-Markovian dynamics and unifies the treatment of linear and nonlinear system-bath interactions.
We extend a standard stochastic theory to study open quantum systems coupled to a generic quantum environment. We exemplify the general framework by studying a two-level quantum system coupled bilinearly to the three fundamental classes of non-interacting particles: bosons, fermions, and spins. In this unified stochastic approach, the generalized stochastic Liouville equation (SLE) formally captures the exact quantum dissipations when noise variables with appropriate statistics for different bath models are applied. Anharmonic effects of a non-Gaussian bath are precisely encoded in the bath multi-time correlation functions that noise variables have to satisfy. Starting from the SLE, we devise a family of generalized hierarchical equations by averaging out the noise variables and expand bath multi-time correlation functions in a complete basis of orthonormal functions. The general hierarchical equations constitute systems of linear equations that provide numerically exact simulations of quantum dynamics. For bosonic bath models, our general hierarchical equation of motion reduces exactly to an extended version of hierarchical equation of motion which allows efficient simulation for arbitrary spectral densities and temperature regimes. Similar efficiency and flexibility can be achieved for the fermionic bath models within our formalism. The spin bath models can be simulated with two complementary approaches in the present formalism. (I) They can be viewed as an example of non-Gaussian bath models and be directly handled with the general hierarchical equation approach given their multi-time correlation functions. (II) Alternatively, each bath spin can be first mapped onto a pair of fermions and be treated as fermionic environments within the present formalism.