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

Generalization of the hierarchical equations of motion theory for efficient calculations with arbitrary correlation functions

Tatsushi Ikeda, Gregory D. Scholes

Published May 22, 2020• 85 citations
arXivHTMLPDF

Abstract

The hierarchical equations of motion (HEOM) theory is one of the standard methods to rigorously describe open quantum dynamics coupled to harmonic environments. Such a model is used to capture non-Markovian and non-perturbative effects of environments appearing in ultrafast phenomena. In the regular framework of the HEOM theory, the environment correlation functions are restricted to linear combinations of exponential functions. In this article, we present a new formulation of the HEOM theory including treatment of non-exponential correlation functions, which enables us to describe general environmental effects more efficiently and stably than the original theory and other generalizations. The library and its Python binding we developed to perform simulations based on our approach, named LibHEOM and PyHEOM, respectively, are provided as the supplementary material.

DOI
10.1063/5.0007327

Related papers

Collective bath coordinate mapping of “hierarchy” in hierarchical equations of motion

Tatsushi Ikeda, Akira Nakayama

The theory of hierarchical equations of motion (HEOM) is one of the standard methods to give exact evaluations of the dynamics as coupled to harmonic oscillator environments. However, the theory is numerically demanding due to its hierarchy, which is the set of auxiliary elements introduced to capture the non-Markovian and non-perturbative effects of environments. When system-bath coupling becomes relatively strong, the required computational resources and precision move beyond the regime that can be currently handled. This article presents a new representation of HEOM theory in which the hierarchy is mapped into a continuous space of a collective bath coordinate and several auxiliary coordinates as the form of the quantum Fokker-Planck equation. This representation gives a rigorous time evolution of the bath coordinate distribution and is more stable and efficient than the original HEOM theory, particularly when there is a strong system-bath coupling. We demonstrate the suitability of this approach to treat vibronic system models coupled to environments.

arXiv:2112.09861•The Journal of Chemical Physics•Feb 9, 2022•14 citations

Input-output hierarchical equations of motion

Mauro Cirio, Pengfei Liang, Neill Lambert

We derive an extended version of the hierarchical equations of motion (HEOM) to compute output physical properties of a bosonic environment, which is allowed to be initially prepared at an earlier time in a non-Gaussian input state and then non-perturbatively interact with a quantum system with a linear environmental operator. While spectral assumptions analogous to the ones used in the regular HEOM are imposed to compute dynamical output bath observables, they are not required to model input states or output observables at a fixed time, in this case leading to time-dependent contributions to the equations. In the Markovian limit, we use this formalism to derive an input-output Lindblad equation which can be used to extend the applicability of the regular version. For a given desired input state and output observable, the range of the indexes extending the regular HEOM is, by construction, bounded. Overall, the aim of this formalism is to take advantage of the efficiency and generality of the HEOM framework to model non-Gaussian input states and the dynamics of environmental observables in bosonic, non-Markovian open quantum systems.

arXiv:2408.12221•Physical Review A•Jun 16, 2025•3 citations

Numerically “exact” approach to open quantum dynamics: The hierarchical equations of motion (HEOM)

Yoshitaka Tanimura

An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.

arXiv:2006.05501•The Journal of Chemical Physics•Jul 8, 2020•536 citations

Theoretical method for open quantum systems: Progresses and perspectives

YiJing Yan

The hierarchical equations of motion (HEOM) formalism developed recently is becoming an important method for investigating the properties and dynamical processes of quantum open systems. A comprehensive account on its establishment and recent development was given, covering theoretical formulations, numerical algorithm and practical applications. By construction HEOM resolves in a nonperturbative manner the combined effects of many-body interaction, system-environment coupling and non-Markovian memory. While it is mathematically equivalent to an exact path integral formalism, HEOM is also very convenient in computation and versatile in various systems of study. Moreover, it is shown recently that HEOM is actually a correlated system-environment dynamics theory. The practicality of HEOM to address physical and chemical problems was exemplified with examples relevant to quantum dissipation and quantum transport.

•Jan 1, 2013•0 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

One-to-one correspondence between Hierarchical Equations of Motion and Pseudomodes for Open Quantum System Dynamics

Kai Müller, Walter T. Strunz

We unite two of the most widely used approaches for strongly damped, non-Markovian open quantum dynamics, the Hierarchical Equations of Motion (HEOM) and the pseudomode method by proving two statements: First, every physical bath correlation function (BCF) that can be written as a sum of $N$ exponential terms can be obtained from a physical model with $N$ interacting pseudomodes which are damped in Lindblad form. Second, for every such BCF there exists a non-unitary, linear transformation which mirrors the evolution of the system-pseudomode state onto the HEOM hierarchy, and vice versa. Our proofs are constructive and we give explicit expressions for the mirror transformation as well as for the pseudomode Lindbladian corresponding to a given exponential BCF. This approach also gives insight and provides elegant derivations of the corresponding Hierarchy of stochastic Pure States (HOPS) method and its nearly-unitary version, nuHOPS. Our result opens several avenues for further optimization of non-Markovian open quantum system dynamics methods.

