Quantum Thermodynamics at Strong Coupling: Operator Thermodynamic Functions and Relations

Identifying or constructing a fine-grained microscopic theory that will emerge under specific conditions to a known macroscopic theory is always a formidable challenge. Thermodynamics is perhaps one of the most powerful theories and best understood examples of emergence in physical sciences, which can be used for understanding the characteristics and mechanisms of emergent processes, both in terms of emergent structures and the emergent laws governing the effective or collective variables. Viewing quantum mechanics as an emergent theory requires a better understanding of all this. In this work we aim at a very modest goal, not quantum mechanics as thermodynamics, not yet, but the thermodynamics of quantum systems, or quantum thermodynamics. We will show why even with this minimal demand, there are many new issues which need be addressed and new rules formulated. The thermodynamics of small quantum many-body systems strongly coupled to a heat bath at low temperatures with non-Markovian behavior contains elements, such as quantum coherence, correlations, entanglement and fluctuations, that are not well recognized in traditional thermodynamics, built on large systems vanishingly weakly coupled to a non-dynamical reservoir. For quantum thermodynamics at strong coupling, one needs to reexamine the meaning of the thermodynamic functions, the viability of the thermodynamic relations and the validity of the thermodynamic laws anew. After a brief motivation, this paper starts with a short overview of the quantum formulation based on Gelin & Thoss and Seifert. We then provide a quantum formulation of Jarzynski’s two representations. We show how to construct the operator thermodynamic potentials, the expectation values of which provide the familiar thermodynamic variables. Constructing the operator thermodynamic functions and verifying or modifying their relations is a necessary first step in the establishment of a viable thermodynamics theory for quantum systems. We mention noteworthy subtleties for quantum thermodynamics at strong coupling, such as in issues related to energy and entropy, and possible ambiguities of their operator forms. We end by indicating some fruitful pathways for further developments.

Why is a paper on this subject matter appearing in this special issue of Entropy? The short answer is the same reason why EmQM appears in this journal Entropy, which is generally considered as treating topics in statistical mechanics: emergence. The long answer, serving as a justification for our dwelling on quantum thermodynamics in the realm of emergent quantum mechanics, is as follows. One way to see quantum mechanics as emergent is by analogy with hydrodynamics and thermodynamics, probably the two best known emergent theories because we know exactly what the collective variables are (thermodynamic functions and relations), the laws they obey (the four laws), how they are related to the basic constituents (molecular dynamics) and the mediating theory from which we can derive both the fundamental and the emergent theories (kinetic theory).
In analogy to emergent gravity [1][2][3][4][5], one of the present authors has championed the thesis of "general relativity as geometro-hydrodynamics" [6,7]. Since Verlinde's "gravity is entropic force" popularization of Jacobson's "Einstein equation of state" thesis [8][9][10], gravity as thermodynamics has caught a wider attention in the quantum gravity community. However, a new challenge arises. The question one of us posed for enthusiasts of this theme is the following: Since both gravity and thermodynamics are old subjects established centuries before the advent of quantum mechanics, and both can make sense and stand alone at the classical level without quantum theory, well, what exactly is quantum doing here?-What is the role of quantum in emergent gravity? Do we really need quantum if we view "gravity as thermodynamics"?

Quantum in "Gravity as Thermodynamics"
This question, "Wither the Quantum?", as Hu calls it, puts the spotlight on quantum, in how it contributes to the emergent phenomena which gives us both gravity and thermodynamics. One answer to this is to also consider quantum mechanics as emergent. For example, in the emergent theories of Adler and 't Hooft [11,12] probability theory, stochastic and statistical mechanics as a slate play a pivotal role, just as they do for thermodynamics and hydrodynamics arising from molecular dynamics. The same applies to emergent gravity: classical gravity captured by general relativity is an effective theory emergent from some fundamental theories of the basic constituents of spacetime functioning at the sub-Planckian scale. How these basic constituents interact, how their interaction strength varies with energy, how at some specific scale(s) some set(s) of collective variables and the law(s) governing them emerge, and in succession, leading to the physics at the lowest energy as we know it in today's universe is perhaps just as interesting as the manifestation of the relevant physics at the different scales familiar to us-from molecules to atoms to nucleons to quarks and below. Putting aside gravity for now in this investigation, we wish to see a deeper connection between micro and macro, quantum and thermo.

David Bohm: Quantum in Classical Terms
Here, Bohm's philosophical influence is evident. His pilot wave theory may not offer a better description or explanation of quantum phenomena, but the view that quantum mechanics is not a fundamental theory any more than a classical wave theory is, provides an inspiration for asking a deeper layer of questions: If we view quantum theory functioning in the capacity as thermodynamics, we should ask: What are the fundamental constituents, the laws governing them, and how quantum mechanics emerges from the sub-structures and theories depicting them (Long before we get to this point, many readers may have raised this objection: This is obviously nonsense: The second law of thermodynamics ostensibly shows the effects of an arrow of time, while quantum mechanics is time-reversal invariant. Well, if the mechanical processes which we can observe are in the underdamped regime where the dissipative effects are not strong enough, they would appear to obey time-reversal symmetry. This is not an outlandish explanation: For most physical systems, in the open system perspective, quantum phenomena in the system appears within the decoherence time which is many orders of magnitude shorter than the relaxation time, as is the case in many well-controlled environments (e.g., cavity QED). Or, if the system is near a nonequilibrium critical point. On this issue, cosmology, despite its seemingly remote bearing, may actually enter in a basic way, in terms of the origin of the arrow of time, and the mere fact that nonequilibrium conditions prevail in an expanding universe.).

