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:2202.04059 · Physical Review Letters

Taming Quantum Noise for Efficient Low Temperature Simulations of Open Quantum Systems

Meng Xu, Yaming Yan, Qiang Shi, Joachim Ankerhold, Jürgen T. Stockburger

Published Nov 30, 2022• 110 citations
arXivHTMLPDF

Abstract

The hierarchical equations of motion (HEOM), derived from the exact Feynman-Vernon path integral, is one of the most powerful numerical methods to simulate the dynamics of open quantum systems. Its applicability has so far been limited to specific forms of spectral reservoir distributions and relatively elevated temperatures. Here we solve this problem and introduce an effective treatment of quantum noise in frequency space by systematically clustering higher order Matsubara poles, equivalent to an optimized rational decomposition. This leads to an elegant extension of the HEOM to arbitrary temperatures and very general reservoirs in combination with efficiency, high accuracy, and long-time stability. Moreover, the technique can directly be implemented in other approaches such as Green's function, stochastic, and pseudomode formulations. As one highly nontrivial application, for the subohmic spin-boson model at vanishing temperature the Shiba relation is quantitatively verified which predicts the long-time decay of correlation functions.

DOI
10.1103/physrevlett.129.230601

Related papers

Efficient low-temperature simulations for fermionic reservoirs with the hierarchical equations of motion method: Application to the Anderson impurity model

Xiaohan Dan, Meng Xu, Jürgen T. Stockburger +2 more

The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent effects of the bath, the reservoir correlation functions are usually decomposed into summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. This method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultralow temperature and general band structures. We demonstrate the efficiency, accuracy, and long-time stability of this decomposition scheme by applying it to the Anderson impurity model in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.

arXiv:2211.04089•Physical review. B./Physical review. B•May 22, 2023•45 citations

Hierarchical equations of motion method based on Fano spectrum decomposition for low temperature environments

Hou-Dao Zhang, Lei Cui, Hong Gong +3 more

The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.

The Journal of Chemical Physics•Feb 11, 2020•66 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

About the performance of perturbative treatments of the spin-boson dynamics within the hierarchical equations of motion approach

Meng Xu, Joachim Ankerhold

Abstract The hierarchical equations of motion (HEOM) provide a numerically exact approach for simulating the dynamics of open quantum systems coupled to a harmonic bath. However, its applicability has traditionally been limited to specific spectral forms and relatively high temperatures. Recently, an extended version called Free-Pole HEOM (FP-HEOM) has been developed to overcome these limitations. In this study, we demonstrate that the FP-HEOM method can be systematically employed to investigate higher order master equations by truncating the FP-HEOM hierarchy at a desired tier. We focus on the challenging scenario of the spin-boson problem with a sub-Ohmic spectral distribution at zero temperature and analyze the performance of the corresponding master equations. Furthermore, we compare the memory kernel for population dynamics obtained from the exact FP-HEOM dynamics with that of the approximate Non-Interacting-Blip Approximation (NIBA).

arXiv:2308.02010•The European Physical Journal Special Topics•Oct 11, 2023•5 citations

Exploiting the path-integral radius of gyration in open quantum dynamics

Andrew C. Hunt, Stuart C. Althorpe

A major challenge in open quantum dynamics is the inclusion of Matsubara-decay terms in the memory kernel, which arise from the quantum-Boltzmann delocalisation of the bath modes. This delocalisation can be quantified by the radius of gyration squared ${\mathcal R}^2(ω)$ of the imaginary-time Feynman paths of the bath modes as a function of the frequency $ω$. In a Hierarchical Equations of Motion (HEOM) calculation with a Debye--Drude spectral density, ${\mathcal R}^2(ω)$ is the only quantity that is treated approximately (assuming convergence with respect to hierarchy depth). Here, we show that the well-known Ishizaki--Tanimura correction is equivalent to separating smooth from `Brownian' contributions to ${\mathcal R}^2(ω)$, and that modifying the correction leads to a more efficient HEOM in the case of fast baths. We also develop a simple `A4' adaptation of the `AAA' (Adaptive Antoulas--Anderson) algorithm in order to fit ${\mathcal R}^2(ω)$ to a sum over poles, which results in an extremely efficient implementation of the standard HEOM method at low temperatures.

quant-pharXiv:2602.14647•Feb 17, 2026•0 citations

Chebyshev hierarchical equations of motion for systems with arbitrary spectral densities and temperatures

Hasan Rahman, Ulrich Kleinekathöfer

The time evolution in open quantum systems, such as a molecular aggregate in contact with a thermal bath, still poses a complex and challenging problem. The influence of the thermal noise can be treated using a plethora of schemes, several of which decompose the corresponding correlation functions in terms of weighted sums of exponential functions. One such scheme is based on the hierarchical equations of motion (HEOM), which is built using only certain forms of bath correlation functions. In the case where the environment is described by a complex spectral density or is at a very low temperature, approaches utilizing the exponential decomposition become very inefficient. Here, we utilize an alternative decomposition scheme for the bath correlation function based on Chebyshev polynomials and Bessel functions to derive a HEOM approach up to an arbitrary order in the environmental coupling. These hierarchical equations are similar in structure to the popular exponential HEOM scheme, but are formulated using the derivatives of the Bessel functions. The proposed scheme is tested up to the fourth order in perturbation theory for a two-level system and compared to benchmark calculations for the case of zero-temperature quantum Ohmic and super-Ohmic noise. Furthermore, the benefits and shortcomings of the present Chebyshev-based hierarchical equations are discussed.

