Hierarchical Equations of Motion Approach to Quantum Thermodynamics

Quantum thermal transport beyond second order with the reaction coordinate mapping

Standard quantum master equation techniques, such as the Redfield or Lindblad equations, are perturbative to second order in the microscopic system–reservoir coupling parameter λ. As a result, the characteristics of dissipative systems, which are beyond second order in λ, are not captured by such tools. Moreover, if the leading order in the studied effect is higher-than-quadratic in λ, a second-order description fundamentally fails even at weak coupling. Here, using the reaction coordinate (RC) quantum master equation framework, we are able to investigate and classify higher-than-second-order transport mechanisms. This technique, which relies on the redefinition of the system–environment boundary, allows for the effects of system–bath coupling to be included to high orders. We study steady-state heat current beyond second-order in two models: The generalized spin-boson model with non-commuting system–bath operators and a three-level ladder system. In the latter model, heat enters in one transition and is extracted from a different one. Crucially, we identify two transport pathways: (i) System’s current, where heat conduction is mediated by transitions in the system, with the heat current scaling as j q ∝ λ2 to the lowest order in λ. (ii) Inter-bath current, with the thermal baths directly exchanging energy between them, facilitated by the bridging quantum system. To the lowest order in λ, this current scales as j q ∝ λ4. These mechanisms are uncovered and examined using numerical and analytical tools. We contend that the RC mapping brings, already at the level of the mapped Hamiltonian, much insight into transport characteristics.

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.

Exact description of excitonic dynamics in molecular aggregates weakly driven by light

We present a rigorous theoretical description of excitonic dynamics in molecular light-harvesting aggregates photoexcited by weak-intensity radiation of arbitrary properties. While the interaction with light is included up to the second order, the treatment of the excitation-environment coupling is exact and results in an exact expression for the reduced excitonic density matrix that is manifestly related to the spectroscopic picture of the photoexcitation process. This expression takes fully into account the environmental reorganization processes triggered by the two interactions with light. This is particularly important for slow environments and/or strong excitation-environment coupling. Within the exponential decomposition scheme, we demonstrate how our result can be recast as the hierarchy of equations of motion (HEOM) that explicitly and consistently includes the photoexcitation step. We analytically describe the environmental reorganization dynamics triggered by a delta-like excitation of a single chromophore and demonstrate how our HEOM, in appropriate limits, reduces to the Redfield equations comprising a pulsed photoexcitation and the nonequilibrium Förster theory. We also discuss the relation of our formalism to the combined Born-Markov-HEOM approaches in the case of excitation by thermal light.

Numerically "exact" simulations of entropy production in the fully quantum regime: Boltzmann entropy vs von Neumann entropy

We present a scheme to evaluate thermodynamic variables for a system coupled to a heat bath under a time-dependent external force using the quasi-static Helmholtz energy from the numerically "exact" hierarchical equations of motion (HEOM). We computed the entropy produced by a spin system strongly coupled to a non-Markovian heat bath for various temperatures. We showed that when changes to the external perturbation occurred sufficiently slowly, the system always reached thermal equilibrium. Thus, we calculated the Boltzmann entropy and the von Neumann entropy for an isothermal process, as well as various thermodynamic variables, such as changes in internal energies, heat, and work, for a system in quasi-static equilibrium based on the HEOM. We found that although the characteristic features of the system entropies in the Boltzmann and von Neumann cases as a function of the system-bath coupling strength are similar, those for the total entropy production are completely different. The total entropy production in the Boltzmann case is always positive, whereas that in the von Neumann case becomes negative if we chose a thermal equilibrium state of the total system (an unfactorized thermal equilibrium state) as the initial state. This is because the total entropy production in the von Neumann case does not properly take into account the contribution of the entropy from the system-bath interaction. Thus, the Boltzmann entropy must be used to investigate entropy production in the fully quantum regime. Finally, we examined the applicability of the Jarzynski equality.