Quantum as Thermodynamics?
In quantum thermodynamics, we may not see much in terms of what fundamental theory quantum mechanics emerges from (Adler and 't Hooft may have their answers: trace dynamics and cellular automaton, for instance, respectively), but even the juxtaposition or crossing of what is traditionally considered as governing the two opposite ends in the macro/micro and classical/quantum spectrum may reveal some deeper meaning in both. Macroscopic quantum phenomena is another such arena. For example, in small quantum systems, at low temperatures, or when the system is strongly coupled to its environment, is there a lower limit to the validity of the laws of thermodynamics, which play such an important role in our understanding of the macro world? Under what conditions will macroscopic entities show quantum phenomena? Is there an upper limit to quantum mechanics governing the meso domain? Is there a limit to quantum commanding the macro world? The above explains the philosophical issues which motivated us to take up a study of quantum thermodynamics. We are also of the opinion that useful philosophical discourses of any subject matter should be based on the hard-core scientific knowledge of that subject, down to all the nitty-gritty details of each important topic that makes up that body of knowledge. Thus we start with the basic demands in the formulation of quantum thermodynamics and try to meet them in a rigorous, no-nonsense way. The specific goal of this paper is to define the operator thermodynamic quantities and spell out their relations for quantum many-body systems in thermal equilibrium.

Quantum Thermodynamics
Quantum thermodynamics is a fast developing field, emergent from quantum many body physics and nonequilibrium statistical mechanics. Simply, it is the study of the thermodynamic properties of quantum many-body systems. Quantum now refers not just to the particle spin-statistics (boson vs. fermion) aspects in traditional quantum statistical mechanics, but also includes in the present era the quantum phase aspects, such as quantum coherence, quantum correlations, and quantum entanglement, where quantum information enters. The new challenges arise from several directions not falling under the assumptions of traditional classical thermodynamics: finding the quantum properties of small systems, at zero or very low temperature, strongly coupled to an environment, which could have non-ohmic spectral densities and colored noise, while the system evolves following a non-Markovian dynamics (with memory).

This Work
In this paper we discuss the issues and the technical challenges encountered in the first stage in the construction of a viable theory of quantum thermodynamics, where the system is strongly coupled with a heat bath. We wish to present in a systematic way how to introduce the operator thermodynamic functions and construct their relations. Here, we treat this problem in an equilibrium setting. There are other ways to construct such a theory, such as pursued in the so-called "eigenvalue thermalization" program [13][14][15][16][17][18][19][20][21][22], which treats the system and the environment as a closed system, or along the lines expounded in the fluctuation theorems [23][24][25] where the system under an external drive is allowed to evolve in a nonequilibrium albeit controlled manner [26,27], or in a fully nonequilibrium open quantum system dynamics program (see, e.g., [28,29] and references therein).

New Challenges in Quantum Thermodynamics
Small quantum many-body systems strongly coupled to a heat reservoir at low temperatures are the new focuses of interest for quantum thermodynamics [30][31][32]. Under these hitherto lesser explored conditions, one needs to re-examine the meaning of the thermodynamic functions, the viability of the thermodynamic relations and the validity of the thermodynamic laws anew. Traditional thermodynamics is built on large systems weakly coupled to a reservoir [33], and for quantum systems, only the spin-statistics aspect is studied, as in quantum statistical mechanics, leaving the important factors of quantum coherence, correlations, entanglement and fluctuations as new challenges for quantum thermodynamics. To see the difference strong coupling makes, we take as an example the definition of heat, the energy transferred between the system and the reservoir, for systems strongly coupled to a bath. Esposito et al. [34], for example, show that any heat definition expressed as an energy change in the reservoir energy plus any fraction of the system-reservoir interaction is not an exact differential when evaluated along reversible isothermal transformations, except when that fraction is zero. Even in that latter case the reversible heat divided by temperature, namely entropy, does not satisfy the third law of thermodynamics and diverges in the low temperature limit. For quantum systems, as pointed out by Ankerhold and Pekola [35], in actual measurements, especially for solid state structures, quantum correlations between system and reservoir may be of relevance not only far from but also close to and in thermal equilibrium. Even in the weak coupling regime, this heat flow is substantial at low temperatures and may become comparable to typical predictions for the work based on conventional weak coupling approaches. It further depends sensitively on the non-Markovian features of the reservoir. These observations exemplify the intricacies involved in defining heat for strongly-coupled systems and added complexities for quantum systems, especially at low temperatures.
This incertitude regarding heat translates to ambiguity in the definition of thermodynamic functions and the thermodynamic relations. For example, it was shown [36][37][38][39][40][41] that the expressions for the specific heat derived from the internal energies of a quantum-mechanical harmonic oscillator bilinearly coupled to a harmonic bath calculated by two different approaches can have dramatically different behavior in the low temperature regime. To illustrate this point, Gelin and Thoss [42] compared these two approaches of calculating the internal energy of the system, which give identical results if the system-bath coupling is negligible, but predict significantly differently for finite system-bath coupling. In the first approach, the mean energy of the system given by the expectation value of the system Hamiltonian is evaluated with respect to the total (system + bath) canonical equilibrium distribution. The second approach is based on the partition function of the system, Z s , which is postulated to be given as the ratio of the total (system + bath) and the bath partition functions. Gelin and Thoss [42] introduce a bath-induced interaction operator∆ s , which would account for the effects of finite system-bath coupling and analyze the two approaches for several different systems including several quantum and classical point particles and nonlinear system bath coupling. They found that Approach II leads to very different results from Approach I, their differences exist already within classical mechanics, provided the system-bath interaction is not bilinear and/or the system of interest consists of more than a single particle. Similar ambiguity appears in the entropy of the quantum system in the same setup. In the first approach, the von Neumann entropy is chosen to be the entropy of the system, while in the second approach, the entropy is given by the derivative of the system's free energy with respect to the inverse temperature. Both definitions are equivalent in the limit of weak system-bath coupling. However, at finite coupling strength, it has been noticed that [36,37,43,44] the von Neumann entropy may not vanish when the bath temperature is close to zero, even for a simple quantum system that consists of harmonic oscillators. This nonvanishing behavior of the von Neumann entropy is related to the quantum entanglement between the system and the bath [45][46][47][48]. On the other hand, the entropy defined in the second approach of the same system gives an expected vanishing result, consistent with the physical picture described by [49]. These disparities between different choices of the same thermodynamic function of the system can be traced to the degree in which the system-bath interaction enters into or counted as the system.
In a small quantum system, quantum fluctuations of a physical observable of the system can grow to the same order of magnitude as its average value. Thus a description of the thermodynamics of such a system based on the mean value is insufficient. To include the effects of quantum fluctuations of the thermodynamic variables, we need their corresponding quantum operators, from which we can calculate its higher moments. Likewise for quantities related to the quantum correlations and quantum coherence of such a system. Constructing the operator form of a thermodynamic variable is nontrivial and ambiguous at times. Take as an example the internal energy of the system. As is mentioned earlier, it can be given either as the expectation value of the system's Hamiltonian operator, or as the partial derivative of the (effective) partition function of the system. In the former case, the operator of the internal energy is already chosen but with arbitrariness, while in the latter case, as will be discussed in Section 7.3, our capabilities of identifying the corresponding operator is limited by our knowledge of the thermodynamic variables based on either its expectation or its classical counterpart. Either way there exists non-uniqueness in the introduction of the thermodynamic operator form of the system variables.
Even with a successful extraction and identification of the operator form from the expectation value or the partial derivative of the system's partition, such a thermodynamic operator for the reduced system in general acts on the Hilbert space of the full composite (system + environment). This is not always a desirable feature because most of the time we are interested in what happens to the system under the influence of the environment, whose details are of no particular concern. A description of the full composite does not easily translate to useful information about the system. Thus it is more preferable to derive operators of the system which solely act on the Hilbert space of the system, without reference to the environment. This is what we try to accomplish here; it is different from what has been reported in the literature.

Goal and Findings of Present Work
In traditional weak coupling thermodynamics, knowing the thermodynamic potentials and their relations enables one to construct a theory of thermodynamics. We believe they are also necessary for the establishment of a theory of thermodynamics for quantum systems. The new challenge is twofold: constructing operator thermodynamic potentials and treating quantum systems strong coupled to their baths. This is our goal in this paper. We shall provide a quantum formulation of Jarzynski's [27] classical thermodynamics at strong coupling. In contrast to the study of the nonequilibrium dynamics of an open quantum system (called ONEq in [29,50,51]), the combined system + environment studied here is a closed system, called the composite, which is assumed to stay in a global thermal state. In such a configuration, even though the interaction between the system and the bath is non-negligible, the partition function of the composite can be defined. This facilitates the introduction of thermodynamic potentials in a way similar to the traditional vanishing-coupling thermodynamics. We will focus on thermodynamic functions such as internal energy, enthalpy, entropy, and free energies, but exclude the consideration of heat and work in this paper, as quantum work is not well understood and requires a separate treatment.
This paper is organized as follows: In Section 3, we give a quick summary of the familiar thermodynamic relations, which we call traditional or weak-coupling thermodynamics, if only to establish notations. We then consider interacting quantum systems with the help of the Hamiltonian of mean force [52][53][54][55][56][57][58] and discuss some conceptual and technical difficulties we may face when the coupling strength becomes strong. In Sections 4-6, we overview the equilibrium quantum formulations of thermodynamics at strong coupling, based on Gelin & Thoss' and Seifert's work. In Section 7, we present a equilibrium quantum formulation of Jarzynski's thermodynamics in the "bare" and "partial molar" representations, respectively in the same setting. We conclude in Section 8 with a short summary, a brief discussion of two issues involved and a suggestion for further developments.

Thermodynamic Functions, Hamiltonian of Mean Force
We first summarize the familiar traditional thermodynamic relations, if only to establish notations. We then consider interacting quantum systems with the help of the Hamiltonian of mean force (Hamiltonian of mean force is a useful yet not indispensable concept for this purpose. It is useful because in the same representation, the formal expressions associated with it resemble the counterparts in the traditional weak coupling thermodynamics). We outline two quantum formulations of thermodynamic functions and relations; one based on Gelin and Thoss [42] and the other on Seifert [26]. With the abundance of thermodynamic quantities, a word about notations is helpful: quantum expectation values or classical ensemble averages are denoted by math calligraphic, quantum operators associated with the variable O will carry an overhatÔ.

Traditional (Weak-Coupling) Equilibrium Thermodynamic Relations
The pre-conditions for the traditional weak-coupling thermodynamic theory to be well-defined and operative for a classical or quantum system are very specific despite their wide ranging applicability: (a) A system S of relatively few degrees of freedom is in contact with a thermal bath of a large number or infinite degrees of freedom (We shall consider only heat but no particle transfer here and thus the thermodynamics refers only to canonical, not grand canonical ensembles); (b) the coupling between the system and the bath is vanishingly small; and (c) the system is eternally in a thermal equilibrium state by proxy with the bath which is impervious to any change in the system. In weak-coupling thermodynamics, the bath variables are not dynamical variables (Dynamical variables are those which are determined consistently by the interplay between the system and the bath through their coupled equations of motion); they only provide weak-coupling thermodynamic parameters such as a temperature in canonical ensemble, or, in addition, a chemical potential, in grand canonical ensemble.
The classical thermodynamic relations among the internal energy U , enthalpy H, Helmholtz free energy F and Gibbs free energy G in conjunction with the temperature T, entropy S, pressure P and volume V are well-known. From the first law, By virtue of the differential form of the internal energy, the enthalpy H = U + PV obeys dH = T dS + V dP. Since it is a function of the entropy and pressure, we can identify Likewise, for the Helmholtz free energy F = U − TS, we have so the Helmholtz free energy is a function of the temperature and the volume, F = F (T, V ). Finally, the Gibbs free energy, defined by G = H − TS, obeys Thus G = G(T, P). Many more relations can be derived from these three basic relations. These relations are mutually compatible based on differential calculus.
Next we turn to the weak-coupling thermodynamics of quantum systems (Hereafter, we will choose the units such that the Boltzmann constant k B = 1 and the reduced Planck constanth = 1. In addition, to distinguish them from strong-coupling thermodynamics, all quantities defined in the context of traditional (weak-coupling) thermodynamics are identified with a subscript Θ). The state of a quantum system in contact with a heat bath at temperature T = β −1 with vanishing coupling is described by the density matrixρ s = e −βĤ s /Z Θ where Z Θ = Tr s e −βĤ s is the canonical partition function. HereĤ s is the Hamiltonian of the system and is assumed to be independent of the inverse temperature β = T −1 . The notation Tr with a subscript s or b represents the sum over the states of the system or the bath respectively. The density matrixρ s is a time-independent Hermitian operator and is normalized to unity, i.e., Tr sρs = 1 to ensure unitarity.
The free energy F Θ of a quantum system in a canonical distribution is F Θ = −β −1 ln Z Θ . The quantum expectation value Ĥ s is identified with the internal energy U Θ of the quantum system, and can be found by If we define the entropy S Θ of the system by then it will be connected with the internal energy by the relation F Θ = U Θ − T S Θ . These two expressions imply that the entropy of the quantum system can be expressed in terms of the density matrix S Θ = − Tr s ρ s lnρ s , which is the von Neumann entropy. The von Neumann entropy plays an important role in quantum information as a measure of quantum entanglement, and can be used to measure the non-classical correlation in a pure-state system. (Beware of issues at zero temperature as discussed in Section 6.) Here we note that both internal energy and the entropy of the system can be equivalently defined in terms of the expectation values of the quantum operators, or as the derivative of the free energy.
The heat capacity It is always semi-positive. Up to this point, under the vanishing system-bath coupling assumption, all the quantum thermodynamic potentials and relations still resemble their classical counterparts.

Quantum System in a Heat Bath with Nonvanishing Coupling
In formulating the quantum thermodynamics at strong coupling, we immediately face some conceptual and technical issues. At strong coupling, the interaction energy between the system and the bath is not negligible, so the total energy cannot be simply divided as the sum of the energies of the system and the bath. This introduces an ambiguity in the definition of, for example, internal energy. We may have more than one way to distribute the interaction energy between the system and the bath. The same ambiguity also arises in the other thermodynamic functions such as enthalpy and entropy, thus affecting the relations among the thermodynamic functions. On the technical side, in the course of formulating quantum thermodynamics, the non-commutative natures of the quantum operators make formidable what used to be straightforward algebraic manipulations in the classical thermodynamics.
Let us illustrate the previous points by an example. Consider in general an interacting quantum system C whose evolution is described by the HamiltonianĤ c =Ĥ s +Ĥ i +Ĥ b , whereĤ s ,Ĥ b are the Hamiltonians of the system S and the bath B, respectively andĤ i accounts for the interaction between them. Suppose initially the composite C = S + B is in a global thermal equilibrium state which is stationary, and thus has reversible dynamics, described by the density matrixρ c = e −βĤ c /Z c at the inverse temperature β −1 . The quantity Z c = Tr sb e −βĤ c is the partition function of the composite for the global thermal state.
In the case of vanishing coupling between the system and the bath, we may approximate the total HamiltonianĤ c to the leading order byĤ c Ĥ s +Ĥ b . Since [Ĥ s ,Ĥ b ] = 0, we notice that with the partition function of the free bath being given by Z b = Tr b e −βĤ b . Equation (8) implies that the reduced stateρ r = Tr bρc , which is also stationary, will assume a canonical form that is, Z c Z s Z b in the limit of vanishing system-bath coupling. In addition, (9) ensures the proper normalization condition Tr s ρ r = 1. Thus in the weak limit of the system-bath interaction, the reduced density matrix of the interacting composite system in the global thermal state will take the canonical form, hence it to some degree justifies the choice of the system state in the context of quantum thermodynamics in the textbooks [59,60]. Hereafter we will denote the reduced density matrix of the system byρ s . When the interaction between the system and the bath cannot be neglected, the righthand side of (8) no longer holds. In addition, non-commutating nature among the operatorsĤ s ,Ĥ i andĤ b prevents us from writing In fact, according to the Baker-Campbell-Haussdorff (BCH) formula, the previous decomposition will have the form The exponent on the righthand side typically contains an infinite number of terms. This makes algebraic manipulation of the strongly interacting system rather complicated, in contrast to its classical or quantum weak-coupling counterpart.

Hamiltonian of Mean Force
To account for non-vanishing interactions (in this paper, we apply the Hamiltonian of mean force to a quantum system in the global thermal state setup without any time-dependent driving force. See [53][54][55] for its use in nonequilibrium systems at strong coupling), one can introduce the Hamiltonian of mean force H * s for the system defined by [56][57][58] It depends only on the system operator but has included all the influences from the bath. In the limitĤ i is negligibleĤ * s Ĥ s ; otherwise, in generalĤ * s =Ĥ s . The corresponding partition function Z * is then given by Z * = Tr s e −βĤ * s = Z −1 b Tr sb e −βĤ c = Z c /Z b . If one followed the procedure of traditional weak-coupling thermodynamics to define the free energy as F = −β −1 ln Z, then the total free energy F c of the composite system can be given by a simple additive expression In addition, one can write the reduced density matrixρ s in a form similar to (9), with the replacement ofĤ s byĤ * s , in the hope that the conventional procedures of weak-coupling thermodynamics will follow. However, in Section 7 or in [42], we see that at strong coupling even though we already have the (reduced) density matrix and the free energy of the system, we can introduce two different sets of thermodynamic potentials for the system. The thermodynamic potentials in each set are mathematically self-consistent, but they are not compatible with their counterparts in the other set, in contrast to the weak coupling limit, where both definitions are equivalent.
Two earlier approaches to introduce the thermodynamic potentials in a strongly interacting system had been proposed by Gelin and Thoss [42], and by Seifert [26]. We shall summarize the approach I in Gelin and Thoss' work below and present a more detailed quantum formulation of Seifert's approach following, both for a configuration that the composite is initially in a global thermal state without any external force. A recent proposal by Jarzynski [27] for classical systems will be formulated quantum-mechanically in the same setting in Section 7.

Quantum Formulation of Gelin and Thoss' Thermodynamics at Strong Coupling
The first approach, based on Approach I of Gelin & Thoss [42], is rather intuitive, because their definitions of the internal energy and the entropy are the familiar ones in traditional thermodynamics. They define the internal energy U s of the (reduced) system by the quantum expectation value of the system Hamiltonian operator alone, U s = Tr s ρ sĤs , and choose the entropy to be the von Neumann (vN) entropy S s = S vN = − Tr s ρ s lnρ s . These are borrowed from the corresponding definitions in weak-coupling thermodynamics.
They write the same reduced density matrix (9) in a slightly different representation to highlight the difference from the weak-coupling thermodynamics case, where∆ s depends only on the system variables but includes all of the influence from the bath from their interaction. Comparing this with (11), we note that∆ s is formally related to the Hamiltonian of mean force by∆ s =Ĥ * s −Ĥ s + F b . Finally, they let the partition function of the system take on the value Z c , which is distinct from Z * . Thus the corresponding free energy will be given by F c which contains all the contributions from the composite C.
Although in this approach the definitions of internal energy and entropy of the system are quite intuitive, these two thermodynamic quantities do not enjoy simple relations with the partition function Z c , as in (5). From (13), we can show (Here some discretion is advised in taking the derivative with respect to β because in general an operator will not commute with its own derivative. See details in Appendix A) with the corresponding free energy F c = U s + ∆ s + β ∂ β∆s − β ∂ β F c . Here · · · represents the expectation value taken with respect to the density matrixρ c of the composite. For a system operator O s , this definition yields an expectation value equal to that with respect to the reduced density matrix ρ s , namely, Ô s s = Tr s ρ sÔs = Tr sb ρ cÔs = Ô s .
Likewise, the von Neumann entropy S vN can be expressed in terms of the free energy F c by which does not resemble the traditional form in (3). Additionally, we observe the entropy so defined is not additive, that is, Here S c and S b are the von Neumann entropies of the composite and the free bath, respectively. Note that theρ b in this formulation is the density matrix of the free bath, not the reduced density matrix of the bath, namely,ρ b = Tr sρc . The reduced density matrix of the bath will contain an additional overlap with the system owing to their coupling.
When the internal energy of the system given by the expectation value of the system Hamiltonian operatorĤ s , the specific heat C s will take the form, In general ∂ β∆s = ∂ β ∆ s since the reduced density matrixρ s also has a temperature dependence. We thus see in this case the heat capacity cannot be directly given as the derivative of the (von Neumann) entropy with respect to β, as in (7).
In short, in this formulation, the thermodynamic potentials of the system are defined in a direct and intuitive way. They introduce an operator∆ s to highlight the foreseen ambiguity when the system is strongly coupled with the bath. Formally we see that∆ s = −β −1 ln Tr b e −β(Ĥ s +Ĥ i +Ĥ b ) −Ĥ s . In the limit of weak coupling,Ĥ i ≈ 0, the operator∆ s reduces tô Hence in this limit,∆ s reduces to a c-number and plays the role of the free energy F b of the free bath. Observe∆ s ≈ F b in the weak coupling limit annuls the expression in the square brackets in (18) and restores the traditional relation (7) between the heat capacity and the entropy. However, even in the weak coupling limit, the internal energy still cannot be given by (5). The disparity lies in the identification of Z c as the partition function of the system. As is clearly seen from (14), in the weak coupling limit, we have This implies that Z c is not a good candidate for the partition function of the system. A more suitable option would be Z c /Z b .

Quantum Formulation of Seifert's Thermodynamics at Strong Coupling
If we literally follow (11) and identifyĤ * s as the effective Hamiltonian operator of the (reduced) system, we will nominally interpret that the reduced system assumes a canonical distribution. Thus it is natural to identify Z * as the partition function associated with the reduced state of the system.
Suppose we maintain the thermodynamic relations regardless of the coupling strength between the system and the bath. From (5) to (6), we will arrive at expressions of the internal energy and entropy of the system. This is essentially Seifert's approach [26] to the thermodynamics at strong coupling. Here we will present the quantum-mechanical version of it for a equilibrium configuration without the external drive, that is, λ = 0 in Seifert's notion.
First, from (11), we have the explicit form of the Hamiltonian of mean forceĤ * This is the operator form of H(ξ s , λ) in Equation (5) of [26]. Noting the non-commutative characters of the operators. Since [Ĥ s ,Ĥ i ] = 0, If one prefers factoring out e −βĤ s from e −β(Ĥ s +Ĥ i +Ĥ b ) , one can use the Baker-Campbell-Haussdorff formula, outlined in Appendix A, to expand out the operator products to a certain order commensurate with a specified degree of accuracy.
Second, it is readily seen that p eq (ξ s |λ= 0) in Equation (4) of [26] is the reduced density matrixρ s of the system (12). The (Helmholtz) free energy F in Seifert's Equation (7) is exactly the free energy of the reduced system F * in Section 3.3.
With these identifications, it is easier to find the rest of the physical quantities in Seifert's strong coupling thermodynamics for the equilibrium configuration. We now proceed to derive the entropy and the internal energy, i.e., Equations (8) and (9) of [26], for quantum systems in his framework in the equilibrium setting.
In general, an operator does not commute with its derivative, so taking the derivative of an operator-valued function or performing integration by parts on an operator-valued function can be nontrivial. Their subtleties are discussed in Appendix A, where we show that the derivative of an operator function is in general realized by its Taylor's expansion in a symmetrized form (A5). However, when such a form appears in the trace, the cyclic property of the trace allow us to manipulate the derivative of a operator-valued function as that of a c-number function thus sidestepping the symmetrized ordering challenge. Hence from the thermodynamic relation (6), we have Here we recall that even though the operatorsĤ * s and ∂ βĤ * s in general do not commute, the trace operation allowing for cyclic permutations of the operator products eases the difficulties in their manipulation. Since (12) implies the operator identity βĤ * s = βF * − lnρ s , we can recast (23) to This is the quantum counterpart of Seifert's entropy, Equation (8) of [26]. This entropy is often called the "thermodynamic" entropy in the literature. Note that it is not equal to the von Neumann ("statistical") entropy S vN of the system. The internal energy can be given by the thermodynamic relation Thus from (23), we obtain, This deviation results from the fact thatĤ * s , introduced in (11) may depend on β. When we take this into consideration, we can verify that the internal energy can also be consistently given by In fact, we can show, by recognizing with Ĥ b b ≡ Tr b ρ bĤb , Ĥ b ≡ Tr sb ρ cĤb , and Ĥ s ≡ Tr sb ρ cĤs = Ĥ s s . Equation (28) implies that the internal energy, defined by (27), accommodates more than mere Ĥ s s . The additional pieces contain contributions from the bath and the interaction. In particular, when the coupling between the system and the bath is not negligible, we have Ĥ b = Ĥ b b in general. As a matter of fact, even the internal energy defined by the expectation value of the system Hamiltonian operator in approach I of Gelin & Thoss' work also encompasses influence from the bath because the reduced density matrixρ s includes all the effects of the bath on the system. Hitherto, we have encountered three possible definitions of internal energies, namely, Ĥ s , Ĥ * s , and U s . As can be clearly seen from (26) and (28), they essentially differ by the amount of the bath and the interaction energy which are counted toward the system energy. This ambiguity arises from strong coupling between the system and the bath. When the system-bath interaction is negligibly small, we have Ĥ i ≈ 0, and since in this limit, the full density matrix of the composite is approximately given by the product of that of the system and the bath, we arrive at Ĥ b ≈ Ĥ b b , and these three energies become equivalent.
To explicate the physical content of Ĥ * s , from (12) we can write Ĥ * s as This offers an interesting comparison with (25), where F * = U s − β −1 S s . It may appear that we can replace the pair (U s , S s ) by another pair ( Ĥ * s , S vN ), leaving F * unchanged, thus suggesting an alternative definition of internal energy by Ĥ * s and that of entropy by S vN . However, in so doing, the new energy and entropy will not satisfy a simple thermodynamic relation like (5) and (6). This is a good sign, as it is an indication that certain internal consistency exists in the choice of the thermodynamic variables.
We now investigate the differences between the two definitions of entropy. From (24), we obtain The factor ∂ βρs can be written as Thus, part of the difference between the two entropies result from the correlation between the full HamiltonianĤ c and the Hamiltonian of mean forceĤ * s . This correlation will disappear in the vanishing coupling limit because there is no interaction to bridge the system and the bath. We also note that in the same limit, Ĥ * s ≈ Ĥ s becomes temperature-independent, and both definitions of the entropy turn synonymous.
Since the von Neumann entropy S vN can be used as a measure of entanglement between the system and the bath at zero temperature, we often introduce the quantum mutual information I sb to quantify how they are correlated, where S b is the von Neumann entropy associated with the reduced density matrixˆ b of the bath, in contrast to S b we have met earlier. This mutual information can be related to the quantum relative entropy S(ρ c ρ s ⊗ˆ b ) by becauseˆ b = Tr sρc . On the other hand, Equation (23) imply that the thermodynamic entropy S s is additive S s + S b = S c , from which we find This and (31) provide different perspectives on how the difference between the two system entropies is related to the system-bath entanglement, and how the system-bath coupling has a role in establishing such correlations.
Following the definitions of the internal energy (27) and the entropy (23), the heat capacity of the system still satisfies a familiar relation C s = −β 2 ∂ β U s = β 2 ∂ 2 β ln Z * = −β ∂ β S s . Compared with (18), with the help of (28), we clearly see their difference, caused by different definitions of internal energy, is given by

Issues of These Two Approaches: Entropy and Internal Energy
Both equilibrium quantum formulations for thermodynamics at strong coupling are based on plausible assumptions and are mathematically sound. In Approach I outlined in Section 4, one starts with intuitive definitions of the thermodynamics quantities, inspired by traditional thermodynamics for classical systems premised on vanishingly weak coupling between the system and the bath. This leads to modifications in the thermodynamic relations of the relevant thermodynamics quantities. In Approach II delineated in Section 5, one opts to maintain the familiar thermodynamic relations but is compelled to deal with a rather obscure interpretation of the thermodynamic potentials. Although both approaches in the vanishing system-bath coupling limit are compatible, as shown in Section 3.1, they in general entail distinct definitions of the thermodynamic functions. This disparity is amplified with strong system-bath coupling in the deep quantum regime, where quantum coherence plays an increasingly significant role. Thus, even though both approaches possess the same correct classical thermodynamic limit, they are not guaranteed to give unique physical results in the deep quantum regime, even for simple quantum systems, which are areas for interesting further investigations.
To highlight the issues more explicitly, we can apply these two methods to a simple and completely solvable model, namely, a Brownian oscillator linearly but strongly coupled with a large (or infinitely large, as modeled by a scalar field) bath. We will see both approaches at some point, or others that produce ambiguous or paradoxical results. We make a few observations in the following section.

Entropy
(1) It has been discussed in [36,43,44] that the von Neumann entropy S vN will not approach to zero for the finite system-bath coupling in the limit of zero temperature, but the thermodynamic entropy S s , defined in Approach II, behaves nicely in the same limit. (2) It has been shown [49] that if the composite is in a global thermal state, the discrete energy spectrum of the undamped oscillator will become a continuous one with a unique ground level. This supports physics described by the thermodynamic entropy S s .
It has been argued [44][45][46] that the entanglement between the system and the bath prevents the von Neumann entropy from approaching zero at zero temperature. Without quantum entanglement between the system and the bath, the lowest energy level of the composite system will be given by the tensor product of the ground state of the unperturbed system and bath, that is, a pure state. In this case, the von Neumann entropy will go to zero as expected, and this is the scenario that occurred in traditional quantum/classical thermodynamics in the vanishing system-bath coupling limit.

Internal Energy
• It has been discussed [37][38][39][40] that the internal energy defined in Approach II can lead to anomalous behavior of the heat capacity in the low temperature limit. When the system, consisting of a quantum oscillator [40] or a free particle [37][38][39] is coupled to a heat bath modeled by a large number of quantum harmonic oscillators, the heat capacity of the system can become negative if the temperature of the bath is sufficiently low. If the internal energy defined in Approach I is used to compute the heat capacity, then it has been shown that the heat capacity remains positive for all nonzero temperatures but vanishes in the zero bath temperature limit, for a system with one harmonic oscillator [37], or a finite number of coupled harmonic oscillators [29]. This discrepancy may result from the fact that the internal energy defined in Approach II contains contributions from the interaction and the bath Hamiltonian.
It seems to imply that in the low-temperature, strong coupling regime, it remains an open question how to properly define the thermodynamic functions; being able to show the well-known behaviors in the classical thermodynamic limit is a necessary but not sufficient condition.

Quantum Formulation of Jarzynski's Strong Coupling Thermodynamics
We now provide a quantum formulation of Jarzynski's classical results [27] but for composite system C kept under thermal equilibrium. The Hamiltonian operator of the composite C = S + B is assumed to take the formĤ Here in this paper, J will be some external, but constant c-number drive acting on the bath via a bath operatorÂ b . It can be a constant pressure, as in Jarzynski's classical formulation, and thenÂ b will be an operator corresponding to V b , conjugated to P. However, in general,Â b will be the operator of the quantity conjugated to J. This analogy, though formal, provides an alternative route to introduce the operator conjugated to J.
If the composite system is in thermal equilibrium at the temperature β −1 , its state is described by the density matrix operatorρ c = e −βĤ c /Z c , where Z c = Tr sb e −βĤ c , a c-number, is the partition function of the composite. For later convenience, we also define the corresponding quantities for the bath B when it is coupled to the system S, . We introduce the Hamiltonian operator of mean forceĤ * s by such that the reduced density matrix of the system S takes the form The quantity Z s can be viewed as an effective partition function of the system S. This is motivated by the observation that, in the absence of coupling between S and B, or in the weak coupling limit, the composite is additive so its partition function is the product of those of the subsystems, i.e., Z c = Z s Z b . The differenceĤ * s −Ĥ s modifies the dynamics of the system S due to its interaction with the bath B. In fact, by the construction, e −β H * s , once sandwiched by the appropriate states of the system S and expressed in the imaginary-time path integral formalism (for further details regarding the connection with the influence action, please refer to [28,61,62]), is formally e −S CG , where S CG is the coarse-grained effective action of the system S, wick-rotated to the imaginary time. Thus formally β(Ĥ * s −Ĥ s ) is equivalent to the influence action in the imaginary time formalism.
Similar to the classical formulations, we may have two different representations of the operator A s of the system.

"Bare" Representation
In the bare representation, we may defineÂ s = (Ĥ * s −Ĥ s )/J, and the internal energy operatorÛ s and the enthalpy operatorĤ s , respectively, byÛ s =Ĥ s andĤ s =Ĥ * s , with expectation values given by U s = Tr s ρ sÛs and H s = Tr s ρ sĤs = U s + J · A s , corresponding to the internal energy and the enthalpy we are familiar with, respectively. Here A s = Tr s ρ sÂs . The entropy is chosen to be the von Neumann entropy of the system S s = Tr s ρ s lnρ s = β H s − G s .
The Gibbs free energy G s is defined as G s = −β −1 ln Z s . These definitions are in exact parallel to those in the classical formulation [27].

"Partial Molar" Representation
In contrast to the bare representation, we can alternatively define the operatorÂ s of the system S that corresponds toÂ b of the bath B byÂ s = ∂ Ĥ * s −Ĥ s /∂J = ∂Ĥ * s /∂J. The last equality results from the fact thatĤ s has no dependence on the external parameter J. Owing to the non-commutativity of operators, the micro-physics interpretation of the operatorÂ s is not so transparent. We first focus on its quantum expectation value A s As stressed earlier, since the operator does not commute with its derivative, care must be taken when we move the derivative around in an operator expression. However, from (A7), we learn that the righthand side of (40) can be identified as and thus we have A s = −β −1 ∂ ln Z s /∂J. The advantage of this expression is that the observation of In particular we can check that A b indeed is the expectation value of the operatorÂ b , that is, A b = Tr b ρ bÂb = Tr sb ρ bÂb . The latter expression can nicely bridge with A c for the composite, A c = Tr sb ρÂ b . Thus the expectation value A is additive. Its value for the combined systems is equal to the sum of those of the subsystems, A c = A s + A b . In fact, this additive property holds for all the thermodynamics potentials introduced afterwards. This is a nice feature in Jarzynski's partial molar representation or in Seifert's approach. From this aspect, we can interpret A s as the change of A b due to the intervention of the system S. For example, consider a photon gas inside a cavity box, one side of which is a movable classical mirror and is exerted by a constant pressure. Assume originally the photon gas and the mirror are in thermal equilibrium. In this cavity we now place a Brownian charged oscillator and maintain the new composite system in thermal equilibrium at the same temperature and the same pressure (The equilibration process in this example can be awfully complicated if we mind the subtleties regarding whether the photon gas can ever reach thermal equilibrium in a cavity whose walls are perfectly reflective and so on. For the present argument, we assume equilibration is possible and there is no leakage of the photons). Then we should note that there is a minute change in the mean position of the mirror before and after the Brownian charged oscillator is placed into the cavity. This change can also be translated to an effective or dynamical size of the charged oscillator due to its interaction with the photon gas, and thus is accounted for in A s when J is identified as the constant pressure applied to the wall.
From this example, it is tempting to identify J ·Â s as some quantum work operator (its value depends on the interaction between the system and the bath and when this interaction is switched on. It is thus path-dependent in the parameter space of the coupling constant). Alternatively we may view it or its expectation as some additional "energy content" of the system S due to its interaction with the bath when the composite is acted upon by an external agent J, sinceÂ s is related toĤ * s −Ĥ s [63,64]. Inspired by this observation and taking the hint from the definition of A s , we introduce the enthalpy of the system S by where we have identified the enthalpies of the composite C = S + B and the bath B as H c = −∂ ln Z c /∂β and H b = −∂ ln Z b /∂β. We may rewrite them as It implies that (1) the system enthalpy can be decomposed as and (2) the internal energy U s of the system S can be consistently inferred as This is exactly the same internal energy obtained in Seifert's approach in the equilibrium setting. We can define the internal energy of the composite system and of the bath by U c = Ĥ s + Ĥ i + Ĥ b and U b = Ĥ b b , and then we may conclude U s = U c − U b . Thus the internal energy U s also includes contributions that naïvely we will not ordinarily attribute to the system, such as Ĥ b − Ĥ b b . Doing so will complicate the physical connotation of the internal energy of the system.
Up to this moment, we essentially write the thermodynamic quantities by the quantum expectation value and in terms of the partition functions. Thus it is appropriate to introduce the Gibbs free energies of the composite C, the system S, and the bath B by G a = −β −1 ln Z a , where a = c, s and b, and they obey the additive property of the Gibbs energy, G c = G s + G b . Furthermore, in the composite, we note that From (45), we can consistently define the entropy S of the composite by S c = β H c − G c and, similarly, the entropy S b of the bath: The additive property of the free energy and the enthalpy implies that the entropy S s of the system in this representation is also additive, S s = S c − S b , and is given by Note it is not equal to the von Neumann entropy, which is defined as the entropy of the system in the "bare" representation.

Operator Forms of the Thermodynamic Functions
In trying to formulate a set of laws to describe the thermodynamics of a quantum system (even the existence of such a theory, under certain appropriate conditions, is not a matter of presumption or prescription, but by demonstration and proof) it would be most convenient if we could define operators of the thermodynamic functions in such a way that the quantum expectation values of those operators give the familiar expressions for the thermodynamic functions. As we see it, this is the paramount challenge in the formal establishment of quantum thermodynamics as a viable theory. The laws of thermodynamics have been understood in terms of the mean values of the relevant operator quantities. For a system where the fluctuations of the thermodynamic functions become comparable to the corresponding mean values, the thermodynamic laws governing the mean values need be supplanted by laws governing their quantum fluctuations or higher order quantum correlations. A case in point for classical systems where fluctuations are as important as the mean values is near the critical point of the system. The truly quantum properties would impact on the quantum thermodynamics for small quantum systems in the regimes of strong couplings to its environment, and at low temperatures, where quantum coherence effects take center stage. Having the operator forms of these thermodynamic potentials allows us to calculate the higher-order quantum correlations of those quantities existent in larger fluctuations.
In the following sections, we will attempt to identify the operator form of the thermodynamic function for the reduced system.
This result in (48) is nonsensical because (1) the righthand side still explicitly depends on the bath degree of freedom; (2) we cannot take its trace with respect to the state of the system,ρ s ; and thus (3) the expectation value will not be H s . This is because the operators defined this way act on Hilbert spaces different from that ofρ s ;Ĥ c is an operator in the Hilbert space of the composite whileĤ b is an operator in the Hilbert space of the bath. Neither operator acts exclusively in the Hilbert space of the system. Thus, extreme care is needed when manipulating the operator forms of the thermodynamical potentials. What one needs to do is to seek the local forms of these operators, i.e., operators which act only on the Hilbert space of the system. This can be done in parallel to Jarzynski's classical formulation.

System Enthalpy Operator: Approved
We first inspect the internal energy operator. Since the averaged internal energy of the composite system is given by U c = Tr sb e −βĤ c Ĥ s +Ĥ i +Ĥ b /Z c , we can rewrite the expressions inside the trace into in a way analogous to Jarzynksi's classical formulation. Here we have used the fact that Z c = Z s Z b and the identity for the operatorẐ î If we define an internal energy operatorÛ i bŷ then we obtain a new representation of U c U c = Tr s ρ sÛi .
Equation (50) is an operator expression of the internal energy of the composite system, on account of the non-commutativity of the operators, but its expectation value is taken with respect to the system's density matrixρ s . With the help of (A9), this is equivalent to Equation (28) of [26] in the J = 0 case. In addition, we note thatẐ i is an operator, not a c number. Since Tr sρs = 1, we may define the operatorÛ s byÛ such that Tr s ρ sÛs = Tr s ρ sÛi − Tr s ρ s U b = U c − U b = U s . The advantage of (50), (52) is that, unlike those introduced in the previous subsection, they are all operators in the Hilbert space of the system S. Indeed, using the identity operatorÎ s in the Hilbert space of the system S we can also definê In the same fashion, we may rewrite A c by Thus we can defineÂ so that A c = Tr s ρ sÂi . We then can have a local form for theÂ s given byÂ s =Â i −Â b in close resemblance to its classical expression in [27], if we re-defineÂ b asÂ b = A bÎs . The expectation value ofÂ s is then given by Tr s ρ sÂs = Tr s ρ sÂi − A b = A c − A b = A s . Now we proceed with constructing a local form of the enthalpy operator of the system. From (52) and the definition of the operatorÂ s , we claim that the local formĤ s iŝ We can straightforwardly show that Tr s ρ sĤs = Tr s ρ sÛs + J · Tr s ρ sÂs = U s + J · A s . Thus we have succeeded in constructing the operators that correspond to A s , U s , H s in forms local in the Hilbert space of the system S. However, as can be seen from their expressions, their meanings are not transparent a priori. They are determined a posteriori because we would like their expectation values to take certain forms. This can pose a question about the uniqueness of these operators (A similar issue is also raised in [58] for the classical formulation. However, in this context it is not clear whether this ambiguity can be fixed by calculating the cumulants associated with these operators. If there exist physical, measurable observables that correspond to the expectation values of the moments of these operators, then one may entertain the possibility of using them to uniquely determine these operators.). At least for a given reduced density matrixρ s of the system, we can always attach a system operator Λ s that satisfies Tr sρsΛs = 0 to the definitions of those local operators, that is, any system operator that has a zero mean. The choice ofΛ s may not be unique in the sense that in the basis {|n } that diagonalizesρ S , we can write Tr s {ρ SΛS } = 0 as It says that the vectors that are respectively composed of the diagonal elements ofρ S andΛ S are orthogonal, but it does not place any restriction on the off-diagonal elements ofΛ S on this basis.

Summary
In this paper we provide quantum formulations of equilibrium thermodynamic functions and their relations for Jarzynski's classical thermodynamics at strong coupling [27] without the consideration of heat and work. The combined system + environment, called the composite, is assumed initially to be in a global thermal state. In such a configuration, even though the interaction between the system and the bath is non-negligible, the partition function of the composite is well defined. This facilitates the introduction of thermodynamic potentials in a way similar to the traditional vanishing-coupling thermodynamics. Such a configuration allows for the introduction of enthalpy by introducing an external agent (keeping at constant pressure in Jarzynski's case) acting on the conjugate bath operator. The effect can be represented by an equivalent effect on the system, which then appears in the expression of the enthalpy of the system. There are two ways to capture this effect, called the "bare" and "partial molar" representations by Jarzynski. We have worked out a quantum formulation for each of these two representations of Jarzynski's classical thermodynamics. In addition, we attempt to identify the operator forms of the thermodynamic functions, which can be potentially useful in directly addressing the effects of quantum fluctuations on thermodynamics at strong coupling, where the fluctuations can reach the same order of magnitude as the corresponding mean values.

