Hierarchical Equations Of Motion Theory Research Articles

Collective bath coordinate mapping of "hierarchy" in hierarchical equations of motion.

The theory of hierarchical equations of motion (HEOM) is one of the standard methods to give exact evaluations of the dynamics as coupled to harmonic oscillator environments. However, the theory is numerically demanding due to its hierarchy, which is the set of auxiliary elements introduced to capture the non-Markovian and non-perturbative effects of environments. When system-bath coupling becomes relatively strong, the required computational resources and precision move beyond the regime that can be currently handled. This article presents a new representation of HEOM theory in which the hierarchy is mapped into a continuous space of a collective bath coordinate and several auxiliary coordinates as the form of the quantum Fokker-Planck equation. This representation gives a rigorous time evolution of the bath coordinate distribution and is more stable and efficient than the original HEOM theory, particularly when there is a strong system-bath coupling. We demonstrate the suitability of this approach to treat vibronic system models coupled to environments.

Numerically "exact" approach to open quantum dynamics: The hierarchical equations of motion (HEOM).

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.

Stochastic equation of motion approach to fermionic dissipative dynamics. I. Formalism.

In this work, we establish formally exact stochastic equation of motion (SEOM) theory to describe the dissipative dynamics of fermionic open systems. The construction of the SEOM is based on a stochastic decoupling of the dissipative interaction between the system and fermionic environment, and the influence of environmental fluctuations on the reduced system dynamics is characterized by stochastic Grassmann fields. Meanwhile, numerical realization of the time-dependent Grassmann fields has remained a long-standing challenge. To solve this problem, we propose a minimal auxiliary space (MAS) mapping scheme with which the stochastic Grassmann fields are represented by conventional c-number fields along with a set of pseudo-levels. This eventually leads to a numerically feasible MAS-SEOM method. The important properties of the MAS-SEOM are analyzed by making connection to the well-established time-dependent perturbation theory and the hierarchical equations of motion theory. The MAS-SEOM method provides a potentially promising approach for the accurate and efficient simulation of fermionic open systems at ultra-low temperatures.

Generalization of the hierarchical equations of motion theory for efficient calculations with arbitrary correlation functions.

The hierarchical equations of motion (HEOM) theory is one of the standard methods to rigorously describe open quantum dynamics coupled to harmonic environments. Such a model is used to capture non-Markovian and non-perturbative effects of environments appearing in ultrafast phenomena. In the regular framework of the HEOM theory, the environment correlation functions are restricted to linear combinations of exponential functions. In this article, we present a new formulation of the HEOM theory including treatment of non-exponential correlation functions, which enables us to describe general environmental effects more efficiently and stably than the original theory and other generalizations. The library and its Python binding we developed to perform simulations based on our approach, named LibHEOM and PyHEOM, respectively, are provided as the supplementary material.

Precise simulation of strongly correlated quantum impurity systems

In recent years, nano-sized systems involving local charge or spin states have received wide interests because of their potential application in emerging fields such as quantum information and quantum computation. For instance, organometallic molecular complexes may serve as building blocks of quantum storage devices, because the spin-unpaired d or f electrons at transition metal centers may be employed to construct spin qubits. Moreover, if the system contains strong electron-electron interaction, the involving local quantum states are subject to prominent electron correlation effects (such as the Kondo effect). Theoretically, spatially confined nano-systems are often described by quantum impurity models. Thus the accurate prediction of the intrinsic properties of quantum impurity systems and the deep understanding on the response and evolution of local quantum states under external fields or in dissipative environment are fundamentally important for the design and fabrication of quantum devices. The accurate characterization of quantum coherence, correlation, and entanglement in quantum impurity systems remains a great challenge. Enormous efforts have been made to achieve this goal. A variety of theoretical methods have been developed, including the numerical renormalization group method, the quantum Monte Carlo method, and many others. However, all the existing methods are subject to certain limitations regarding accuracy or efficiency. Therefore, we choose to view this problem from a new perspective—the perspective of open quantum systems. Over the past decade, we have developed a formally exact quantum dissipation theory, the hierarchical equations of motion (HEOM) theory, for fermionic open systems. The HEOM theory captures the combined effects of system-environment dissipation, many-body interaction, and non-Markovian memory in a nonperturbative manner. It is capable of addressing static and dynamic responses of system observables in both equilibrium and nonequilibrium situations. We have implemented the HEOM method in our self-designed computer program HEOM-QUICK. We have also devised a series of advanced algorithms which substantially enhance the numerical efficiency of HEOM. The HEOM-QUICK program thus provides an accurate, efficient, and versatile theoretical tool for the investigation of strongly correlated quantum impurity systems. We have applied the HEOM method to study a variety of problems associated with quantum impurity systems. These include the intrinsic properties of quantum impurity models and lattice models, the transient current response of quantum dots to ac voltages, the thermopower and local heating effect in nonequilibrium quantum dots, and the tuning of local spin states in adsorbed molecular magnets. In particular, the HEOM method has been combined with density functional theory (DFT) method, and thus allows for first-principles-based simulation on the precise tuning of local spin states in adsorbed magnetic molecules. Numerical simulations have discovered or reproduced many important quantum phenomena that are essential for relevant experiments, including the Kondo effect, Mott metal-insulator transition, quantum memristive effect, etc. This paper gives a brief overview of our development of the HEOM method for fermionic open systems, as well as its applications to various strongly correlated quantum impurity systems.

