1. Quantum and Thermodynamics—Why?
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”?
1.1. 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.
1.2. 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.).
1.3. 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).
1.4. 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).
3. 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
.
3.1. 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 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
, enthalpy
, Helmholtz free energy
and Gibbs free energy
in conjunction with the temperature
T, entropy
, pressure
P and volume
are well-known. From the first law,
. With
, we have
By virtue of the differential form of the internal energy, the enthalpy
obeys
. Since it is a function of the entropy and pressure, we can identify
Likewise, for the Helmholtz free energy
, we have
so the Helmholtz free energy is a function of the temperature and the volume,
. Finally, the Gibbs free energy, defined by
, obeys
Thus . 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 and the reduced Planck constant . 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 with vanishing coupling is described by the density matrix where is the canonical partition function. Here is the Hamiltonian of the system and is assumed to be independent of the inverse temperature . 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 is a time-independent Hermitian operator and is normalized to unity, i.e., to ensure unitarity.
The free energy
of a quantum system in a canonical distribution is
. The quantum expectation value
is identified with the internal energy
of the quantum system, and can be found by
If we define the entropy
of the system by
then it will be connected with the internal energy by the relation
. These two expressions imply that the entropy of the quantum system can be expressed in terms of the density matrix
, 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
is given by
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.
3.2. 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 whose evolution is described by the Hamiltonian , where are the Hamiltonians of the system and the bath , respectively and accounts for the interaction between them. Suppose initially the composite is in a global thermal equilibrium state which is stationary, and thus has reversible dynamics, described by the density matrix at the inverse temperature . The quantity 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
to the leading order by
. Since
, we notice that
with the partition function of the free bath being given by
. Equation (
8) implies that the reduced state
, which is also stationary, will assume a canonical form
that is,
in the limit of vanishing system-bath coupling. In addition, (
9) ensures the proper normalization condition
. 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
.
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
,
and
prevents us from writing
, due to
and
in general. 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.
3.3. 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
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 is negligible ; otherwise, in general . The corresponding partition function is then given by .
If one followed the procedure of traditional weak-coupling thermodynamics to define the free energy as
, then the total free energy
of the composite system can be given by a simple additive expression
, with
and
. In addition, one can write the reduced density matrix
in a form similar to (
9), with the replacement of
by
,
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.
4. 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
of the (reduced) system by the quantum expectation value of the system Hamiltonian operator alone,
, and choose the entropy to be the von Neumann (vN) entropy
. 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
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
is formally related to the Hamiltonian of mean force by
. Finally, they let the partition function of the system take on the value
, which is distinct from
. Thus the corresponding free energy will be given by
which contains all the contributions from the composite
.
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
, 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
. Here
represents the expectation value taken with respect to the density matrix
of the composite. For a system operator
, this definition yields an expectation value equal to that with respect to the reduced density matrix
, namely,
Likewise, the von Neumann entropy
can be expressed in terms of the free energy
by
which does not resemble the traditional form in (
3). Additionally, we observe the entropy so defined is not additive, that is,
Here and are the von Neumann entropies of the composite and the free bath, respectively. Note that the in this formulation is the density matrix of the free bath, not the reduced density matrix of the bath, namely, . 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
, the specific heat
will take the form,
In general
since the reduced density matrix
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
to highlight the foreseen ambiguity when the system is strongly coupled with the bath. Formally we see that
. In the limit of weak coupling,
, the operator
reduces to
Hence in this limit,
reduces to a
c-number and plays the role of the free energy
of the free bath. Observe
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
as the partition function of the system. As is clearly seen from (
14), in the weak coupling limit, we have
This implies that is not a good candidate for the partition function of the system. A more suitable option would be .
5. Quantum Formulation of Seifert’s Thermodynamics at Strong Coupling
If we literally follow (
11) and identify
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
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,
in Seifert’s notion.
First, from (
11), we have the explicit form of the Hamiltonian of mean force
This is the operator form of
in Equation (5) of [
26]. Noting the non-commutative characters of the operators. Since
,
If one prefers factoring out
from
, 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
in Equation (4) of [
26] is the reduced density matrix
of the system (
12). The (Helmholtz) free energy
in Seifert’s Equation (
7) is exactly the free energy of the reduced system
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
and
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
, 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
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
, 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 Equation (
5)
In fact, we can show, by recognizing
, that
with
,
, and
. Equation (
28) implies that the internal energy, defined by (
27), accommodates more than mere
. 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
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
includes all the effects of the bath on the system.
Hitherto, we have encountered three possible definitions of internal energies, namely,
,
, and
. 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
, 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
, and these three energies become equivalent.
To explicate the physical content of
, from (
12) we can write
as
This offers an interesting comparison with (
25), where
. It may appear that we can replace the pair
by another pair
, leaving
unchanged, thus suggesting an alternative definition of internal energy by
and that of entropy by
. 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
can be written as
with
. We then obtain
Thus, part of the difference between the two entropies result from the correlation between the full Hamiltonian and the Hamiltonian of mean force . 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, becomes temperature-independent, and both definitions of the entropy turn synonymous.
Since the von Neumann entropy
can be used as a measure of entanglement between the system and the bath at zero temperature, we often introduce the quantum mutual information
to quantify how they are correlated,
where
is the von Neumann entropy associated with the reduced density matrix
of the bath, in contrast to
we have met earlier. This mutual information can be related to the quantum relative entropy
by
because
. On the other hand, Equation (
23) imply that the thermodynamic entropy
is additive
, 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
. Compared with (
18), with the help of (
28), we clearly see their difference, caused by different definitions of internal energy, is given by
6. 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.
6.1. Entropy
- (1)
It has been discussed in [
36,
43,
44] that the von Neumann entropy
will not approach to zero for the finite system-bath coupling in the limit of zero temperature, but the thermodynamic entropy
, 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
.
- (3)
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.
6.2. 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.
7. Quantum Formulation of Jarzynski’s Strong Coupling Thermodynamics
We now provide a quantum formulation of Jarzynski’s classical results [
27] but for composite system
kept under thermal equilibrium. The Hamiltonian operator of the composite
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 . It can be a constant pressure, as in Jarzynski’s classical formulation, and then will be an operator corresponding to , conjugated to P. However, in general, 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
, its state is described by the density matrix operator
, where
, a
c-number, is the partition function of the composite. For later convenience, we also define the corresponding quantities for the bath
when it is coupled to the system
,
with the bath partition function
. We introduce the Hamiltonian operator of mean force
by
such that the reduced density matrix of the system
takes the form
The quantity can be viewed as an effective partition function of the system . This is motivated by the observation that, in the absence of coupling between and , or in the weak coupling limit, the composite is additive so its partition function is the product of those of the subsystems, i.e., . The difference modifies the dynamics of the system due to its interaction with the bath .
In fact, by the construction,
, once sandwiched by the appropriate states of the system
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
, where
is the coarse-grained effective action of the system
, wick-rotated to the imaginary time. Thus formally
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 of the system.
7.1. “Bare” Representation
In the bare representation, we may define
, and the internal energy operator
and the enthalpy operator
, respectively, by
and
, with expectation values given by
and
, corresponding to the internal energy and the enthalpy we are familiar with, respectively. Here
. The entropy is chosen to be the von Neumann entropy of the system
The Gibbs free energy
is defined as
. These definitions are in exact parallel to those in the classical formulation [
27].
7.2. “Partial Molar” Representation
In contrast to the bare representation, we can alternatively define the operator
of the system
that corresponds to
of the bath
by
. The last equality results from the fact that
has no dependence on the external parameter
J. Owing to the non-commutativity of operators, the micro-physics interpretation of the operator
is not so transparent. We first focus on its quantum expectation value
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
. The advantage of this expression is that the observation of
enables us to write
as
, if we have defined the corresponding expectation values for the composite
and the bath
by
and
. In particular we can check that
indeed is the expectation value of the operator
, that is,
. The latter expression can nicely bridge with
for the composite,
. Thus the expectation value
is additive. Its value for the combined systems is equal to the sum of those of the subsystems,
. 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 as the change of due to the intervention of the system . 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 when J is identified as the constant pressure applied to the wall.
From this example, it is tempting to identify
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
due to its interaction with the bath when the composite is acted upon by an external agent
J, since
is related to
[
63,
64]. Inspired by this observation and taking the hint from the definition of
, we introduce the enthalpy of the system
by
where we have identified the enthalpies of the composite
and the bath
as
and
. We may rewrite them as
and
. It implies that (1) the system enthalpy can be decomposed as
and (2) the internal energy
of the system
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 and , and then we may conclude . Thus the internal energy also includes contributions that naïvely we will not ordinarily attribute to the system, such as . 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
, the system
, and the bath
by
, where
,
s and
b, and they obey the additive property of the Gibbs energy,
. Furthermore, in the composite, we note that
From (
45), we can consistently define the entropy
of the composite by
and, similarly, the entropy
of the bath:
The additive property of the free energy and the enthalpy implies that the entropy
of the system in this representation is also additive,
, 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.
7.3. 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.
7.3.1. Enthalpy and Energy Operators: Caution
In fact, we may deduce the operator form of the quantities introduced earlier. For example, we may intuitively define the enthalpy operator of the composite by , and then it is clear to see that the expectation value is related to this operator by . Likewise, the enthalpy operator of the bath can be defined by , and its expectation value gives . Moreover, the internal energy operator of the composite system and the expectation value can be chosen such that such that . For the bath, the internal energy operator is, intuitively, with expectation values that is consistent with the expressions of the internal energy discussed earlier.
Despite their intuitively appealing appearances, these operator forms of the enthalpies and internal energies are not very useful. Inadvertent use of them may result in errors. For example, we cannot define the enthalpy operator of system
simply by the difference of
and
, since
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,
; and thus (3) the expectation value will not be
. This is because the operators defined this way act on Hilbert spaces different from that of
;
is an operator in the Hilbert space of the composite while
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.
7.3.2. System Enthalpy Operator: Approved
We first inspect the internal energy operator. Since the averaged internal energy of the composite system is given by
, 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
and the identity for the operator
If we define an internal energy operator
by
then we obtain a new representation of
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
. With the help of (
A9), this is equivalent to Equation (28) of [
26] in the
case. In addition, we note that
is an operator, not a
c number. Since
, we may define the operator
by
such that
. 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
. Indeed, using the identity operator
in the Hilbert space of the system
we can also define
as
.
In the same fashion, we may rewrite
by
Thus we can define
so that
. We then can have a local form for the
given by
in close resemblance to its classical expression in [
27], if we re-define
as
. The expectation value of
is then given by
.
Now we proceed with constructing a local form of the enthalpy operator of the system. From (
52) and the definition of the operator
, we claim that the local form
is
We can straightforwardly show that
. Thus we have succeeded in constructing the operators that correspond to
,
,
in forms local in the Hilbert space of the system
. 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
of the system, we can always attach a system operator
that satisfies
to the definitions of those local operators, that is, any system operator that has a zero mean. The choice of
may not be unique in the sense that in the basis
that diagonalizes
, we can write
as
It says that the vectors that are respectively composed of the diagonal elements of
and
are orthogonal, but it does not place any restriction on the off-diagonal elements of
on this basis.