Issues
We mention two outstanding issues of the two representations or approaches, namely, internal energy and entropy of the system at strong coupling in the global thermal state. When the quantum versions of these two approaches are applied to a small quantum system that strongly couples with a low-temperature bath, some nonintuitive results have been reported in the literature [29,[36][37][38][39][40]43,44]. If the von Neumann entropy is adopted as the system entropy, then when the system and the bath are entangled, this entropy will not approach zero for a simple system such as a harmonic oscillator in the zero temperature bath, contradicting the result in [49], where it has been shown that the ground state of such a composite system is non-degenerate in general, thus implying vanishing entropy at zero temperature. On the other hand, if the system's internal energy is defined as the derivative of the partition function of the system, the heat capacity derived therefrom can take on negative values in the low temperature regime when the system consists of free particles or coupled harmonic oscillators. This anomalous behavior, not seen when the internal energy defined as the expectation value of the system Hamiltonian, may be traced to an excessive inclusion of the interaction and the bath contributions, as shown in (44), in the definition of the internal energy of the system.

Further Developments
To complete a theory for equilibrium quantum thermodynamics, we need to include the considerations of heat and work. For understanding quantum work, the physical meaning of the operatorÂ s conjugate to J is a key issue. At the next level of investigation, reaching out to nonequilibrium conditions, also at strong coupling, we posit that the operator∆ s , introduced by Gelin and Thoss, and the Hamiltonian operator of mean force can be related to the influence action or the coarse-grained effective action [28,61,62] of the system when they are sandwiched by the states of the system and formulated in the finite temperature imaginary-time path integral method. This and an earlier observation we made for the partition function point out a way to extend the present equilibrium quantum thermodynamics at strong coupling for a closed system in a global thermal state to a nonequilibrium setting, by employing the real-time closed-time-path formalism used in [28,29,61,62]. This is the goal of our next paper in this series [65].

