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 Physical Chemistry B

Expanding the Range of Hierarchical Equations of Motion by Tensor-Train Implementation

Raffaele Borrelli, Sergey Dolgov

Published May 13, 2021• 38 citations
Publisher

Abstract

The non-equilibrium thermo-field dynamics formulation of the hierarchical equations of motion combined with the tensor-train representation of the density matrix is discussed, and a new numerical integration scheme is introduced. The numerical methodology is based on an adaptive low-rank Galerkin reduction scheme and can preserve linear invariants (such as the trace of the density matrix). The method is applied to the study of the charge transfer dynamics in model pentacene molecular aggregates. The combined effect of a discrete set of molecular vibrational modes and a thermal bath is investigated, paying special attention to the coherent-incoherent transition of the charge transport. The new computational framework is shown to be a very promising methodology for the study of the quantum dynamics of complex molecular systems in the condensed phase.

DOI
10.1021/acs.jpcb.1c02724

Related papers

Tensor-Train Format Hierarchical Equations of Motion Formalism: Charge Transfer in Organic Semiconductors via Dissipative Holstein Models

H. Takahashi, Raffaele Borrelli

Hierarchical Equations of Motion (HEOM) in the Tensor-Train (TT) representation is applied to study the charge-transfer dynamics in organic semiconductors (OSCs). The theoretical formulation as well as the basic computational aspects of HEOM-TT are discussed in detail. Charge transfer in OSCs is modeled using dissipative polaronic models that incorporate the effects of both high- and low-frequency molecular vibrations, and it is simulated in a fully quantum and nonperturbative manner, which has not been studied intensively. The capability of treating complex electron-vibrational systems is examined by analyzing and comparing the numerical behavior of the time-dependent variational approach and the time-Alternating Minimal Energy methods and by calculating the current autocorrelation function and diffusivity across various models. Our results indicate that the HEOM-TT framework offers a robust tool for the detailed analysis of complex polaronic systems, suggesting its potential for broader applications.

Journal of Chemical Theory and Computation•Aug 17, 2024•13 citations

Finite temperature quantum dynamics of complex systems: Integrating thermo‐field theories and tensor‐train methods

Raffaele Borrelli, Maxim F. Gelin

Abstract This review provides the fundamental theoretical tools for the development of a complete wave‐function formalism for the study of time‐evolution of chemico‐physical systems at finite temperature. The methodology is based on the non‐equilibrium thermo‐field dynamics (NE‐TFD) representation of quantum mechanics, which is alternative to the commonly used density matrix representation. TFD concepts are extended and integrated with the tensor‐train (TT) numerical tools leading to a novel and powerful theoretical and computational framework for the study of complex quantum dynamical problems. In addition, NE‐TFD techniques are extended to enable the study of dissipative open systems via a new formulation of the hierarchical equations of motion (HEOM) fully integrated with TT methodologies. We demonstrate that the combination of the TFD machinery with computational advantages of TTs results in a powerful theoretical and computational framework for scrutinizing dynamics of complex multidimensional electron‐vibrational systems. We illustrate the validity and the computational advantages of the developed methodologies by applying them to the study of quantum coherence effects in the energy‐transfer processes in antenna systems, to the analysis of fingerprints of vibrational modes in electron‐transfer and charge‐transfer processes in various model and realistic multidimensional molecular systems, as well as to simulation of other fundamental models of physical chemistry. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and Kinetics Theoretical and Physical Chemistry > Statistical Mechanics

Wiley Interdisciplinary Reviews Computational Molecular Science•May 10, 2021•95 citations

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 approach to hybrid fermionic and bosonic environments: Matrix product state formulation in twin space

Yaling Ke, Raffaele Borrelli, Michael Thoss

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.

arXiv:2202.10273•The Journal of Chemical Physics•Apr 28, 2022•49 citations

Quantum electron-vibrational dynamics at finite temperature: Thermo field dynamics approach

Raffaele Borrelli, Maxim F. Gelin

Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on the thermo field dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. A comparison with the theoretically equivalent density matrix formulation shows the key numerical advantages of the present approach. The solution of thermo field dynamics equations with a novel technique for the propagation of tensor trains (matrix product states) is discussed. Numerical applications to model spin-boson systems show that the present approach is a promising tool for the description of quantum dynamics of complex molecular systems at finite temperature.

