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.10273 · The Journal of Chemical Physics

Hierarchical equations of motion approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space

Yaling Ke, Raffaele Borrelli, Michael Thoss

Published Apr 28, 2022• 49 citations
arXivHTMLPDF

Abstract

We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation [R. Borrelli, J. Chem. Phys. 150, 234102 (2019)] to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.

DOI
10.1063/5.0088947

Related papers

Nonequilibrium open quantum systems with multiple bosonic and fermionic environments: A hierarchical equations of motion approach

Jakob Bätge, Yaling Ke, Christoph Kaspar +1 more

We present a hierarchical equations of motion approach, which allows a numerically exact simulation of nonequilibrium transport in general open quantum systems involving multiple macroscopic bosonic and fermionic environments. The performance of the method is demonstrated for a model of a nanosystem, which involves interacting electronic and vibrational degrees of freedom and is coupled to fermionic and bosonic baths. The results show the intricate interplay of electronic and vibrational degrees of freedom in this nonequilibrium transport scenario for both voltage and thermally driven transport processes. Furthermore, the use of importance criteria to improve the efficiency of the method is discussed.

arXiv:2102.09484•Physical review. B./Physical review. B•Jun 8, 2021•30 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

An efficient Julia framework for hierarchical equations of motion in open quantum systems

Yi-Te Huang, Po-Chen Kuo, Neill Lambert +5 more

Abstract The hierarchical equations of motion (HEOM) approach can describe the reduced dynamics of a system simultaneously coupled to multiple bosonic and fermionic environments. The complexity of exactly describing the system-environment interaction with the HEOM method usually results in time-consuming calculations and a large memory cost. Here, we introduce an open-source software package called HierarchicalEOM.jl: a Julia framework integrating the HEOM approach. HierarchicalEOM.jl features a collection of methods to compute bosonic and fermionic spectra, stationary states, and the full dynamics in the extended space of all auxiliary density operators (ADOs). The required handling of the ADOs multi-indexes is achieved through a user-friendly interface. We exemplify the functionalities of the package by analyzing a single impurity Anderson model, and an ultra-strongly coupled charge-cavity system interacting with bosonic and fermionic reservoirs. HierarchicalEOM.jl achieves a significant speedup with respect to the corresponding method in the Quantum Toolbox in Python (QuTiP), upon which this package is founded.

arXiv:2306.07522•Communications Physics•Oct 23, 2023•26 citations

Hierarchical Equations of Motion in Matrix Product States: Formalism and Applications for Charge Transport

Hengrui Yang, Zirui Sheng, Liqi Zhou +1 more

We present a Hierarchical Equations of Motion (HEOM) approach in the Matrix Product State (MPS) formalism to simulate carrier transport in molecular aggregates described by an electron-phonon Hamiltonian with bosonic dissipation. Transport properties are evaluated through time-dependent population analysis and mobility calculations. The method's validity is rigorously established through benchmarking against conventional HEOM. Comparative analysis with Thermo-Field Dynamics combined with MPS (TFD + MPS) reveals fundamental similarities and differences in their effective Hamiltonians and demonstrates the superior accuracy and computational efficiency of our HEOM + MPS framework. For single-electron systems, we introduce state-vector space configurations that enhance performance beyond traditional Fock space approaches. Results confirm that our method provides a robust, nearly exact, and efficient numerical quantum dynamic approach for carrier transport in dissipative bosonic environments.

Journal of Chemical Theory and Computation•Oct 16, 2025•1 citation

Density matrix dynamics in twin-formulation: An efficient methodology based on tensor-train representation of reduced equations of motion

Raffaele Borrelli

The twin-formulation of quantum statistical mechanics is employed to describe a new methodology for the solution of the equations of motion of the reduced density matrix in their hierarchical formulation. It is shown that the introduction of tilde operators and of their algebra in the dual space greatly simplifies the application of numerical techniques for the propagation of the density matrix. The application of tensor-train representation of a vector to solve complex quantum dynamical problems within the framework of the twin-formulation is discussed. Next, applications of the hierarchical equations of motion to a dissipative polaron model are presented showing the validity and accuracy of the new approach.