arXiv:1904.06982•The Journal of Chemical Physics•Jun 24, 2019•39 citations
The hierarchical equations of motion (HEOM), derived from the exact Feynman-Vernon path integral, is one of the most powerful numerical methods to simulate the dynamics of open quantum systems. Its applicability has so far been limited to specific forms of spectral reservoir distributions and relatively elevated temperatures. Here we solve this problem and introduce an effective treatment of quantum noise in frequency space by systematically clustering higher order Matsubara poles, equivalent to an optimized rational decomposition. This leads to an elegant extension of the HEOM to arbitrary temperatures and very general reservoirs in combination with efficiency, high accuracy, and long-time stability. Moreover, the technique can directly be implemented in other approaches such as Green's function, stochastic, and pseudomode formulations. As one highly nontrivial application, for the subohmic spin-boson model at vanishing temperature the Shiba relation is quantitatively verified which predicts the long-time decay of correlation functions.
The hierarchical equations of motion (HEOM) approach is an accurate method to simulate open system quantum dynamics, which allows for systematic convergence to numerically exact results. To represent effects of the bath, the reservoir correlation functions are usually decomposed into summation of multiple exponential terms in the HEOM method. Since the reservoir correlation functions become highly non-Markovian at low temperatures or when the bath has complex band structures, a present challenge is to obtain accurate exponential decompositions that allow efficient simulation with the HEOM. In this work, we employ the barycentric representation to approximate the Fermi function and hybridization functions in the frequency domain. This method, by approximating these functions with optimized rational decomposition, greatly reduces the number of basis functions in decomposing the reservoir correlation functions, which further allows the HEOM method to be applied to ultralow temperature and general band structures. We demonstrate the efficiency, accuracy, and long-time stability of this decomposition scheme by applying it to the Anderson impurity model in the low-temperature regime with the Lorentzian and tight-binding hybridization functions.
The hierarchical equations of motion (HEOM) method has become one of the most popular methods for the studies of the open quantum system. However, its applicability to systems at ultra-low temperatures is largely restrained by the enormous computational cost, which is caused by the numerous exponential functions required to accurately characterize the non-Markovian memory of the reservoir environment. To overcome this problem, a Fano spectrum decomposition (FSD) scheme has been proposed recently [Cui et al., J. Chem. Phys. 151, 024110 (2019)], which expands the reservoir correlation functions using polynomial-exponential functions and hence greatly reduces the size of the memory basis set. In this work, we explicitly establish the FSD-based HEOM formalisms for both bosonic and fermionic environments. The accuracy and efficiency of the FSD-based HEOM are exemplified by the calculated low-temperature dissipative dynamics of a spin-boson model and the dynamic and static properties of a single-orbital Anderson impurity model in the Kondo regime. The encouraging numerical results highlight the practicality and usefulness of the FSD-based HEOM method for general open systems at ultra-low temperatures.
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.
Abstract The hierarchical equations of motion (HEOM) provide a numerically exact approach for simulating the dynamics of open quantum systems coupled to a harmonic bath. However, its applicability has traditionally been limited to specific spectral forms and relatively high temperatures. Recently, an extended version called Free-Pole HEOM (FP-HEOM) has been developed to overcome these limitations. In this study, we demonstrate that the FP-HEOM method can be systematically employed to investigate higher order master equations by truncating the FP-HEOM hierarchy at a desired tier. We focus on the challenging scenario of the spin-boson problem with a sub-Ohmic spectral distribution at zero temperature and analyze the performance of the corresponding master equations. Furthermore, we compare the memory kernel for population dynamics obtained from the exact FP-HEOM dynamics with that of the approximate Non-Interacting-Blip Approximation (NIBA).
A major challenge in open quantum dynamics is the inclusion of Matsubara-decay terms in the memory kernel, which arise from the quantum-Boltzmann delocalisation of the bath modes. This delocalisation can be quantified by the radius of gyration squared R2(ω){\mathcal R}^2(ω)R2(ω) of the imaginary-time Feynman paths of the bath modes as a function of the frequency ωωω. In a Hierarchical Equations of Motion (HEOM) calculation with a Debye--Drude spectral density, R2(ω){\mathcal R}^2(ω)R2(ω) is the only quantity that is treated approximately (assuming convergence with respect to hierarchy depth). Here, we show that the well-known Ishizaki--Tanimura correction is equivalent to separating smooth from `Brownian' contributions to R2(ω){\mathcal R}^2(ω)R2(ω), and that modifying the correction leads to a more efficient HEOM in the case of fast baths. We also develop a simple `A4' adaptation of the `AAA' (Adaptive Antoulas--Anderson) algorithm in order to fit R2(ω){\mathcal R}^2(ω)R2(ω) to a sum over poles, which results in an extremely efficient implementation of the standard HEOM method at low temperatures.
The time evolution in open quantum systems, such as a molecular aggregate in contact with a thermal bath, still poses a complex and challenging problem. The influence of the thermal noise can be treated using a plethora of schemes, several of which decompose the corresponding correlation functions in terms of weighted sums of exponential functions. One such scheme is based on the hierarchical equations of motion (HEOM), which is built using only certain forms of bath correlation functions. In the case where the environment is described by a complex spectral density or is at a very low temperature, approaches utilizing the exponential decomposition become very inefficient. Here, we utilize an alternative decomposition scheme for the bath correlation function based on Chebyshev polynomials and Bessel functions to derive a HEOM approach up to an arbitrary order in the environmental coupling. These hierarchical equations are similar in structure to the popular exponential HEOM scheme, but are formulated using the derivatives of the Bessel functions. The proposed scheme is tested up to the fourth order in perturbation theory for a two-level system and compared to benchmark calculations for the case of zero-temperature quantum Ohmic and super-Ohmic noise. Furthermore, the benefits and shortcomings of the present Chebyshev-based hierarchical equations are discussed.