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

The Journal of Chemical Physics

Non-Markovian reduced dynamics based upon a hierarchical effective-mode representation

Irène Burghardt, Rocco Martinazzo, Keith H. Hughes

Published Oct 12, 2012• 26 citations
Publisher

Abstract

A reduced dynamics representation is introduced which is tailored to a hierarchical, Mori-chain type representation of a bath of harmonic oscillators which are linearly coupled to a subsystem. We consider a spin-boson system where a single effective mode is constructed so as to absorb all system-environment interactions, while the residual bath modes are coupled bilinearly to the primary mode and among each other. Using a cumulant expansion of the memory kernel, correlation functions for the primary mode are obtained, which can be suitably approximated by truncated chains representing the primary-residual mode interactions. A series of reduced-dimensional bath correlation functions is thus obtained, which can be expressed as Fourier-Laplace transforms of spectral densities that are given in truncated continued-fraction form. For a master equation which is second order in the system-bath coupling, the memory kernel is re-expressed in terms of local-in-time equations involving auxiliary densities and auxiliary operators.

DOI
10.1063/1.4752078

Related papers

Effective-mode representation of non-Markovian dynamics: A hierarchical approximation of the spectral density. I. Application to single surface dynamics

Keith H. Hughes, Clara D. Christ, Irène Burghardt

An approach to non-Markovian system-environment dynamics is described which is based on the construction of a hierarchy of coupled effective environmental modes that is terminated by coupling the final member of the hierarchy to a Markovian bath. For an arbitrary environment, which is linearly coupled to the subsystem, the discretized spectral density is replaced by a series of approximate spectral densities involving an increasing number of effective modes. This series of approximants, which are constructed analytically in this paper, guarantees the accurate representation of the overall system-plus-bath dynamics up to increasing times. The hierarchical structure is manifested in the approximate spectral densities in the form of the imaginary part of a continued fraction similar to Mori theory. The results are described for cases where the hierarchy is truncated at the first-, second-, and third-order level. It is demonstrated that the results generated from a reduced density matrix equation of motion and large dimensional system-plus-bath wavepacket calculations are in excellent agreement. For the reduced density matrix calculations, the system and hierarchy of effective modes are treated explicitly and the effects of the bath on the final member of the hierarchy are described by the Caldeira-Leggett equation and its generalization to zero temperature.

The Journal of Chemical Physics•Jul 9, 2009•134 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

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

Rui‐Xue Xu, YiJing Yan

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.

arXiv:cond-mat/0612649•Physical Review E•Mar 15, 2007•177 citations

Entanglement dynamics of non-linearly coupled spin-boson model using hierarchical equation of motion approach

V A Wasesatama, Jusak Sali Kosasih

Abstract One of the most intensively studied problems in the field of open quantum systems is the spin-boson model. In this paper, we consider an extension of the spin-boson model with non-linear coupling. We solve the entanglement dynamics of two independent spins quadratically coupled to a bath of harmonic oscillators, using numerically ”exact” approach of hierarchical equations of motion (HEOM). Without Markovian approximation, generalized influence functional formulation by Hu et.al is used to describe the environmental effect of noisy quadratic bath in high-temperature regime. Wotters concurrence is then applied to numerically solved reduced density matrix to visualize the entanglement dynamics. By adjusting the environmental parameters, we vary the degree of non-markovianity and found the elongated time of entanglement.

Journal of Physics Conference Series•Jun 1, 2022•1 citation

Projection-Based Memory Kernel Coupling Theory for Quantum Dynamics: A Stable Framework for Non-Markovian Simulations

Wei Liu, Rui-Hao Bi, Yu Su +4 more

We present a projection-based, stability-preserving methodology for computing time correlation functions in open quantum systems governed by generalized quantum master equations with non-Markovian effects. Building upon the memory kernel coupling theory framework, our approach transforms the memory kernel hierarchy into a system of coupled linear differential equations through Mori-Zwanzig projection, followed by spectral projection onto stable eigenmodes to ensure numerical stability. By systematically eliminating unstable modes while preserving the physically relevant dynamics, our method guaranties long-time convergence without introducing artificial damping or ad hoc modifications. The theoretical framework maintains mathematical rigor through orthogonal projection operators and spectral decomposition. Benchmark calculations on the spin-boson model show excellent agreement with exact hierarchical equations of motion results while achieving significant computational efficiency. This approach provides a versatile and reliable framework for simulating non-Markovian dynamics in complex systems.

quant-pharXiv:2602.10629•Feb 12, 2026•0 citations

Reduced dynamics in spin-boson models: A method for both slow and fast bath

Andrei A. Golosov, Richard A. Friesner, Philip Pechukas

