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:cond-mat/0612649 · Physical Review E

Dynamics of quantum dissipation systems interacting with bosonic canonical bath: Hierarchical equations of motion approach

Rui‐Xue Xu, YiJing Yan

Published Mar 15, 2007• 177 citations
arXivHTMLPDF

Abstract

A nonperturbative theory is developed, aiming at an exact and efficient evaluation of a general quantum system interacting with arbitrary bath environment at any temperature and in the presence of arbitrary time-dependent external fields. An exact hierarchical equations of motion formalism is constructed on the basis of a calculus-on-path-integral algorithm, via the auxiliary influence generating functionals related to the interaction bath correlation functions in a parametrization expansion form. The corresponding continued-fraction Green's functions formalism for quantum dissipation is also presented. Proposed further is the principle of residue correction, not just for truncating the infinite hierarchy, but also for incorporating the small residue dissipation that may arise from the practical difference between the true and parametrized bath correlation functions. The final residue-corrected hierarchical equations of motion can therefore be used practically for the evaluation of arbitrary dissipative quantum systems.

DOI
10.1103/physreve.75.031107

Related papers

Dynamics of quantum dissipation systems interacting with Fermion and Boson grand canonical bath ensembles: Path integral versus equations of motion formalism

Jinshuang Jin, Sven Welack, JunYan Luo +4 more

An exact hierarchical equations of motion formalism for a quantum dissipation system in a grand canonical bath ensemble surrounding is constructed, on the basis of the calculus-on-path-integral algorithm, together with the parametrization of arbitrary non-Markovin bath that satisfies fluctuation-dissipation theorem. The influence functionals for both the Fermion or Boson bath interaction are found to be of the same path-integral expression as the canonical bath, assuming they all satisfy the Gaussian statistics. However, the equation of motion formalism are different, due to the fluctuation-dissipation theories that are distinct and used explicitly. The implications of the present work to quantum transport through molecular wires and electron transfer in complex molecular systems are discussed.

•Dec 6, 2006•0 citations

Dynamics of dissipative electronic systems and quantum transport: Hierarchical equations of motion approach

YiJing Yan

A quantum dissipation theory is formulated in terms of hierarchically coupled equations of motion for an arbitrary fermionic system coupled with grand canonical fermionic bath ensembles.1 The theoretical construction starts with the second-quantization influence functional in path integral formalism, in which the fermionic creation and annihilation operators are represented by Grassmann variables. Temporal derivatives on influence functionals are then performed in a hierarchical manner, on the basis of calculus-on-path-integral algorithm. Both the multiple-frequency-dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting formalism is in principle exact, applicable to interacting systems, with arbitrary time-dependent external fields. It renders an exact tool to evaluate various transient and stationary quantum transport properties of many-electron systems. At the second-tier truncation level the present theory recovers the real-time diagrammatic formalism developed by Schön and coworkers.2 For a single-particle system, the hierarchical formalism terminates at the second tier exactly, and the Landuer-Büttiker's transport current expression is readily recovered. Numerical studies will be presented to highlight the richness of transient current through both interacting and noninteracting model systems.

OSA Workshop on Entanglement and Quantum Decoherence•Jan 1, 2008•0 citations

Exact dynamics of dissipative electronic systems and quantum transport: Hierarchical equations of motion approach

Jinshuang Jin, Xiao Hong Zheng, YiJing Yan

A generalized quantum master equation theory that governs the exact, nonperturbative quantum dissipation and quantum transport is formulated in terms of hierarchically coupled equations of motion for an arbitrary electronic system in contact with electrodes under either a stationary or a nonstationary electrochemical potential bias. The theoretical construction starts with the influence functional in path integral, in which the electron creation and annihilation operators are Grassmann variables. Time derivatives on the influence functionals are then performed in a hierarchical manner. Both the multiple-frequency dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting hierarchical equations of motion formalism is in principle exact and applicable to arbitrary electronic systems, including Coulomb interactions, under the influence of arbitrary time-dependent applied bias voltage and external fields. Both the conventional quantum master equation and the real-time diagrammatic formalism of Schon and co-workers can be readily obtained at well defined limits of the present theory. We also show that for a noninteracting electron system, the present hierarchical equations of motion formalism terminates at the second tier exactly, and the Landuer-Buttiker transport current expression is recovered. The present theory renders an exact and numerically tractable tool to evaluate various transient and stationary quantum transport properties of many-electron systems, together with the involving nonperturbative dissipative dynamics.

