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:1805.09723 · Physical Review A

Hierarchical Schrödinger equations of motion for open quantum dynamics

Kiyoto Nakamura, Yoshitaka Tanimura

Published Jul 6, 2018• 45 citations
arXivHTMLPDF

Abstract

We rigorously investigate the quantum non-Markovian dissipative dynamics of a system coupled to a harmonic-oscillator bath by deriving hierarchical Schr\"odinger equations of motion (HSEOM) and studying their dynamics. The HSEOM are the equations for wave functions derived on the basis of the Feynman-Vernon influence functional formalism for the density operator, $\ensuremath{\langle}q|\ensuremath{\rho}(t)|{q}^{\ensuremath{'}}\ensuremath{\rangle}$, where $\ensuremath{\langle}q|$ and $|{q}^{\ensuremath{'}}\ensuremath{\rangle}$ are the left- and right-hand elements. The time evolution of $\ensuremath{\langle}q|$ is computed from time 0 to $t$, and, subsequently, the time evolution of $|{q}^{\ensuremath{'}}\ensuremath{\rangle}$ is computed from time $t$ to 0 along a contour in the complex time plane. By appropriately choosing functions for the bath correlation function and the spectral density, we can take advantage of an HSEOM method to carry out simulations without the need for a great amount of computational memory. As a demonstration, quantum annealing simulation for a ferromagnetic $p$-spin model is studied.

DOI
10.1103/physreva.98.012109

Related papers

Hierarchical Equations of Motion for Quantum Dissipation and Quantum Transport

Zheng Xiao, Rui‐Xue Xu, Jian Xu +3 more

In this review we give a comprehensive account of a hierarchical equations of motion (HEOM) approach to the characterization of stationary and dynamic properties of open quantum systems. This approach is rooted at the Feynman-Vernon influence functional path integral formalism, but much more implementable numerically and operationally for the study of various complex molecular dynamics and quantum transport in strongly correlated electronic systems. By construction, HEOM resolves nonperturbatively the combined effects of many-particle interaction, system-bath coupling,and non-Markovian memory. Finally the practicality of HEOM to address physical and chemical problems is exemplified with a model simulation of coherent two-dimensional spectroscopy signals of a biological light-harvesting system and a time-dependent quantum transport system involving dynamic Kondo transition.級聯方程已成為研究量子開放系統的穩態性質和動力學過程的重要方法。本文旨在系統綜述量子耗散和量子輸運的級聯方程組方法的建立、發展以及在理論、算法和應用方面的一些最新進展。級聯方程形式理論的建立以影響泛函路徑積分為基礎,并具有數值上的高效性和應用上的靈活性,可用于研究分子體系的復雜動力學過程以及強關聯電子體系中的量子輸運。其級聯耦合結構以非微擾的方式揭示了多體相互作用、體系-環境耦合、非馬爾可夫記憶等的綜合效應。作為應用示例,我們采用級聯方程模擬了生物光富集體系的二維相干動力學光譜以及含時電子輸運過程中的動態近藤效應。

Huaxue jinzhan•Jun 24, 2012•41 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

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

Reduced quantum dynamics with arbitrary bath spectral densities: Hierarchical equations of motion based on several different bath decomposition schemes

Hao Liu, Lili Zhu, Shuming Bai +1 more

We investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.

The Journal of Chemical Physics•Apr 4, 2014•151 citations

Markovian Embeddings of Non-Markovian Open System Dynamics

Meng Xu, J. T. Stockburger, J. Ankerhold

Embedding non-Markovian open quantum dynamics into an enlarged Markovian space offers a powerful route to nonperturbative simulations, where the dynamics of the extended space can be governed by multiple distinct Markovian equations. We show that these distinct embeddings arise from different unravelings of Gaussian bath self-energies, generating a family of deterministic, time-local equations for the extended system. Using the Brownian-oscillator spectral density as an illustrative example, we clarify the relationships among existing approaches, including the Hierarchical Equations of Motion (HEOM) and the Lindblad--pseudomode formalism, and demonstrate how this framework enables numerically stable and efficient simulations. This work provides both a transparent theoretical foundation for embedding techniques and a flexible platform for developing new methods to simulate non-Markovian quantum dynamics.