We study a model for treating dissipative systems, a one dimensional quantum system coupled to a harmonic bath. The dynamics of such a system can be described by Feynman’s path integral expression for the reduced density matrix. In this formulation the interaction of the system with the environment is stored in the influence functional. Recently we showed that fast environmental modes that give rise to correlations in the influence functional which are short range in time can be treated efficiently by a memory equation algorithm, which is a discretized version of a master equation. In this work we extend this approach to treat slow environmental modes as well, thereby efficiently linking adiabatic and nonadiabatic regimes. In this extended method the long range correlations in the influence functional arising from slow bath modes are taken into account through Stock’s semiclassical self-consistent-field approach.

The Journal of Chemical Physics•Feb 1, 2000•46 citations
A reduced dynamics representation is introduced which is tailored to a hierarchical, Mori-chain type representation of a bath of harmonic oscillators which are linearly coupled to a subsystem. We consider a spin-boson system where a single effective mode is constructed so as to absorb all system-environment interactions, while the residual bath modes are coupled bilinearly to the primary mode and among each other. Using a cumulant expansion of the memory kernel, correlation functions for the primary mode are obtained, which can be suitably approximated by truncated chains representing the primary-residual mode interactions. A series of reduced-dimensional bath correlation functions is thus obtained, which can be expressed as Fourier-Laplace transforms of spectral densities that are given in truncated continued-fraction form. For a master equation which is second order in the system-bath coupling, the memory kernel is re-expressed in terms of local-in-time equations involving auxiliary densities and auxiliary operators.
An approach to non-Markovian system-environment dynamics is described which is based on the construction of a hierarchy of coupled effective environmental modes that is terminated by coupling the final member of the hierarchy to a Markovian bath. For an arbitrary environment, which is linearly coupled to the subsystem, the discretized spectral density is replaced by a series of approximate spectral densities involving an increasing number of effective modes. This series of approximants, which are constructed analytically in this paper, guarantees the accurate representation of the overall system-plus-bath dynamics up to increasing times. The hierarchical structure is manifested in the approximate spectral densities in the form of the imaginary part of a continued fraction similar to Mori theory. The results are described for cases where the hierarchy is truncated at the first-, second-, and third-order level. It is demonstrated that the results generated from a reduced density matrix equation of motion and large dimensional system-plus-bath wavepacket calculations are in excellent agreement. For the reduced density matrix calculations, the system and hierarchy of effective modes are treated explicitly and the effects of the bath on the final member of the hierarchy are described by the Caldeira-Leggett equation and its generalization to zero temperature.
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.
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.
Abstract One of the most intensively studied problems in the field of open quantum systems is the spin-boson model. In this paper, we consider an extension of the spin-boson model with non-linear coupling. We solve the entanglement dynamics of two independent spins quadratically coupled to a bath of harmonic oscillators, using numerically ”exact” approach of hierarchical equations of motion (HEOM). Without Markovian approximation, generalized influence functional formulation by Hu et.al is used to describe the environmental effect of noisy quadratic bath in high-temperature regime. Wotters concurrence is then applied to numerically solved reduced density matrix to visualize the entanglement dynamics. By adjusting the environmental parameters, we vary the degree of non-markovianity and found the elongated time of entanglement.
We present a projection-based, stability-preserving methodology for computing time correlation functions in open quantum systems governed by generalized quantum master equations with non-Markovian effects. Building upon the memory kernel coupling theory framework, our approach transforms the memory kernel hierarchy into a system of coupled linear differential equations through Mori-Zwanzig projection, followed by spectral projection onto stable eigenmodes to ensure numerical stability. By systematically eliminating unstable modes while preserving the physically relevant dynamics, our method guaranties long-time convergence without introducing artificial damping or ad hoc modifications. The theoretical framework maintains mathematical rigor through orthogonal projection operators and spectral decomposition. Benchmark calculations on the spin-boson model show excellent agreement with exact hierarchical equations of motion results while achieving significant computational efficiency. This approach provides a versatile and reliable framework for simulating non-Markovian dynamics in complex systems.
We study a model for treating dissipative systems, a one dimensional quantum system coupled to a harmonic bath. The dynamics of such a system can be described by Feynman’s path integral expression for the reduced density matrix. In this formulation the interaction of the system with the environment is stored in the influence functional. Recently we showed that fast environmental modes that give rise to correlations in the influence functional which are short range in time can be treated efficiently by a memory equation algorithm, which is a discretized version of a master equation. In this work we extend this approach to treat slow environmental modes as well, thereby efficiently linking adiabatic and nonadiabatic regimes. In this extended method the long range correlations in the influence functional arising from slow bath modes are taken into account through Stock’s semiclassical self-consistent-field approach.