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

Physical Review A

Extracting bath information from open-quantum-system dynamics with the hierarchical equations-of-motion method

Chang Woo Kim

Published Oct 31, 2022• 5 citations
Publisher

Abstract

We present a simple and effective approach for extracting the bath information in open-quantum-system dynamics under the framework of the numerically exact hierarchical equations-of-motion method, for which it is usually not straightforward to adjust the boundary between the system and bath. The key component of our approach is an additional bath mode that weakly interacts with the system through the same coupling channel as the mode of interest, which we call the probe mode. By combining the probe mode with the system and calculating its expectation value, we can deduce the corresponding expectation value of the mode of interest via simple relations. The effectiveness of our approach is demonstrated by exact decomposition of dissipation by a continuous spectral density.

DOI
10.1103/physreva.106.042223

Related papers

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

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

RC-HEOM Hybrid Method for Non-Perturbative Open System Dynamics

Po-Rong Lai, Jhen-Dong Lin, Yi-Te Huang +3 more

The Hierarchical equations of motion (HEOM) method is an important non-perturbative technique, allowing numerically exact treatment of open quantum systems with strong coupling and non-Markovian memory. However, its encoding of bath memory into auxiliary density operators often limits direct access to detailed bath information. In contrast, the reaction-coordinate (RC) mapping allows direct and transparent access to the dominant collective bath mode, but its perturbative and often Markovian treatment of the residual bath restricts its reliability. In this work, we introduce RC-HEOM, a hybrid method that unifies the strengths of both approaches by combining RC mapping with a fully non-perturbative HEOM description of the residual bath. RC-HEOM simultaneously retains exact non-Markovian memory and access to the RC mode, which enables analysis of system-RC information. Applying this method to the Anderson impurity models, we directly track the emergence of the Kondo singlet from the growth of the Kondo resonance and uncover a nontrivial RC-mediated coherence revival. These results demonstrate that RC-HEOM is a promising method for characterizing open quantum systems in regimes that are difficult to access with conventional master-equation methods.

quant-pharXiv:2603.22833•Mar 25, 2026•0 citations

Hierarchical Equations-of-Motion Method for Momentum System–Bath Coupling

Maxim F. Gelin, Raffaele Borrelli, Lipeng Chen

For a broad class of quantum models of practical interest, we demonstrate that the Hamiltonian of the system nonlinearly coupled to a harmonic bath through the system and bath coordinates can be equivalently mapped into the Hamiltonian of the system bilinearly coupled to the bath through the system and bath momenta. We show that the Hamiltonian with bilinear system-bath momentum coupling can be treated by the hierarchical equations-of-motion (HEOM) method and present the corresponding proof-of-principle simulations. The developed methodology creates the opportunity to scrutinize a new family of nonlinear quantum systems by the numerically accurate HEOM method.

The Journal of Physical Chemistry B•Apr 30, 2021•5 citations

Tree tensor network state approach for solving hierarchical equations of motion

Yaling Ke

The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. The method is rooted in an exponential expansion of the bath correlation function, which in essence strategically reshapes a continuous environment into a set of effective bath modes that allow for more efficient cutoff at finite temperatures. Based on this understanding, one can map the HEOM method into a Schrödinger-like equation, with a non-Hermitian super-Hamiltonian for an extended wave function being the tensor product of the central system wave function and the Fock state of these effective bath modes. In this work, we explore the possibility of representing the extended wave function as a tree tensor network state (TTNS) and the super-Hamiltonian as a tree tensor network operator of the same structure as the TTNS, as well as the application of a time propagation algorithm using the time-dependent variational principle. Our benchmark calculations based on the spin-boson model with a slow-relaxing bath show that the proposed HEOM+TTNS approach yields consistent results with those of the conventional HEOM method, while the computation is considerably sped up. In addition, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme.

arXiv:2304.05151•The Journal of Chemical Physics•Jun 1, 2023•46 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 present a simple and effective approach for extracting the bath information in open-quantum-system dynamics under the framework of the numerically exact hierarchical equations-of-motion method, for which it is usually not straightforward to adjust the boundary between the system and bath. The key component of our approach is an additional bath mode that weakly interacts with the system through the same coupling channel as the mode of interest, which we call the probe mode. By combining the probe mode with the system and calculating its expectation value, we can deduce the corresponding expectation value of the mode of interest via simple relations. The effectiveness of our approach is demonstrated by exact decomposition of dissipation by a continuous spectral density.
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.
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) method is an important non-perturbative technique, allowing numerically exact treatment of open quantum systems with strong coupling and non-Markovian memory. However, its encoding of bath memory into auxiliary density operators often limits direct access to detailed bath information. In contrast, the reaction-coordinate (RC) mapping allows direct and transparent access to the dominant collective bath mode, but its perturbative and often Markovian treatment of the residual bath restricts its reliability. In this work, we introduce RC-HEOM, a hybrid method that unifies the strengths of both approaches by combining RC mapping with a fully non-perturbative HEOM description of the residual bath. RC-HEOM simultaneously retains exact non-Markovian memory and access to the RC mode, which enables analysis of system-RC information. Applying this method to the Anderson impurity models, we directly track the emergence of the Kondo singlet from the growth of the Kondo resonance and uncover a nontrivial RC-mediated coherence revival. These results demonstrate that RC-HEOM is a promising method for characterizing open quantum systems in regimes that are difficult to access with conventional master-equation methods.
For a broad class of quantum models of practical interest, we demonstrate that the Hamiltonian of the system nonlinearly coupled to a harmonic bath through the system and bath coordinates can be equivalently mapped into the Hamiltonian of the system bilinearly coupled to the bath through the system and bath momenta. We show that the Hamiltonian with bilinear system-bath momentum coupling can be treated by the hierarchical equations-of-motion (HEOM) method and present the corresponding proof-of-principle simulations. The developed methodology creates the opportunity to scrutinize a new family of nonlinear quantum systems by the numerically accurate HEOM method.
The hierarchical equations of motion (HEOM) method is a numerically exact open quantum system dynamics approach. The method is rooted in an exponential expansion of the bath correlation function, which in essence strategically reshapes a continuous environment into a set of effective bath modes that allow for more efficient cutoff at finite temperatures. Based on this understanding, one can map the HEOM method into a Schrödinger-like equation, with a non-Hermitian super-Hamiltonian for an extended wave function being the tensor product of the central system wave function and the Fock state of these effective bath modes. In this work, we explore the possibility of representing the extended wave function as a tree tensor network state (TTNS) and the super-Hamiltonian as a tree tensor network operator of the same structure as the TTNS, as well as the application of a time propagation algorithm using the time-dependent variational principle. Our benchmark calculations based on the spin-boson model with a slow-relaxing bath show that the proposed HEOM+TTNS approach yields consistent results with those of the conventional HEOM method, while the computation is considerably sped up. In addition, the simulation with a genuine TTNS is four times faster than a one-dimensional matrix product state decomposition scheme.
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.