quant-pharXiv:2602.21430•The Journal of Chemical Physics•Feb 26, 2026•3 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
We rigorously investigate the quantum non-Markovian dissipative dynamics of a system coupled to a harmonic-oscillator bath by deriving hierarchical Schr"odinger equations of motion (HSEOM) and studying their dynamics. The HSEOM are the equations for wave functions derived on the basis of the Feynman-Vernon influence functional formalism for the density operator, \ensuremath⟨q∣\ensuremathρ(t)∣q\ensuremath′\ensuremath⟩\ensuremath{\langle}q|\ensuremath{\rho}(t)|{q}^{\ensuremath{'}}\ensuremath{\rangle}\ensuremath⟨q∣\ensuremathρ(t)∣q\ensuremath′\ensuremath⟩, where \ensuremath⟨q∣\ensuremath{\langle}q|\ensuremath⟨q∣ and ∣q\ensuremath′\ensuremath⟩|{q}^{\ensuremath{'}}\ensuremath{\rangle}∣q\ensuremath′\ensuremath⟩ are the left- and right-hand elements. The time evolution of \ensuremath⟨q∣\ensuremath{\langle}q|\ensuremath⟨q∣ is computed from time 0 to ttt, and, subsequently, the time evolution of ∣q\ensuremath′\ensuremath⟩|{q}^{\ensuremath{'}}\ensuremath{\rangle}∣q\ensuremath′\ensuremath⟩ is computed from time ttt to 0 along a contour in the complex time plane. By appropriately choosing functions for the bath correlation function and the spectral density, we can take advantage of an HSEOM method to carry out simulations without the need for a great amount of computational memory. As a demonstration, quantum annealing simulation for a ferromagnetic ppp-spin model is studied.
In this review we give a comprehensive account of a hierarchical equations of motion (HEOM) approach to the characterization of stationary and dynamic properties of open quantum systems. This approach is rooted at the Feynman-Vernon influence functional path integral formalism, but much more implementable numerically and operationally for the study of various complex molecular dynamics and quantum transport in strongly correlated electronic systems. By construction, HEOM resolves nonperturbatively the combined effects of many-particle interaction, system-bath coupling,and non-Markovian memory. Finally the practicality of HEOM to address physical and chemical problems is exemplified with a model simulation of coherent two-dimensional spectroscopy signals of a biological light-harvesting system and a time-dependent quantum transport system involving dynamic Kondo transition.級聯方程已成為研究量子開放系統的穩態性質和動力學過程的重要方法。本文旨在系統綜述量子耗散和量子輸運的級聯方程組方法的建立、發展以及在理論、算法和應用方面的一些最新進展。級聯方程形式理論的建立以影響泛函路徑積分為基礎,并具有數值上的高效性和應用上的靈活性,可用于研究分子體系的復雜動力學過程以及強關聯電子體系中的量子輸運。其級聯耦合結構以非微擾的方式揭示了多體相互作用、體系-環境耦合、非馬爾可夫記憶等的綜合效應。作為應用示例,我們采用級聯方程模擬了生物光富集體系的二維相干動力學光譜以及含時電子輸運過程中的動態近藤效應。
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.
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 investigated applications of the hierarchical equation of motion (HEOM) method to perform high order perturbation calculations of reduced quantum dynamics for a harmonic bath with arbitrary spectral densities. Three different schemes are used to decompose the bath spectral density into analytical forms that are suitable to the HEOM treatment: (1) The multiple Lorentzian mode model that can be obtained by numerically fitting the model spectral density. (2) The combined Debye and oscillatory Debye modes model that can be constructed by fitting the corresponding classical bath correlation function. (3) A new method that uses undamped harmonic oscillator modes explicitly in the HEOM formalism. Methods to extract system-bath correlations were investigated for the above bath decomposition schemes. We also show that HEOM in the undamped harmonic oscillator modes can give detailed information on the partial Wigner transform of the total density operator. Theoretical analysis and numerical simulations of the spin-Boson dynamics and the absorption line shape of molecular dimers show that the HEOM formalism for high order perturbations can serve as an important tool in studying the quantum dissipative dynamics in the intermediate coupling regime.
Embedding non-Markovian open quantum dynamics into an enlarged Markovian space offers a powerful route to nonperturbative simulations, where the dynamics of the extended space can be governed by multiple distinct Markovian equations. We show that these distinct embeddings arise from different unravelings of Gaussian bath self-energies, generating a family of deterministic, time-local equations for the extended system. Using the Brownian-oscillator spectral density as an illustrative example, we clarify the relationships among existing approaches, including the Hierarchical Equations of Motion (HEOM) and the Lindblad--pseudomode formalism, and demonstrate how this framework enables numerically stable and efficient simulations. This work provides both a transparent theoretical foundation for embedding techniques and a flexible platform for developing new methods to simulate non-Markovian quantum dynamics.
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.