quant-pharXiv:2604.06466•Apr 9, 2026•0 citations
The hierarchical equations of motion (HEOM) theory is one of the standard methods to rigorously describe open quantum dynamics coupled to harmonic environments. Such a model is used to capture non-Markovian and non-perturbative effects of environments appearing in ultrafast phenomena. In the regular framework of the HEOM theory, the environment correlation functions are restricted to linear combinations of exponential functions. In this article, we present a new formulation of the HEOM theory including treatment of non-exponential correlation functions, which enables us to describe general environmental effects more efficiently and stably than the original theory and other generalizations. The library and its Python binding we developed to perform simulations based on our approach, named LibHEOM and PyHEOM, respectively, are provided as the supplementary material.
The theory of hierarchical equations of motion (HEOM) is one of the standard methods to give exact evaluations of the dynamics as coupled to harmonic oscillator environments. However, the theory is numerically demanding due to its hierarchy, which is the set of auxiliary elements introduced to capture the non-Markovian and non-perturbative effects of environments. When system-bath coupling becomes relatively strong, the required computational resources and precision move beyond the regime that can be currently handled. This article presents a new representation of HEOM theory in which the hierarchy is mapped into a continuous space of a collective bath coordinate and several auxiliary coordinates as the form of the quantum Fokker-Planck equation. This representation gives a rigorous time evolution of the bath coordinate distribution and is more stable and efficient than the original HEOM theory, particularly when there is a strong system-bath coupling. We demonstrate the suitability of this approach to treat vibronic system models coupled to environments.
We derive an extended version of the hierarchical equations of motion (HEOM) to compute output physical properties of a bosonic environment, which is allowed to be initially prepared at an earlier time in a non-Gaussian input state and then non-perturbatively interact with a quantum system with a linear environmental operator. While spectral assumptions analogous to the ones used in the regular HEOM are imposed to compute dynamical output bath observables, they are not required to model input states or output observables at a fixed time, in this case leading to time-dependent contributions to the equations. In the Markovian limit, we use this formalism to derive an input-output Lindblad equation which can be used to extend the applicability of the regular version. For a given desired input state and output observable, the range of the indexes extending the regular HEOM is, by construction, bounded. Overall, the aim of this formalism is to take advantage of the efficiency and generality of the HEOM framework to model non-Gaussian input states and the dynamics of environmental observables in bosonic, non-Markovian open quantum systems.
An open quantum system refers to a system that is further coupled to a bath system consisting of surrounding radiation fields, atoms, molecules, or proteins. The bath system is typically modeled by an infinite number of harmonic oscillators. This system-bath model can describe the time-irreversible dynamics through which the system evolves toward a thermal equilibrium state at finite temperature. In nuclear magnetic resonance and atomic spectroscopy, dynamics can be studied easily by using simple quantum master equations under the assumption that the system-bath interaction is weak (perturbative approximation) and the bath fluctuations are very fast (Markovian approximation). However, such approximations cannot be applied in chemical physics and biochemical physics problems, where environmental materials are complex and strongly coupled with environments. The hierarchical equations of motion (HEOM) can describe the numerically "exact" dynamics of a reduced system under nonperturbative and non-Markovian system-bath interactions, which has been verified on the basis of exact analytical solutions (non-Markovian tests) with any desired numerical accuracy. The HEOM theory has been used to treat systems of practical interest, in particular, to account for various linear and nonlinear spectra in molecular and solid state materials, to evaluate charge and exciton transfer rates in biological systems, to simulate resonant tunneling and quantum ratchet processes in nanodevices, and to explore quantum entanglement states in quantum information theories. This article presents an overview of the HEOM theory, focusing on its theoretical background and applications, to help further the development of the study of open quantum dynamics.
The hierarchical equations of motion (HEOM) formalism developed recently is becoming an important method for investigating the properties and dynamical processes of quantum open systems. A comprehensive account on its establishment and recent development was given, covering theoretical formulations, numerical algorithm and practical applications. By construction HEOM resolves in a nonperturbative manner the combined effects of many-body interaction, system-environment coupling and non-Markovian memory. While it is mathematically equivalent to an exact path integral formalism, HEOM is also very convenient in computation and versatile in various systems of study. Moreover, it is shown recently that HEOM is actually a correlated system-environment dynamics theory. The practicality of HEOM to address physical and chemical problems was exemplified with examples relevant to quantum dissipation and quantum transport.
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 unite two of the most widely used approaches for strongly damped, non-Markovian open quantum dynamics, the Hierarchical Equations of Motion (HEOM) and the pseudomode method by proving two statements: First, every physical bath correlation function (BCF) that can be written as a sum of NNN exponential terms can be obtained from a physical model with NNN interacting pseudomodes which are damped in Lindblad form. Second, for every such BCF there exists a non-unitary, linear transformation which mirrors the evolution of the system-pseudomode state onto the HEOM hierarchy, and vice versa. Our proofs are constructive and we give explicit expressions for the mirror transformation as well as for the pseudomode Lindbladian corresponding to a given exponential BCF. This approach also gives insight and provides elegant derivations of the corresponding Hierarchy of stochastic Pure States (HOPS) method and its nearly-unitary version, nuHOPS. Our result opens several avenues for further optimization of non-Markovian open quantum system dynamics methods.