arXiv:0710.5367•The Journal of Chemical Physics•Jun 17, 2008•452 citations

Exact quantum master equation via the calculus on path integrals

Rui‐Xue Xu, Ping Cui, Xin-Qi Li +2 more

An exact quantum master equation formalism is constructed for the efficient evaluation of quantum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. A novel truncation scheme is further proposed and compared with other approaches to close the resulting hierarchically coupled equations of motion. The interplay between system-bath interaction strength, non-Markovian property, and required level of hierarchy is also demonstrated with the aid of simple spin-boson systems.

arXiv:quant-ph/0409007•The Journal of Chemical Physics•Jan 10, 2005•276 citations

Quantum-Classical Hierarchical Equations of Motion

Amartya Bose

We develop a quantum-classical hierarchical equations of motion (QC-HEOM) approach for simulating non-Markovian open quantum systems. The method combines the ensemble-averaged classical path reference of the quantum-classical path integral formalism with a hierarchy of auxiliary quantum influence functionals. By incorporating thermal fluctuations through an ensemble average over reference trajectories, the hierarchy is required to represent only the residual quantum memory associated with the imaginary part of the bath response function. Consequently, unlike conventional hierarchical equations of motion, QC-HEOM does not require Matsubara or Padé expansions of the thermal kernel and exhibits only weak temperature dependence of the hierarchy size. Furthermore, because thermal fluctuations are supplied through reference classical trajectories, the framework naturally extends beyond harmonic baths and enables the incorporation of anharmonic and molecular environments through externally generated trajectories. We derive the formalism and demonstrate its exactness for a harmonic bath. Applications to an asymmetric spin-boson model and the seven-site Fenna--Matthews--Olson complex illustrate the accuracy of QC-HEOM. It reproduces benchmark quasi-adiabatic path integral and hierarchical equations of motion results while requiring substantially fewer auxiliary objects, particularly at low temperatures. These results establish QC-HEOM as an efficient framework for treating residual quantum memory in quantum-classical descriptions of open-system dynamics. The separation of thermal fluctuations from residual quantum memory through the use of Wigner trajectories provides an approximate route toward hierarchical treatments of complex anharmonic environments that are inaccessible to conventional HEOM approaches.

