In the previous paper [C. W. Kim and I. Franco, J. Chem. Phys. 160, 214111-1–214111-13 (2024)], we developed a theory called MQME-D, which allows us to decompose the overall energy dissipation process in open quantum system dynamics into contributions by individual components of the bath when the subsystem dynamics is governed by a Markovian quantum master equation (MQME). Here, we contrast the predictions of MQME-D against the numerically exact results obtained by combining hierarchical equations of motion (HEOM) with a recently reported protocol for monitoring the statistics of the bath. Overall, MQME-D accurately captures the contributions of specific bath components to the overall dissipation while greatly reducing the computational cost compared to exact computations using HEOM. The computations show that MQME-D exhibits errors originating from its inherent Markov approximation. We demonstrate that its accuracy can be significantly increased by incorporating non-Markovianity by exploiting time scale separations (TSS) in different components of the bath. Our work demonstrates that MQME-D combined with TSS can be reliably used to understand how energy is dissipated in realistic open quantum system dynamics.
more »
« less
General framework for quantifying dissipation pathways in open quantum systems. I. Theoretical formulation
We present a general and practical theoretical framework to investigate how energy is dissipated in open quantum system dynamics. This is performed by quantifying the contributions of individual bath components to the overall dissipation of the system. The framework is based on the Nakajima–Zwanzig projection operator technique, which allows us to express the rate of energy dissipation into a specific bath degree of freedom by using traces of operator products. The approach captures system-bath interactions to all orders, but is based on second-order perturbation theory on the off-diagonal subsystem's couplings and a Markovian description of the bath. The usefulness of our theory is demonstrated by applying it to various models of open quantum systems involving harmonic oscillators or spin baths and connecting the outcomes to existing results such as our previously reported formula derived for locally coupled harmonic baths [Kim and Franco, J. Chem. Phys. 154, 084109 (2021)]. We also prove that the dissipation calculated by our theory rigorously satisfies thermodynamic principles such as energy conservation and detailed balance. Overall, the strategy can be used to develop the theory and simulation of dissipation pathways to interpret and engineer the dynamics of open quantum systems.
more »
« less
- PAR ID:
- 10516898
- Publisher / Repository:
- AIP
- Date Published:
- Journal Name:
- The Journal of Chemical Physics
- Volume:
- 160
- Issue:
- 21
- ISSN:
- 0021-9606
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
An efficient, fully quantum mechanical real-time path integral method for including the effects of static disorder in the dynamics of systems coupled to common or local harmonic baths is presented. Rather than performing a large number of demanding calculations for different realizations of the system Hamiltonian, the influence of the bath is captured through a single evaluation of the path sum by grouping the system paths into equivalence classes of fixed system amplitudes. The method is illustrated with several analytical and numerical examples that show a variety of nontrivial effects arising from the interplay among coherence, dissipation, thermal fluctuations and geometric phases.more » « less
-
Abstract Conical intersections in two-state systems require a coordinate-dependent coupling. This paper identifies and investigates conical intersections in cyclic tight-binding system-bath Hamiltonians with an odd number of sites and a constant site-to-site coupling. In the absence of bath degrees of freedom, such tight-binding systems with a positive coupling parameter exhibit electronic frustration and a doubly-degenerate ground state. When these systems interact with a harmonic bath, the degeneracy becomes a conical intersection between the adiabatic ground and first excited states. Under weak system-bath coupling, overlapping wavefunctions associated with different sites give rise to distinct pathways with interfering geometric phases, which lead to considerably slower transfer dynamics. The effect is most pronounced in the presence of low-temperature dissipative baths characterized by a continuous spectral density. It is found that the transfer dynamics and equilibration time of a cyclic dissipative three-site system with a positive coupling exceeds that of a similar three-site system with a negative coupling, as well as that of cyclic four-site systems, by an order of magnitude.more » « less
-
null (Ed.)Extending our previous work, quantum dynamic simulations are performed to study low temperature heat transport in a spin-boson model where a two-level subsystem is coupled to two independent harmonic baths. Multilayer multiconfiguration time-dependent Hartree theory is used to numerically evaluate the thermal flux, for which the bath is represented by hundreds to thousands of modes. The simulation results are compared with the approximate Redfield theory approach, and the physics is analyzed versus different physical parameters.more » « less
-
We develop a quasiparticle approach to capture the dynamics of open quantum systems coupled to bosonic thermal baths of arbitrary complexity based on the Hierarchical Equations of Motion (HEOM). This is done by generalizing the HEOM dynamics and mapping it into that of the system in interaction with a few bosonic fictitious quasiparticles that we call bexcitons. Bexcitons arise from a decomposition of the bath correlation function into discrete features. Specifically, bexciton creation and annihilation couple the auxiliary density matrices in the HEOM. The approach provides a systematic strategy to construct exact quantum master equations that include the system–bath coupling to all orders even for non-Markovian environments. Specifically, by introducing different metrics and representations for the bexcitons it is possible to straightforwardly generate different variants of the HEOM, demonstrating that all these variants share a common underlying quasiparticle picture. Bexcitonic properties, while unphysical, offer a coarse-grained view of the correlated system–bath dynamics and its numerical convergence. For instance, we use it to analyze the instability of the HEOM when the bath is composed of underdamped oscillators and show that it leads to the creation of highly excited bexcitons. The bexcitonic picture can also be used to develop more efficient approaches to propagate the HEOM. As an example, we use the particle-like nature of the bexcitons to introduce mode-combination of bexcitons in both number and coordinate representation that uses the multi-configuration time-dependent Hartree to efficiently propagate the HEOM dynamics.more » « less