Efficient steady-state solver for hierarchical quantum master equations

Theory and quantum dynamics simulations of exciton-polariton motional narrowing

The motional narrowing effect has been extensively studied for cavity exciton-polariton systems in recent decades both experimentally and theoretically, which is featured by (1) the subaverage behavior and (2) the asymmetric linewidths for the upper polariton and the lower polariton. However, a minimal theoretical model that is clear and adequate to address all these effects as well as the linewidth scaling relations remains missing. In this work, based on the single mode 1D Holstein-Tavis-Cummings (HTC) model, we studied the motional narrowing effect of the polariton linear absorption spectra via both semi-analytic derivations and numerically exact quantum dynamics simulations using the hierarchical equations of motion approach. The results reveal that under collective light-matter coupling between a cavity mode and N molecules, the polariton linewidth scales as 1/N under the slow limit, while scales as 1/N under the fast limit, due to the polaron decoupling effect. Furthermore, by varying the detunings, the polariton linewidths exhibit significant motional narrowing, covering both characters mentioned above. Our analytic linewidth expressions [Eqs. (34) and (35)] agree well with the numerical exact simulations in all the parameter regimes we explored. These results indicate that the physics of motional narrowing is adequately accounted for by the single-mode 1D HTC model. We envision that both the numerical results and the analytic polariton linewidths expression presented in this work will offer great theoretical value for providing a better understanding of the exciton-polariton motional narrowing based on the HTC model.

Application of the imaginary time hierarchical equations of motion method to calculate real time correlation functions

We investigate the application of the imaginary time hierarchical equations of motion method to calculate real time quantum correlation functions. By starting from the path integral expression for the correlated system–bath equilibrium state, we first derive a new set of equations that decouple the imaginary time propagation and the calculation of auxiliary density operators. The new equations, thus, greatly simplify the calculation of the equilibrium correlated initial state that is subsequently used in the real time propagation to obtain the quantum correlation functions. It is also shown that a periodic decomposition of the bath imaginary time correlation function is no longer necessary in the new equations such that different decomposition schemes can be explored. The applicability of the new method is demonstrated in several numerical examples, including the spin-Boson model, the Holstein model, and the double-well model for proton transfer reaction.

Nonequilibrium reaction rate theory: Formulation and implementation within the hierarchical equations of motion approach

The study of chemical reactions in environments under nonequilibrium conditions has been of interest recently in a variety of contexts, including current-induced reactions in molecular junctions and scanning tunneling microscopy experiments. In this work, we outline a fully quantum mechanical, numerically exact approach to describe chemical reaction rates in such nonequilibrium situations. The approach is based on an extension of the flux correlation function formalism to nonequilibrium conditions and uses a mixed real and imaginary time hierarchical equations of motion approach for the calculation of rate constants. As a specific example, we investigate current-induced intramolecular proton transfer reactions in a molecular junction for different applied bias voltages and molecule-lead coupling strengths.

Strong Exciton-Vibrational Coupling in Molecular Assemblies. Dynamics Using the Polaron Transformation in HEOM Space

In Frenkel exciton dynamics of aggregated molecules, the polaron transformation (PT) technique leads to decoupling of diagonal elements in the subspace of excited electronic states from vibrations. In this article we describe for the first time how PT becomes applicable in the framework of the "Hierarchical Equations of Motion" (HEOM) approach for treatment of open quantum systems. We extend the concept of formulating operators in HEOM space by deriving hierarchical equations of PT which lead to a shift in the excited state potential energy surface to compensate its displacement. While the assumption of thermal equilibration of the vibrational oscillators, introduced by PT, results in a stationary state in a monomer, in a dimer under the same assumption nonequilibrium dynamics appears because of the interplay of the transfer process and vibrational equilibration. Both vertical transitions generating a vibrationally hot state and initially equilibrated vibrational oscillators evolve toward the same stationary asymptotic state associated with polaron formation. The effect of PT on the dynamics of this process depends on initial excitation and basis representation of the electronic system. The developed approach facilitates a generic formulation of quantum master equations involving perturbative treatment of polaron dynamics.

Optical response of laser-driven charge-transfer complex described by Holstein-Hubbard model coupled to heat baths: Hierarchical equations of motion approach

We investigate the optical response of a charge-transfer complex in a condensed phase driven by an external laser field. Our model includes an instantaneous short-range Coulomb interaction and a local optical vibrational mode described by the Holstein-Hubbard (HH) model. Although characterization of the HH model for a bulk system has typically been conducted using a complex phase diagram, this approach is not sufficient for investigations of dynamical behavior at finite temperature, in particular for studies of nonlinear optical properties, where the time irreversibility of the dynamics that arises from the environment becomes significant. We therefore include heat baths with infinite heat capacity in the model to introduce thermal effects characterized by fluctuation and dissipation to the system dynamics. By reducing the number of degrees of freedom of the heat baths, we derive numerically "exact" hierarchical equations of motion for the reduced density matrix of the HH system. As demonstrations, we calculate the optical response of the system in two- and four-site cases under external electric fields. The results indicate that the effective strength of the system-bath coupling becomes large as the number of sites increases. Excitation of electrons promotes the conductivity when the Coulomb repulsion is equivalent to or dominates the electron-phonon coupling, whereas excitation of optical vibrations always suppresses the conductivity.

