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.
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.
By using the quasi-equilibrium Helmholtz energy, which is defined as the thermodynamic work in a quasi-static process, we investigate the thermal properties of both an isothermal process and a transition process between the adiabatic and isothermal states (adiabatic transition). Here, the work is defined by the change in energy from a steady state to another state under a time-dependent perturbation. In particular, the work for a quasi-static change is regarded as thermodynamic work. We employ a system-bath model that involves time-dependent perturbations in both the system and the system-bath interaction. We conduct numerical experiments for a three-stroke heat machine (a Kelvin-Planck cycle). For this purpose, we employ the hierarchical equations of motion (HEOM) approach. These experiments involve an adiabatic transition field that describes the operation of an adiabatic wall between the system and the bath. Thermodynamic-work diagrams for external fields and their conjugate variables, similar to the P-V diagram, are introduced to analyze the work done for the system in the cycle. We find that the thermodynamic efficiency of this machine is zero because the field for the isothermal processes acts as a refrigerator, whereas that for the adiabatic wall acts as a heat engine. This is a numerical manifestation of the Kelvin-Planck statement, which states that it is impossible to derive the mechanical effects from a single heat source. These HEOM simulations serve as a rigorous test of thermodynamic formulations because the second law of thermodynamics is only valid when the work involved in the operation of the adiabatic wall is treated accurately.
A system--bath (SB) model is considered to examine the Jarzynski equality in the fully quantum regime. In our previous paper [J. Chem. Phys. 153 (2020) 234107], we carried out "exact" numerical experiments using hierarchical equations of motion (HEOM) in which we demonstrated that the SB model describes behavior that is consistent with the first and second laws of thermodynamics and that the dynamics of the total system are time irreversible. The distinctive quantity in the Jarzynski equality is the "work characteristic function (WCF)", ⟨exp[−βW]⟩, where W is the work performed on the system and β is the inverse temperature. In the present investigation, we consider the definitions based on the partition function (PF) and on the path, and numerically evaluate the WCF using the HEOM to determine a method for extending the Jarzynski equality to the fully quantum regime. We show that using the PF-based definition of the WCF, we obtain a result that is entirely inconsistent with the Jarzynski equality, while if we use the path-based definition, we obtain a result that approximates the Jarzynski equality, but may not be consistent with it.
We investigate the efficiency of a quantum Carnot engine based on open quantum dynamics theory. The model includes time-dependent external fields for the subsystems controlling the isothermal and isentropic processes and for the system-bath (SB) interactions controlling the transition between these processes. Numerical simulations are conducted in a nonperturbative and non-Markovian SB coupling regime by using the hierarchical equations of motion under these fields at different cycle frequencies. The work applied to the total system and the heat exchanged with the baths are rigorously evaluated. In addition, by regarding quasi-static work as free energy, we compute the quantum thermodynamic variables and analyze the simulation results by using thermodynamic work diagrams for the first time. Analysis of these diagrams indicates that, in the strong SB coupling region, the fields for the SB interactions are major sources of work, while in other regions, the field for the subsystem is a source of work. We find that the maximum efficiency is achieved in the quasi-static case and is determined solely by the bath temperatures, regardless of the SB coupling strength, which is a numerical manifestation of Carnot's theorem.
We investigate the dynamics of a two-qubit open quantum system, in particular the two-qubit spin-boson model in the strong coupling regime, coupled to two thermal bosonic baths under non-Markovian and non-equilibrium conditions. Two complementary approaches, the Hierarchical Equations of Motion (HEOM) and Reaction Coordinate Mapping (RCM), are employed to examine various coupling regimes between the qubits and their respective baths. The dynamical features of the model and the impact of the tunneling amplitude on quantum coherence of the system are probed using the l1-norm of coherence. The model is further shown to have non-Markovian evolution. The nontrivial task of calculating entropy production in the strong-coupling regime is performed using auxiliary density operators in HEOM. Motivated by the realization of a quantum thermal device in the strong-coupling regime, the non-equilibrium steady-state behavior of the system is investigated. Furthermore, the relationship between the heat and spin currents and the tunneling amplitude is probed.
A central open question in quantum thermodynamics is the identification of thermodynamic entropy away from equilibrium. Although thermodynamic entropy coincides with the von Neumann entropy for a Gibbs state at equilibrium, the thermodynamic role of quantum coherence generated during nonequilibrium evolution remains unclear. Here, we show that von Neumann entropy does not, in general, represent the thermodynamic entropy governing heat exchange in coherent nonequilibrium quantum systems. By integrating the resource theory of quantum coherence with thermodynamic principles, we develop a general framework for nonequilibrium quantum thermodynamics. We identify thermodynamic entropy with the energy entropy derived from energy-measurement statistics and introduce a dynamical temperature extending the equilibrium concept to nonequilibrium states. Consequently, we obtain an exact entropy-balance relation separating the entropy flux associated with heat exchange from the coherence-loss contribution to entropy production, establishing an operational connection between quantum coherence, heat flow, and total entropy change. Unlike approaches relying on equilibrium reference states or stationary Gibbs ensembles, our framework uses instantaneous dynamical quantities and remains valid throughout transient nonequilibrium evolution. It provides consistent dynamical definitions of temperature, work, heat, and free energy. Applying the formalism to an exactly solvable open quantum system, we show that the first and second laws hold throughout transient evolution and that equilibrium thermodynamics emerges dynamically in the weak-coupling long-time limit. Our results establish a unified framework for coherent nonequilibrium quantum thermodynamics and suggest opportunities for controlling heat dissipation through coherence engineering in quantum thermal technologies.
Thermodynamics is a phenomenological theory describing the energy conversion of work and heat. At its origins thermodynamics was developed in order to understand and improve heat engines. In conventional thermodynamics, however, only such processes are completely describable which are slow enough to keep the system of interest in an equilibrium state with its thermal environment at all times. On the contrary, all real physical processes are accompanied by non-equilibrium phenomena. These are mathematically described with the help of the irreversible entropy production, i.e. the dissipated work divided by the temperature of the bath. The second law of thermodynamics merely states that the irreversible entropy production is on average always positive. The modern trend of miniaturization, however, leads to smaller and smaller devices. On short length scales thermal noise as well as quantum fluctuations become important. Thus, usual thermodynamic quantities as work and heat acquire stochastic nature. Moreover, in the quantum regime a completely new theory had to be invented, since classical notions of work and heat are no longer valid. The present dissertation contributes to this prevailing field by the derivation of analytical expressions for the entropy production in open and closed quantum system far from thermal equilibrium. The theoretical treatment is motivated by an experimental point of view. The obtained results were derived from the reduced dynamics of the quantum mechanical system, which may be coupled to a thermal bath. To this end, it was dealt with methods and approaches of statistical physics, conventional thermodynamics, quantum information theory and the theory of open quantum systems. In this dissertation the most important physical systems are treated. First, an isolated system is considered before later on a heat bath is included. The limit of weak as well as of strong coupling are analyzed in detail. For isolated systems it turns out that a geometric approach is able to capture the thermodynamic properties, completely. To this end, the irreversible entropy production is sharply estimated from below by the angle between the real non-equilibrium state and the corresponding equilibrium one. Further, it is shown that the Heisenberg uncertainty relation for energy and time can be formulated more precisely in terms of that angle. As experimental system cold ion traps are introduced, in which the theoretical predictions can be verified. For weakly coupled systems, then, the irreversible entropy production is identified and a fluctuation theorem is derived by merely thermodynamic arguments. As experimental system we, again, propose the cold ion traps, where the heat bath can be simulated with the help of LASER light. In the limit of strong coupling a semiclassical description of the reduced quantum system is discussed. Again, an analytically exact expression for the entropy production and the corresponding fluctuation theorem is derived. A physical system, in which strong coupling to the environment can be analyzed, is given by Josephson junctions. This dissertation discusses the entropy production for isolated, weakly coupled, and overdamped quantum systems. For these physical regimes analytical relation are derived and experimental systems proposed, with which the predictions can be verified.