Conflicts of Interest:
The authors declare no conflict of interest.

Appendix A. Handling Operator Products in Quantum Thermodynamics
In deriving various thermodynamic relations in the context of quantum thermodynamics, we often come up with expressions involving exponential of the sum of two operators, sayλ andμ. Unlike its c-number counterpart, in general such an exponential cannot be written as a product of exponential of the respective operators, that is, because these two operatorsλ,μ may not commute. From Baker-Campbell-Haussdorff formulas, the righthand side of (A1) in fact is given by for any two operatorsλ,μ. However, in the special case that [λ,μ] = 0, the equality in (A1) indeed is valid. The other useful expression in the Baker-Campbell-Hausdorff formulas is eλμ e −λ =μ + λ ,μ + 1 2 λ , λ ,μ + · · · . (A3) This is particularly useful in deriving the unitary transformation ofμ by the operator eλ. We also often come to a situation that we need to take a derivative of an exponential of the operator. This is less straightforward than is expected due to the fact that the operator in the exponent may not commute with its own derivative. For example, consider an operatorÔ(χ) of the form O(χ) = α(χ)X + β(χ)P, with α, β being functions of χ, but the operatorsX,P of the canonical variables have no explicit χ dependence. We immediately see the trivial result [Ô(χ),Ô(χ)] = 0, but Ô (χ), ∂ χÔ (χ) = αβ −αβ X ,P = 0, where the overhead dot represents the derivative with respect to χ. This introduces complications in taking the derivative of, say, e −Ô(χ) with respect to χ. If we realize an operator function in terms of its Taylor's expansion, then Taking the derivative of (A4) with respective to χ yields k ∂ χÔ Ô · · · · · ·Ô (k−1) terms +Ô ∂ χÔ Ô · · · · · ·Ô (k−2) terms + · · · · · · +Ô · · · · · ·Ô (k−1) terms ∂ χÔ where we define the symmetrized product (Ô 1Ô2 · · ·Ô k ) sym as a generalization of the anti-commutator by (Ô 1Ô2 · · ·Ô k ) sym = 1 # of perm. ∑ # of perm.Ô σ 1Ô σ 2 · · ·Ô σ k , with σ being the permutations of 1, 2, · · · , k. Thus the expression, say, (∂ βĤ * s ) e −βĤ * s , and similar derivative expressions will be understood in this manner as a symmetrized product of ∂ βĤ * s and the Taylor-expanded e −βĤ * s , shown in (A5). However, if the derivative like (A5) is taken within a trace, then the complicated expression (A5) will reduce to a simple form due to the cyclic property of the trace formula. Hence in general we have ∂ χ e −Ô = − ∂ λÔ e −Ô sym , but its trace gives as if the operatorÔ is a c-number. Note here we have assumed the traces applied in (A7) and (A8) are not a partial trace; otherwise the same symmetrization procedure is still needed. A special case of (A7) is ∂ χ e −χÔ , in whichÔ has no explicit dependence on χ. Then it is straightforward to perform the differentiation, and since [Ô, e −χÔ ] = 0, we obtain ∂ χ e −χÔ = −Ô e −χÔ . (A9) Next we give an explicit application of (A7). In particular, we focus on the expression Tr s e −β(Ĥ s +∆ s ) , with Z c = Tr s e −β(Ĥ s +∆ s ) .