Nonequilibrium Kondo regime current noise spectrum of quantum dot systems with the single impurity Anderson model

We investigate the nonequilibrium current noise spectrum of single impurity Anderson model quantum dot systems on the basis of the accurate dissipation equation of motion evaluations. By comparing between the equilibrium and nonequilibrium cases and between the non-Kondo and Kondo regimes, we identify the current noise spectrum of the nonequilibrium Kondo features that actually appear in the entire region of ω ∈ [-eV, eV]. It is well known that the primary Kondo characteristics at ω = ±eV = ±(μ_L - μ_R) display asymmetrical upturns and remarkable peaks in S(ω) and dS(ω)/dω, respectively. These features are originated from the Rabi interference of the transport current dynamics, with the Kondo oscillation frequency of |eV|. Moreover, we also identify the minor but very distinguishable inflections, crossing over from ω = -eV to ω = +eV. This uncovered feature would be related to the interference between two Kondo resonance channels.

Correlated vibration-solvent effects on the non-Condon exciton spectroscopy

Excitation energy transfer is crucially involved in a variety of systems. During the process, the non-Condon vibronic coupling and the surrounding solvent interaction may synergetically play important roles. In this work, we study the correlated vibration-solvent influences on the non-Condon exciton spectroscopy. Statistical analysis is elaborated for the overall vibration-plus-solvent environmental effects. Analytic solutions are derived for the linear absorption of monomer systems. General simulations are accurately carried out via the dissipaton-equation-of-motion approach. The resulted spectra in either the linear absorption or strong field regime clearly demonstrate the coherence enhancement due to the synergetic vibration-solvent correlation.

Efficient Steady-State Solver for the Hierarchical Equations of Motion Approach: Formulation and Application to Charge Transport through Nanosystems

An iterative approach is introduced, which allows the efficient solution of the hierarchical equations of motion (HEOM) for the steady-state of open quantum systems. The approach combines the method of matrix equations with an efficient preconditioning technique to reduce the numerical effort of solving the HEOM. Illustrative applications to simulate nonequilibrium charge transport in single-molecule junctions demonstrate the performance of the method.

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.

Nonequilibrium steady-state picture of incoherent light-induced excitation harvesting

We formulate a comprehensive theoretical description of excitation harvesting in molecular aggregates photoexcited by weak incoherent radiation. An efficient numerical scheme that respects the continuity equation for excitation fluxes is developed to compute the nonequilibrium steady state (NESS) arising from the interplay between excitation generation, excitation relaxation, dephasing, trapping at the load, and recombination. The NESS is most conveniently described in the so-called preferred basis in which the steady-state excitonic density matrix is diagonal. The NESS properties are examined by relating the preferred-basis description to the descriptions in the site or excitonic bases. Focusing on a model photosynthetic dimer, we find that the NESS in the limit of long trapping time is quite similar to the excited-state equilibrium in which the stationary coherences originate from the excitation-environment entanglement. For shorter trapping times, we demonstrate how the properties of the NESS can be extracted from the time-dependent description of an incoherently driven but unloaded dimer. This relation between stationary and time-dependent pictures is valid, provided that the trapping time is longer than the decay time of dynamic coherences accessible in femtosecond spectroscopy experiments.

Equilibrium and transient thermodynamics: A unified dissipaton-space approach

This work presents a unified dissipaton-equation-of-motion (DEOM) theory and its evaluations on the Helmholtz free energy change due to the isotherm mixing of two isolated subsystems. One is a local impurity, and the other is a nonlocal Gaussian bath. DEOM constitutes a fundamental theory for such open quantum mixtures. To complete the theory, we also construct the imaginary-time DEOM formalism via an analytical continuation of dissipaton algebra, which would be limited to equilibrium thermodynamics. On the other hand, the real-time DEOM deals with both equilibrium structural and nonequilibrium dynamic properties. Its combination with the thermodynamic integral formalism would be a viable and accurate means to both equilibrium and transient thermodynamics. As illustrations, we report the numerical results on a spin-boson system, with elaborations on the underlying anharmonic features, the thermodynamic entropy vs the von Neumann entropy, and an indication of "solvent-cage" formation. Beside the required asymptotic equilibrium properties, the proposed transient thermodynamics also supports the basic spontaneity criterion.

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

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.