quant-pharXiv:2606.14363•Jun 15, 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
A nonperturbative theory is developed, aiming at an exact and efficient evaluation of a general quantum system interacting with arbitrary bath environment at any temperature and in the presence of arbitrary time-dependent external fields. An exact hierarchical equations of motion formalism is constructed on the basis of a calculus-on-path-integral algorithm, via the auxiliary influence generating functionals related to the interaction bath correlation functions in a parametrization expansion form. The corresponding continued-fraction Green's functions formalism for quantum dissipation is also presented. Proposed further is the principle of residue correction, not just for truncating the infinite hierarchy, but also for incorporating the small residue dissipation that may arise from the practical difference between the true and parametrized bath correlation functions. The final residue-corrected hierarchical equations of motion can therefore be used practically for the evaluation of arbitrary dissipative quantum systems.
An exact hierarchical equations of motion formalism for a quantum dissipation system in a grand canonical bath ensemble surrounding is constructed, on the basis of the calculus-on-path-integral algorithm, together with the parametrization of arbitrary non-Markovin bath that satisfies fluctuation-dissipation theorem. The influence functionals for both the Fermion or Boson bath interaction are found to be of the same path-integral expression as the canonical bath, assuming they all satisfy the Gaussian statistics. However, the equation of motion formalism are different, due to the fluctuation-dissipation theories that are distinct and used explicitly. The implications of the present work to quantum transport through molecular wires and electron transfer in complex molecular systems are discussed.
A quantum dissipation theory is formulated in terms of hierarchically coupled equations of motion for an arbitrary fermionic system coupled with grand canonical fermionic bath ensembles.1 The theoretical construction starts with the second-quantization influence functional in path integral formalism, in which the fermionic creation and annihilation operators are represented by Grassmann variables. Temporal derivatives on influence functionals are then performed in a hierarchical manner, on the basis of calculus-on-path-integral algorithm. Both the multiple-frequency-dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting formalism is in principle exact, applicable to interacting systems, with arbitrary time-dependent external fields. It renders an exact tool to evaluate various transient and stationary quantum transport properties of many-electron systems. At the second-tier truncation level the present theory recovers the real-time diagrammatic formalism developed by Schön and coworkers.2 For a single-particle system, the hierarchical formalism terminates at the second tier exactly, and the Landuer-Büttiker's transport current expression is readily recovered. Numerical studies will be presented to highlight the richness of transient current through both interacting and noninteracting model systems.
A generalized quantum master equation theory that governs the exact, nonperturbative quantum dissipation and quantum transport is formulated in terms of hierarchically coupled equations of motion for an arbitrary electronic system in contact with electrodes under either a stationary or a nonstationary electrochemical potential bias. The theoretical construction starts with the influence functional in path integral, in which the electron creation and annihilation operators are Grassmann variables. Time derivatives on the influence functionals are then performed in a hierarchical manner. Both the multiple-frequency dispersion and the non-Markovian reservoir parametrization schemes are considered for the desired hierarchy construction. The resulting hierarchical equations of motion formalism is in principle exact and applicable to arbitrary electronic systems, including Coulomb interactions, under the influence of arbitrary time-dependent applied bias voltage and external fields. Both the conventional quantum master equation and the real-time diagrammatic formalism of Schon and co-workers can be readily obtained at well defined limits of the present theory. We also show that for a noninteracting electron system, the present hierarchical equations of motion formalism terminates at the second tier exactly, and the Landuer-Buttiker transport current expression is recovered. The present theory renders an exact and numerically tractable tool to evaluate various transient and stationary quantum transport properties of many-electron systems, together with the involving nonperturbative dissipative dynamics.
An exact quantum master equation formalism is constructed for the efficient evaluation of quantum non-Markovian dissipation beyond the weak system-bath interaction regime in the presence of time-dependent external field. A novel truncation scheme is further proposed and compared with other approaches to close the resulting hierarchically coupled equations of motion. The interplay between system-bath interaction strength, non-Markovian property, and required level of hierarchy is also demonstrated with the aid of simple spin-boson systems.
We develop a quantum-classical hierarchical equations of motion (QC-HEOM) approach for simulating non-Markovian open quantum systems. The method combines the ensemble-averaged classical path reference of the quantum-classical path integral formalism with a hierarchy of auxiliary quantum influence functionals. By incorporating thermal fluctuations through an ensemble average over reference trajectories, the hierarchy is required to represent only the residual quantum memory associated with the imaginary part of the bath response function. Consequently, unlike conventional hierarchical equations of motion, QC-HEOM does not require Matsubara or Padé expansions of the thermal kernel and exhibits only weak temperature dependence of the hierarchy size. Furthermore, because thermal fluctuations are supplied through reference classical trajectories, the framework naturally extends beyond harmonic baths and enables the incorporation of anharmonic and molecular environments through externally generated trajectories. We derive the formalism and demonstrate its exactness for a harmonic bath. Applications to an asymmetric spin-boson model and the seven-site Fenna--Matthews--Olson complex illustrate the accuracy of QC-HEOM. It reproduces benchmark quasi-adiabatic path integral and hierarchical equations of motion results while requiring substantially fewer auxiliary objects, particularly at low temperatures. These results establish QC-HEOM as an efficient framework for treating residual quantum memory in quantum-classical descriptions of open-system dynamics. The separation of thermal fluctuations from residual quantum memory through the use of Wigner trajectories provides an approximate route toward hierarchical treatments of complex anharmonic environments that are inaccessible to conventional HEOM approaches.
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.