The Journal of Chemical Physics•Dec 8, 2016•91 citations

Efficient quantum dynamics simulations of complex molecular systems: A unified treatment of dynamic and static disorder

Maxim F. Gelin, Amalia Velardo, Raffaele Borrelli

We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.

The Journal of Chemical Physics•Oct 1, 2021•32 citations
The non-equilibrium thermo-field dynamics formulation of the hierarchical equations of motion combined with the tensor-train representation of the density matrix is discussed, and a new numerical integration scheme is introduced. The numerical methodology is based on an adaptive low-rank Galerkin reduction scheme and can preserve linear invariants (such as the trace of the density matrix). The method is applied to the study of the charge transfer dynamics in model pentacene molecular aggregates. The combined effect of a discrete set of molecular vibrational modes and a thermal bath is investigated, paying special attention to the coherent-incoherent transition of the charge transport. The new computational framework is shown to be a very promising methodology for the study of the quantum dynamics of complex molecular systems in the condensed phase.
Hierarchical Equations of Motion (HEOM) in the Tensor-Train (TT) representation is applied to study the charge-transfer dynamics in organic semiconductors (OSCs). The theoretical formulation as well as the basic computational aspects of HEOM-TT are discussed in detail. Charge transfer in OSCs is modeled using dissipative polaronic models that incorporate the effects of both high- and low-frequency molecular vibrations, and it is simulated in a fully quantum and nonperturbative manner, which has not been studied intensively. The capability of treating complex electron-vibrational systems is examined by analyzing and comparing the numerical behavior of the time-dependent variational approach and the time-Alternating Minimal Energy methods and by calculating the current autocorrelation function and diffusivity across various models. Our results indicate that the HEOM-TT framework offers a robust tool for the detailed analysis of complex polaronic systems, suggesting its potential for broader applications.
Abstract This review provides the fundamental theoretical tools for the development of a complete wave‐function formalism for the study of time‐evolution of chemico‐physical systems at finite temperature. The methodology is based on the non‐equilibrium thermo‐field dynamics (NE‐TFD) representation of quantum mechanics, which is alternative to the commonly used density matrix representation. TFD concepts are extended and integrated with the tensor‐train (TT) numerical tools leading to a novel and powerful theoretical and computational framework for the study of complex quantum dynamical problems. In addition, NE‐TFD techniques are extended to enable the study of dissipative open systems via a new formulation of the hierarchical equations of motion (HEOM) fully integrated with TT methodologies. We demonstrate that the combination of the TFD machinery with computational advantages of TTs results in a powerful theoretical and computational framework for scrutinizing dynamics of complex multidimensional electron‐vibrational systems. We illustrate the validity and the computational advantages of the developed methodologies by applying them to the study of quantum coherence effects in the energy‐transfer processes in antenna systems, to the analysis of fingerprints of vibrational modes in electron‐transfer and charge‐transfer processes in various model and realistic multidimensional molecular systems, as well as to simulation of other fundamental models of physical chemistry. This article is categorized under: Theoretical and Physical Chemistry > Reaction Dynamics and Kinetics Theoretical and Physical Chemistry > Statistical Mechanics
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.
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.
Quantum electron-vibrational dynamics in molecular systems at finite temperature is described using an approach based on the thermo field dynamics theory. This formulation treats temperature effects in the Hilbert space without introducing the Liouville space. A comparison with the theoretically equivalent density matrix formulation shows the key numerical advantages of the present approach. The solution of thermo field dynamics equations with a novel technique for the propagation of tensor trains (matrix product states) is discussed. Numerical applications to model spin-boson systems show that the present approach is a promising tool for the description of quantum dynamics of complex molecular systems at finite temperature.
We present a unified and highly numerically efficient formalism for the simulation of quantum dynamics of complex molecular systems, which takes into account both temperature effects and static disorder. The methodology is based on the thermo-field dynamics formalism, and Gaussian static disorder is included into simulations via auxiliary bosonic operators. This approach, combined with the tensor-train/matrix-product state representation of the thermalized stochastic wave function, is applied to study the effect of dynamic and static disorders in charge-transfer processes in model organic semiconductor chains employing the Su-Schrieffer-Heeger (Holstein-Peierls) model Hamiltonian.