On the exact truncation tier of fermionic hierarchical equations of motion.

The hierarchical equations of motion (HEOM) theory is in principle exact for describing the dissipative dynamics of quantum systems linearly coupled to Gaussian environments. In practice, the hierarchy needs to be truncated at a finite tier. We demonstrate that, for general systems described by the fermionic HEOM, the (n+L̃)th-tier truncation with L̃=2NσNν yields the exact density operators up to the nth tier. Here, Nσ = 2 for fermionic systems and Nν is the system degrees of freedom. For noninteracting systems, L̃ is further reduced by half. Such an exact termination pattern originates from the Pauli exclusion principle for fermions, and it holds true regardless of the system-environment coupling strength, the number of coupling reservoirs, or the specific scheme employed to unravel the environment memory contents. The relatively small L̃ emphasizes the nonperturbative nature of the HEOM theory. We also propose a simplified HEOM approach to further reduce the memory cost for practical calculations.

Solvent-induced polarization dynamics and coherent two-dimensional spectroscopy: Dissipaton equation of motion approach

This work exploits the dissipaton equation of motion (DEOM) approach to study the solvent-induced non-Condon polarization dynamics and its two-dimensional coherent spectroscopy for model excitonic systems. DEOM is a second quantization generation of the celebrated hierarchical equations of motion theory, with the capability of accurate evaluation for both reduced system and hybrid bath dynamics. The underlying dissipaton algebra includes the generalized Wick’s theorem and DEOM-space techniques for hybrid bath dynamics. Methods discussed also involve an efficient derivative–resum level truncation scheme that preserves the prescription invariance, and the mixed Heisenberg–Schrödinger picture approach for efficient simulation of nonlinear response functions. Our model simulations show clearly the correlated system-and-bath interference and the resulting Fano entanglement spectroscopy profiles.

Optimizing hierarchical equations of motion for quantum dissipation and quantifying quantum bath effects on quantum transfer mechanisms

We present an optimized hierarchical equations of motion theory for quantum dissipation in multiple Brownian oscillators bath environment, followed by a mechanistic study on a model donor-bridge-acceptor system. We show that the optimal hierarchy construction, via the memory-frequency decomposition for any specified Brownian oscillators bath, is generally achievable through a universal pre-screening search. The algorithm goes by identifying the candidates for the best be just some selected Padé spectrum decomposition based schemes, together with a priori accuracy control criterions on the sole approximation, the white-noise residue ansatz, involved in the hierarchical construction. Beside the universal screening search, we also analytically identify the best for the case of Drude dissipation and that for the Brownian oscillators environment without strongly underdamped bath vibrations. For the mechanistic study, we quantify the quantum nature of bath influence and further address the issue of localization versus delocalization. Proposed are a reduced system entropy measure and a state-resolved constructive versus destructive interference measure. Their performances on quantifying the correlated system-environment coherence are exemplified in conjunction with the optimized hierarchical equations of motion evaluation of the model system dynamics, at some representing bath parameters and temperatures. Analysis also reveals the localization to delocalization transition as temperature decreases.

Optimized hierarchical equations of motion theory for Drude dissipation and efficient implementation to nonlinear spectroscopies

Hierarchical equations of motion theory for Drude dissipation is optimized, with a convenient convergence criterion proposed in advance of numerical propagations. The theoretical construction is on the basis of a Padé spectrum decomposition that has been qualified to be the best sum-over-poles scheme for quantum distribution function. The resulting hierarchical dynamics under the a priori convergence criterion are exemplified with a benchmark spin-boson system, and also the transient absorption and related coherent two-dimensional spectroscopy of a model exciton dimer system. We combine the present theory with several advanced techniques such as the block hierarchical dynamics in mixed Heisenberg-Schrödinger picture and the on-the-fly filtering algorithm for the efficient evaluation of third-order optical response functions.