The Journal of Chemical Physics•Jun 17, 2019•75 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
We extend the twin-space formulation of the hierarchical equations of motion approach in combination with the matrix product state representation [R. Borrelli, J. Chem. Phys. 150, 234102 (2019)] to nonequilibrium scenarios where the open quantum system is coupled to a hybrid fermionic and bosonic environment. The key ideas used in the extension are a reformulation of the hierarchical equations of motion for the auxiliary density matrices into a time-dependent Schrödinger-like equation for an augmented multi-dimensional wave function as well as a tensor decomposition into a product of low-rank matrices. The new approach facilitates accurate simulations of non-equilibrium quantum dynamics in larger and more complex open quantum systems. The performance of the method is demonstrated for a model of a molecular junction exhibiting current-induced mode-selective vibrational excitation.
We present a hierarchical equations of motion approach, which allows a numerically exact simulation of nonequilibrium transport in general open quantum systems involving multiple macroscopic bosonic and fermionic environments. The performance of the method is demonstrated for a model of a nanosystem, which involves interacting electronic and vibrational degrees of freedom and is coupled to fermionic and bosonic baths. The results show the intricate interplay of electronic and vibrational degrees of freedom in this nonequilibrium transport scenario for both voltage and thermally driven transport processes. Furthermore, the use of importance criteria to improve the efficiency of the method is discussed.
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.
Abstract The hierarchical equations of motion (HEOM) approach can describe the reduced dynamics of a system simultaneously coupled to multiple bosonic and fermionic environments. The complexity of exactly describing the system-environment interaction with the HEOM method usually results in time-consuming calculations and a large memory cost. Here, we introduce an open-source software package called HierarchicalEOM.jl: a Julia framework integrating the HEOM approach. HierarchicalEOM.jl features a collection of methods to compute bosonic and fermionic spectra, stationary states, and the full dynamics in the extended space of all auxiliary density operators (ADOs). The required handling of the ADOs multi-indexes is achieved through a user-friendly interface. We exemplify the functionalities of the package by analyzing a single impurity Anderson model, and an ultra-strongly coupled charge-cavity system interacting with bosonic and fermionic reservoirs. HierarchicalEOM.jl achieves a significant speedup with respect to the corresponding method in the Quantum Toolbox in Python (QuTiP), upon which this package is founded.
We present a Hierarchical Equations of Motion (HEOM) approach in the Matrix Product State (MPS) formalism to simulate carrier transport in molecular aggregates described by an electron-phonon Hamiltonian with bosonic dissipation. Transport properties are evaluated through time-dependent population analysis and mobility calculations. The method's validity is rigorously established through benchmarking against conventional HEOM. Comparative analysis with Thermo-Field Dynamics combined with MPS (TFD + MPS) reveals fundamental similarities and differences in their effective Hamiltonians and demonstrates the superior accuracy and computational efficiency of our HEOM + MPS framework. For single-electron systems, we introduce state-vector space configurations that enhance performance beyond traditional Fock space approaches. Results confirm that our method provides a robust, nearly exact, and efficient numerical quantum dynamic approach for carrier transport in dissipative bosonic environments.
The twin-formulation of quantum statistical mechanics is employed to describe a new methodology for the solution of the equations of motion of the reduced density matrix in their hierarchical formulation. It is shown that the introduction of tilde operators and of their algebra in the dual space greatly simplifies the application of numerical techniques for the propagation of the density matrix. The application of tensor-train representation of a vector to solve complex quantum dynamical problems within the framework of the twin-formulation is discussed. Next, applications of the hierarchical equations of motion to a dissipative polaron model are presented showing the validity and accuracy of the new approach.
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.